R help - May 2009 - Kumaraswamy distribution

Dear R users,

Does anyone know how to write function for Kumaraswamy distribution in R? Since
I cannot write dkumar, pkumar, etc. in R.

 

Please help.

 

Thanks a lot,

Debbie

_________________________________________________________________
[[elided Hotmail spam]]

	[[alternative HTML version deleted]]

Dear Debbie,

On Tue, May 12, 2009 at 5:37 AM, Debbie Zhang <debbie0621 at hotmail.com>
wrote:
>
> Dear R users,
>
> Does anyone know how to write function for Kumaraswamy distribution in R?
Since I cannot write dkumar, pkumar, etc. in R.
>
>
>
> Please help.
>
>
>
> Thanks a lot,
>
> Debbie
>
Check the CRAN Task View on Probability distributions:

http://cran.r-project.org/web/views/Distributions.html

and in particular, check out package VGAM.

HTH,
Jay

dkumar = function(x,a,b) a*b*(x^(a-1))*((1-x^a)^(b-1))
pkumar = function(x,a,b) 1-(1-x^a)^b

(I've based this entirely on the wikipedia entry on the Kumaraswamy
distribution [http://en.wikipedia.org/wiki/Kumaraswamy_distribution],
so best to check both my replication of the formula there and the
accuracy of the wikipedia entry itself)

On Tue, May 12, 2009 at 6:37 AM, Debbie Zhang <debbie0621 at hotmail.com>
wrote:
>
> Dear R users,
>
> Does anyone know how to write function for Kumaraswamy distribution in R?
Since I cannot write dkumar, pkumar, etc. in R.
>
>
>
> Please help.
>
>
>
> Thanks a lot,
>
> Debbie
>
> _________________________________________________________________
> [[elided Hotmail spam]]
>
> ? ? ? ?[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at 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.
>


-- 
Mike Lawrence
Graduate Student
Department of Psychology
Dalhousie University

Looking to arrange a meeting? Check my public calendar:
http://tr.im/mikes_public_calendar

~ Certainty is folly... I think. ~