ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
# requires lm.beta
if (requireNamespace("lm.beta", quietly=TRUE))
{
lm.D9 <- lm.beta::lm.beta(lm(weight ~ group))
pprint(summary(lm.D9),
format = "summary.lm.beta.coeff",
coeff = 2L)
}
Run the code above in your browser using DataLab