Hi,
I am a newbie at R and was thinking of ways to vectorize this loop but can't
quite figure out if anything more can be done :
Here the dimensions of the arrays are (N=2000, K=2, D=1000):
responsibilities : N * K
pi: K
corp: N*D
theta: K*D
for (n in 1:N)
{
for (k in 1:K)
{
responsibilities[n, k] <- (pi[k]) +
(log1p(dmultinom(x=corp[n, ], prob=exp(theta[k, ]))))
}
denom <- log.sum(responsibilities[n, ])
responsibilities[n,] <- responsibilities[n, ] - denom
}
Any help is much appreciated!
thanks.
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
