[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