## library(metaSEM) ## ## Clear all objects ## rm(list=ls()) ## ## Set seed for replication ## set.seed(1000000) ## ## Copy Bornmann07 to my.df ## my.df <- Bornmann07 ## ## "Fellowship": 1; "Grant": 0 ## my.df$Type_MCAR <- ifelse(Bornmann07$Type=="Fellowship", yes=1, no=0) ## ## Create 17 out of 66 missingness with MCAR ## my.df$Type_MCAR[sample(1:66, 17)] <- NA ## metaSEM.fit <- meta3X(y=logOR, v=v, cluster=Cluster, x2=Type_MCAR, data=my.df) ## ## Get the OpenMx model ## OpenMx.model <- metaSEM.fit$mx.model ## save(metaSEM.fit, OpenMx.model, my.df, file="test5b.txt") ## Version "2.0.0.3578" library(OpenMx) sessionInfo() ## Clear all objects rm(list=ls()) ## Load the workspace load("test5b.txt") ## Got an error OpenMx.fit <- mxRun(OpenMx.model)