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 ______________________________________________ R-help@r-project.org mailing list 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.
