you could help, Now I
>>> want to pass this curve to a function that can generate random numbers
>>> distributed according to DPLN ( right curve).
>>>
>>> I found the package Runuran can do that
>>> http://cran.r-project.org/web/packages/Runuran/Runuran.pdf
e
>>> http://www.math.uvic.ca/faculty/reed/dPlN.3.pdf i.e. the original paper on
>>> Double Pareto Log Normal Distribution ). > > can you please see the
>>> correction in this
>>> linkhttp://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0
>> ******************
>> Bander
>> *
>>
>>
>>
>> From: d...@vims.edu
>> To: cs_2...@hotmail.com
>> Subject: Re: [R] Double Pareto Log Normal Distribution DPLN
>> Da
>
> norm2<-pnorm((log(x)-v+(b*t^2))/t)
> expo2<- -b*t+(b^2*t^2/2)
>
> y<- x^(b-1)*exp(expo2)*(1-norm2) # 1-norm is the complementary CDF of N(0,1)
> j*(z+y)
> }
> ******
> Bander Alzahrani, Teacher
ander Alzahrani, Teacher Assistant
> Information Systems Department
> Faculty of Computing & Information Technology
> King Abdulaziz University
>
> *********
>
>
>
>> From: d...@vims.edu
>> To: cs_2...@hotmail.com
>> CC: r-h...@stat.math.ethz
***
> From: d...@vims.edu
> To: cs_2...@hotmail.com
> CC: r-h...@stat.math.ethz.ch
> Subject: Re: [R] Double Pareto Log Normal Distribution
> Date: Tue, 12 Nov 2013 16:51:22 +
>
>
> http://www.math.uvic.ca/faculty/reed/dPlN.3.pdf is the original ref and
http://www.math.uvic.ca/faculty/reed/dPlN.3.pdf is the original ref and has the
equations.
library(VGAM) for *pareto() and library(stats) for *lnorm() should get you most
of the way there.
On Nov 12, 2013, at 10:47 AM, "b. alzahrani"
wrote:
> Hi guys
> I would like to generate random number
There's a recipe involving two exponentials and a normal deviate as formula (6)
in (1st hit on Google for "dpln distribution").
http://www.math.uvic.ca/faculty/reed/dPlN.3.pdf
It should be a no-brainer to code it up in R.
-pd
On 12 Nov 2013, at 16:47 , b. alzahrani wrote:
> Hi guys
> I woul
Hi guys
I would like to generate random number Double Pareto Log Normal Distribution
(DPLN). does anyone know how to do this in R or if there is any built-in
function.
Thanks
**
Bander
*
9 matches
Mail list logo
