Fitting Mixed Distributions in the fitdistrplus package
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:
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
[[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.