Question

我正在尝试绘制模型的点和置信区间。我已经掌握了要点，但是我对c.i的了解却很迷失。我应该对信心区间使用polygon（）。我在发布之前查看了建议的欺骗，但由于我是r的新手，这有点令人困惑。我尝试在该解决方案之后添加多边形，但是它不起作用。这是我的位置：

head(bill)

   month year kwh days est cost  avgT dT.yr kWhd.1
93     4 2011 294   28   a 37.19   39   -13 10.50000
57     4 2008 367   29   a 32.01   39    -2 12.65517
73     8 2009 415   30   a 39.41   71    -3 13.83333
95     6 2011 461   33   a 52.19   66     1 13.96970
83     6 2010 502   33   a 54.12   65     4 15.21212
99    10 2011 416   27   a 48.03   64     0 15.40741

bill <- bill[order(bill$avgT), ]

lm1 <- lm(kWhd.1~avgT+I(avgT^2), data=bill)
summary(lm1)

ci <- predict(lm1, bill, interval="confidence")

par(mar=c(4.3,4.3,1,1))
with(bill, plot(avgT, kWhd.1, xlab="Average Temperature (F)", ylab="Daily Energy Use (KWH/d)", type="n", las=1))

newx <- seq(min(bill$avgT), max(bill$avgT), length.out=101)

polygon(c(rev(newx), newx), c(rev(ci[ ,2]), ci[ ,3]), col="lightblue1", border = NA)

with(bill, curve(lm1$coef[1] + lm1$coef[2]*x + lm1$coef[3]*x^2, from=min(avgT), to=max(avgT), add=T))

在基数r中绘制置信区间

0 个答案: