[R] double gaussian with mixdist: what's wrong?
Sepp Tannhuber
sepp.tannhuber at yahoo.de
Tue Nov 20 11:54:05 CET 2012
Dear all,
I am trying to fit a double gaussian to some data using the mixdist package:
--- begin code ---
library(mixdist)
time <- seq(673,723)
counts <- c(3,12,8,12,18,24,39,48,64,88,101,132,198,253,331,419,563,781,1134,1423,1842,2505,374,6099,9343,13009,
15097,13712,9969,6785,4742,3626,3794,4737,5494,5656,4806,3474,2165,1290,799,431,213,137,66,57,41,35,27,27,27)
data.df <- data.frame(time=time, counts=counts)
data.mix <- as.mixdata(data.df)
plot(data.mix)
startparam <- mixparam(c(699,707),1
data.fit <- mix(data.mix, startparam, "norm")
--- end code ---
This results in:
Error in nlm(mixlike, lmixdat = mixdat, lmixpar = fitpar, ldist = dist, :
missing value in parameter
I have used this fit method successfully with other data sets. Why does it not work here?
Best regards
Joseph
More information about the R-help
mailing list