#naming convention:
#  function.FEATURE/argument.number, e.g.
#  tex.usedims.1

library(vcov)

test = vcov:::test

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)
reg_lm = lm(weight ~ group)
test('vcov.lm.VALUE.1', Vcov(reg_lm), vcov(reg_lm),
	 #stats:::summary.lm uses sqrt(rss/rdf)^2 to calculate vcov;
	 #  this package instead just uses rss/rdf, leading to
	 #  numerical discrepancies between the two
	 approx = TRUE)

counts = c(18,17,15,20,10,20,25,13,12)
outcome = gl(3,1,9)
treatment = gl(3,3)
reg_glm = glm(counts ~ outcome + treatment, family = poisson)
test('vcov.glm.VALUE.1', Vcov(reg_glm), vcov(reg_glm))
