使用ggplot汇总图表中的Mann Kendall测试结果

Question

我一直在使用代码来总结图形中的一些线性回归结果（均附有代码和图形）。 我采用了从这里引用各种来源的代码。

现在，我想为Mann-Kendall测试结果制作完全相同的图形。 并代替我图表的“斜率”部分，我要插入“ Sen斜率估计”的结果。 所以我想从两个不同的测试结果显示在此图中。

我认为问题出在“ geom（label）”的“ Paste”部分。 因为我必须确定要粘贴在图形中的特定数字结果，所以我不知道该怎么做。

如果有任何机构可以解决此问题。 我将不胜感激。

ggplotRegressionxpl <- function (fit) {
 require(ggplot2)
ggplot(fit$model, aes_string(x = names(fit$model)[2], y = names(fit$model)[1])) + 
geom_point() + geom_line() + geom_label(aes(2000, 33, hjust = 0, vjust = 0, label = paste("R^2 = ",signif(summary(fit)$adj.r.squared, 3),"\n",
                                                                                          "Slope =",signif(fit$coef[[2]], 3),"\n",
                                                                                          "p-value =",signif(summary(fit)$coef[2,4], 3)))) +
stat_smooth(method = "lm", col = "red") + xlab("Year") + ylab("Total Precipitation") +
 labs(title = "Pullman (1941–2018)") + scale_y_continuous(limits = c(0, 40)) +
theme(plot.title = element_text(hjust = 0.5))



}

Answer 1

几年前，我曾经使用过Mann-Kendall / Sens。 我使用的软件包是wql 。

您可以对数据result <- mannKen(data$Precipitation)运行mannKen函数。 完成此操作后，您可以使用result $sen.slope和result$p.value result $sen.slope访问Sens Slope和p值。

为了绘制趋势线，我通常选择一个点来固定该线-例如，时间序列中的中值浓度与时间序列的中点相结合。 在上面的图中，这为我提供了（middleYear，中位数）-（1978，中位数降水）的参考点。

然后，我使用sens斜率，时间序列的起点（在上面的示例中为1941），数据的终点（在示例中为2019）和一些线性代数来计算直线的起点（1941年的降水）如果您将sens斜率从参考点向后延伸）和直线的终点（如果您将sens斜率从参考点向前延伸则为2019年的降水）。

我很想知道是否还有其他软件包可以简化此过程。