我想将参考分布与使用 Kolmogorov-Smirnov 距离按比例绘制d_1
的样本进行比较。d_2
w_1
鉴于这d_2
是加权的,我正在考虑使用 R 中的加权经验累积分布函数(使用ewcdf {spatstat}
)来解决这个问题。
下面的例子表明我可能没有指定权重,因为当lenght(d_1) == lenght(d_2)
Kolmogorov-Smirnov 没有给出 0 值时。
有人可以帮我弄这个吗?为清楚起见,请参阅下面的可重现示例。
#loop for testing sample sizes 1:length(d_1)
d_stat <- data.frame(1:1000, rep(NA, 1000))
names(d_stat) <- c("sample_size", "ks_distance")
for (i in 1:1000) {
#reference distribution
d_1 <- rpois(1000, 500)
w_1 <- d_1/sum(d_1)
m_1 <- data.frame(d_1, w_1)
#sample from the reference distribution
m_2 <-m_1[(sample(nrow(m_1), size=i, prob=w_1, replace=F)),]
d_2 <- m_2$d_1
w_2 <- m_2$w_1
#ewcdf for the reference distribution and the sample
f_d_1 <- ewcdf(d_1)
f_d_2 <- ewcdf(d_2, 1/w_2, normalise=F, adjust=1/length(d_2))
#kolmogorov-smirnov distance
d_stat[i,2] <- max(abs(f_d_1(d_2) - f_d_2(d_2)))
}
d_stat[1000,2]