set more off
set linesize 150
global filesave "C:\Dropbox\PilesOfVariance\Chapter3b\STATA"
use "$filesave\STATA_Chapter3b.dta", clear
log using $filesave\STATA_Chapter3b_Output, replace name(STATA_Chapter3b)
display as result "Chapter 3b Example: Means by session for RT outcome"
tabulate session, summarize(rt)
display as result "Eq 3b.10: Between-Person Independent ANOVA"
mixed rt i.session,
|| personid: , noconstant variance reml covariance(unstructured)
residuals(independent,t(session)),
estat ic, n(101),
estat wcorrelation, covariance,
estat wcorrelation,
contrast i.session,
margins i.session,
margins i.session, pwcompare(pveffects)
estimates store FitBPANOVA,
display as result "Eq 3b.10: Univariate Repeated Measures ANOVA"
mixed rt i.session,
|| personid: , noconstant variance reml covariance(unstructured)
residuals(exchangeable,t(session)),
estat ic, n(101),
estat wcorrelation, covariance,
estat wcorrelation,
contrast i.session,
margins i.session,
margins i.session, pwcompare(pveffects)
estimates store FitUnivANOVA,
lrtest FitUnivANOVA FitBPANOVA,
display as result "Eq 3b.10: Multivariate Repeated Measures ANOVA"
mixed rt i.session,
|| personid: , noconstant variance reml covariance(unstructured)
residuals(unstructured,t(session)),
estat ic, n(101),
estat wcorrelation, covariance,
estat wcorrelation,
contrast i.session,
margins i.session,
margins i.session, pwcompare(pveffects)
estimates store FitMultivANOVA,
lrtest FitMultivANOVA FitUnivANOVA,
log close STATA_Chapter3b
log2html $filesave\STATA_Chapter3b_Output, replace