异常值的敏感性分析

Question

假设我有两个向量及其相关性：

set.seed(123)    
x<-c(rnorm(19), 99.1)
y<-rnorm(20, 5, 8)
cor(x,y)
#[1] -0.0721101
cor(x[1:19], y[1:19])
#[1] -0.0331376

我的问题是如何根据潜在异常值衡量sensitivity（例如使用p.value）x和y之间的相关性。

Answer 1

您可以使用包outliers来检测异常值：

outliers::grubbs.test(c(x,y))



Grubbs test for one outlier

data:  c(x, y)
G = 5.84030, U = 0.10299, p-value < 2.2e-16
alternative hypothesis: highest value 99.1 is an outlier