我試圖用nls()函數來擬合數據,其中數據的性質使得我可以限制一個係數和兩個係數的總和。讓我來介紹一個簡短的例子,看看問題出在哪裏。我希望參數b1介於0和1之間,我希望參數b1和b2之和也介於0和1之間。 set.seed(123)
# example where everything is OK
x <- 1:200
g <- rbinom(200, 1, 0.5)
y <- 3 +
我只是用NLS 我試圖適應 y ~ exp(-x^a) :
m<- nls(y ~ exp(-x^a),start=list(a=0))
它seemes做工精細:) 我的問題是,我想開始爲了將值寫入文件而跳出「a」。 我可以看到該值是3.612,但如何獲得該值? m
Nonlinear regression model
model: y ~ exp(-x^a)
data