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!
Running ACE with 4 parameters

Beginning initial fit attempt
Running ACE with 4 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: 8184596658079.31, 17189.0834202889, 5900.46398143576, -24553789907417.1

Beginning fit attempt 1 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: 8184596658111.92, 17189.1014608634, 5900.48044220316, -24553789907514.9

Beginning fit attempt 2 of at maximum 10 extra tries
Running ACE with 4 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: 3845246209940.38, 17189.0926590304, 5900.57813446911, -11535738563000.2

Beginning fit attempt 3 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: 8117843533230.31, 17189.1297998051, 5900.50923171273, -24353530532870.2

Beginning fit attempt 4 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: 8117843533228.87, 17189.1295814238, 5900.53353039109, -24353530532865.8

Beginning fit attempt 5 of at maximum 10 extra tries
Running ACE with 4 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: 8117843533235.14, 17189.0436551939, 5900.46701615534, -24353530532884.6

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: 8117644863216.61, 17189.0921289886, 5900.49956483185, -24352934522829

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: 8117843533235.54, 17189.1062075198, 5900.59182422924, -24353530532885.9

Beginning fit attempt 8 of at maximum 10 extra tries
Running ACE with 4 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: 8117843533233.47, 17189.1079907032, 5900.52280473856, -24353530532879.6

Beginning fit attempt 9 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: 8117843533232.04, 17189.1245123274, 5900.48796819019, -24353530532875.3

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

 Fit attempt worse than current best:  22567.6187804623 vs 22092.8159384325

Retry limit reached



 Retry limit reached; Best fit=22092.816 (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:
-2.18706521430023e-08,2.66348091472849e-07,0.612219287908221,0.790371746276854
Running drop_a with 3 parameters

Beginning initial fit attempt
Running drop_a with 3 parameters

 Fit attempt worse than current best:  22699.0155309412 vs 22092.8159384325

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

 Fit attempt worse than current best:  35457.0190220442 vs 22092.8159384325

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

 Fit attempt worse than current best:  36856.1954822885 vs 22092.8159384325

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

 Fit attempt worse than current best:  36569.9370457271 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:  22699.0155309414 vs 22092.8159384325

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

 Fit attempt worse than current best:  40491.1477378527 vs 22092.8159384325

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

 Fit attempt worse than current best:  37733.789401899 vs 22092.8159384325

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

 Fit attempt worse than current best:  33771.2043221764 vs 22092.8159384325

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

 Fit attempt generated errors

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

 Fit attempt worse than current best:  22699.0155309358 vs 22092.8159384325

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

 Fit attempt worse than current best:  35500.5244890284 vs 22092.8159384325

Retry limit reached



 All fit attempts resulted in errors - check starting values or model specification

For 80% power, you need 2 MZ and 2 DZ pairs
   method = ncp
        n = 3
    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!
Running ACE with 4 parameters

Beginning initial fit attempt
Running ACE with 4 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: 3909810819396.12, 17189.1500951842, 5900.61154327801, -11729432391367.5

Beginning fit attempt 1 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: 8117644863215.48, 17189.2505989158, 5900.58959644677, -24352934522825.7

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: 8117843533229.72, 17189.1074900007, 5900.51621589512, -24353530532868.5

Beginning fit attempt 3 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: 8117843533228.59, 17189.089609189, 5900.5164603159, -24353530532864.9

Beginning fit attempt 4 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: 8184596658037.09, 17189.0738862581, 5900.55262343337, -24553789907290.5

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

 Fit attempt worse than current best:  22608.9216113309 vs 22092.8159384325

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

 Fit attempt worse than current best:  22638.9796193259 vs 22092.8159384325

Beginning fit attempt 7 of at maximum 10 extra tries
Running ACE with 4 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: 3845643562568.02, 17189.2159987051, 5900.66501325504, -11536930620883.3

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: 4144446093539.65, 17189.1375743847, 5900.56134181345, -12433338213798.2

Beginning fit attempt 9 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: 8184596658108.76, 17189.0371451878, 5900.47926283943, -24553789907505.3

Beginning fit attempt 10 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: 8117644863215.72, 17189.0421627775, 5900.51913914135, -24352934522826.3

Retry limit reached



 Retry limit reached; Best fit=22092.816 (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:
-2.31173443581308e-08,2.93932089271245e-10,0.612219284149609,0.790371760529176
Running drop_a with 3 parameters

Beginning initial fit attempt
Running drop_a with 3 parameters

 Fit attempt worse than current best:  22699.0155309413 vs 22092.8159384325

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

 Fit attempt worse than current best:  36461.0059684024 vs 22092.8159384325

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

 Fit attempt worse than current best:  38503.6516393036 vs 22092.8159384325

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

 Fit attempt worse than current best:  42451.603812271 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:  22699.0155309412 vs 22092.8159384325

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

 Fit attempt worse than current best:  34542.2145377764 vs 22092.8159384325

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

 Fit attempt worse than current best:  22699.0155273146 vs 22092.8159384325

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

 Fit attempt worse than current best:  33946.2456274839 vs 22092.8159384325

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

 Fit attempt worse than current best:  22699.0155307462 vs 22092.8159384325

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

 Fit attempt worse than current best:  22699.0155309414 vs 22092.8159384325

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

 Fit attempt worse than current best:  37558.3374983428 vs 22092.8159384325

Retry limit reached



 All fit attempts resulted in errors - check starting values or model specification

   N     power lower upper
1  0 0.0500000    NA    NA
2  1 0.5025754    NA    NA
3  1 0.5025754    NA    NA
4  1 0.5025754    NA    NA
5  1 0.5025754    NA    NA
6  1 0.5025754    NA    NA
7  2 0.7941267    NA    NA
8  2 0.7941267    NA    NA
9  2 0.7941267    NA    NA
10 2 0.7941267    NA    NA
11 2 0.7941267    NA    NA
12 2 0.7941267    NA    NA
13 3 0.9258816    NA    NA
14 3 0.9258816    NA    NA
15 3 0.9258816    NA    NA
16 3 0.9258816    NA    NA
17 3 0.9258816    NA    NA
18 4 0.9757328    NA    NA
19 4 0.9757328    NA    NA
20 4 0.9757328    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" (the default)
> power.ACE.test(AA= .5, CC= 0, update = "a")
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:  22092.8159384325
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 8184016040015.48, 17189.06857331, 5900.51289833048, -24552048053225.7

Beginning fit attempt 1 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: 8117843533230.94, 17189.1320127365, 5900.55857225606, -24353530532871.9

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: 8181422485143.22, 17189.1424329021, 5900.56302930492, -24544267388608.8

Beginning fit attempt 3 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: 8181613459313.77, 17189.1838405882, 5900.57833593374, -24544840311120.6

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

 Lowest minimum so far:  22092.8159384324
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 8184008344061.19, 17189.1317672277, 5900.51894261381, -24552024965362.6

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

 Lowest minimum so far:  22092.8159384324
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 8182010799334.65, 17189.1275921421, 5900.47813739044, -24546032331183.1

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

 Lowest minimum so far:  22092.8159384321
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 8117064245193.74, 17189.1530883146, 5900.4742763464, -24351192668760.4

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: 8248763924060.82, 17189.188339479, 5900.60789404757, -24746291705361.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: 8117064245188.35, 17189.1645306452, 5900.60997860241, -24351192668744.3

Beginning fit attempt 9 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: 4077494281687.22, 17189.1675170386, 5900.56282387897, -12232482778240.9

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

 Fit attempt generated errors

Retry limit reached



 Retry limit reached; Best fit=22092.816 (started at 22414.561)  (11 attempt(s): 10 valid, 1 errors)


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

 Start values from best fit:
1.77249876686774e-09,9.8086045619646e-07,0.612219332008635,0.790371752134642
Running drop_a with 3 parameters

Beginning initial fit attempt
Running drop_a with 3 parameters

 Fit attempt worse than current best:  22699.0155309412 vs 22092.8159384325

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

 Fit attempt worse than current best:  35954.1636788864 vs 22092.8159384325

Beginning fit attempt 2 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 3 of at maximum 10 extra tries
Running drop_a with 3 parameters

 Fit attempt worse than current best:  33975.5661309145 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:  22699.0155309413 vs 22092.8159384325

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

 Fit attempt worse than current best:  22699.0155309414 vs 22092.8159384325

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

 Fit attempt worse than current best:  41383.396033901 vs 22092.8159384325

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

 Fit attempt worse than current best:  40422.1750369855 vs 22092.8159384325

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

 Fit attempt generated errors

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

 Fit attempt worse than current best:  22699.0155309334 vs 22092.8159384325

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

 Fit attempt worse than current best:  32769.6176501438 vs 22092.8159384325

Retry limit reached



 All fit attempts resulted in errors - check starting values or model specification

For 80% power, you need 2 MZ and 2 DZ pairs
   method = ncp
        n = 3
    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")
This is beta code: I likely will alter the interface!
Running ACE with 4 parameters

Beginning initial fit attempt
Running ACE with 4 parameters

 Fit attempt worse than current best:  20266.5478053732 vs 20174.8414520721

Beginning fit attempt 1 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: 7324342936431.27, 30688.0831126268, 2475.14345035572, -21973028703592.5

Beginning fit attempt 2 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 3 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 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: 7324858657999.59, 30688.076508023, 2475.10876661832, -21974575868297.5

Beginning fit attempt 6 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: 7324858657999.61, 30688.0639070513, 2475.1764565557, -21974575868297.6

Beginning fit attempt 7 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 8 of at maximum 10 extra tries
Running ACE with 4 parameters

 Fit attempt worse than current best:  20266.5478054398 vs 19935.2075053042

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

 Fit attempt worse than current best:  20266.5478053815 vs 19935.2075053042

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

 Fit attempt generated errors

Retry limit reached



 Retry limit reached; Best fit=19935.208 (started at 20174.841)  (11 attempt(s): 10 valid, 1 errors)


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

 Start values from best fit:
-5.17737093618877e-08,0.706930021125827,0.547585513665838,0.44710173602757
Running drop_c with 3 parameters

Beginning initial fit attempt
Running drop_c with 3 parameters

 Fit attempt generated errors

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

 Fit attempt worse than current best:  27232.6689200668 vs 20408.3386140284

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

 Fit attempt worse than current best:  35695.8429454294 vs 20408.3386140284

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

 Fit attempt worse than current best:  22699.0155309414 vs 20408.3386140284

Beginning fit attempt 4 of at maximum 10 extra tries
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: 5430296603419.19, 27764.7864286332, -20266142670099.8

Beginning fit attempt 5 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: 5424117498984.74, 27764.8833777839, -20243081938405.2

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

 Fit attempt worse than current best:  31251.9168970895 vs 20015.8223706487

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

 Fit attempt worse than current best:  26489.9098639422 vs 20015.8223706487

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

 Fit attempt worse than current best:  22699.0155305454 vs 20015.8223706487

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

 Fit attempt worse than current best:  22699.0155243958 vs 20015.8223706487

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

 Fit attempt worse than current best:  30164.6509641401 vs 20015.8223706487

Retry limit reached



 Retry limit reached; Best fit=20015.822 (started at 20408.339)  (11 attempt(s): 10 valid, 1 errors)


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

 Start values from best fit:
4.7555041121107e-09,0.881661981846575,0.436315101622322
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!
Running ACE with 4 parameters

Beginning initial fit attempt
Running ACE with 4 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: 8117644863216.83, 17189.1893212518, 5900.64522303143, -24352934522829.7

Beginning fit attempt 1 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: 8184397988087.44, 17189.1454039541, 5900.51131229613, -24553193897441.5

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: 8117255219193.28, 17189.1366136606, 5900.56324362081, -24351765590759

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

 Fit attempt generated errors

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

 Fit attempt generated errors

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: 8184405684107.23, 17189.0584396035, 5900.54949770303, -24553216985500.9

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

 Fit attempt worse than current best:  22668.5868015314 vs 22092.8159384325

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

 Fit attempt worse than current best:  22234.1774255541 vs 22092.8159384325

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: 8117843533228.59, 17189.1509243982, 5900.56547826697, -24353530532865

Beginning fit attempt 9 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: 8117644863215.32, 17189.1167936064, 5900.49677714936, -24352934522825.2

Beginning fit attempt 10 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: 3845643543333.34, 17189.0903647111, 5900.54319539537, -11536930563179.1

Retry limit reached



 Retry limit reached; Best fit=22092.816 (started at 22414.561)  (11 attempt(s): 9 valid, 2 errors)


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

 Start values from best fit:
-3.48894908706943e-08,1.08186242109658e-07,0.612219366263576,0.79037170826496
Running drop_a with 3 parameters

Beginning initial fit attempt
Running drop_a with 3 parameters

 Fit attempt worse than current best:  22258.8947184954 vs 22146.7875197418

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

 Fit attempt worse than current best:  34028.6469404217 vs 22146.7875197418

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

 Fit attempt worse than current best:  22258.8947184967 vs 22146.7875197418

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

 Fit attempt worse than current best:  34729.3423069054 vs 22146.7875197418

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

 Fit attempt worse than current best:  22258.8947184959 vs 22146.7875197418

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

 Fit attempt generated errors

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

 Fit attempt worse than current best:  35929.8391152505 vs 22146.7875197418

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

 Fit attempt worse than current best:  36471.1322089844 vs 22146.7875197418

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

 Fit attempt generated errors

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

 Fit attempt worse than current best:  22258.8947184945 vs 22146.7875197418

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

 Fit attempt worse than current best:  22258.8947184967 vs 22146.7875197418

Retry limit reached



 All fit attempts resulted in errors - check starting values or model specification

For 80% power, you need 98 MZ and 98 DZ pairs
   method = ncp
        n = 197
    power = 0.8
sig.level = 0.05
statistic = LRT
> # TODO get power.ACE.test to print the value of A in the null model.
> 
> # ========================================
> # = Drop More than one parameter (A & C) =
> # ========================================
> # E vs AE: the hypothesis that twins show no familial similarity.
> power.ACE.test(update = "a_after_dropping_c", AA= .5, CC= .3)
This is beta code: I likely will alter the interface!
Running ACE with 3 parameters

Beginning initial fit attempt
Running ACE with 3 parameters

 Fit attempt worse than current best:  28892.8035183118 vs 21070.5635414025

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

 Fit attempt worse than current best:  25283.0407039872 vs 21070.5635414025

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

 Fit attempt worse than current best:  24678.59500601 vs 21070.5635414025

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

 Fit attempt worse than current best:  22699.0155309414 vs 21070.5635414025

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

 Fit attempt worse than current best:  22699.0155309297 vs 21070.5635414025

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

 Fit attempt worse than current best:  34843.0704249433 vs 21070.5635414025

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

 Fit attempt worse than current best:  22699.0155309409 vs 21070.5635414025

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

 Fit attempt worse than current best:  33769.5571014792 vs 21070.5635414025

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

 Fit attempt worse than current best:  22699.0155309353 vs 21070.5635414025

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

 Fit attempt worse than current best:  22699.015530941 vs 21070.5635414025

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

 Fit attempt worse than current best:  35817.8430731261 vs 21070.5635414025

Retry limit reached



 All fit attempts resulted in errors - check starting values or model specification

Running drop_a_after_dropping_c with 2 parameters

Beginning initial fit attempt
Running drop_a_after_dropping_c with 2 parameters

 Lowest minimum so far:  22699.0155309414
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3454618508307.89, -20134992028180.2

Beginning fit attempt 1 of at maximum 10 extra tries
Running drop_a_after_dropping_c with 2 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3454618508308.14, -20134992028181.7

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

 Lowest minimum so far:  22699.0155309414
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3510335486642.28, -20459734376013.1

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

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

Beginning fit attempt 5 of at maximum 10 extra tries
Running drop_a_after_dropping_c with 2 parameters

 Lowest minimum so far:  22699.0155309413
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3454701420481.47, -20135475275741.7

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

Beginning fit attempt 7 of at maximum 10 extra tries
Running drop_a_after_dropping_c with 2 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3482394085305.88, -20296879954560.8

Beginning fit attempt 8 of at maximum 10 extra tries
Running drop_a_after_dropping_c with 2 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3482394085304.88, -20296879954555.3

Beginning fit attempt 9 of at maximum 10 extra tries
Running drop_a_after_dropping_c with 2 parameters
 OpenMx status code 6 not in list of acceptable status codes, (0,0)
 Not all eigenvalues of the Hessian are positive: 3454618508311.1, -20134992028198.8

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

Retry limit reached



 Retry limit reached; Best fit=22699.016 (started at 35817.843)  (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:
-3.38122423634325e-10,0.99974994727963
Error in mxPowerSearch(trueModel, falseModel, probes = probes, gdFun = gdFun,  : 
  falseModel fit better than trueModel?
Calls: power.ACE.test -> mxPower -> mxPowerSearch
Execution halted