如何对R中的“用户定义”公式进行线性回归?

时间:2013-09-05 16:31:49

标签: r regression

我有一个包含5个独立变量的数据框,我希望线性方程式采用以下形式:

y = A (a + pA + qB + rC + sD + tE) 

其中ABCDE是我的自变量,pqrst是我需要找到的系数。

2 个答案:

答案 0 :(得分:2)

使用lm()中的R命令,您可以执行以下操作:

#Pseudo Data
y = rnorm(100) 

A = rnorm(100)
B = rnorm(100)
C = rnorm(100)
D = rnorm(100)
E = rnorm(100)

AB = A*B
AC = A*C
AD = A*D
AE = A*E

model = lm(y~-1+A+AB+AC+AD+AE)

产生:

> model

Call:
lm(formula = y ~ -1 + A + AB + AC + AD + AE)

Coefficients:
         A          AB          AC          AD          AE  
 0.1896753  -0.0835971  -0.0183475  -0.0007795  -0.0174815  

> 
> summary(model)

Call:
lm(formula = y ~ -1 + A + AB + AC + AD + AE)

Residuals:
     Min       1Q   Median       3Q      Max 
-2.05531 -0.58641  0.08847  0.73281  2.86074 

Coefficients:
     Estimate Std. Error t value Pr(>|t|)  
A   0.1896753  0.1084157   1.750   0.0834 .
AB -0.0835971  0.1088133  -0.768   0.4442  
AC -0.0183475  0.1264781  -0.145   0.8850  
AD -0.0007795  0.0930502  -0.008   0.9933  
AE -0.0174815  0.1140712  -0.153   0.8785  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.957 on 95 degrees of freedom
Multiple R-squared:  0.03374,   Adjusted R-squared:  -0.01712 
F-statistic: 0.6634 on 5 and 95 DF,  p-value: 0.6521

答案 1 :(得分:1)

:中的变量之间使用formula命令获取没有主效应的乘法交互。例如:

y = A (a + pA + qB + rC + sD + tE)

我不知道你是否试图压制拦截,但我假设你没有将a放在所需的估算列表中。

y ~ 0 + A + A:B + A:C + A:D + A:E