Gamma distribution in R

[This article was first published on Data Science Tutorials, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

The post Gamma distribution in R appeared first on Data Science Tutorials

What do you have to lose?. Check out Data Science tutorials here Data Science Tutorials.

Gamma distribution in R, This guide demonstrates how to use R to fit a gamma distribution to a dataset.

Gamma Distribution Fitting in R

Let’s say you have a dataset z that was produced using the following method:

Create 30 random values that are distributed according to a gamma distribution with shape parameters of 4 and 9 along with some gaussian noise.

Gama <- rgamma(30, 4, 9) + rnorm(30, 0, .03)

Now we can view the first 6 values

head(Gama)
[1] 0.1362240 0.5979568 0.4930604 0.2808689 0.4361617

With the help of the fitdistrplus package in R, we can examine how well a gamma distribution fits this dataset Gama.

if not already installed, install the ‘fitdistrplus’ package.

install.packages('fitdistrplus')

Yes, now we can load the package

library(fitdistrplus)

To fit a distribution using this package, the following general syntax should be used:

fitdist(dataset, distr = “your distribution choice”, method = “your method of fitting the data”)

In this instance, we’ll use the gamma distribution and maximum likelihood estimation approach to suit the dataset z that we created earlier:

Using mle, we fitted our dataset to a gamma distribution.

fit <- fitdist(Gama, distr = "gamma", method = "mle")

Let’s view the summary of the fit

summary(fit)

Fitting of the distribution ‘ gamma ‘ by maximum likelihood

Parameters :
      estimate Std. Error
shape 3.551416   0.647940
rate  7.019582   1.375659
Loglikelihood:  -0.1783264   AIC:  4.356653   BIC:  8.371319
Correlation matrix:
          shape      rate
shape 1.0000000 0.9309661
rate  0.9309661 1.0000000

The following syntax can then be used to create several charts that demonstrate how well the gamma distribution fits the dataset.

par(mar=c(1, 1, 1, 1))
plot(fit)

Further Resources:-
The following tutorials provide guidance on using R to manipulate various probability distributions:

The Multinomial Distribution in R – Data Science Tutorials

The Uniform Distribution in R – Data Science Tutorials

Test for Normal Distribution in R-Quick Guide – Data Science Tutorials

The post Gamma distribution in R appeared first on Data Science Tutorials

Learn how to expert in the Data Science field with Data Science Tutorials.

To leave a comment for the author, please follow the link and comment on their blog: Data Science Tutorials.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)