cox-stuart实现方法

1、R语言代码# 处理成为符号数据 cox_test <- function(dat = customers) { n <- length(dat) c <- (n - 1)/2 d <- rep(0, c) if (c %% 2 == 0) { # dat长度为偶数 for (i in 1:c) {d[i] <- dat[c + i] - dat[i]} } else { # dat长度为奇数时要去掉中间数 dat <- dat[-(c + 1)] for (i in 1:c) {d[i] <- dat[c + i] - dat[i]} } d <- d[d != 0] binom.test(sum(d > 0), length(d), alt = 'less') } customers <- c(5, 9, 12, 18, 17, 16, 19, 20, 4, 3, 18, 16, 17, 15) cox_test(dat = customers)