我有一个名为 的数据集carcass，其中包含过去 5 年中挪威出现的牲畜尸体（每个原始尸体对应一个尸体，其中包含、ID number以及和年份）。我还有一个名为region的 shapefile ，它是挪威不同地区的矢量文件（每个地区也有一个 ID 号，在 column 中指定）。longitudelatitudenumber_reg
我设法在 R 中绘制了区域地图，并在地图上绘制了与每个胴体相对应的点。
regions<-readShapeSpatial("\\\\homer.uit.no/fma023/Desktop/Filippo Marolla/Datasets/rbd2012_siida.shp")
plot(regions)
carcass<-readShapePoints("\\\\homer.uit.no/fma023/Desktop/Filippo Marolla/Datasets/carcass_data_edit.shp")
plot(carcass, add=T, pch=20, cex=.5, col="red")
然后，我设法使用 over 函数计算了每个区域内有多少点：
res<-table(over(carcass, regions)$number_reg
我还设法绘制了整个国家尸体发生的时间趋势：
plot(table(carcass$YEAR), type="o", ylab="Carcass occurrences", main="Temporal trend of carcass occurrences in Norway")
现在我需要绘制每个区域随时间推移的胴体发生情况，即我想要一个图表，其中 y 轴上的胴体数量和区域 1 的 x 轴上的时间（年份），区域 2 的相同图表，区域的相同图表3，以此类推。
但是，我遇到了麻烦，因为 R 中的over函数“（...）在对象 x 的空间位置从空间对象 y 中检索索引或属性”，因此在我的res对象中，我没有包含year列（因为年份在carcass数据集中，而不是在区域数据集中）。