## 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.fit <- metaSEM.fit$mx.fit

## save(metaSEM.fit, OpenMx.fit, my.df, file="test5a.txt")

## Version "1.4-3532"
library(OpenMx)

sessionInfo()

## Clear all objects
rm(list=ls())

## Load the workspace
load("test5a.txt")

summary(OpenMx.fit)