Hi all, I wrote the following code for the function that is attached. I want to simply the code. Can some one give me some help?
The function is attached and the following is the code: #################################################################### ## f1 is the function calculate the value outside the sum f1 <- function(t, a, b) {t/beta(a, b)} ## f2 is the function that calculates the summand f2 <- function(j, k, t, a, b){choose(j-1,k-1)*(1-t)^(j-1) * beta(a+k, b+j- k)} ## f3 calculates the sum f3 <- function(n, t, a, b){ ss <- matrix(0, nrow=n, ncol=n) for (j in 1:n){ for (i in 1:j) { ss[i,j]<- f2(j,i, t, a, b) } } sum <- sum(ss) return(sum) } ## f is the final function f <- function(n,t, a, b) {f1(t,a ,b)*f3(n,t,a,b)} ########################################################################## Many thanks! Hannah
question to R-heip.pdf
Description: Adobe PDF document
______________________________________________ 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.