[R] Help with function "fitdistr" in "MASS"
Peter Ehlers
ehlers at ucalgary.ca
Mon Jan 4 18:24:47 CET 2010
Saji Ren wrote:
> Thank you,man. the problem solved.
> Plus. when I got the parameters of the data.
> And I used the "truehist(mydata)" to get a histogram of the data,
> How can I draw a line of the distribution of the estimated parameters in the
> histogram plot?
>
> for example:
>> fitdistr(na.exclude(mydata),"normal")
> mean sd
> 4.052594 75.620350
> ( 1.420743) ( 1.004617)
>> truehist(mydata)
>
> After that, I got the histogram of mydata.
> And I want plot a extra line of the density of a normal distribution of
> mean=4.052594 and sd=75.620350 in the histogram plot?
add this:
curve(dnorm(x, mean=4.052594, sd=75.620350), add=TRUE)
-Peter Ehlers
>
> thank you
>
>
>
> Bernardo Rangel tura wrote:
>> On Sat, 2010-01-02 at 23:20 -0800, Saji Ren wrote:
>>> Hi, R users:
>>>
>>> I want to fit my data into a normal distribution by using the command
>>> "fitdistr" in "MASS".
>>> I changed my data class from "ts" to "numeric" by
>>>
>>>> class(mydata)="numeric"
>>> but after using "fitdistr", I got the result below
>>>
>>>> fitdistr(mydata,"normal")
>>> mean sd
>>> NA NA
>>> (NA) (NA)
>>>
>>> the help doc of "fitdistr" does not mention anything about that, thus I
>>> need
>>> your help.
>>>
>>> Thank you in advanced,
>>> Saji from Shanghai
>> Hi Sajj,
>>
>> You hava NA in your data
>>
>> try: fitdistr(na.exclude(mydata),"normal")
>>
>> --
>> Bernardo Rangel Tura, M.D,MPH,Ph.D
>> National Institute of Cardiology
>> Brazil
>>
>> ______________________________________________
>> 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.
>>
>>
>
--
Peter Ehlers
University of Calgary
403.202.3921
More information about the R-help
mailing list