1

使用EnvStats软件包中的 Kendall 季节性趋势测试来评估水质数据，但我无法将代码调整为：

  1. 返回单个站点/站点的趋势统计数据（到目前为止唯一有效的方法是对每个站点/变量组合进行子集化）

  2. 将输出合并为列表或其他紧凑形式

测试本身很简单：

kendallSeasonalTrendTest(WQconc ~ Season + Year, data = ____)

由于要评估的变量和监测站点组合众多，因此必须有更好的方法。看起来 apply 家族在这里可能是相关的，但到目前为止还没有奏效：

 data.frame':   761 obs. of  13 variables:
 $ Year        : int  2007 2007 2007 2007 2008 2008 2008 2008 2009 2009 ...
 $ Quarter     : Factor w/ 4 levels "Fall","Spring",..: 4 2 3 1 4 2 3 1 4 2 ...
 $ Date        : Factor w/ 319 levels "1/10/2007","1/10/2013",..: 1 177 285 40 5 183 261 43 6 184 ...
 $ Site_ID     : Factor w/ 24 levels "s01541000","s01541500",..: 17 17 17 17 17 17 17 17 17 17 ...
 $ DO          : num  7.57 6.44 6.44 7.79 8.68 8.5 8.21 9.91 9.55 11.7 ...
 $ Flow        : num  33259 8912 2002 803 27680 ...
 $ Iron        : num  196 108 103 126 136 77 17 214 90 8 ...
 $ Magnesium   : num  205 309 75 92 190 285 71 98 320 263 ...
 $ pH          : num  6.75 7.4 8 7.8 6.1 7.4 7.7 7.8 7.6 6.8 ...
 $ Temperature : num  3.5 5.7 27.7 20.1 5.5 9.7 23.8 17.6 0.1 8.7 ...
 $ TDS         : num  163 15 80 92 178 23 77 93 35 33 ...
 $ TP          : num  1 1 1 8 12 7 16 2 8 63 ...

 kendallSeasonalTrendTest(DO ~ Quarter + Year, data =water1)

但这并不完全有效，它返回整个数据集的测试统计信息，而不是按站点 ID。如果结果以列表形式出现，那就太好了

我尝试了 apply 系列的转换，但我也真的不明白它们。欣赏任何想法！

4

1 回答 1

1

只需使用lapplywhich 返回kendallSeasonalTrendTest()等于输入向量/列表长度的结果列表，特别是使用unique(). 两种子集类型如下所示。

然后lapply遍历每个变量。此外，用于setNames()命名每个列表元素以便于参考（例如，DO_results$s01541000）

# DO Results
DO_results <- lapply(unique(water1$Site_ID), function(s) 
                kendallSeasonalTrendTest(DO ~ Quarter + Year, 
                                         data = subset(water1, Site_ID == s)))
# NAMING EACH LIST ELEMENT BY CORRESPONDING SITE_ID
DO_results <- setNames(DO_results, unique(water1$Site_ID))

# Flow Results
Flow_results <- lapply(unique(water1$Site_ID), function(s) 
                  kendallSeasonalTrendTest(Flow ~ Quarter + Year, 
                                           data = water1[water1$Site_ID == s,]))  
# NAMING EACH LIST ELEMENT BY CORRESPONDING SITE_ID
Flow_results <- setNames(Flow_results, unique(water1$Site_ID))
...
于 2017-04-19T15:43:32.940 回答