va(formula, data=list(hospital, community), nsymp=16, n.subset=300, method=``quadOpt'', fix=NA, bound=NA, prob.wt=1, boot.se=FALSE, nboot=300, printit=TRUE, print.reg.size=FALSE)