rollapply

    2热度

    1回答

    在动物园对象的滚动窗口上应用函数通常很顺利,移动平均线: z <- zoo(1:10, as.Date(31:40)) rollapply(z, 4, mean, align="right") 现在我想要做同样的事情用一个统计检验,即适用于每个窗口上Augmented Dickey-Fuller test并得到检验统计量就像我在上面的例子中得到的平均值。 所以基本上我找下面的一段代码相当于(

    0热度

    3回答

    如何总结两列完整案例的总数? 随着c等于: a b [1,] NA NA [2,] 1 1 [3,] 1 1 [4,] NA 1 应用类似 rollapply(c, 2, function(x) sum(complete.cases(x)),fill=NA) 我想拿回单号,2在这种情况下。这将是一个有很多列的大型数据集,所以我想在整个集合中使用rollapply,而不是简单地执行

    0热度

    1回答

    我有一个矩阵(V),每天的股票收益在一定的时间段内,我试图获得之前 1个月(= 21天)和6个月(= 126天)。 这就是我将如何根据未来的1个月和6个月的数据来做到这一点。 Variance.x <- rollapply(data=V, width=21, var, by.column=F, align="right") Variance.x <- rollapply(data=V, widt

    1热度

    1回答

    我有时间系列数据作为data.table类和每列(观察点)有值,我想在滑动窗口(30宽度)内计数它们。 我试图用rle(sort(x))来计算rollapply中的每个值,但它不起作用。 例如,如果我有表像下面, dt <- data.frame(v1=c(1,0,1,4,4,4,4,4),v2=c(1,1,1,4,3,3,3,3), v3=c(0,1,1,3,3,3,3,2),v4

    0热度

    1回答

    我对R很陌生,正试图获得基因表达数据集的手段。这些文件是一个csv文件,有185行(基因名称)和30列(时间点重复三次)。 我使用以下命令导入csv。文件为R read.csv( '〜/桌面/ WNV Spleen.csv',row.names = 1) CM1 CM2 CM3 D1M1 D1M2 D1M3 D2M1 D2M2 etc..... Atf2 3152.9 6178.12 4