Nothing
ZFnlm <-
function(Ulstv,erg_estep,startOBJ,reshOBJ,quads)
{
SKEL <- startOBJ$stwm1
Q <- reshOBJ$Qmat
if(all(!is.na(startOBJ$setC)))
{
bigv <- vector(mode="numeric",length=ncol(Q))
bigv[-startOBJ$setC$whichetas] <- Ulstv
bigv[startOBJ$setC$whichetas] <- startOBJ$setC$whichconstant
Ulstv <- bigv
}
opp <- as.vector(Q %*% Ulstv)
relstv <- relist(opp,SKEL)
ALL_G <- mapply(function(stvl,ql,RI,FI0)
{
#riq_quer_mat <- sapply(RI$riq_querG,function(x)x)
#f_iq <- sapply(RI$riqv_querG,colSums) + riq_quer_mat # nodes x items
ZQsternl <- lapply(stvl,function(pitem)
{
# 2PL Part
# ----------------------
Km <- matrix(c(rep(1,length(ql$nodes)),ql$nodes),ncol=2)
woba <- grep("(beta|alpha)",names(pitem))
abpar <- pitem[woba]
solit <- twoplpart(Km=Km, abpar=abpar)
dosolit <- 1-solit
# NRM Part
# ----------------------
pitemNRM <- pitem[-woba]
ZQstern <- coP_nrm(pitemNRM,Km)
ZQstern_all <- cbind(solit,ZQstern)
ZQstern_all
})
# extract and compute
logZQstern_M <- log(t(do.call("cbind",ZQsternl)))
EmPe <- 1-sapply(ZQsternl,function(getfirst)getfirst[,1])
# sum to loglikelihood
T1uT3 <- sum(logZQstern_M * RI)
T2 <- sum(FI0 * (log(EmPe)))
allin <- T1uT3 + T2
# logPiXq <- sapply(ZQsternl,function(x)log(x[,1]))
# T1 <- sum(logPiXq * riq_quer_mat)
#
# sumriqv_quer <- sapply(RI$riqv_querG,function(x) colSums(x))
# log_EMINUS_PiXq <- sapply(ZQsternl,function(x)log(1 - x[,1]))
# T2 <- sum(log_EMINUS_PiXq * sumriqv_quer)
#
# logPivuis0 <- lapply(ZQsternl,function(x)log(x[,-1]))
# T3 <- sum(mapply(function(A,B){sum(t(A) * B)},A=RI$riqv_querG,B=logPivuis0))
#
# allin <- T1 + T2 + T3
return(allin)
# },stvl=relstv,ql=quads,RI=erg_estep,SIMPLIFY = FALSE)
},stvl=relstv,ql=quads,RI=erg_estep$riqv_querG,FI0=erg_estep$fique0G,SIMPLIFY = FALSE)
# stvl=relstv[[1]]
# ql=quads[[1]]
# FI = erg_estep$fiqG[[1]]
# RI=erg_estep$riqv_querG[[1]]
# FI0 = erg_estep$fique0G[[1]]
endsum <- do.call("sum",ALL_G)
#endsum <- sapply(ALL_G,sum)
return(endsum)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.