Cox-Staut 趋势存在性检验year - seq1971 2002rain - c206 223 235 264 229 217 188 204 182 230 223 227 242 238 207 208 216 233 233 274 234 227 221 214 226 228 235 237 243 240 231 210自编函数r语言
cox_stuart <- function(x){ n <- length(x) S <- numeric(n) for(i in 1:n){ for(j in 1:(i-1)){ S[i] <- S[i] + sign(x[i]-x[j]) } } T <- sqrt(n*(n-1)(2n+5)/18) * (sum(S)/n) p_value <- 2*pt(-abs(T), df=n-2) result <- list(statistic=T, p_value=p_value) return(result) }
使用函数进行趋势检验
cox_stuart(rain) # 结果表明存在显著的正向趋势
原文地址: https://www.cveoy.top/t/topic/fUH4 著作权归作者所有。请勿转载和采集!