У меня есть данные о годовом количестве осадков с 1086 станций в отдельных текстовых файлах для каждой станции. Имена файлов выглядят как 1.txt, 2.txt, 36.txt, 156.txt, 1098.txt .... в качестве идентификатора станции. Каждая станция имеет данные за разное количество лет (от 48 до 119 лет). Каждый файл данных содержит только значение осадков, без заголовка и информации о году. Пример данных выглядит следующим образом:
198
117
104
114
120
.
.
.
Теперь я хочу выполнить тест тренда Манна Кендалла, используя R для данных каждой станции, и хочу получить вывод в одном кадре данных с идентификатором станции (имя файла), чтобы я мог подготовитьсякарта, показывающая растущие или убывающие тренды.
Любая помощь?
Спасибо!
Привет, Parfait и All, Спасибо за комментарий Parfait. В качестве вашего предложения я попытался написать следующие синтаксисы, но результат содержит только информацию из последних наборов данных.
library(Kendall)
datas <- list.files(pattern="*.txt", full.names=TRUE, recursive=TRUE)
for (i in 1:length(datas)){
data<-read.csv(datas[i],header=F)
mk.test <- MannKendall(data$V1)
values<-data.frame(tau = mk.test$tau[1],p_value = mk.test$sl[1])
}
Можете ли вы помочь сохранить информацию из всех данных и связать их в одном фрейме данных с идентификатором станции?
Заранее спасибо.