我正在尝试运行我在R中创建的名为自相关的函数,并希望对保存在列表中的所有数据框运行该函数
这是我当前的代码
autocorrelation <- function(x,y){
mean_x <- mean(x)
Denominador <- (x[] - mean(x))^2 #Denominator
k=1; n=length(x)-k; u=1
for (u in 1:n) {
Numerador[u] <- (x[u] - mean_x)*(x[u+k] - mean_x) #Numerador
}
Lag1 <- ifelse(sum(Denominador)!=0,((1/(length(x)-1))*sum(Numerador))/((1/length(x))*sum(Denominador)),0) #Lag 1
low <- (-1 - (1.96*sqrt(length(x)-2)))/(length(x)-1) #5% significance level
up <- (-1 + (1.96*sqrt(length(x)-2)))/(length(x)-1) #5% significance level
Significance <- ifelse( up <= Lag1 & Lag1 >= low,"Independent","No independent") #Significance
Autocorrelation_results <- structure(list(Name=y,Lag1=Lag1,Low=low,Up=up,Significance=Significance)) #save the results
}
for (i in seq_along(listDF_Ann)){
x<- listDF_Ann[[i]]$Qmean
y <- names(listDF)[i]
autocorrelation(x,y)
}
这是我的输入
> dput(listDF_Ann[[1]]$Qmean)
c(80.74, 64.152602739726, 68.7273224043716, 132.860273972603,
337.284931506849, 78.3835616438356, 204.416120218579, 58.7824657534247,
315.110684931507, 146.238356164384, 127.573770491803, 286.835616438356,
353.802739726027, 88.3093150684931, 304.866666666667, 172.287123287671,
90.4186301369863, 51.1808219178082, 131.634153005464, 608.193150684932,
49.1172602739726, 257.268493150685, 112.969945355191, 50.0246575342466,
275.166849315068, 548.178082191781, 405.07650273224, 97.0849315068493,
202.31095890411, 452.841095890411, 206.071038251366, 518.246575342466,
182.787945205479, 309.183561643836, 233.625683060109, 83.5260273972603,
51.4739726027397, 340.13698630137, 308.68306010929, 616.715068493151,
107.912328767123, 342.956164383562, 37.9808743169399, 151.323287671233,
231.164383561644, 90.786301369863, 38.0136612021858, 185.421917808219,
38.0167123287671, 142.664109589041, 114.443442622951, 281.054246575342,
92.0813698630137, 96.5923287671233, 289.146721311475, 149.043287671233,
251.462739726027, 54.4093150684932, 81.118306010929, 168.744931506849,
73.9879452054794, 124.941917808219, 151.301366120219, 258.280547945205,
29.5304109589041, 263.763561643836, 218.915027322404)
当前代码可以工作,但我不知道如何保存结果(结构?)当我运行函数时,对于每个数据帧,x
。在我想导出所有x的结果之后。你知道吗?
转载请注明出处:http://www.lechuangzk.com/article/20230526/2602498.html