demo/Ex_07_MultiReg.R

data(AMI)
summary(AMI)
mlm1 <- lm(cbind(TOT, AMI) ~ GEN + AMT + PR + DIAP + QRS, data = AMI)
summary(mlm1)

m1 <- lm(TOT ~ GEN + AMT + PR + DIAP + QRS, data =AMI)
summary(m1)
m2 <- lm(AMI ~ GEN + AMT + PR + DIAP + QRS, data = AMI)
summary(m2)

head(resid(mlm1))
head(fitted(mlm1))
coef(mlm1)
sigma(mlm1)
vcov(mlm1)

Try the glmbayes package in your browser

Any scripts or data that you put into this service are public.

glmbayes documentation built on Aug. 5, 2026, 1:07 a.m.