過去 5 年間のノルウェーでの家畜の死骸の発生を含むというデータセットがあります (各 raw は、 、、および年carcassを持つ単一の死体に対応します)。ノルウェーのさまざまな地域のベクター ファイルであるregionというシェープファイルもあります (すべての地域には、列で指定された ID 番号もあります)。ID numberlongitudelatitudenumber_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")
今、私はすべての単一地域の経時的な枝肉の発生をプロットする必要があります。つまり、地域1のy軸に枝肉の数、x軸に時間（年）のグラフ、地域2の同じグラフ、地域の同じグラフが必要です3など。
ただし、R のover関数「(...) オブジェクト x の空間位置で、空間オブジェクト y からインデックスまたは属性を取得する」ため、resオブジェクトに年列が含まれていないため、問題が発生しています。 (年は地域のデータセットではなく、枝肉のデータセットにあるため)。