Full_Name: Marek Ancukiewicz Version: 1.8.1 OS: Linux Submission from: (NULL) (132.183.12.87) Function cor() incorrectly handles missing observation with method="spearman":
x <- c(1,2,3,NA,5,6) y <- c(4,NA,2,5,1,3) cor(x,y,use="complete.obs",method="s")
[1] -0.1428571
cor(x[!is.na(x)&!is.na(y)],y[!is.na(x)&!is.na(y)],method="s")
[1] -0.4 These two results should be the same.