Description Usage Arguments Details Value Author(s) See Also Examples
This function makes it easy to compare different models based on different design matrices.
1 | model_comp(object, ...)
|
object |
An object of class |
... |
Insert different design matrices here (Input of at least one design-matrix is required!). Use the help function |
The main advantage of this function is first, that different models with different design-matrices are fitted in one step and second that reasonable starting values are used, based on the values of the submitted model object. This should accelerate the model fitting process (by reducing the number of iterations until convergence).
The loop which fits the models always try()
s to fit the model, which means that an output is returned even if one or more models fail to converge.
The resulting object can be used to apply likelihood-ratio tests by using the lrt()
function.
A list with 2 elements
object |
the submitted object - which is the reference object |
modres |
a list of fitted objects, which used the same data like the reference object but different designs |
Manuel Reif
designTemp
lrt
reshMG
nrm
nelm
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 | ## Not run:
# Don't run this example because it takes more than 5sec.
NUMBI <- 5
ParList <- lapply(1:NUMBI,function(x)
{
Item1 <- c(c(-2,-1,1,2),c(-1.2,0.3,0.2,0.7))
names(Item1) <- c(paste("zeta",1:4,sep=""),paste("lamb",1:4,sep=""))
Item1
})
names(ParList) <- paste("item",1:NUMBI,sep="")
set.seed(111)
perp1 <- rnorm(5000,0,1)
perp2 <- rnorm(5000,1,1)
simdat1 <- NRM.sim(ParList,perp1)
simdat2 <- NRM.sim(ParList,perp2)
simdat1 <- data.frame(ID=1:5000,simdat1)
simdat2 <- data.frame(ID=5001:10000,simdat2)
simdatalla <- merge(simdat1,simdat2,all=TRUE)
simdatall <- simdatalla[,-1]
head(simdatall)
gruAB <- factor(rep(c("A","B"),each=5000))
DAT1 <- data.frame(simdatall,ABgroup = gruAB)
head(DAT1)
reshOBJ <- reshMG(DAT1,items=1:NUMBI,groups=NUMBI+1,correct=rep(0,NUMBI))
# ---------------------------------------------------------------------------------
ergNRM <- nrm(reshOBJ)
mydes <- designTemp(ngru=2,nit=NUMBI,TYPE="NRM")
mydes[[1]][2,1] <- 2
erg_modc <- model_comp(ergNRM,mydes)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.