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

> ### Name: power.ACE.test
> ### Title: Test the power of an ACE model to detect paths of interest.
> ### Aliases: power.ACE.test
> 
> ### ** Examples
> 
> 
> # ===============================================
> # = Power to detect a^2=.5 with equal MZ and DZ =
> # ===============================================
> power.ACE.test(AA = .5, CC = 0, update = "a") 
This is beta code: I likely will alter the interface!
Warning: In model 'drop_a' Optimizer returned a non-zero status code 6. The model does not satisfy the first-order optimality conditions to the required accuracy, and no improved point for the merit function could be found during the final linesearch (Mx status RED)
For 80% power, you need 166 MZ and 166 DZ pairs
   method = ncp
        n = 333
    power = 0.8
sig.level = 0.05
statistic = LRT
> # Suggests n = 84 MZ and 94 DZ pairs.
> 
> # ================================
> # = Show power across range of N =
> # ================================
> power.ACE.test(AA= .5, CC= 0, update = "a", search = TRUE)
This is beta code: I likely will alter the interface!
Warning: In model 'drop_a' Optimizer returned a non-zero status code 6. The model does not satisfy the first-order optimality conditions to the required accuracy, and no improved point for the merit function could be found during the final linesearch (Mx status RED)
     N     power lower upper
1   68 0.2527548    NA    NA
2   97 0.3387609    NA    NA
3  127 0.4232452    NA    NA
4  156 0.4990226    NA    NA
5  186 0.5704155    NA    NA
6  215 0.6323162    NA    NA
7  245 0.6890179    NA    NA
8  274 0.7370188    NA    NA
9  304 0.7800923    NA    NA
10 333 0.8159042    NA    NA
11 363 0.8475326    NA    NA
12 392 0.8734572    NA    NA
13 421 0.8953717    NA    NA
14 451 0.9143768    NA    NA
15 480 0.9296981    NA    NA
16 510 0.9428569    NA    NA
17 539 0.9533703    NA    NA
18 569 0.9623254    NA    NA
19 598 0.9694253    NA    NA
20 628 0.9754296    NA    NA
> 
> # Salutary note: You need well fitting models with correct betas in the data
> # for power to be valid.
> # tryHard helps ensure this, as does the default nSim= 4000 pair data.
> # Power is important to get right, so I recommend using tryHard = "yes"
> power.ACE.test(AA= .5, CC= 0, update = "a", tryHard= "yes")
This is beta code: I likely will alter the interface!
Running ACE with 4 parameters

Beginning initial fit attempt
Running ACE with 4 parameters

 Lowest minimum so far:  21994.5743438305
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 8211695686692.22, 14072.55996328, 6628.72814397701, -24635086988923.7

Beginning fit attempt 1 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22092.8159384325 vs 21994.5743438305

Beginning fit attempt 2 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22092.8159384325 vs 21994.5743438305

Beginning fit attempt 3 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22092.8159384326 vs 21994.5743438305

Beginning fit attempt 4 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22092.8159384326 vs 21994.5743438305

Beginning fit attempt 5 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22092.8159384325 vs 21994.5743438305

Beginning fit attempt 6 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22699.0155309381 vs 21994.5743438305

Beginning fit attempt 7 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22129.9139739754 vs 21994.5743438305

Beginning fit attempt 8 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22092.8159384325 vs 21994.5743438305

Beginning fit attempt 9 of at maximum 10 extra tries
Running ACE with 4 parameters

 Lowest minimum so far:  21994.5743437627
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 8081555359387.54, 14094.8001975453, 6678.40654827431, -24244666007189.2

Beginning fit attempt 10 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  22092.8159384325 vs 21994.5743437627

Retry limit reached



 Retry limit reached; Best fit=21994.574 (started at 22414.561)  (11 attempt(s): 11 valid, 0 errors)


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

 Start values from best fit:
-4.90431716631184e-08,0.706929964374927,8.3089179480229e-05,0.706929936521183
Running drop_a with 3 parameters

Beginning initial fit attempt
Running drop_a with 3 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 6874518275547.46, 51992.2145600415, -25656051267170

Beginning fit attempt 1 of at maximum 10 extra tries
Running drop_a with 3 parameters

 Lowest minimum so far:  22092.8159384326
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 2912654592054.57, 13604.4434944163, -10870174809439.1

Beginning fit attempt 2 of at maximum 10 extra tries
Running drop_a with 3 parameters

 Lowest minimum so far:  22092.8159384325
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 5942534124694.73, 13604.4856648443, -22177839166065.1

Beginning fit attempt 3 of at maximum 10 extra tries
Running drop_a with 3 parameters

 Fit attempt worse than current best:  22699.0155309413 vs 22092.8159384325

Beginning fit attempt 4 of at maximum 10 extra tries
Running drop_a with 3 parameters

 Fit attempt worse than current best:  26637.4179904811 vs 22092.8159384325

Beginning fit attempt 5 of at maximum 10 extra tries
Running drop_a with 3 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 2863927717993.05, 13604.5558204752, -10688323639747.5

Beginning fit attempt 6 of at maximum 10 extra tries
Running drop_a with 3 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 2863787913230.55, 13604.4926378934, -10687801881270.6

Beginning fit attempt 7 of at maximum 10 extra tries
Running drop_a with 3 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3033520080685.57, 13604.4938899784, -11321250953891.7

Beginning fit attempt 8 of at maximum 10 extra tries
Running drop_a with 3 parameters

 Fit attempt worse than current best:  44209.9308566266 vs 22092.8159384325

Beginning fit attempt 9 of at maximum 10 extra tries
Running drop_a with 3 parameters

Solution found




 Solution found!  Final fit=22092.816 (started at 25153.839)  (10 attempt(s): 10 valid, 0 errors)

 Start values from best fit:
-5.0229286891743e-08,0.612219300590885,0.790371703463739
For 80% power, you need 166 MZ and 166 DZ pairs
   method = ncp
        n = 333
    power = 0.8
sig.level = 0.05
statistic = LRT
> 
> # =====================
> # = Power to detect C =
> # =====================
> 
> # 102 of each of MZ and DZ pairs for 80% power.
> power.ACE.test(AA= .5, CC= .3, update = "c", tryHard= "yes")
This is beta code: I likely will alter the interface!
Running ACE with 4 parameters

Beginning initial fit attempt
Running ACE with 4 parameters

 Lowest minimum so far:  19935.2075053042
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 7443016550817.3, 30687.9970771025, 2475.11328581816, -22329049546750.6

Beginning fit attempt 1 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  20266.5478054396 vs 19935.2075053042

Beginning fit attempt 2 of at maximum 10 extra tries
Running ACE with 4 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 7380194001473.96, 30687.9582222784, 2475.07635142732, -22140581898720.5

Beginning fit attempt 3 of at maximum 10 extra tries
Running ACE with 4 parameters

 Lowest minimum so far:  19935.2075053042
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 7324881987755.5, 30688.037349933, 2475.12278253207, -21974645857565.2

Beginning fit attempt 4 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  20266.5478054399 vs 19935.2075053042

Beginning fit attempt 5 of at maximum 10 extra tries
Running ACE with 4 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 354403705533.19, 30687.9951559667, 2475.16991383769, -118134533277.287

Beginning fit attempt 6 of at maximum 10 extra tries
Running ACE with 4 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3466641998338.61, 30687.9015681464, 2475.08243838464, -10399925889314.5

Beginning fit attempt 7 of at maximum 10 extra tries
Running ACE with 4 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3411329984616.64, 30688.0067038529, 2475.10309361035, -10233989848148.6

Beginning fit attempt 8 of at maximum 10 extra tries
Running ACE with 4 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 353278943389.318, 30688.0677513596, 2475.17769867212, -117759612562.677

Beginning fit attempt 9 of at maximum 10 extra tries
Running ACE with 4 parameters

Solution found




 Solution found!  Final fit=19935.208 (started at 20174.841)  (10 attempt(s): 10 valid, 0 errors)

 Start values from best fit:
-4.52288982667131e-08,0.706929981581117,0.547585574730543,0.447101727454412
Running drop_c with 3 parameters

Beginning initial fit attempt
Running drop_c with 3 parameters

 Lowest minimum so far:  20015.8223706487
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 5383852040912.34, 27764.8393275464, -20092809203088.9

Beginning fit attempt 1 of at maximum 10 extra tries
Running drop_c with 3 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 5429917984717.38, 27764.8201523739, -20264729645868

Beginning fit attempt 2 of at maximum 10 extra tries
Running drop_c with 3 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 5430158375509.94, 27764.7669768666, -20265626796519.5

Beginning fit attempt 3 of at maximum 10 extra tries
Running drop_c with 3 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 5470544028986.41, 27764.7528097517, -20416348107190.5

Beginning fit attempt 4 of at maximum 10 extra tries
Running drop_c with 3 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 149241005976.932, 27764.8537689327, -39988965757.0855

Beginning fit attempt 5 of at maximum 10 extra tries
Running drop_c with 3 parameters

Solution found




 Solution found!  Final fit=20015.822 (started at 20408.339)  (6 attempt(s): 6 valid, 0 errors)

 Start values from best fit:
-4.51399141700716e-08,0.881661999578546,0.436315077515404
For 80% power, you need 203 MZ and 203 DZ pairs
   method = ncp
        n = 406
    power = 0.8
sig.level = 0.05
statistic = LRT
> 
> # ========================================
> # = Set 'a' to a fixed, but non-zero value =
> # ========================================
> 
> power.ACE.test(update= "a", value= sqrt(.2), AA= .5, CC= 0)
This is beta code: I likely will alter the interface!
Warning: In model 'ACE' Optimizer returned a non-zero status code 5. The Hessian at the solution does not appear to be convex. See ?mxCheckIdentification for possible diagnosis (Mx status RED).
Error in mxPowerSearch(trueModel, falseModel, probes = probes, gdFun = gdFun,  : 
  falseModel fit better than trueModel?
Calls: power.ACE.test -> mxPower -> mxPowerSearch
Execution halted