if(requireNamespace("GDINA")){
library(GDINA)
dat <- sim10GDINA$simdat
Q <- matrix(c(1,0,0,
0,1,0,
0,0,1,
1,0,1,
0,1,1,
1,1,0,
0,0,1,
1,0,0,
1,1,1,
1,0,1),byrow = TRUE,ncol = 3)
fit.object <- GDINA(dat = dat, Q = Q, model = "GDINA", att.dist = "independent",verbose = FALSE)
CEP <- CEP_1t(fit.object)
}
Run the code above in your browser using DataLab