Figure 3.3b

Title:

Estimated prevalence of past-year hypertension and corresponding 95% confidence interval across age, separately by male (solid green line) and female (dashed purple line).

Syntax:

For generating a plot of the estimated odds of hypertension among male adults only:

Copy to Clipboard

%WeightedTVEM(
data = nesarc_chapter3a,
dist =binary,
cluster = caseid,
time = nage,
dv = highbp,
tvary_effect = int,
knots = 1,
weight=audweight, 
domain = male,
which = 1);

 
%WeightedTVEM(
data = nesarc_chapter3a,
dist =binary,
cluster = caseid,
time = nage,
dv = highbp,
tvary_effect = int,
knots = 1,
weight=audweight, 
domain = male,
which = 1);

For generating a plot of the estimated odds of hypertension among female adults only:

Copy to Clipboard

%WeightedTVEM(
data = nesarc_chapter3a,
dist =binary,
cluster = caseid,
time = nage,
dv = highbp,
tvary_effect = int,
knots = 1,
weight=audweight, 
domain = male,
which = 0);

 
%WeightedTVEM(
data = nesarc_chapter3a,
dist =binary,
cluster = caseid,
time = nage,
dv = highbp,
tvary_effect = int,
knots = 1,
weight=audweight, 
domain = male,
which = 0);

Equation:

\ln \left(\frac{p\left(H I G H B P_{i}\right)}{1+p\left(H I G H B P_{i}\right.}\right)=\beta_{0}(t)

Notes:

To export the male dataset into Excel for plotting purposes:

Copy to Clipboard

data highbp_male;
set tvem_1_plot_data_or;
mprop=exp_int/(1+exp_int);
mprop_l=exp_int_l/(1+exp_int_l);
mprop_u=exp_int_u/(1+exp_int_u); 
run;
 
PROC EXPORT DATA=highbp_male DBMS=xlsx
OUTFILE="C:\FolderName\highbp_male.xlsx";
run;

 
data highbp_male;
set tvem_1_plot_data_or;
mprop=exp_int/(1+exp_int);
mprop_l=exp_int_l/(1+exp_int_l);
mprop_u=exp_int_u/(1+exp_int_u); 
run;
 
PROC EXPORT DATA=highbp_male DBMS=xlsx
OUTFILE="C:\FolderName\highbp_male.xlsx";
run;

To export the female dataset into Excel:

Copy to Clipboard

data highbp_female
set tvem_0_plot_data_or;
fprop=exp_int/(1+exp_int);
fprop_l=exp_int_l/(1+exp_int_l);
fprop_u=exp_int_u/(1+exp_int_u);
run;
 
PROC EXPORT DATA=highbp_female DBMS=xlsx 
OUTFILE=”C:\FolderName\highbp_female.xlsx”
run;

 
data highbp_female
set tvem_0_plot_data_or;
fprop=exp_int/(1+exp_int);
fprop_l=exp_int_l/(1+exp_int_l);
fprop_u=exp_int_u/(1+exp_int_u);
run;
 
PROC EXPORT DATA=highbp_female DBMS=xlsx 
OUTFILE=”C:\FolderName\highbp_female.xlsx”
run;

Return to all Figures HERE.

Time-Varying Effect Modeling for the Behavioral, Social, and Health Sciences

Authors: Stephanie Lanza and Ashley Linden-Carmichael

Purchase This Book

Reprinted by permission from Springer Nature: Generalized Time-Varying Effect Models for Binary and Count Outcomes by Stephanie T. Lanza and Ashley N. Linden-Carmichael, 2021.