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]]

______________________________________________
R-help@r-project.org 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.

Reply via email to