Re: [R] code for double sum

On Wed, Jun 03, 2009 at 03:02:26AM -0700, Roslina Zakaria wrote: > Hi R-users, > ? > I wrote a code to evaluate double sum as follows: > ? > ff2 <- function(bb,eta,z,k) > { r <- length(z) > for (i in 1:r) > { sm1 <- > sum((z[i]*bb/2)*(psigamma((0:k)+eta+1,deriv=0)/(factorial(0:k)*gamma((0:k)+eta

Re: [R] code for double sum

Roslina Zakaria wrote: Hi R-users, I wrote a code to evaluate double sum as follows: ff2 <- function(bb,eta,z,k) { r <- length(z) for (i in 1:r) { sm1 <- sum((z[i]*bb/2)*(psigamma((0:k)+eta+1,deriv=0)/(factorial(0:k)*gamma((0:k)+eta+1 sm2 <- sum((besselI(z[i]*bb,eta)*log(z[i]*bb/2) -