Compute goodness of fit as described in AIC. The number of parameters used correspond to the number of variables free to vary in the general model., u, method = c("AIC", "BIC"), k = 2)



A list of parameters as defined in rtheta. For t this function, it will usually be the output of fit.full.GMCM.


An n by d matrix of marginally uniform observations. Rows corresponds to observations and columns to the dimensions of the variables. I.e. these are often ranked and scaled test statistics or other observations.


A character of length 1 which specifies the goodness of fit to compute. Default is "AIC". "BIC" is also a option.


A integer specifying the default used constant "k" in AIC. See AIC.


A single number giving the goodness of fit as requested.


set.seed(2) data(u133VsExon) u <- Uhat(u133VsExon[sample(19577, 500), ]) # Subset for faster fitting theta1 <- fit.full.GMCM(u, m = 2, method = "L-BFGS")
#> iter 10 value -110.257470 #> iter 20 value -110.738198 #> final value -110.874731 #> converged, u) # AIC
#> [1] -205.7495, u, method = "BIC")
#> [1] -172.0326
if (FALSE) { theta2 <- fit.full.GMCM(u, m = 3, method = "L-BFGS"), u), u, method = "BIC") }