[問題] 迴歸分析中的交互作用項消失

看板R_Language作者時間6年前 (2017/09/03 21:55), 編輯推噓0(004)
留言4則, 2人參與, 最新討論串1/1
先附上程式碼: ############## x1 <- c(750, 1000, 1250, 1500, 1750, 2000, 1250, 1250, 1250, 1250, 1250, 1250, 1250, 1250, 1250, 1250) x2 <- c(11, 11, 11, 11, 11, 11, 7, 9, 11, 13, 15, 11, 11, 11, 11, 11) y <- c(0.9 ,1.3, 1.3, 1.4, 1.4, 1.4, 1.9, 1.6, 1.4, 1.2, 1.1, 1.1, 1.1, 1.3, 1.7, 1.7) df <- data.frame(x1, x2, y) model <- lm(y ~ x1 + x2 + x1:x2, data = df) summary(model) ############## 手邊有一些data想做迴歸 上面的資料 是從real data 擷取出問題的部分下來 我的問題是 執行完lm之後 不知道為什麼 交互作用項 x1*x2 一直是NA值 summary之後 還會出現: ####### Coefficients: (1 not defined because of singularities) ####### 但 x1.x2 的值 如果自己隨便亂給 又不會有這個問題 顯然這筆data有點問題?? 上網爬過文 這人好像有和我一樣的錯誤 https://goo.gl/37a3yc 沒有線性獨立??? 看起來應該有啊 而且如果只有x1 , x2時 是沒有問題的 應該不會因為加入交互作用項 就變得線性相依吧??? 求指點 謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.235.148.245 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1504446926.A.026.html

09/03 22:04, , 1F
kappa(model.matrix(y ~ x1 * x2, df))
09/03 22:04, 1F

09/03 22:05, , 2F
check the condition number of model matrix
09/03 22:05, 2F

09/03 22:05, , 3F
算出來是10的20次方.... 應該是variable相依
09/03 22:05, 3F

09/03 22:49, , 4F
試試看先標準化?
09/03 22:49, 4F
文章代碼(AID): #1Ph0dE0c (R_Language)
文章代碼(AID): #1Ph0dE0c (R_Language)