该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
函数以便它可以接受公式中的参数向量参数而不是许多单独的参数?