va.gcv(formula, data=list(hospital, community), nsymp.vec, n.subset=300, prob.wt=1, boot.se=FALSE, nboot=1, printit=FALSE, print.reg.size=TRUE)