我有数百个文本文件,每个文件中包含以下信息:
*****Auto-Corelation Results******
1 .09 -.19 .18 non-Significant
*****STATISTICS FOR MANN-KENDELL TEST******
S= 609
VAR(S)= 162409.70
Z= 1.51
Random : No trend at 95%
*****SENs STATISTICS ******
SEN SLOPE = .24
现在,我想读取所有这些文件,并从每个文件(例如)中“收集” Sen 的统计.24
信息,并与相应的文件名一起编译成一个文件。我必须在 R 中做到这一点。
我使用过 CSV 文件,但不确定如何使用文本文件。
这是我现在使用的代码:
require(gtools)
GG <- grep("*.txt", list.files(), value = TRUE)
GG<-mixedsort(GG)
S <- sapply(seq(GG), function(i){
X <- readLines(GG[i])
grep("SEN SLOPE", X, value = TRUE)
})
spl <- unlist(strsplit(S, ".*[^.0-9]"))
SenStat <- as.numeric(spl[nzchar(spl)])
SenStat<-data.frame( SenStat,file = GG)
write.table(SenStat, "sen.csv",sep = ", ",row.names = FALSE)
当前代码无法正确读取所有值并给出此错误:
Warning message:
NAs introduced by coercion
我也没有得到输出的另一列的文件名。请帮忙!
诊断 1
代码也在读取 = 符号。这是 print(spl) 的输出
[1] "" "5.55" "" "-.18" "" "3.08" "" "3.05" "" "1.19" "" "-.32"
[13] "" ".22" "" "-.22" "" ".65" "" "1.64" "" "2.68" "" ".10"
[25] "" ".42" "" "-.44" "" ".49" "" "1.44" "" "=-1.07" "" ".38"
[37] "" ".14" "" "=-2.33" "" "4.76" "" ".45" "" ".02" "" "-.11"
[49] "" "=-2.64" "" "-.63" "" "=-3.44" "" "2.77" "" "2.35" "" "6.29"
[61] "" "1.20" "" "=-1.80" "" "-.63" "" "5.83" "" "6.33" "" "5.42"
[73] "" ".72" "" "-.57" "" "3.52" "" "=-2.44" "" "3.92" "" "1.99"
[85] "" ".77" "" "3.01"
诊断 2
发现了我认为的问题。负号有点棘手。在某些文件中是
SEN SLOPE =-1.07
SEN SLOPE = -.11
由于 = 之后的差距,我得到了第一个的 NA,但代码正在读取第二个。如何修改正则表达式来解决这个问题?谢谢!