[R] Plot gamma distribution with relative frequency histogram
David Winsemius
dwinsemius at comcast.net
Fri Dec 5 09:32:24 CET 2014
On Dec 4, 2014, at 6:55 PM, dila radi wrote:
> Dear R users,
>
> I am looking to fit a gamma curve onto a histogram of the data.
>
> Consider dt1 as my data:
> c(203.9, 91.5, 24.5, 34.5, 164, 144, 160.5, 195, 191.5, 189,
> 133, 110.5, 155, 80.5, 250.5, 116, 145, 118.5, 406, 183.5, 142.5,
> 197, 367, 134, 153.5)
>
> I am using this code in R
> x<-dt1
> h<-hist(x,breaks=seq(0,500,25),col=rgb(0,0,1), plot=F)
> h$counts <- h$counts / sum(h$counts)
> plot(h, freq=TRUE, ylab="Relative Frequency", xlab="Rainfall (mm)",
> main = 'Histogram KBU-ML Gamma Disribution')
> est <- fitdistr(dt1,"gamma")$estimate
> curve(dgamma(x, rate=0.02220921, shape=3.63421746),from=0, to=500,
> main="Gamma distribution", add = T, col = "red")
>
> but the gamma fit is not fitted onto histogram as much as it supposed to
> be. How should I overcome this?
>
First read ?hist. Then set freq=FALSE. Then you will be comparing "like to like".
Otherwise you cna just hack it by multiplying the dgamma expression by 25.
> sum(h$counts) #before the standardization
[1] 25
> Thank you in advanced.
>
> Dila
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.
David Winsemius
Alameda, CA, USA
More information about the R-help
mailing list