*read in antisocial data in wide format; use stata_MLM_data, clear summarize *reshaping data from long to wide; reshape long age anti, i(id) j(year) replace age = age-6 drop if age==. & anti==. tabstat anti, by(age) statistics(n mean sd var) nototal *random intercept multilevel growth model; *homescedastic time-specific errors; mixed anti || id: , reml /// dfmethod(repeated) covariance(unstructured) estat recovariance, correlation *random intercept and slope multilevel growth model; *homescedastic time-specific errors; mixed anti age || id: age, reml /// dfmethod(repeated) covariance(unstructured) estat recovariance, correlation gen age_male = age*male *random intercept and slope multilevel growth model; *homescedastic time-specific errors; *regressed on MALE; mixed anti age male age_male || id: age, reml /// dfmethod(repeated) covariance(unstructured) estat recovariance, correlation