[R] Likelihood of multiple random processes

Dan Bebber danbebber at forestecology.co.uk
Thu Nov 2 12:34:12 CET 2006


> Do you know which observations came from which process (you do in your 
> example, but it whatever it is supposed to emulate)?

No- I need to determine whether more than one process could have been 
involved.

> If not, you want to estimate a mixture distribution.  This is covered in 
> chapter 16 of the book fitdistr() supports as well as in several CRAN 
> packages.

Many thanks- I didn't know they were called mixture models.
Page 437 is in front of me now.

>
> On Thu, 2 Nov 2006, Dan Bebber wrote:
>
>> Is there a method in R for estimating the likelihood that two (or more)
>> random variables generated a univariate sample?
>>
>> For a single random variable there is:
>> x <- c(rnorm(20,10,3),rnorm(20,20,5))
>> plot(density(x))
>> logLik(fitdistr(x,"normal"))
>>
>> Is there a way of of specifying that more than one normal distribution
>> should be fitted?
>> If so, then I suppose that the AIC of the two models could be calculated.
>>
>> Many thanks,
>> Dan Bebber
>>
>> Department of Plant Sciences
>> University of Oxford
>>
>> ______________________________________________
>> 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
>> and provide commented, minimal, self-contained, reproducible code.
>>
>
> -- 
> Brian D. Ripley,                  ripley at stats.ox.ac.uk
> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
> University of Oxford,             Tel:  +44 1865 272861 (self)
> 1 South Parks Road,                     +44 1865 272866 (PA)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595
>
>
> -- 
> No virus found in this incoming message.
> Checked by AVG Free Edition.
> Version: 7.1.409 / Virus Database: 268.13.23/513 - Release Date: 
> 02/11/2006
>
>



More information about the R-help mailing list