[英]Best fit quadratic regression
我遇到一个奇怪的问题; 在这里得到我的数据集: 数据集
我只需要一个简单的图表,显示rao和obs_richness之间的最佳拟合回归(二次回归); 但是相反,我得到了截然不同的多项式模型。 对于如何解决这个问题,有任何的建议吗?
#read in data
F_Div<-read.csv('F_Div.csv', header=T)
str(F_Div)
pairs(F_Div[2:12], pch=16)
#richness vs functional diversity
par(mfrow=c(1,1))
lm1<-lm ( rao~Obs_Richness, data=F_Div)
summary (lm1)
plot (rao~Obs_Richness, data=F_Div, pch=16, xlab="Species Richness", ylab="Rao's Q")
abline(lm1, lty=3)
lines (lowess (F_Div$rao~F_Div$Obs_Richness))
poly.mod<- lm (F_Div$rao ~ poly (F_Div$Obs_Richness, 2, raw=T))
summary (poly.mod)
lines (F_Div$Obs_Richness, predict(poly.mod))
我需要最接近最低线(一条简单的曲线)的线,而不是这弯弯曲曲的曲线。
我也试过了,但没有什么需要:
xx <- seq(0,30, length=67)
plot (rao~Obs_Richness, data=F_Div, pch=16, xlab="Species Richness", ylab="Rao's Q")
lines(xx, predict(poly.mod, data.frame(x=xx)), col="blue")
因为line(...)
以数据的原始顺序在连续点之间绘制线,所以造成了混乱的混乱。 最后尝试一下。
p <- data.frame(x=F_Div$Obs_Richness,y=predict(poly.mod))
p <- p[order(p$x),]
lines(p)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.