Dear R users,

I have been trying to compute the following function and need it to work
with n=15000, but it would only compute for smaller ns, such as n=1000 and
not above. I was wondering if anyone would have a solution for this problem!
Thank you very much for your kind support!

Sincerely,

Nan

------

Wi <- function(n)    {
    fun <- function(w,i){
        w*(-log(w))*((w-w*log(w))^(i-1))*((1-w+w*log(w))^(n-i));
    }
    out <- as.double(n);

    for(i in 1:n){
        out[i] <- integrate(fun,lower=0,upper=1,i=i)$value;
        out[i] <- out[i]*n*choose(n-1,i-1);
    }
    out
}

	[[alternative HTML version deleted]]

