我正在使用split()
和lapply
函数批量运行Mann Kendall趋势测试。在下面的代码中,split()
通过分析器(水质参数)分隔结果(ConcLow)。然后,lapply分别运行MannKendall
和summary
。输出将发送到控制台(示例代码如下所示),但我希望将其放入Excel或cvs文档中,以便我可以使用它。理想情况下,Excel文档的第一列中应该有分析物(例如TOC),然后是结束列= tau value,第三列= pvalue。然后,下一个选项卡或后面的列将显示summary
函数的结果。如果您能提供任何帮助,我们将非常感谢!我对R很陌生。
mk.analyte <- split(BarkTop$ConcLow, BarkTop$Analyte) lapply(mk.analyte, MannKendall) lapply(mk.analyte, summary)
每个分析物的输出如下所示(在这里缩写,但它是一个很长的列表):
$TOC
tau = 0.0108, 2-sided pvalue =0.8081
$TOC
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.378 2.054 2.255 2.434 2.600 4.530
数据如下所示:
Date Location Analyte ConcLow Units
5/8/2000 Barker Res. Hardness 3.34 mg/L (as CaCO3)
11/24/2000 Barker Res. Hardness 9.47 mg/L (as CaCO3)
6/12/2001 Barker Res. Hardness 1.4 mg/L (as CaCO3)
12/29/2001 Barker Res. Hardness 21.9 mg/L (as CaCO3)
7/17/2002 Barker Res. Fe (diss 81 ug/L
2/2/2003 Barker Res. Fe (diss 90 ug/L
8/21/2003 Barker Res. Fe (diss 0.08 ug/L
3/8/2004 Barker Res. Fe (diss 15.748 ug/L
9/24/2004 Barker Res. TSS 6.2 mg/L
4/12/2005 Barker Res. TSS 8 mg/L
10/29/2005 Barker Res. TSS 10 mg/L
转载请注明出处:http://www.bqysc.com/article/20230526/2080637.html