[R] Applying a function to a vector

From: Florent Bresson <f_bresson_at_yahoo.fr>
Date: Tue 01 Nov 2005 - 01:40:15 EST

I just have a try with sapply. The problem is that my function pbeta2 has two parameters z and p (wich is a vector of two parameters). If I use sapply , R returns a message incicating that parameter p is missing. It is a problem since both z and p are varying along my data.frame.

>Florent
>have a look at:
> help(sapply)

>On Mon, 31 Oct 2005, Florent Bresson wrote:

>> 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 ?
>>
>> Here'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

---

R-help@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 Received on Tue Nov 01 01:47:59 2005