悬赏 pairwise interaction in multiple linear regression

这是我写的一个ｒｅｇｒｅｓｓｉｏｎ　ｍｏｄｅｌ(mod4.lm)，下面两个是同一个模型的不同种写法，因为原模型中有自变量没有意义。　　
housing.value <- read.table("http://archive.ics.uci.edu/ml/machine-learning-databases/housing/housing.data", header=F)
colnames(housing.value) <- c("CRIM", "ZN", "INDUS", "CHAS", "NOX", "RM", "AGE", "DIS", "RAD", "TAX", "PTRAIO", "B", "LSTAT", "MEDV")
colnames(housing.value)

mod4.lm <- lm(MEDV~CRIM + ZN + CHAS + NOX + RM + DIS + PTRAIO + B + LSTAT, data=housing.value)
mod4.lm <- lm(MEDV~.-TAX-RAD-AGE-INDUS, data=housing.value)

我想把里面加入ｐａｉｒｗｉｓｅ　ｉｎｔｅｒａｃｔｉｏｎ　ｆａｃｔｏｒｓ　为了方便，重新定义了自变量。　但是，为什么我写的ｆｏｒ　ｌｏｏｐ　不对啊？

x1 <- housing.value$CRIM
x2 <- housing.value$ZN
x3 <- housing.value$CHAS
x4 <- housing.value$NOX
x5 <- housing.value$RM
x6 <- housing.value$DIS
x7 <- housing.value$PTRAIO
x8 <- housing.value$B
x9 <- housing.value$LSTAT

for(w in 1:9)
  for(v in 1:9)
{
   {
    a <- sum(x[w]*x[v])  #这是为了把那些成对的interaction factors 加和。
   }
}
   mod5.lm <- lm(MEDV~.-TAX-RAD-AGE-INDUS + a, data=housing.value)

特此悬赏解答，谢谢关注。

用 (x1+x2….)^2 就可以得到up to 2的interaction

如果只需要interaction

(x1+x2….)^2 - x1 - x2 -...

for(w in 1:9)
  for(v in 1:9)
{
   {
    a <- sum(x[w]*x[v])  #这是为了把那些成对的interaction factors 加和。
    mod5.lm <- lm(MEDV~.-TAX-RAD-AGE-INDUS + a, data=housing.value)
   }
}
不是太清楚你的a是否每对interaction factors后就进行一次回归分析，如果不是而是所有interaction factors和之后进行回归分析应该把回归方程放在所有循环在外，这是我的一点理解，不知道对不对
for(w in 1:9)
  for(v in 1:9)
{
   {
    a <- sum(x[w]*x[v])  #这是为了把那些成对的interaction factors 加和。
      }
}
mod5.lm <- lm(MEDV~.-TAX-RAD-AGE-INDUS + a, data=housing.value)

统计R浪人 发表于 2013-10-29 13:49
for(w in 1:9)
  for(v in 1:9)
{

Error: object 'x' not found 这个是需要解决的， 然后在算那些成对的interaction factors 的和。

你这个a是 通过    两两相互乘积后 组建个矩阵， 然后 矩阵以行sum，变为一个向量，是吗？
如果这么的话，看看有用不？

1. x <- data.frame(x1=housing.vaule$...,x2=....)

复制代码

1. dat<-NULL
   
2. dat<-as.data.frame(dat)
   
3. for(w in 1:9)
   
4.   {
   
5.   for(v in 1:9)
   
6.   {
   
7.   a1 <- x[w]*x[v]
   
8.   dat<-rbind(dat,a1)
   
9. }
   
10.   }
    
11. a <- apply(t(dat),1,sum)

复制代码

formula里面的加号不是加减乘除的加号

ntsean 发表于 2013-10-29 21:26
用 (x1+x2….)^2 就可以得到up to 2的interaction

如果只需要interaction

length 不符合。

jmpamao 发表于 2013-10-29 15:05

还是有那个 error object 'x' not found.

reduce_fat 发表于 2013-10-30 00:12
length 不符合。

lm(MEDV~(CRIM + ZN + CHAS + NOX + RM + DIS + PTRAIO + B + LSTAT)^2, data=housing.value)

这个不行吗?