2016-10-17 231 views
3

我已經存儲爲XTS對象的時間序列數據XTS對象和虛設交互項。當迴歸對獨立變量和一個虛設的相互作用術語因變量,輸出原來是與獨立變量,相互作用術語和虛設本身自動迴歸。下面是我做過的一個例子:線性迴歸用

x <- xts(rnorm(100,0,1), Sys.Date()-100:1) 
y <- xts(rnorm(100, 1, 1), Sys.Date()-100:1) 
d <- xts(order.by = index(x)) 
d <- merge(d, dummy = 1) 
d["/2016-09-06"] <- 0 

Call: 
lm(formula = y ~ x + x * d) 

> Coefficients: 
    (Intercept)   x   d   x:d 
     0.95559  0.07350  0.29469  -0.09851 

這對我來說看起來有點奇怪......這是正確的還是我做錯了什麼?

謝謝! (和請原諒我奧赫的一個基本問題。)

回答

2

這就是*意味着配方食品中。如果你只是想交互作用項,使用:代替。來自?formula

這些術語本身包含由':'運算符分隔的變量和因子名稱。這個術語被解釋爲術語中出現的所有變量和因素的相互作用。

的「*」算子表示因子交叉:「A * B」解釋爲「A + B + A:B」。

所以,你會希望使用lm(y ~ x + x:d)(和你最初的嘗試可以減少到lm(y ~ x*d) - 首x是多餘的)。