我正在尝试通过循环创建点模式ppp.object {spatstat}的列表。 我的数据集如下:
> names(OT1);head(OT1);dim(OT1)
[1] "EID" "latitude" "longitude" "month" "year" "CPUE" "TSUM"
[8] "fTSUM"
EID latitude longitude month year CPUE TSUM fTSUM
1 167-1-1996-1135 67.70000 -61.81667 9 1996 0 0 F
2 167-10-1996-1135 67.71667 -59.18333 9 1996 0 0 F
3 167-100-1996-1135 67.86667 -59.43333 10 1996 0 0 F
4 167-101-1996-1135 67.95000 -59.58333 10 1996 0 0 F
5 167-102-1996-1135 68.10000 -59.76667 10 1996 0 0 F
6 167-103-1996-1135 67.81667 -59.38333 10 1996 0 0 F
[1] 2707 8
我想做的是选择每个月的数据并创建一个 ppp.object。
> sort(unique(OT1$month))
[1] 7 8 9 10 11 12
以下循环可以正常运行,并且我可以看到每一个图像:
for(i in sort(unique(OT1$month))){
a<-OT1[OT1$month==i,]
b<-ppp(a$longitude,a$latitude,marks=a$fTSUM,window=newW)
plot(b,main=i)
}
我想创建一个所有ppp.object的列表,方便我单独访问。我已经尝试在循环命令中添加list(),但没有成功...任何帮助将不胜感激!
谢谢!