如何拟合线性模型

Question

我想创建一个这样的函数：

在 R 中，可以将其用作一般示例，我发现：

fit <- lm(y ~ x1 + x2, data=mydata)

我的数据是：

Ay Ty year
20 11 2009
48 27 2010
35 58 2011

怎么可能为我的案件制作它？

Answer 1

显然你想拟合一个线性模型。这是通过lm函数完成的。通常，您有更多的数据行而不是要拟合的系数，但它会起作用：

> d <- read.table(header=TRUE, text="Ay Ty year
+ 20 11 2009
+ 48 27 2010
+ 35 58 2011")
> d 
  Ay Ty year
1 20 11 2009
2 48 27 2010
3 35 58 2011
> model <- lm(Ay ~ Ty + year, data = d)
> summary(model)
Call:
lm(formula = Ay ~ Ty + year, data = d)
Residuals:
ALL 3 residuals are 0: no residual degrees of freedom!
Coefficients:
              Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.441e+05         NA      NA       NA
Ty          -2.733e+00         NA      NA       NA
year         7.173e+01         NA      NA       NA
Residual standard error: NaN on 0 degrees of freedom
Multiple R-squared:      1, Adjusted R-squared:    NaN 
F-statistic:   NaN on 2 and 0 DF,  p-value: NA

这给你留下了

Ay = -144100 - 2.733 * Ty + 71.73 * year + error

结果具有最小平方残差。