setwd("")
mydata<-read.csv(file="",header=T)
n <- length(unique(mydata$sp))
sp1 <- unique(mydata$sp)
for(i in 1:n ){
name=sp1
c=mydata[mydata$sp==name, ]
plot(c(22,30), c(28,400), type = "n", xlab = "x", ylab = "y")
z<-lm(c$SLA~c$T1,data=c)
abline(z)
}
数据格式
| sp(20个物种) | traits |


雷达卡



京公网安备 11010802022788号







