从一个数据帧创建时间序列的子集，以便按站点运行SeasonalMannKendall趋势测试

Question

我有大约13个不同电台的以下数据：

> dput(RawData[1:10,])
structure(list(Station = c(469, 469, 469, 469, 469, 469, 469, 
469, 469, 469), Classification = structure(c(1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L), .Label = c("  Approved  ", "  Conditionally Approved  ", 
"  Prohibited  "), class = "factor"), SampleDate = structure(c(8504, 
8504, 8505, 8505, 8532, 8532, 8533, 8533, 8561, 8561), class = "Date"), 
Year = c(1993, 1993, 1993, 1993, 1993, 1993, 1993, 1993, 
1993, 1993), SWTemp = c(10, 11, 10, 10, 14, 15, 12, 14, 15, 
16), Salinity = c(26, 28, 28, 30, NA, NA, 30, 30, 28, 18), 
FecalColiform = c(1.8, 2, 2, 1.8, 2, 2, 1.8, 1.8, 4.5, 2)), .Names = c("Station", 
"Classification", "SampleDate", "Year", "SWTemp", "Salinity", 
"FecalColiform"), row.names = c(NA, 10L), class = "data.frame")

我想对每个站的粪大肠菌群数据分别运行一个SeasonalMannKendal。 我知道这必须是一个时间序列。 如何使每个电台进入自己的时间序列，以便可以运行这些测试？

我试图重塑数据以按采样日期列出测站结果，但这会为某些日期创建NA，因此我也无法以这种方式运行测试。

我最好的方法是什么？

先感谢您！

Answer 1

您可以通过几种方式来执行此操作，但是这种方式将为您提供一个数据帧列表，其中每个数据帧仅具有一个Station的数据。 我还添加了as.ts以使用TSA库将数据转换为时间序列。 您可能需要对日期进行调整，才能获得所需的日期。

library("TSA")
mylist = list()
for(i in 1:unique(RawData$Station)){
  mylist[[i]] = as.ts(RawData[RawData$Station == i,])
}

names(mylist) = unique(RawData$Station)

现在，您可以通过执行mylist[["429"]]类的操作来选择所需的站点。

如果您只想要不在列表中的单独数据框，则可以执行以下操作：

Station1 = RawData[RawData$Station == 1,]
Station2 = RawData[RawData$Station == 2,]

...等等。

PS，将来使用dput(yourdata)会很有帮助，以便某人可以轻松地复制粘贴您的问题并进行处理。 不过，感谢您给出一个明确的例子！