Hi, I have a dataset of outcome measures at difference time points and an anchor score. I want to calculate the 95% CI for difference in means between two groups who score a certain value on an external anchor criterion. the code I am using is below, but the ci command just brings up an error message. Can soemone please assist. PS. I'm a stata novice.
quietly sum pem_change_3mths if groc_3mths==1 | groc_3mths<0
local mean_change = r(mean)
dis `mean_change'
quietly sum pem_change_3mths if groc_3mths==0
local mean_unchanged = r(mean)
dis `mean_unchanged'
local betweenp = `mean_change' - `mean_unchanged'
dis `betweenp'
ci means `betweenp'
quietly sum pem_change_3mths if groc_3mths==1 | groc_3mths<0
local mean_change = r(mean)
dis `mean_change'
quietly sum pem_change_3mths if groc_3mths==0
local mean_unchanged = r(mean)
dis `mean_unchanged'
local betweenp = `mean_change' - `mean_unchanged'
dis `betweenp'
ci means `betweenp'
Comment