Skip to content

Fitting Mixed Distributions in the fitdistrplus package

2 messages · Charles Thuo, Abby Spurdle

Charles Thuo
# 
Dear Sirs,

The below listed code fits a gamma and a pareto distribution to a data set
danishuni. However the distributions are not appropriate to fit both  tails
of the data set hence a mixed distribution is required  which has ben
defined as "mixgampar"
as shown below.

library(fitdistrplus)
x<- danishuni$Loss
 fgam<- fitdist(x,"gamma",lower=0)
 fpar<- fitdist(x,"pareto",start = list(shape=2,scale=2),lower=0)
fmixgampar<- fitdist(x,"mixgampar",start =
list(prob=1/2,nu=1,lambda=1,alpha=1,theta=1),lower=0)
Error in fitdist(x, "mixgampar", start = list(prob = 1/2, nu = 1, lambda =
1,  :

 The  dmixgampar  function must be defined

Kindly assist to define the dmixgampar

Charkes
Abby Spurdle
# 
dmixgampar <- function (x, param1, param2, ...)
{
        #compute density at x
        <your code>
}
On Wed, Oct 21, 2020 at 8:03 PM Charles Thuo <tcmuigai at gmail.com> wrote: