Displaying 1 result from an estimated 1 matches for "mweib".
Did you mean:
mweb
2003 Jul 28
1
Optimization failed in fitting mixture 3-parameter Weibull distri bution using fitdistr()
...#
Distribution 2
> fitdistr(x4, w3den, start=list(a = 5.0, b = 90, c = 3))
a b c
-67.572244 159.020786 5.835588
(105.975870) (107.842846) ( 4.288019)
Warning messages:
1: NaNs produced in: log(x)
2: NaNs produced in: log(x)
> mweib <- function(x, p, a,b,c,a1, b1, c1) {p*(c/b*((x -
a)/b)^(c-1)*exp(-((x-a)/b)^c))+
+ +(1-p)*(c1/b1*((x -a1)/b1)^(c1-1)*exp(-((x-a1)/b1)^c1))}
# define mixture distribution
> x5 <- c(x3, x4)
> fitdistr(x5, mweib, start=list(p = 0.7, a = 8.9, b=90.77, c = 3.68, a1 =
-67.57, b1...