我刚刚发现了使用Betareg软件包的Beta回归,并且我想看到一个使用predict()并绘制拟合的回归曲线的完整示例,例如betareg小插图文档中第2页第9页(使用汽油产量)来自Prater的数据)。作者的示例(如下)仅提供了回归的摘要。http://cran.r-project.org/web/packages/betareg/vignettes/betareg.pdf。
data("GasolineYield", package = "betareg")
gy_logit <- betareg(yield ~ batch + temp, data = GasolineYield)
summary(gy_logit)
...
使用以下方法ggplot2
:
首先,将模型拟合为数据batch == 6
:
library(betareg)
data("GasolineYield", package = "betareg")
gy_logit <- betareg(yield ~ temp, data = GasolineYield, subset = batch == 6)
gy_loglog <- betareg(yield ~ temp, data = GasolineYield, subset = batch == 6,
link = "loglog")
其次,绘制数据:
library(ggplot2)
ggplot(GasolineYield, aes(x = temp, y = yield)) +
geom_point(size = 4, aes(fill = batch), shape = 21) +
scale_fill_grey() +
geom_line(aes(y = predict(gy_loglog, GasolineYield),
colour = "log-log", linetype = "log-log")) +
geom_line(aes(y = predict(gy_logit, GasolineYield),
colour = "logit", linetype = "logit")) +
scale_colour_manual("", values = c("red", "blue")) +
scale_linetype_manual("", values = c("solid", "dashed")) +
theme_bw()
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句