library(OpenMx)
library(merror)
data(redshift)
merror.pairs(redshift)
# estimates of accuracy and precision
# parameters for a one-factor
# measurement error model
head(redshift)
merror.pairs(redshift)
if (FALSE) {
red <- omx(redshift[,-1],reps=200) # Drop the spectroscopic measurements
summary(red$fit)
red$ci
red$q.boot
cplot(redshift[,-1],1,2,alpha.beta.sigma=red$abs)
}
Run the code above in your browser using DataLab