具有2級因子誤差的R中的線性回歸

Question

我正在嘗試使用2級性別作為分類變量在R中進行線性回歸，但是R會導致錯誤，請問我該如何解決？ 您可以在這里下載我的數據https://www.dropbox.com/s/jbuyearzlw5zf9q/covar.csv?dl=0

非常感謝您的幫助！

> covar = read.csv("covar.csv")
> yv = as.numeric(covar$yv)
> sex = as.factor(covar$sex)
> model = lm(yv~sex)
Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : 
  contrasts can be applied only to factors with 2 or more levels
> head(sex)
[1] 0 0 1 1 0 0
Levels: 0 1
> head(yv)
[1] 23 NA NA NA NA 23

Answer 1

不幸的是，您的性別為1的數據中100％的點具有yv = NA。 您的數據實際上僅包含sex = 0數據。

summary(covar[covar$sex==1,])
       yv              sex   
 Min.   : NA      Min.   :1  
 1st Qu.: NA      1st Qu.:1  
 Median : NA      Median :1  
 Mean   :NaN      Mean   :1  
 3rd Qu.: NA      3rd Qu.:1  
 Max.   : NA      Max.   :1  
 NA's   :187533          

table(covar$sex)
     0      1 
220685 187533 

sum(is.na(covar[covar$sex ==1,1]))
[1] 187533