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
