Печать модели (nls) характеристик на участке - PullRequest
Новая
       21

Печать модели (nls) характеристик на участке

2 голосов
/ 21 июня 2011

У меня есть экспоненциальная модель, рассчитанная с

mod <- nls(y ~ exp(- (x/a)^b), data = DF, start = list(a = 200, b = 1.4))

который я построил используя:

plot(x,y)<br> lines(sort(DF$x),predict(mod, list(x=sort(DF$x))), lwd=2, col="red")

пишу mod Я получаю: 

Nonlinear regression model 
  model:  y ~ exp(-(x/a)^b) 
   data:  DF 
       a        b 
    211.7098   0.3908 
 residual sum-of-squares: 17.69

Number of iterations to convergence: 5 
Achieved convergence tolerance: 1.477e-07

С summary(mod) Я получаю: 

Formula: y ~ exp(-(x/a)^b)
Parameters:
   Estimate Std. Error t value Pr(>|t|)       
a 2.117e+02  2.799e+00   75.64   <2e-16 ***
b 3.908e-01  9.154e-03   42.69   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
Residual standard error: 0.08832 on 2268 degrees of freedom
Number of iterations to convergence: 5 
Achieved convergence tolerance: 1.477e-07 
  (8 observations deleted due to missingness)

Мой вопрос: есть ли способ добавить формулу, параметры и остаточную сумму квадратов к графику? Единственный способ, которым я мог придумать, - это использовать text(), копировать и вставлять информацию, которую я хотел, но мне придется делать много разных моделей, поэтому мне бы хотелось что-то автоматическое. Спасибо!

1 Ответ

2 голосов
/ 21 июня 2011

Попробуйте это, когда первая строка создаст график, который мы можем использовать в этом примере, и Следующая строка добавляет желаемый текст: 

example(nls)
legend("topleft", legend = capture.output(summary(nlmod)), cex = .5)
...