Running examples in ‘metaSEM-Ex.R’ failed
The error most likely occurred in:

> ### Name: Berkey98
> ### Title: Five Published Trails from Berkey et al. (1998)
> ### Aliases: Berkey98
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(Berkey98)
> 
> #### ML estimation method
> ## Multivariate meta-analysis
> x <- meta(y=cbind(PD, AL), v=cbind(var_PD, cov_PD_AL, var_AL), data=Berkey98)
> x <- rerun(x)
Running Meta analysis with ML with 5 parameters

Beginning initial fit attempt
Running Meta analysis with ML with 5 parameters

 Lowest minimum so far:  -11.6813137695501
 OpenMx status code 6 not in list of acceptable status codes, (0,1)
 Not all eigenvalues of the Hessian are positive: 25565313.331146, 125323.603271849, 18950.5806256828, 1916.14742049567, -66914467.0311619

Beginning fit attempt 1 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,1)
 Not all eigenvalues of the Hessian are positive: 391490.284254166, 125114.536157842, 19101.0227339667, 1946.52078064902, -143415.288753153

Beginning fit attempt 2 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters

 Fit attempt generated errors

Beginning fit attempt 3 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,1)
 Not all eigenvalues of the Hessian are positive: 46146975.0028634, 125323.225688536, 18951.0237237051, 1916.22500711413, -120797958.094035

Beginning fit attempt 4 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters

 Fit attempt generated errors

Beginning fit attempt 5 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters

 Fit attempt generated errors

Beginning fit attempt 6 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters

 Fit attempt generated errors

Beginning fit attempt 7 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters

 Fit attempt generated errors

Beginning fit attempt 8 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters

 Fit attempt generated errors

Beginning fit attempt 9 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters

 Fit attempt generated errors

Beginning fit attempt 10 of at maximum 10 extra tries
Running Meta analysis with ML with 5 parameters

 Fit attempt generated errors

Retry limit reached



 Retry limit reached; Best fit=-11.681314 (started at -11.681314)  (11 attempt(s): 3 valid, 8 errors)


 Uncertain solution found - consider parameter validity, try again, increase extraTries, change inits, change model, or check data!

> summary(x)

Call:
meta(y = cbind(PD, AL), v = cbind(var_PD, cov_PD_AL, var_AL), 
    data = Berkey98)

95% confidence intervals: z statistic approximation (robust=FALSE)
Coefficients:
             Estimate  Std.Error     lbound     ubound  z value Pr(>|z|)    
Intercept1  0.3448392         NA         NA         NA       NA       NA    
Intercept2 -0.3379381  0.0273638 -0.3915702 -0.2843060 -12.3498   <2e-16 ***
Tau2_1_1    0.0070020  0.0085440 -0.0097439  0.0237479   0.8195   0.4125    
Tau2_2_1    0.0094607  0.0084524 -0.0071057  0.0260270   1.1193   0.2630    
Tau2_2_2    0.0261445  0.0164950 -0.0061850  0.0584741   1.5850   0.1130    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Q statistic on the homogeneity of effect sizes: 128.2267
Degrees of freedom of the Q statistic: 8
P value of the Q statistic: 0

Heterogeneity indices (based on the estimated Tau2):
                             Estimate
Intercept1: I2 (Q statistic)   0.6021
Intercept2: I2 (Q statistic)   0.9250

Number of studies (or clusters): 5
Number of observed statistics: 10
Number of estimated parameters: 5
Degrees of freedom: 5
-2 log likelihood: -11.68131 
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.

> plot(x)
Warning in .solve(x = object$mx.fit@output$calculatedHessian, parameters = my.name) :
  Error in solving the Hessian matrix. Generalized inverse is used. The standard errors may not be trustworthy.

Warning in sqrt(c(x[xind, xind], x[yind, yind])) : NaNs produced
Error in if (scale[1] > 0) r <- r/scale[1] : 
  missing value where TRUE/FALSE needed
Calls: plot -> plot.meta -> points -> ellipse -> ellipse.default
Execution halted