基於R適合Y〜EXP（-x ^一）

Question

我只是用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: parent.frame() 
    a 
3.612 
residual sum-of-squares: 0.06132 

Number of iterations to convergence: 9 
Achieved convergence tolerance: 3.654e-06 

Answer 1

您可以通過執行

coef(m)[["a"]] 

示例（型號從help(nls)）獲得的a值：

x <- 1:10 
y <- 2*x + 3        
yeps <- y + rnorm(length(y), sd = 0.01) 
(m <- nls(yeps ~ a + b*x, start = list(a = 0.12345, b = 0.54321))) 
# Nonlinear regression model 
# model: yeps ~ a + b * x 
# data: parent.frame() 
# a b 
# 3 2 
# residual sum-of-squares: 0.000834 
# 
# Number of iterations to convergence: 2 
# Achieved convergence tolerance: 5.23e-09 
coef(m)[["a"]] 
# [1] 3.0052