1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
# computes a list of expected probabilities for objects of class Rm
# with 0th category included!
plist.internal <- function(object, theta){
X <- object$X
mt_vek <- apply(X, 2L, max, na.rm = TRUE) # number of categories - 1 for each item
mt_ind <- rep(seq_along(mt_vek), mt_vek)
#--------compute list matrix of probabilites for fixed theta)
p.list <- tapply(object$betapar, mt_ind, function(beta.i){
beta.i <- c(0, beta.i)
ind.h <- 0:(length(beta.i)-1)
theta.h <- tcrossprod(ind.h, theta) # ind.h %*% t(theta) # multiply category with
tb <- exp(theta.h + beta.i)
denom <- colSums(tb)
pi.mat <- apply(tb, 1L, function(y){ y/denom })
return(pi.mat)
})
return(p.list)
}
|
