要求每个物种的线性回归,并把这些回归线放在一张图中,写了代码,但每次都只有最后一个物种的回归线
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)
}
数据格式