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.
