tbates / umx

Making Structural Equation Modeling (SEM) in R quick & powerful
https://tbates.github.io/
44 stars 17 forks source link

get power.ACE.test() to print the value of the parameter being dropped #111

Closed tbates closed 4 years ago

tbates commented 4 years ago

Be nice to print the original value of the parameter being tested in ACE power test. (current output)

power.ACE.test(update= "a", value= sqrt(.2), AA= .5, CC= 0)

For 80% power you need 463 MZ and 463 DZ pairs
   method = ncp
        n = 926
    power = 0.8
sig.level = 0.05
statistic = LRT
tbates commented 4 years ago

now says:

power.ACE.test(update= "a", value= sqrt(.2), AA= .5, CC= 0, try="no")

For 80% power to detect 'a_r1c1' of size 0.5, you need 463 MZ and 463 DZ pairs
   method = ncp
        n = 926
    power = 0.8
sig.level = 0.05
statistic = LRT