我正在对我的数据使用adf.test来检查其平稳性。我希望将测试的pvalue存储在一个新的变量中,以便我可以在进一步的处理中使用它。基本上,我想做这样的事情: x <- adf.test(Timeseries_1)$pvalue 但这不管用!有什么帮助吗?为了继续上面的内容,我试图从准确性检查中提取MAPE的值,但得到了以下错误。
我正在使用split()和lapply函数批量运行Mann Kendall趋势测试。在下面的代码中,split()通过分析器(水质参数)分隔结果(ConcLow)。然后,lapply分别运行MannKendall和summary。输出将发送到控制台(示例代码如下所示),但我希望将其放入Excel或cvs文档中,以便我可以使用它。理想情况下,Excel文档的第一列中应该有分析物(例如TOC),然后是结束列= tau value,第三列= pval
我有两个系列的列表,从相同的长度开始。在执行以下代码之后,第二个系列比第一个系列少一个元素。是否有一种通用的方法只删除包含n+1元素的系列的最后一个元素,以便我的列表中的所有系列都有n个元素?如果我的列表中包含n、n+1和n+2元素的序列组合如何?下面是一个最小的可重复的例子。<- lapply(tseries, function(x) tseries::adf.test(x)$p.value)
# index only series that need
我是个新手 我在R中执行置换测试,以确定某些SNP的出现百分比是否与我的数据集52K值的向量无关。我正通过Loop原样执行测试: R1_H <- R1[,12] #extract the vector from a dataframe
out <- rep(0,length(R1_H))
out = o