set more off
set linesize 150
global filesave "C:\Dropbox\PilesOfVariance\Chapter8\STATA"
use "$filesave\STATA_Chapter8.dta", clear
egen pmmood = mean(mood), by(personid)
egen pmstressor = mean(stressor), by(personid)
label variable pmstressor "pmstressor: Person Mean Stressors"
label variable pmmood "pmmood: Person Mean Negative Mood"
gen age80 = baseage - 80
gen mood2 = mood - 2
gen wpmood = mood - pmmood
gen pmmood2 = pmmood - 2
gen wpstressor = stressor - pmstressor
gen pmstressor40 = pmstressor - .40
label variable age80 "age80: Baseline Age (0=80)"
label variable mood2 "mood2: Daily Negative Mood (0=2)"
label variable wpmood "wpmood: Within-Person Negative Mood (0=PM)"
label variable pmmood2 "pmmood2: Person Mean Negative Mood (0=2)"
label variable wpstressor "wpstressor: Within-Person Stressors (0=PM)"
label variable pmstressor40 "pmstressor40: Mean #Days with Stressor (0=0.40)"
egen nummiss = rowmiss(women baseage symptoms mood stressor mood)
drop if nummiss>0
drop if studyday>14
log using $filesave\STATA_Chapter8_Output, replace name(STATA_Chapter8)
display as result "Chapter 8: Descriptive Statistics for Time-Invariant Variables"
preserve
collapse women baseage mood stressor, by(personid)
summarize women baseage mood stressor
restore
display as result "Chapter 8: Descriptive Statistics for Time-Varying Variables"
summarize symptoms mood stressor wpmood wpstressor
display as result "Eq 8.1: Conditional Baseline Model"
display as result "Just Sex and Age in the Model for the Means"
mixed symptoms c.women c.age80 c.women#c.age80,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
estimates store FitSexAge,
predict PredSexAge, xb,
display as result "Ch 8: Empty Means, Random Intercept Model for Negative Mood Predictor"
mixed mood ,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
estat icc,
estat wcorrelation, covariance,
estat wcorrelation,
display as result "Eq 8.3: Adding BP and WP Effects of Negative Mood to the Model for the Means"
display as result "Using Person-Mean-Centering for Negative Mood"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
test (c.pmmood2=0) (c.wpmood=0)
lincom c.age80*1 + c.women#c.age80*0
lincom c.age80*1 + c.women#c.age80*1
lincom c.wpmood*-1 + c.pmmood2*1
lincom _cons*1 + c.pmmood2*0 + c.wpmood*-1
lincom _cons*1 + c.pmmood2*0 + c.wpmood*0
lincom _cons*1 + c.pmmood2*0 + c.wpmood*1
lincom _cons*1 + c.pmmood2*1 + c.wpmood*-1
lincom _cons*1 + c.pmmood2*1 + c.wpmood*0
lincom _cons*1 + c.pmmood2*1 + c.wpmood*1
lincom _cons*1 + c.pmmood2*-1 + c.wpmood*0
lincom _cons*1 + c.pmmood2*0 + c.wpmood*0
lincom _cons*1 + c.pmmood2*1 + c.wpmood*0
lincom _cons*1 + c.pmmood2*2 + c.wpmood*0
estimates store FitBPWPMood,
lrtest FitBPWPMood FitSexAge,
predict PredBPWPMood, xb,
corr symptoms PredBPWPMood
display as result "Eq 8.5: Adding Smushed Effect of Negative Mood to the Model for the Means"
display as result "Using Grand-Mean-Centering for Negative Mood"
mixed symptoms c.women c.age80 c.women#c.age80
c.mood2,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
lincom c.age80*1 + c.women#c.age80*0
lincom c.age80*1 + c.women#c.age80*1
display as result "Eq 8.7: Adding a Contextual Effect of Negative Mood to the Model for the Means"
display as result "Using Person-Mean-Centering for Negative Mood"
mixed symptoms c.women c.age80 c.women#c.age80
c.mood2 c.pmmood2,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
lincom c.age80*1 + c.women#c.age80*0
lincom c.age80*1 + c.women#c.age80*1
lincom c.mood2*1 + c.pmmood2*1
estimates store FitBPWPMood,
predict PredBPWPMood2, xb,
corr symptoms PredBPWPMood2
display as result "Eq 8.9: Adding a Random Effect of Negative Mood to the Model for the Variance"
display as result "Using Person-Mean-Centering for Negative Mood"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood,
|| personid: wpmood, variance mle covariance(unstructured),
estat ic, n(105),
estat recovariance, relevel(personid),
estat recovariance, relevel(personid) correlation,
estimates store FitRandPMC,
lrtest FitRandPMC FitBPWPMood,
display as result "Eq 8.9: Adding a Random Effect of Negative Mood to the Model for the Variance"
display as result "Using Grand-Mean-Centering for Negative Mood"
mixed symptoms c.women c.age80 c.women#c.age80
c.mood2 c.pmmood2,
|| personid: mood2, variance mle covariance(unstructured),
estat ic, n(105),
estat recovariance, relevel(personid),
estat recovariance, relevel(personid) correlation,
estimates store FitRandGMC,
lrtest FitRandGMC FitBPWPMood,
display as result "Eq 8.10: Adding Contextual and WP Effects of Stressors to the Model for the Means"
display as result "Using Person-Mean-Centering for Negative Mood, Grand-Mean-Centering for Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood c.pmstressor40 c.stressor,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
test (c.pmstressor40=0) (c.stressor=0)
lincom _cons*1 + c.women*1
lincom c.age80*1 + c.women#c.age80*1
lincom c.wpmood*-1 + c.pmmood2*1
lincom c.stressor*1 + c.pmstressor40*1
lincom c.pmmood2*.2
estimates store FitBPWPStressor,
lrtest FitBPWPStressor FitBPWPMood,
lrtest FitBPWPStressor FitSexAge,
predict PredBPWPStressor, xb,
corr symptoms PredBPWPStressor
display as result "Ch 8: Adding Random WP Effect of Stressors to the Model for the Variance"
display as result "Using Person-Mean-Centering for Negative Mood, Grand-Mean-Centering for Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood c.pmstressor40 c.stressor,
|| personid: stressor, variance mle covariance(unstructured),
estat ic, n(105),
estat recovariance, relevel(personid),
estat recovariance, relevel(personid) correlation,
estimates store FitBPWPStressorR,
lrtest FitBPWPStressorR FitBPWPStressor,
display as result "Eq 8.11: Adding 4 Sex*Negative Mood and Sex*Stressors Interactions to the Model for the Means"
display as result "Using Person-Mean-Centering for Negative Mood, Grand-Mean-Centering for Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood c.pmstressor40 c.stressor
c.women#c.pmmood2 c.women#c.wpmood c.women#c.pmstressor40 c.women#c.stressor,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
test (c.women#c.pmmood2=0) (c.women#c.wpmood=0) (c.women#c.pmstressor40=0) (c.women#c.stressor=0)
lincom _cons*1 + c.women*0
lincom _cons*1 + c.women*1
lincom c.women*1
lincom c.age80*1 + c.women#c.age80*0
lincom c.age80*1 + c.women#c.age80*1
lincom c.women#c.age80*1
lincom c.pmmood2*1 + c.women#c.pmmood2*0
lincom c.pmmood2*1 + c.women#c.pmmood2*1
lincom c.women#c.pmmood2*1
lincom c.wpmood*-1 + c.women#c.wpmood *0 + c.pmmood2*1 + c.women#c.pmmood2*0
lincom c.wpmood*-1 + c.women#c.wpmood*-1 + c.pmmood2*1 + c.women#c.pmmood2*1
lincom c.women#c.wpmood*-1 + c.women#c.pmmood2*1
lincom c.wpmood*1 + c.women#c.wpmood*0
lincom c.wpmood*1 + c.women#c.wpmood*1
lincom c.women#c.wpmood*1
lincom c.stressor*1 + c.women#c.stressor*0 + c.pmstressor40*1 + c.women#c.pmstressor40*0
lincom c.stressor*1 + c.women#c.stressor*1 + c.pmstressor40*1 + c.women#c.pmstressor40*1
lincom c.women#c.stressor*1 + c.women#c.pmstressor40*1
lincom c.pmstressor40*1 + c.women#c.pmstressor40*0
lincom c.pmstressor40*1 + c.women#c.pmstressor40*1
lincom c.women#c.pmstressor40*1
lincom c.stressor*1 + c.women#c.stressor*0
lincom c.stressor*1 + c.women#c.stressor*1
lincom c.women#c.stressor*1
lincom c.pmstressor40*.2
estimates store FitSex4,
lrtest FitSex4 FitBPWPStressor,
predict PredSex4, xb,
corr symptoms PredSex4
display as result "Eq 8.11reduced: Adding Only 2 Sex*Stressors Interactions to the Model for the Means"
display as result "Using Person-Mean-Centering for Negative Mood, Grand-Mean-Centering for Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood c.pmstressor40 c.stressor
c.women#c.pmstressor40 c.women#c.stressor,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
test (c.women#c.pmstressor40=0) (c.women#c.stressor=0)
estimates store FitSex2,
lrtest FitSex2 FitBPWPStressor,
predict PredSex2, xb,
corr symptoms PredSex2
display as result "Eq 8.12: Adding 4 Negative Mood*Stressors Interactions to the Model for the Means"
display as result "Using Person-Mean-Centering for Negative Mood, Grand-Mean-Centering for Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood c.pmstressor40 c.stressor
c.women#c.pmstressor40 c.women#c.stressor
c.pmmood2#c.stressor c.pmmood2#c.pmstressor40 c.wpmood#c.stressor c.wpmood#c.pmstressor40,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
test (c.pmmood2#c.stressor=0) (c.pmmood2#c.pmstressor40=0) (c.wpmood#c.stressor=0) (c.wpmood#c.pmstressor40=0)
lincom _cons*1 + c.women*0
lincom _cons*1 + c.women*1
lincom c.women*1
lincom c.age80*1 + c.women#c.age80*0
lincom c.age80*1 + c.women#c.age80*1
lincom c.women#c.age80*1
lincom c.pmmood2*1
lincom c.wpmood*-1 + c.pmmood2*1
lincom c.wpmood*1
lincom c.stressor*1 + c.women#c.stressor*0 + c.pmstressor40*1 + c.women#c.pmstressor40*0
lincom c.stressor*1 + c.women#c.stressor*1 + c.pmstressor40*1 + c.women#c.pmstressor40*1
lincom c.women#c.stressor*1 + c.women#c.pmstressor40*1
lincom c.pmstressor40*1 + c.women#c.pmstressor40*0
lincom c.pmstressor40*1 + c.women#c.pmstressor40*1
lincom c.women#c.pmstressor40*1
lincom c.stressor*1 + c.women#c.stressor*0
lincom c.stressor*1 + c.women#c.stressor*1
lincom c.women#c.stressor*1
lincom c.pmstressor40*.2
lincom c.pmmood2#c.stressor*1 + c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.stressor*1
lincom c.wpmood#c.stressor*-1 + c.pmmood2#c.stressor*1 + c.wpmood#c.pmstressor40*-1 + c.pmmood2#c.pmstressor40*1
lincom c.wpmood#c.pmstressor40*-1 + c.pmmood2#c.pmstressor40*1
lincom c.wpmood#c.stressor*-1 + c.pmmood2#c.stressor*1
lincom c.wpmood#c.stressor*1 + c.wpmood#c.pmstressor40*1
lincom c.wpmood#c.pmstressor40*1
lincom c.wpmood#c.stressor*1
estimates store FitMoodStressor4,
lrtest FitMoodStressor4 FitSex2,
predict PredMoodStressor4, xb,
corr symptoms PredMoodStressor4
display as result "Ch 8: Checking Inter-Variable Interactions via PMC Negative Mood by PMC Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood c.pmstressor40 c.wpstressor
c.women#c.pmstressor40 c.women#c.wpstressor
c.pmmood2#c.wpstressor c.pmmood2#c.pmstressor40 c.wpmood#c.wpstressor c.wpmood#c.pmstressor40,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
lincom c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.wpstressor*-1 + c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.wpstressor*1
lincom c.wpmood#c.pmstressor40*-1 + c.pmmood2#c.pmstressor40*1
lincom c.wpmood#c.pmstressor40*-1 + c.pmmood2#c.pmstressor40*1 + c.pmmood2#c.wpstressor*-1 + c.wpmood#c.wpstressor*1
lincom c.wpmood#c.wpstressor*-1 + c.pmmood2#c.wpstressor*1
lincom c.wpmood#c.pmstressor40*1
lincom c.wpmood#c.wpstressor*-1 + c.wpmood#c.pmstressor40*1
lincom c.wpmood#c.wpstressor*1
display as result "Ch 8: Checking Inter-Variable Interactions via GMC Negative Mood by GMC Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.mood2 c.pmstressor40 c.stressor
c.women#c.pmstressor40 c.women#c.stressor
c.pmmood2#c.stressor c.pmmood2#c.pmstressor40 c.mood2#c.stressor c.mood2#c.pmstressor40,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
lincom c.mood2#c.pmstressor40*1 + c.pmmood2#c.pmstressor40*1 + c.stressor#c.pmmood2*1 + c.mood2#c.stressor*1
lincom c.mood2#c.pmstressor40*1 + c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.stressor*1 + c.mood2#c.stressor*1
lincom c.pmmood2#c.stressor*1 + c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.stressor*1
lincom c.mood2#c.pmstressor40*1 + c.mood2#c.stressor*1
lincom c.mood2#c.pmstressor40*1
lincom c.mood2#c.stressor*1
display as result "Ch 8: Checking Inter-Variable Interactions via GMC Negative Mood by PMC Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.mood2 c.pmstressor40 c.wpstressor
c.women#c.pmstressor40 c.women#c.wpstressor
c.pmmood2#c.wpstressor c.pmmood2#c.pmstressor40 c.mood2#c.wpstressor c.mood2#c.pmstressor40,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
lincom c.mood2#c.pmstressor40*1 + c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.wpstressor*-1 + c.mood2#c.wpstressor*-1 + c.mood2#c.pmstressor40*1 + c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.wpstressor*1 + c.mood2#c.wpstressor*1
lincom c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.wpstressor*-1 + c.pmmood2#c.pmstressor40*1
lincom c.pmmood2#c.wpstressor*1
lincom c.mood2#c.pmstressor40*1
lincom c.mood2#c.wpstressor*-1 + c.mood2#c.pmstressor40*1
lincom c.mood2#c.wpstressor*1
display as result "Eq 8.13: Adding 5 Intra-Variable Interactions to the Model for the Means"
display as result "Using Person-Mean-Centering for Negative Mood, Grand-Mean-Centering for Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood c.pmstressor40 c.stressor
c.women#c.pmstressor40 c.women#c.stressor
c.pmmood2#c.pmmood2 c.wpmood#c.wpmood c.pmmood2#c.wpmood c.pmstressor40#c.pmstressor40 c.pmstressor40#c.stressor,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
test (c.pmmood2#c.pmmood2=0) (c.wpmood#c.wpmood=0) (c.pmmood2#c.wpmood=0) (c.pmstressor40#c.pmstressor40=0) (c.pmstressor40#c.stressor=0)
lincom c.pmmood2*1
lincom c.pmmood2#c.pmmood2*1
lincom c.wpmood*-1 + c.pmmood2*1
lincom c.pmmood2#c.wpmood*-1 + c.pmmood2#c.pmmood2*1 + c.wpmood#c.wpmood*1
lincom c.wpmood*1
lincom c.wpmood#c.wpmood*1
lincom c.pmmood2#c.wpmood*1
lincom c.pmmood2#c.wpmood*1 + c.wpmood#c.wpmood*-2
lincom c.pmstressor40*1 + c.stressor*1
lincom c.pmstressor40#c.pmstressor40*1 + c.pmstressor40#c.stressor*1
lincom c.pmstressor40*1
lincom c.pmstressor40#c.pmstressor40*1
lincom c.stressor*1
lincom c.pmstressor40#c.stressor*1
lincom c.pmmood2*1 + c.pmmood2#c.pmmood2*-2
lincom c.pmmood2*1 + c.pmmood2#c.pmmood2*2
estimates store FitIntra5,
lrtest FitIntra5 FitSex2,
predict PredIntra5, xb,
corr symptoms PredIntra5
display as result "Eq 8.13reduced: Baseline Homogeneous Variance Model with Significant Fixed Effects Only"
display as result "Using Person-Mean-Centering for Negative Mood, Grand-Mean-Centering for Stressors"
mixed symptoms c.women c.age80 c.women#c.age80
c.pmmood2 c.wpmood c.pmstressor40 c.stressor
c.women#c.pmstressor40 c.women#c.stressor
c.pmmood2#c.pmmood2 c.pmmood2#c.wpmood,
|| personid: , variance mle covariance(unstructured),
estat ic, n(105),
log close STATA_Chapter8
log2html $filesave\STATA_Chapter8_Output, replace