Displaying 1 result from an estimated 1 matches for "lifit".
Did you mean:
limit
2011 Nov 20
2
ltm: Simplified approach to bootstrapping 2PL-Models?
...or 2PL-models, I tried to mimic the
bootstrap-approach chosen in the GoF.rasch()-function. Not being a
statistician, I was wondering whether the following simplification
(omit the "chi-squared-expressed model fit-step") would be appropriate:
GoF.ltm <- function(object, B = 50, ...){
liFits <- list()
for(i in 1:B){
rndDat <- rmvlogis(nrow(object$X), coef(object))
liFits[[i]] <- ltm(rndDat ~ z1)
}
distr <- sort(sapply(liFits, function(x)return(x$log.Lik)))
return(max(which(distr <= object$log.Lik))/length(distr))
}
The rationale behind was to directly u...