Displaying 1 result from an estimated 1 matches for "4icp".
2004 May 27
1
Getting the same values of adjusted mean and standard errors as SAS
...les at meanAge using
predict
# set up data frame to get predicted values at meanAge
evalDF <- data.frame(Age = meanAge, Gender = c("F","M"),
Group = c("1NC","1NC", "2UCLP", "2UCLP", "3BCLP", "3BCLP",
"4ICP", "4ICP", "5CLPP", "5CLPP"))
mod <-lm(Measurement ~ Age + Gender + Group, data=dd)
pred <- predict(mod, evalDF, se.fit = TRUE)
adjDF <- data.frame(evalDF,fit = pred$fit, se = pred$se.fit)
adjMeanMale <- mean(adjDF$fit[adjDF$Gender=="M"]); # m...