R/confint.est_multi_poly_between.R

In MLCIRTwithin: Latent Class Item Response Theory (LC-IRT) Models under Within-Item Multidimensionality

confint.est_multi_poly_between <-function(object, parm, level=0.95, ...){
	
# preliminaries
	parm = NULL
	out = object
# print output
	if(is.null(out$seTh)) stop("Use option out_se=TRUE in est_multi_poly_within()")
	za = qnorm(1-(1-level)/2)	
    cat("\nConfidence interval for abilities:\n")
    L1Th = out$Th-za*out$seTh
    L2Th = out$Th+za*out$seTh
    Tmp = cbind(L1Th,L2Th)
    ind = NULL
    for(j in 1:ncol(out$Th)) ind = c(ind,j,j+ncol(out$Th))
	Tmp = matrix(Tmp[,ind],nrow(out$Th),2*ncol(out$Th))
    dimnames(Tmp) = list(dimension=1:nrow(out$Th),class=paste((1:ncol(out$Th))%x%c(1,1),c(1,2),sep="_"))
    print(round(Tmp,4))    
    cat("\nConfidence interval for the item parameters:\n")
    Tmp = out$Bec-za*out$seBec
    for(j in 1:ncol(out$Bec)) colnames(Tmp)[j] = paste("beta",j,"_1",sep="")
    L1Items = cbind(gamma_1=out$gac-za*out$segac,Tmp)
    Tmp = out$Bec+za*out$seBec
    for(j in 1:ncol(out$Bec)) colnames(Tmp)[j] = paste("beta",j,"_2",sep="")
    L2Items = cbind(gamma_2=out$gac+za*out$segac,Tmp)
    Tmp = cbind(L1Items,L2Items)
    ind = NULL
    for(j in 1:ncol(L1Items)) ind = c(ind,j,j+ncol(L1Items))
    Tmp = Tmp[,ind]
    print(round(Tmp,4))              
    cat("\nConfidence interval for the regression coefficients:\n")
    L1De = out$De-za*out$seDe
    L2De = out$De+za*out$seDe
   	Tmp = cbind(L1De,L2De)
    names = dimnames(Tmp)
    ind = NULL
    for(j in 1:ncol(out$De)) ind = c(ind,j,j+ncol(out$De))
    Tmp = Tmp[,ind]
    if(is.vector(Tmp)) Tmp = t(Tmp)
    dimnames(Tmp) = list(names[[1]],logit=paste(names[[2]][ind],c(1,2),sep="_"))
    print(round(Tmp,4))
    cat("\n")
# output
	out = list(L1Th=L1Th,L2Th=L2Th,L1Items=L1Items,L2Items=L2Items,L1De=L1De,L2De=L2De)    
	
}