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)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.