Conditional nls

前端 未结 1 419

I\'m trying to fit the conditional nls with R 2.15.1. The same code was working fine with R 2.13 but now R 2.15.1 throws errors.

x <- seq(from = 17, to =          


        
相关标签:
1条回答
  • 2020-12-22 06:03

    This code works in R 2.15.1.

    x <- seq(from = 17, to = 47, by = 5)
    y <- c(26.2, 173.6, 233.9, 185.9, 115.4, 62.0, 21.7)
    Data <- data.frame(y, x)
    
    Fit <- nls(formula =  y ~ ifelse(test = x <= Mu, ye = c1*exp(-((x-Mu)/Sigma11)^2), no = c1*exp(-((x-Mu)/Sigma12)^2)),
           data = Data, start = list(c1 = 240, Mu = 25, Sigma11 = 5, Sigma12 = 14), algorithm = "port",
           lower = list(Sigma11 = 0, Sigma12 = 0))
    
    0 讨论(0)
提交回复
热议问题