clear all
set more off

**************************************************
*Title: Nuoxintuo

*Authors: Ruolin Zhang
*Date: 2/2024
*Email: rz117@duke.edu
**************************************************
use "/Users/zhangruolin/Desktop/P_/data/Data_english.dta", clear
**************************************************
**************************************************
///screening
///Missing age gender NT-proBNP
keep if !missing(age) & !missing(gender_c) & !missing(NTproBNP_ad) & !missing(NTproBNP_dis)
///ACEI and SV should not be used together
drop if ACEI_c == 1 & SV_c == 1

//////Baseline (Between)
///Demographic(age(t-test), BMI(t-test), hospitalizationtimes(Wilcoxon), gender_c (CHi- square), marriage_c (CHi- square))
swilk age
robvar age, by(SV_c)
ttest age, by(SV_c)

tab SV_c gender_c, chi2
tab SV_c marriage_c1, chi2

swilk BMI
robvar BMI, by(SV_c) 
ttest BMI, by(SV_c)

///Lifestyle (smoking(CHi- square), drinking(CHi- square))
tab SV_c smoking, chi2
tab SV_c drinking, chi2

///Medical History (diabetes(CHi- square), hypertension(CHi- square))
tab SV_c diabetes, chi2
tab SV_c hypertension, chi2

///Complete Blood Count (redblood_ad(t-test), whiteblood_ad(Wilcoxon), hemoglobin_ad(t-test), platelets_ad(t-test))
swilk redblood_ad 
robvar redblood_ad, by(SV_c)
ttest redblood_ad, by(SV_c)

swilk whiteblood_ad
robvar whiteblood_ad, by(SV_c) 
ranksum whiteblood_ad, by(SV_c)
summarize whiteblood_ad if SV_c == 1, detail
summarize whiteblood_ad if SV_c == 0, detail

swilk hemoglobin_ad
robvar hemoglobin_ad, by(SV_c) 
ttest hemoglobin_ad, by(SV_c)

swilk platelets_ad
robvar platelets_ad, by(SV_c) 
ttest platelets_ad, by(SV_c)

///Myocardial Infarction (maindiagnosed_c(chi-squared))
tab SV_c maindiagnosed_c, chi2

///Biochemical Data (triglycerides_ad(ttest), cholesterol_ad(Wilcoxon), HDL_ad(Wilcoxon), LDL_ad(Wilcoxon), creatinine_ad(Wilcoxon), bloodglucose_ad(Wilcoxon))
swilk triglycerides_ad
robvar triglycerides_ad, by(SV_c) 
ranksum triglycerides_ad, by(SV_c)
summarize triglycerides_ad if SV_c == 1, detail
summarize triglycerides_ad if SV_c == 0, detail

swilk cholesterol_ad
robvar cholesterol_ad, by(SV_c) 
ttest cholesterol_ad, by(SV_c)

swilk HDL_ad
robvar HDL_ad, by(SV_c) 
ranksum HDL_ad, by(SV_c)
summarize HDL_ad if SV_c == 1, detail
summarize HDL_ad if SV_c == 0, detail

swilk LDL_ad
robvar LDL_ad, by(SV_c) 
ranksum LDL_ad, by(SV_c)
summarize LDL_ad if SV_c == 1, detail
summarize LDL_ad if SV_c == 0, detail

swilk creatinine_ad
robvar creatinine_ad, by(SV_c) 
ranksum creatinine_ad, by(SV_c)
summarize creatinine_ad if SV_c == 1, detail
summarize creatinine_ad if SV_c == 0, detail

swilk bloodglucose_ad
robvar bloodglucose_ad, by(SV_c) 
ranksum bloodglucose_ad, by(SV_c)
summarize bloodglucose_ad if SV_c == 1, detail
summarize bloodglucose_ad if SV_c == 0, detail

///Cardiac Function-associated Indexes(EF_ad(Wilcoxon), LVEDD_ad(Wilcoxon), NTproBNP_ad(Wilcoxon))
swilk EF_ad
robvar EF_ad, by(SV_c) 
ranksum EF_ad, by(SV_c)
summarize EF_ad if SV_c == 1, detail
summarize EF_ad if SV_c == 0, detail

swilk LVEDD_ad
robvar LVEDD_ad, by(SV_c)
ranksum LVEDD_ad, by(SV_c)
summarize LVEDD_ad if SV_c == 1, detail
summarize LVEDD_ad if SV_c == 0, detail

swilk NTproBNP_ad
robvar NTproBNP_ad, by(SV_c)
ranksum NTproBNP_ad, by(SV_c)
summarize NTproBNP_ad if SV_c == 1, detail
summarize NTproBNP_ad if SV_c == 0, detail

///Blood pressure(SBP_ad(t-test), DBP_ad(t-test))
swilk SBP_ad
robvar SBP_ad, by(SV_c)
ttest SBP_ad, by(SV_c)

swilk DBP_ad
robvar DBP_ad, by(SV_c)
ttest DBP_ad, by(SV_c)

///Hospitalization
swilk hospitalizationtimes
robvar hospitalizationtimes, by(SV_c) 
ranksum hospitalizationtimes, by(SV_c)
summarize hospitalizationtimes if SV_c == 1, detail
summarize hospitalizationtimes if SV_c == 0, detail



//////Comparison of NT-proBNP Levels at Admission, Discharge, and Follow-Up (NTproBNP_ad(Wilcoxon), NTproBNP_dis(Wilcoxon), NTproBNP_drop(Wilcoxon))
replace NTproBNP_drop = NTproBNP_ad - NTproBNP_dis if !missing(NTproBNP_ad) & !missing(NTproBNP_dis) & missing(NTproBNP_drop)

swilk NTproBNP_dis
robvar NTproBNP_dis, by(SV_c)
swilk NTproBNP_drop
robvar NTproBNP_drop, by(SV_c)
///Between
ranksum NTproBNP_ad, by(SV_c)

ranksum NTproBNP_dis, by(SV_c)
summarize NTproBNP_dis if SV_c == 1, detail
summarize NTproBNP_dis if SV_c == 0, detail

ranksum NTproBNP_drop, by(SV_c)
summarize NTproBNP_drop if SV_c == 1, detail
summarize NTproBNP_drop if SV_c == 0, detail
///Within 
signrank NTproBNP_ad = NTproBNP_dis if SV_c == 1
signrank NTproBNP_ad = NTproBNP_dis if SV_c == 0
///follow up
swilk NTproBNP_fu
robvar NTproBNP_fu, by(SV_c)
ranksum NTproBNP_fu, by(SV_c)
summarize NTproBNP_fu if SV_c == 1, detail
summarize NTproBNP_fu if SV_c == 0, detail


//////Comparison of Blood Pressure at Admission, Discharge, and Follow-Up (SBP_ad(t-test), SBP_dis(Wilcoxon), DBP_ad(t-test), DBP_dis(Wilcoxon))
swilk SBP_ad
robvar SBP_ad, by(SV_c)
swilk SBP_dis
robvar SBP_dis, by(SV_c)
swilk DBP_ad
robvar DBP_ad, by(SV_c)
swilk DBP_dis
robvar DBP_dis, by(SV_c)
///Between
ttest SBP_ad, by(SV_c)
ttest SBP_dis, by(SV_c)
ttest DBP_ad, by(SV_c)
ttest DBP_dis, by(SV_c)
///Within
ttest SBP_ad = SBP_dis if SV_c == 1
ttest SBP_ad = SBP_dis if SV_c == 0

ttest DBP_ad = DBP_dis if SV_c == 1
ttest DBP_ad = DBP_dis if SV_c == 0

//////Symptom Improvement in Dyspnea at Hospitalization and Follow-Up (symptomsimprove_ad(chi-squared), symptomsimprove_dis(chi-squared))
///Between
tab SV_c symptomsimprove_ad, chi2
tab SV_c symptomsimprove_dis, chi2

//////One-year Readmission Rate Due to Cardiac Disease for S/V and ACEIs Groups (readmission(chi-squared))
tab SV_c readmission, chi2