Thank you for your reply. I was just trying to make the code less complicated. I got your point. Thanks!
2010/2/11 Petr PIKAL <petr.pi...@precheza.cz> > Hi > > I did not see the answer yet, probably nobody is much interested in such > function, so I will try. > > r-help-boun...@r-project.org napsal dne 10.02.2010 21:57:42: > > > Hi all, > > I wrote the following code for the function that is attached. > > I want to simply the code. > > Why. If it works there is no problem until you encounter performance > issues. If it does not work you shall specify what is wrong with it. > > Many functions are vectorised (including choose and beta) so you could > probably get rid of nested cycles but I am far from elaborating how to > change your code and make it simpler (shorter, quicker ???) > > Regards > Petr > > > 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 > > [pøíloha question to R-heip.pdf odstranìna u¾ivatelem Petr PIKAL/CTCAP] > > ______________________________________________ > > 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<http://www.r-project.org/posting-guide.html> > > and provide commented, minimal, self-contained, reproducible code. > > [[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.