※上記の広告は60日以上更新のないWIKIに表示されています。更新することで広告が下部へ移動します。

plot(dm$year,dm$adductor, ylim=c(0,30), pch=21,xlab="", ylab="8月~9月の平均貝柱重量(g)",type="b",lwd=2,cex=1.5,bg="black",axes=FALSE)
box()
axis(1,at=c(1992:2012), labels=c(1992:2012))
axis(2,at=c(0,10,20,30), labels=c(0,10,20,30))
legend("topright", pch=21,pt.bg=c("black", "white"),  cex=1, lwd=2, lty=c(1,2), legend=c("実測値", "予測値"))
points(dm$year,glm.res$fitted.values,ylim=c(0,30), pch=21, type="b",cex=1.5,lwd=2, lty=2)
|