data(treesf)
## Not run:
# Y <- log(treesf$Volume)
# X1 <- log(treesf$Diameter / 12)
# X2 <- log(treesf$Height)
# trees.lm <- lm(Y ~ X1 + X2)
# summary(trees.lm)
# log(pi/12)
# vcov(trees.lm)## End(Not run)
Run the code above in your browser using DataLab