# Gamma distribution in R

**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.

**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.