如何R中

Question

的nls功能正常工作像下面的參數一長串傳遞給`nls`功能：

x <- 1:10 
y <- 2*x + 3       # perfect fit 
yeps <- y + rnorm(length(y), sd = 0.01) # added noise 
nls(yeps ~ a + b*x, start = list(a = 0.12345, b = 0.54321))# 

因爲我用的模型有很多的參數，或者我不知道事先什麼都會被包含在參數列表，我想是以下

tmp <- function(x,p) { p["a"]+p["b"]*x } 
p0 <- c(a = 0.12345, b = 0.54321) 
nls(yeps ~ tmp(x,p), start = list(p=p0)) 

有誰知道如何修改nls功能，以便它可以接受該公式，而不是許多單獨的參數，在參數向量的說法？

Answer 1

你可以給初始化係數向量是這樣的：

tmp <- function(x, coef){ 
     a <- coef[1] 
     b <- coef[2] 
     a +b*x 
    } 

x <- 1:10 
yeps <- y + rnorm(length(y), sd = 0.01) # added noise 
nls(yeps ~ a + b*x, start = list(a = 0.12345, b = 0.54321))#      
nls(yeps ~ tmp(x,coef), start = list(coef = c(0.12345, 0.54321))) 

Nonlinear regression model 
    model: yeps ~ tmp(x, coef) 
    data: parent.frame() 
coef1 coef2 
    3  2 
residual sum-of-squares: 0.0016 

Number of iterations to convergence: 2 
Achieved convergence tolerance: 3.47e-08 

PS：

example(nls) 

應該是一個良好的開端，以瞭解如何與NLS玩。