Adjusted BIC

4 replies [Last post]
trinewa's picture
Offline
Joined: 11/30/2009
Ryne's picture
Offline
Joined: 07/31/2009
It'll be released in the next

It'll be released in the next binary, along with other added fit indices in 'summary'. For now, you can calculate it by hand with the following formula or function:

sBIC (sample-size adjusted BIC) = -2LL + log((n+2)/24)

or

trinewa <- function(model){
  neg2LL <- model@output$Minus2LogLikelihood
  adjN <- (model@data@numObs+2)/24
  return(neg2LL + log(adjN))
}

Sabha's picture
Offline
Joined: 05/18/2010
Along the same lines, is

Along the same lines, is there any option available at the moment to compute pval and chi square values for ACE model fit. I have pasted an example of the summary output as an example.

> univACESummX

free parameters:
name matrix row col Estimate Std.Error
1 a11 ACE.a 1 1 -4.302609e-01 0.1389135
2 c11 ACE.c 1 1 -3.342676e-08 0.9382380
3 e11 ACE.e 1 1 5.466756e-01 0.1022683
4 mean ACE.Mean 1 1 3.361335e-08 0.1264393

confidence intervals:
lbound estimate ubound
univACE.StdA[1,1] 3.815964e-15 3.825054e-01 0.7105036
univACE.StdC[1,1] 9.630561e-17 2.308673e-15 0.6585438
univACE.StdE[1,1] 2.894965e-01 6.174946e-01 1.0000000

observed statistics: 38
estimated parameters: 4
degrees of freedom: 34
-2 log likelihood: 78.70771
saturated -2 log likelihood: NA
number of observations: 19
chi-square: NA
p: NA
AIC (Mx): 10.70771
BIC (Mx): -10.70161
adjusted BIC:
RMSEA: NA
timestamp: 2011-02-24 14:53:16
frontend time: 0.608 secs
backend time: 0.03200006 secs
independent submodels time: 0 secs
wall clock time: 0.6400001 secs
cpu time: 0.6400001 secs
openmx version number: 1.0.4-1540

Ryne's picture
Offline
Joined: 07/31/2009
Use the mxCompare function to

Use the mxCompare function to compare nested models by the likelihood ratio test.

trinewa's picture
Offline
Joined: 11/30/2009
Great! Thank you!

Great! Thank you!