[R] Applying a function to a vector

Paulo Justiniano Ribeiro Jr paulojus at est.ufpr.br
Mon Oct 31 15:05:41 CET 2005


Florent
have a look at:
> help(sapply)


On Mon, 31 Oct 2005, Florent Bresson wrote:

> I have defined a function to compute the value of a
> beta distribution of the second kind (the existing
> beta distribution of th stats package is the beta
> distribution of the first kind). It works perfectly
> for a single value, but I want to apply it to a vector
> of 22 000 values. I can use a loop for the calculation
> of each value but it runs very very slowly.
> So, what can I change ?
>
> Hers's the function :
> p <-  c(1,1)
> y <-  1
> z <-  1
> truc  <-  function(y)  {y^(p[1]-1)/(1+y)^(p[1]+p[2])}
> pbeta2  <-  function(z,p)
> 1/beta(p[1],p[2])*integrate(truc,0,z)$value
>
> machin <- pbeta2(vector,p) just return a single value
>
> Thanks for your help
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>
>

Paulo Justiniano Ribeiro Jr
LEG (Laboratório de Estatística e Geoinformação)
Departamento de Estatística
Universidade Federal do Paraná
Caixa Postal 19.081
CEP 81.531-990
Curitiba, PR  -  Brasil
Tel: (+55) 41 3361 3573
Fax: (+55) 41 3361 3141
e-mail: paulojus at est.ufpr.br
http://www.est.ufpr.br/~paulojus


More information about the R-help mailing list