library(metaSEM)
## Loading required package: OpenMx
## "SLSQP" is set as the default optimizer in OpenMx.
## mxOption(NULL, "Gradient algorithm") is set at "central".
## mxOption(NULL, "Optimality tolerance") is set at "6.3e-14".
## mxOption(NULL, "Gradient iterations") is set at "2".
library(lavaan)
## This is lavaan 0.6-3
## lavaan is BETA software! Please report any bugs.
##
## Attaching package: 'lavaan'
## The following object is masked from 'package:OpenMx':
##
## vech
library(semPlot)
## Warning: package 'semPlot' was built under R version 3.6.1
## Registered S3 methods overwritten by 'huge':
## method from
## plot.sim BDgraph
## print.sim BDgraph
source("http://www.suzannejak.nl/subgroup.functions.R")
cat("
1.0 -0.27 -0.23 1.0 0.27 1.0\n
1.0 -0.25 -0.29 1.0 0.35 1.0\n
1.0 -0.12 -0.24 1.0 0.15 1.0\n
1.0 -0.17 -0.12 1.0 0.32 1.0\n
1.0 -0.37 -0.19 1.0 0.27 1.0\n
1.0 -0.1 -0.19 1.0 0.33 1.0\n
1.0 -0.13 0 1.0 0.1 1.0\n
1.0 0.1 0.16 1.0 0.17 1.0\n
1.0 -0.38 -0.36 1.0 0.42 1.0\n
1.0 -0.19 -0.13 1.0 0.22 1.0\n
1.0 -0.13 -0.19 1.0 0.27 1.0\n
1.0 -0.47 -0.41 1.0 0.37 1.0\n
1.0 -0.21 -0.25 1.0 0.38 1.0\n
1.0 -0.21 -0.2 1.0 0.39 1.0\n
1.0 -0.2 -0.16 1.0 0.2 1.0\n
1.0 -0.11 -0.23 1.0 0.19 1.0\n
1.0 -0.28 -0.19 1.0 0.29 1.0\n
1.0 -0.18 -0.21 1.0 0.29 1.0\n
1.0 -0.27 -0.23 1.0 0.33 1.0\n
1.0 -0.28 -0.24 1.0 0.59 1.0\n
1.0 -0.31 -0.36 1.0 0.59 1.0\n
1.0 -0.68 -0.23 1.0 0.39 1.0\n
1.0 -0.3 -0.23 1.0 0.26 1.0\n
1.0 -0.18 -0.31 1.0 0.46 1.0\n
1.0 -0.26 -0.2 1.0 0.31 1.0\n
1.0 -0.08 -0.19 1.0 0.37 1.0\n
1.0 -0.38 0.14 1.0 0.28 1.0\n
1.0 -0.2 -0.15 1.0 0.31 1.0\n
1.0 -0.11 -0.07 1.0 0.15 1.0\n
1.0 -0.22 -0.37 1.0 0.25 1.0\n
1.0 0.04 -0.03 1.0 0.24 1.0\n
1.0 -0.22 -0.27 1.0 0.56 1.0\n
1.0 -0.24 -0.19 1.0 0.14 1.0\n
1.0 -0.31 -0.36 1.0 0.27 1.0\n
1.0 -0.28 -0.23 1.0 0.27 1.0\n
1.0 -0.2 -0.3 1.0 0.26 1.0\n
1.0 -0.13 -0.09 1.0 0.09 1.0\n
1.0 -0.28 -0.15 1.0 0.39 1.0\n
1.0 -0.18 -0.2 1.0 0.22 1.0\n
1.0 -0.32 -0.29 1.0 0.51 1.0\n
1.0 -0.61 -0.52 1.0 0.45 1.0\n
1.0 -0.16 -0.07 1.0 0.17 1.0\n
1.0 -0.1 -0.03 1.0 0.08 1.0\n
1.0 -0.14 -0.18 1.0 0.31 1.0\n
1.0 0.04 -0.12 1.0 0.46 1.0\n
1.0 -0.22 -0.1 1.0 0.49 1.0\n
1.0 -0.19 -0.26 1.0 0.38 1.0\n
1.0 -0.45 -0.45 1.0 0.74 1.0\n
1.0 -0.19 -0.1 1.0 0.13 1.0\n
1.0 -0.18 -0.23 1.0 0.39 1.0\n
1.0 -0.11 -0.08 1.0 0.31 1.0\n
1.0 -0.13 -0.16 1.0 0.43 1.0\n
1.0 -0.14 -0.22 1.0 0.38 1.0\n
1.0 -0.14 -0.12 1.0 0.5 1.0\n
1.0 -0.22 -0.31 1.0 0.49 1.0\n
1.0 0 -0.21 1.0 0.27 1.0\n
1.0 -0.1 -0.22 1.0 0.35 1.0\n
1.0 -0.17 -0.2 1.0 0.45 1.0\n
1.0 -0.28 -0.24 1.0 0.59 1.0\n
1.0 -0.03 -0.2 1.0 0.35 1.0\n
1.0 -0.18 -0.18 1.0 0.34 1.0\n
1.0 -0.3 -0.29 1.0 0.45 1.0\n
1.0 -0.23 -0.04 1.0 0.08 1.0\n
1.0 -0.03 0.15 1.0 0.12 1.0\n
1.0 -0.1 -0.01 1.0 0.06 1.0\n
1.0 -0.27 -0.3 1.0 -0.39 1.0\n
1.0 -0.17 -0.09 1.0 0.26 1.0\n
1.0 -0.14 -0.19 1.0 0.27 1.0\n
1.0 -0.35 -0.19 1.0 0.4 1.0\n
1.0 -0.2 -0.23 1.0 0.23 1.0\n
1.0 -0.03 -0.25 1.0 0.02 1.0\n
1.0 -0.26 -0.1 1.0 0.22 1.0\n
1.0 -0.16 -0.13 1.0 0.45 1.0\n
1.0 -0.29 -0.26 1.0 0.45 1.0\n
1.0 -0.17 -0.01 1.0 0.12 1.0\n
1.0 -0.36 -0.07 1.0 0.17 1.0\n
1.0 -0.19 -0.19 1.0 0.31 1.0\n
1.0 -0.15 -0.28 1.0 0.29 1.0\n
1.0 -0.18 -0.31 1.0 0.46 1.0\n
1.0 -0.1 -0.08 1.0 0.16 1.0\n
1.0 -0.12 -0.08 1.0 0.1 1.0\n
1.0 -0.18 -0.3 1.0 0.19 1.0\n
1.0 -0.06 -0.08 1.0 -0.28 1.0\n
1.0 -0.23 -0.15 1.0 0.33 1.0\n
1.0 -0.15 -0.15 1.0 0.2 1.0\n",
file="stackvec.dat", sep="")
my.vec <- readStackVec("stackvec.dat")
my.vec
## $`1`
## x1 x2 x3
## x1 1.00 -0.27 -0.23
## x2 -0.27 1.00 0.27
## x3 -0.23 0.27 1.00
##
## $`2`
## x1 x2 x3
## x1 1.00 -0.25 -0.29
## x2 -0.25 1.00 0.35
## x3 -0.29 0.35 1.00
##
## $`3`
## x1 x2 x3
## x1 1.00 -0.12 -0.24
## x2 -0.12 1.00 0.15
## x3 -0.24 0.15 1.00
##
## $`4`
## x1 x2 x3
## x1 1.00 -0.17 -0.12
## x2 -0.17 1.00 0.32
## x3 -0.12 0.32 1.00
##
## $`5`
## x1 x2 x3
## x1 1.00 -0.37 -0.19
## x2 -0.37 1.00 0.27
## x3 -0.19 0.27 1.00
##
## $`6`
## x1 x2 x3
## x1 1.00 -0.10 -0.19
## x2 -0.10 1.00 0.33
## x3 -0.19 0.33 1.00
##
## $`7`
## x1 x2 x3
## x1 1.00 -0.13 0.0
## x2 -0.13 1.00 0.1
## x3 0.00 0.10 1.0
##
## $`8`
## x1 x2 x3
## x1 1.00 0.10 0.16
## x2 0.10 1.00 0.17
## x3 0.16 0.17 1.00
##
## $`9`
## x1 x2 x3
## x1 1.00 -0.38 -0.36
## x2 -0.38 1.00 0.42
## x3 -0.36 0.42 1.00
##
## $`10`
## x1 x2 x3
## x1 1.00 -0.19 -0.13
## x2 -0.19 1.00 0.22
## x3 -0.13 0.22 1.00
##
## $`11`
## x1 x2 x3
## x1 1.00 -0.13 -0.19
## x2 -0.13 1.00 0.27
## x3 -0.19 0.27 1.00
##
## $`12`
## x1 x2 x3
## x1 1.00 -0.47 -0.41
## x2 -0.47 1.00 0.37
## x3 -0.41 0.37 1.00
##
## $`13`
## x1 x2 x3
## x1 1.00 -0.21 -0.25
## x2 -0.21 1.00 0.38
## x3 -0.25 0.38 1.00
##
## $`14`
## x1 x2 x3
## x1 1.00 -0.21 -0.20
## x2 -0.21 1.00 0.39
## x3 -0.20 0.39 1.00
##
## $`15`
## x1 x2 x3
## x1 1.00 -0.2 -0.16
## x2 -0.20 1.0 0.20
## x3 -0.16 0.2 1.00
##
## $`16`
## x1 x2 x3
## x1 1.00 -0.11 -0.23
## x2 -0.11 1.00 0.19
## x3 -0.23 0.19 1.00
##
## $`17`
## x1 x2 x3
## x1 1.00 -0.28 -0.19
## x2 -0.28 1.00 0.29
## x3 -0.19 0.29 1.00
##
## $`18`
## x1 x2 x3
## x1 1.00 -0.18 -0.21
## x2 -0.18 1.00 0.29
## x3 -0.21 0.29 1.00
##
## $`19`
## x1 x2 x3
## x1 1.00 -0.27 -0.23
## x2 -0.27 1.00 0.33
## x3 -0.23 0.33 1.00
##
## $`20`
## x1 x2 x3
## x1 1.00 -0.28 -0.24
## x2 -0.28 1.00 0.59
## x3 -0.24 0.59 1.00
##
## $`21`
## x1 x2 x3
## x1 1.00 -0.31 -0.36
## x2 -0.31 1.00 0.59
## x3 -0.36 0.59 1.00
##
## $`22`
## x1 x2 x3
## x1 1.00 -0.68 -0.23
## x2 -0.68 1.00 0.39
## x3 -0.23 0.39 1.00
##
## $`23`
## x1 x2 x3
## x1 1.00 -0.30 -0.23
## x2 -0.30 1.00 0.26
## x3 -0.23 0.26 1.00
##
## $`24`
## x1 x2 x3
## x1 1.00 -0.18 -0.31
## x2 -0.18 1.00 0.46
## x3 -0.31 0.46 1.00
##
## $`25`
## x1 x2 x3
## x1 1.00 -0.26 -0.20
## x2 -0.26 1.00 0.31
## x3 -0.20 0.31 1.00
##
## $`26`
## x1 x2 x3
## x1 1.00 -0.08 -0.19
## x2 -0.08 1.00 0.37
## x3 -0.19 0.37 1.00
##
## $`27`
## x1 x2 x3
## x1 1.00 -0.38 0.14
## x2 -0.38 1.00 0.28
## x3 0.14 0.28 1.00
##
## $`28`
## x1 x2 x3
## x1 1.00 -0.20 -0.15
## x2 -0.20 1.00 0.31
## x3 -0.15 0.31 1.00
##
## $`29`
## x1 x2 x3
## x1 1.00 -0.11 -0.07
## x2 -0.11 1.00 0.15
## x3 -0.07 0.15 1.00
##
## $`30`
## x1 x2 x3
## x1 1.00 -0.22 -0.37
## x2 -0.22 1.00 0.25
## x3 -0.37 0.25 1.00
##
## $`31`
## x1 x2 x3
## x1 1.00 0.04 -0.03
## x2 0.04 1.00 0.24
## x3 -0.03 0.24 1.00
##
## $`32`
## x1 x2 x3
## x1 1.00 -0.22 -0.27
## x2 -0.22 1.00 0.56
## x3 -0.27 0.56 1.00
##
## $`33`
## x1 x2 x3
## x1 1.00 -0.24 -0.19
## x2 -0.24 1.00 0.14
## x3 -0.19 0.14 1.00
##
## $`34`
## x1 x2 x3
## x1 1.00 -0.31 -0.36
## x2 -0.31 1.00 0.27
## x3 -0.36 0.27 1.00
##
## $`35`
## x1 x2 x3
## x1 1.00 -0.28 -0.23
## x2 -0.28 1.00 0.27
## x3 -0.23 0.27 1.00
##
## $`36`
## x1 x2 x3
## x1 1.0 -0.20 -0.30
## x2 -0.2 1.00 0.26
## x3 -0.3 0.26 1.00
##
## $`37`
## x1 x2 x3
## x1 1.00 -0.13 -0.09
## x2 -0.13 1.00 0.09
## x3 -0.09 0.09 1.00
##
## $`38`
## x1 x2 x3
## x1 1.00 -0.28 -0.15
## x2 -0.28 1.00 0.39
## x3 -0.15 0.39 1.00
##
## $`39`
## x1 x2 x3
## x1 1.00 -0.18 -0.20
## x2 -0.18 1.00 0.22
## x3 -0.20 0.22 1.00
##
## $`40`
## x1 x2 x3
## x1 1.00 -0.32 -0.29
## x2 -0.32 1.00 0.51
## x3 -0.29 0.51 1.00
##
## $`41`
## x1 x2 x3
## x1 1.00 -0.61 -0.52
## x2 -0.61 1.00 0.45
## x3 -0.52 0.45 1.00
##
## $`42`
## x1 x2 x3
## x1 1.00 -0.16 -0.07
## x2 -0.16 1.00 0.17
## x3 -0.07 0.17 1.00
##
## $`43`
## x1 x2 x3
## x1 1.00 -0.10 -0.03
## x2 -0.10 1.00 0.08
## x3 -0.03 0.08 1.00
##
## $`44`
## x1 x2 x3
## x1 1.00 -0.14 -0.18
## x2 -0.14 1.00 0.31
## x3 -0.18 0.31 1.00
##
## $`45`
## x1 x2 x3
## x1 1.00 0.04 -0.12
## x2 0.04 1.00 0.46
## x3 -0.12 0.46 1.00
##
## $`46`
## x1 x2 x3
## x1 1.00 -0.22 -0.10
## x2 -0.22 1.00 0.49
## x3 -0.10 0.49 1.00
##
## $`47`
## x1 x2 x3
## x1 1.00 -0.19 -0.26
## x2 -0.19 1.00 0.38
## x3 -0.26 0.38 1.00
##
## $`48`
## x1 x2 x3
## x1 1.00 -0.45 -0.45
## x2 -0.45 1.00 0.74
## x3 -0.45 0.74 1.00
##
## $`49`
## x1 x2 x3
## x1 1.00 -0.19 -0.10
## x2 -0.19 1.00 0.13
## x3 -0.10 0.13 1.00
##
## $`50`
## x1 x2 x3
## x1 1.00 -0.18 -0.23
## x2 -0.18 1.00 0.39
## x3 -0.23 0.39 1.00
##
## $`51`
## x1 x2 x3
## x1 1.00 -0.11 -0.08
## x2 -0.11 1.00 0.31
## x3 -0.08 0.31 1.00
##
## $`52`
## x1 x2 x3
## x1 1.00 -0.13 -0.16
## x2 -0.13 1.00 0.43
## x3 -0.16 0.43 1.00
##
## $`53`
## x1 x2 x3
## x1 1.00 -0.14 -0.22
## x2 -0.14 1.00 0.38
## x3 -0.22 0.38 1.00
##
## $`54`
## x1 x2 x3
## x1 1.00 -0.14 -0.12
## x2 -0.14 1.00 0.50
## x3 -0.12 0.50 1.00
##
## $`55`
## x1 x2 x3
## x1 1.00 -0.22 -0.31
## x2 -0.22 1.00 0.49
## x3 -0.31 0.49 1.00
##
## $`56`
## x1 x2 x3
## x1 1.00 0.00 -0.21
## x2 0.00 1.00 0.27
## x3 -0.21 0.27 1.00
##
## $`57`
## x1 x2 x3
## x1 1.00 -0.10 -0.22
## x2 -0.10 1.00 0.35
## x3 -0.22 0.35 1.00
##
## $`58`
## x1 x2 x3
## x1 1.00 -0.17 -0.20
## x2 -0.17 1.00 0.45
## x3 -0.20 0.45 1.00
##
## $`59`
## x1 x2 x3
## x1 1.00 -0.28 -0.24
## x2 -0.28 1.00 0.59
## x3 -0.24 0.59 1.00
##
## $`60`
## x1 x2 x3
## x1 1.00 -0.03 -0.20
## x2 -0.03 1.00 0.35
## x3 -0.20 0.35 1.00
##
## $`61`
## x1 x2 x3
## x1 1.00 -0.18 -0.18
## x2 -0.18 1.00 0.34
## x3 -0.18 0.34 1.00
##
## $`62`
## x1 x2 x3
## x1 1.00 -0.30 -0.29
## x2 -0.30 1.00 0.45
## x3 -0.29 0.45 1.00
##
## $`63`
## x1 x2 x3
## x1 1.00 -0.23 -0.04
## x2 -0.23 1.00 0.08
## x3 -0.04 0.08 1.00
##
## $`64`
## x1 x2 x3
## x1 1.00 -0.03 0.15
## x2 -0.03 1.00 0.12
## x3 0.15 0.12 1.00
##
## $`65`
## x1 x2 x3
## x1 1.00 -0.10 -0.01
## x2 -0.10 1.00 0.06
## x3 -0.01 0.06 1.00
##
## $`66`
## x1 x2 x3
## x1 1.00 -0.27 -0.30
## x2 -0.27 1.00 -0.39
## x3 -0.30 -0.39 1.00
##
## $`67`
## x1 x2 x3
## x1 1.00 -0.17 -0.09
## x2 -0.17 1.00 0.26
## x3 -0.09 0.26 1.00
##
## $`68`
## x1 x2 x3
## x1 1.00 -0.14 -0.19
## x2 -0.14 1.00 0.27
## x3 -0.19 0.27 1.00
##
## $`69`
## x1 x2 x3
## x1 1.00 -0.35 -0.19
## x2 -0.35 1.00 0.40
## x3 -0.19 0.40 1.00
##
## $`70`
## x1 x2 x3
## x1 1.00 -0.20 -0.23
## x2 -0.20 1.00 0.23
## x3 -0.23 0.23 1.00
##
## $`71`
## x1 x2 x3
## x1 1.00 -0.03 -0.25
## x2 -0.03 1.00 0.02
## x3 -0.25 0.02 1.00
##
## $`72`
## x1 x2 x3
## x1 1.00 -0.26 -0.10
## x2 -0.26 1.00 0.22
## x3 -0.10 0.22 1.00
##
## $`73`
## x1 x2 x3
## x1 1.00 -0.16 -0.13
## x2 -0.16 1.00 0.45
## x3 -0.13 0.45 1.00
##
## $`74`
## x1 x2 x3
## x1 1.00 -0.29 -0.26
## x2 -0.29 1.00 0.45
## x3 -0.26 0.45 1.00
##
## $`75`
## x1 x2 x3
## x1 1.00 -0.17 -0.01
## x2 -0.17 1.00 0.12
## x3 -0.01 0.12 1.00
##
## $`76`
## x1 x2 x3
## x1 1.00 -0.36 -0.07
## x2 -0.36 1.00 0.17
## x3 -0.07 0.17 1.00
##
## $`77`
## x1 x2 x3
## x1 1.00 -0.19 -0.19
## x2 -0.19 1.00 0.31
## x3 -0.19 0.31 1.00
##
## $`78`
## x1 x2 x3
## x1 1.00 -0.15 -0.28
## x2 -0.15 1.00 0.29
## x3 -0.28 0.29 1.00
##
## $`79`
## x1 x2 x3
## x1 1.00 -0.18 -0.31
## x2 -0.18 1.00 0.46
## x3 -0.31 0.46 1.00
##
## $`80`
## x1 x2 x3
## x1 1.00 -0.10 -0.08
## x2 -0.10 1.00 0.16
## x3 -0.08 0.16 1.00
##
## $`81`
## x1 x2 x3
## x1 1.00 -0.12 -0.08
## x2 -0.12 1.00 0.10
## x3 -0.08 0.10 1.00
##
## $`82`
## x1 x2 x3
## x1 1.00 -0.18 -0.30
## x2 -0.18 1.00 0.19
## x3 -0.30 0.19 1.00
##
## $`83`
## x1 x2 x3
## x1 1.00 -0.06 -0.08
## x2 -0.06 1.00 -0.28
## x3 -0.08 -0.28 1.00
##
## $`84`
## x1 x2 x3
## x1 1.00 -0.23 -0.15
## x2 -0.23 1.00 0.33
## x3 -0.15 0.33 1.00
##
## $`85`
## x1 x2 x3
## x1 1.00 -0.15 -0.15
## x2 -0.15 1.00 0.20
## x3 -0.15 0.20 1.00
N <- c(239,130,298,1072,4834,214,346,618,754,1015,1015,377,
381,323,1173,1020,197,750,425,985,1675,1101,2472,615,
315,288,407,3449,2491,731,470,180,146,355,359,199,1725,
1409,1034,825,101,1364,763,296,296,1759,784,1054,867,826,
344,1244,3819,1273,802,689,2603,1200,1364,277,555,1037,448,
282,3780,600,2844,1359,917,249,297,2817,546,318,2844,3610,102,
322,603,1197,6504,89,1504,1724,1919)
D<- c(1,1,1,2,2,1,2,1,2,1,1,1,2,2,2,1,1,1,1,2,1,1,1,2,2,2,1,2,1,1,2,2,2,1,1,1,2,1,
2,2,2,2,1,2,2,1,1,2,1,1,2,1,1,1,1,1,1,1,1,1,2,1,2,1,1,1,1,2,1,1,1,2,1,1,1,2,2,
1,1,1,2,1,2,1,2)
D1<- my.vec[D>1]
N1<-N[D>1]
D2<- my.vec[D<2]
N2<-N[D<2]
Lag1 <- c(14.02,14.71,16.23,13.53,16.02,15.89,15.5,14,13,16.01,16,15,
13,13,12.3,12.23,14,13.22,21,14,13.79,15.5,12.386,13.9,13,13,
16,14,11,14,16.9,12.04,16.97,16.97,16,21.5,12.73,12.2,12.07,15,
15.67,14.03,11.93,14,14,13.79,14,16,10.5,16.6,19.8,16.1,18.2,20,16.4,
15.7,18.22,18.33,14,14,40.5,16,15.2,14,14,14,14,12.3,14.38,15.09,15.09,
14,16,23,13,14,31.73,19.72,14,15.04,9.5,15.58,17,17.2,13.5)
is.pd(my.vec)
## 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
## TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
## TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
## TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
## TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75
## TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## 76 77 78 79 80 81 82 83 84 85
## TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
random1 <- tssem1(my.vec, N, method="REM",
RE.type="Diag")
summary(random1)
##
## Call:
## meta(y = ES, v = acovR, RE.constraints = Diag(paste0(RE.startvalues,
## "*Tau2_", 1:no.es, "_", 1:no.es)), RE.lbound = RE.lbound,
## I2 = I2, model.name = model.name, suppressWarnings = TRUE,
## silent = silent, run = run)
##
## 95% confidence intervals: z statistic approximation
## Coefficients:
## Estimate Std.Error lbound ubound z value Pr(>|z|)
## Intercept1 -0.2078397 0.0216136 -0.2502017 -0.1654777 -9.6161 < 2.2e-16
## Intercept2 -0.1759372 0.0194153 -0.2139904 -0.1378840 -9.0618 < 2.2e-16
## Intercept3 0.3038670 0.0195947 0.2654621 0.3422718 15.5076 < 2.2e-16
## Tau2_1_1 0.0381741 NA NA NA NA NA
## Tau2_2_2 0.0305960 NA NA NA NA NA
## Tau2_3_3 0.0305888 0.0053596 0.0200842 0.0410933 5.7073 1.148e-08
##
## Intercept1 ***
## Intercept2 ***
## Intercept3 ***
## Tau2_1_1
## Tau2_2_2
## Tau2_3_3 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Q statistic on the homogeneity of effect sizes: 5030.813
## Degrees of freedom of the Q statistic: 252
## P value of the Q statistic: 0
##
## Heterogeneity indices (based on the estimated Tau2):
## Estimate
## Intercept1: I2 (Q statistic) 0.9786
## Intercept2: I2 (Q statistic) 0.9725
## Intercept3: I2 (Q statistic) 0.9750
##
## Number of studies (or clusters): 85
## Number of observed statistics: 255
## Number of estimated parameters: 6
## Degrees of freedom: 249
## -2 log likelihood: -242.4222
## OpenMx status1: 6 ("0" or "1": The optimization is considered fine.
## Other values may indicate problems.)
## Warning in print.summary.meta(x): OpenMx status1 is neither 0 or 1. You are advised to 'rerun' it again.