我有数百个文本文件,每个文件中都包含以下信息:
*****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
由于等号后面有空格,第一个得到了NAs,但是代码读取了第二个。我该如何修改正则表达式以修复这个问题?谢谢!