我有一个大型数据框，每个plotID包含多个物种代码以及相关的变量来映射。
plotID<-c(rep(c(1),10),rep(c(2),5),rep(c(3),4))
species<-c(1:10,1:5,1:4)
biomass<-sample(1:100,19)
biomass2<-sample(1:100,19)
lat<-c(rep(c(36.7),10),rep(c(38.6),5),rep(c(39.5),4))
long<-c(rep(c(-89.3),10),rep(c(-89.9),5),rep(c(-90.3),4))
df<-data.frame(plotID,species,biomass,biomass2,lat,long)
我想绘制生物量变量的各个物种。因此在这个例子df中，我想绘制总共10种物质图。
我使用split将数据帧拆分为按物种划分的数据帧列表，但是spplot不适用于列表。
请告诉我。
我尝试如下：
library(sp)
library(lattice)
library(latticeExtra)
coordinates(df)<-c("lat","long")
split_df<-split(df,df$species)
spplot(split_df[c("biomass","biomass2")])