[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