我的Mann Kendall R脚本有问题。
基于总和(数量)进行计算时,我得到了错误的P值,但是它基于记录数为我提供了准确的P值。
数据集%>%突变(
MannKendallY = if (length(Total) > 0)
{MK <- MannKendall(TotalQuantity)
pvalue <- as.numeric(MK$s)} else{
pvalue <- 1}
) %>% mutate(
TauY = if (length(Total) > 0)
{MK <- MannKendall(Total)
tau <- as.numeric(MK$tau)} else{
tau <- 0}
)
感谢任何投入。
非常感谢。