You should report the issue to the author/maintainer of the mixtools
package. gammamixEM can get into this situation when the data is not an
obvious mixture so it has a hard time coming up with a good starting point
for the coefficient estimates. E.g.,
> out <- mixtools::gammamixEM(rep(c(0.0001,
Data has no negative values. Values range from 0.001 to 1.01.
Following is the summary, in case that helps:
Min. 1st Qu. MedianMean 3rd Qu.Max.
0.0010 0.8126 0.8536 0.8464 0. 1.0180
SD: 0.07489977
Any clue?
On Thu, Sep 15, 2016 at 10:32 PM, William Dunlap wrote:
> Does t
Does the data contain non-positive values?
> out <- mixtools::gammamixEM(as.numeric(0:100), lambda = c(1, 1, 1)/3,
verb = TRUE)
iteration = 1 log-lik diff = NaN log-lik = NaN
Error in while (diff > epsilon && iter < maxit) { :
missing value where TRUE/FALSE needed
Bill Dunlap
TIBCO Software
I am using a function gammamixEM where it does it by default. I do not have
the option to change it.
Conceptually, what can make the algorithm not able to calculate likelihood
value at all (and hence log-lik=Nan)? Is there sth wrong with the data?
Under what conditions does it happen?
On Wed, Sep
On 14/09/2016 4:46 PM, Aanchal Sharma wrote:
Hi,
I am trying to fit Gamma mixture model to my data (residual values obtained
after fitting Generalized linear Model) using gammamixEM. It is part of the
script which does it for multiple datasets in loop. The code is running
fine for some datasets
Hi,
I am trying to fit Gamma mixture model to my data (residual values obtained
after fitting Generalized linear Model) using gammamixEM. It is part of the
script which does it for multiple datasets in loop. The code is running
fine for some datasets but it terminates for some giving following err
6 matches
Mail list logo
