我需要将响应变量chl,SST和upwelling绘制为解释变量的函数 单个数字上的纬度。因为这些响应变量有不同的范围,我 将不得不在一个单独的面板上绘制每一个。 这是数据:
d1< - read.csv(“http://www.northeastern.edu/synchrony/stats/data/assn4_dataset1.csv”)
此外,我需要使用lm函数来拟合将chl与SST,上升流相关联的多元回归模型, 和纬度(来自相同的数据)。
最后,拟合一个将chl与SST和上升流相关联的多元回归。
到目前为止绘制我拥有的变量:
情节(表格(d1 $ chl,d1 $ SST)
lm.fit< - (d1)
我是初学者而且非常困惑
答案 0 :(得分:0)
par(mfrow=c(2,2)) # Split plot window in 2 row and 2 col
plot(chl ~ latitude, data=d1, type="l") # type="p" for points
plot(sst ~ latitude, data=d1, type="l") # type="b" for points and line
plot(upwelling ~ latitude, data=d1, type="l")
myreg <- lm(upwelling ~ latitude, data=d1)
abline(myreg, col="blue")
summary(myreg) # Sum up results
# See also predict(), coefficients(), confint(), ... ...
plot(myreg) # To plot residuals etc...