处理R中的过函数：如何绘制“A”的时间趋势；超过；对象

我有一个名为

屠体

的数据集，其中包含了过去5年在挪威出现的牲畜屠体（每个生肉对应一个屠体，它有

ID号

，以及

经度

、

纬度

和年份）。我还有一个名为regions的shapefile，它是挪威不同地区的矢量文件（每个地区都有一个ID号，在

number\u 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")

现在，我需要为每个区域绘制一段时间内屠体发生的情况，即，我想要一张图，其中区域1的屠体数量在y轴上，时间（年）在x轴上，区域2的屠体发生率相同，区域3的屠体发生率相同，依此类推

然而，我遇到了麻烦，因为对象x的空间位置的R“（…）中的over函数从空间对象y检索索引或属性”，因此在我的res对象中，我没有包含列year（因为year在carcas数据集中，而不是regions数据集中）

res<-table(over(carcass, regions)$number_reg

plot(table(carcass$YEAR), type="o", ylab="Carcass occurrences", main="Temporal trend of carcass occurrences in Norway")