options(digits = 4)
data("GasolineYield", package = "betareg")
gy <- betareg(yield ~ gravity + pressure + temp10 + temp, data = GasolineYield)
gy_res <- cbind(
"quantile" = residuals(gy, type = "quantile"),
"pearson" = residuals(gy, type = "pearson"),
"deviance" = residuals(gy, type = "deviance"),
"response" = residuals(gy, type = "response"),
"weighted" = residuals(gy, type = "weighted"),
"sweighted" = residuals(gy, type = "sweighted"),
"sweighted2" = residuals(gy, type = "sweighted2")
)
pairs(gy_res)
cor(gy_res)
Run the code above in your browser using DataLab