欠損値をペアワイズで除去
pw.na.rm<-function(x,y){
non.na <- !is.na(x)&!is.na(y)
xx <- x[non.na]
yy <- y[non.na]
return(list(xx,yy))
}
例
x<-c(3,4,NA,NA,4)
y<-c(3,NA,4,NA,1)
pw.na.rm(x,y)
1
[1] 3 4
2
[1] 3 1
# complete.casesを使う
x<-c(3,4,NA,NA,4)
y<-c(3,NA,4,NA,1)
z<-cbind(x,y)
z[complete.cases(x,y),]
