虛擬編碼的參考組

Question

當使用分類變量在R中使用lm函數建模時，是否有任何明確指定將哪個組用作虛擬編碼的參考組的方法？

Answer 1

我知道最簡單的方法就是先重新編碼係數，以便您想用作參考的等級是第一個等級。

你可以用這樣的功能做到這一點：

recodeFactor<-function(f, ref=levels(f)[1]) 
{ 
    lvls<-levels(f) 
    if(ref== lvls[1]) return(f) 

    lvls<-c(ref, setdiff(lvls, ref)) 
    f<-factor(as.character(f), levels=lvls) 
    return(f) 
} 

Answer 2

你可以做這樣的lm電話：

y <- rnorm(100, 0, 1) 
x <- c(rbinom(50,1,.5),(2*rbinom(50,1,.5))) 
lm(y ~ factor(x,c(1,0,2))) # one way 
lm(y ~ factor(x,c(0,1,2))) # another way, etc.