Displaying 1 result from an estimated 1 matches for "smokgp".
Did you mean:
smoke
2010 Feb 12
1
"drop if missing" command?
...'s "drop1" with the following code:
##geno is the dataset
## the dependent variable (casectrln) is dichotomous and coded 0,1
## rs743572_2 is one of the snps (which is coded 0,1,2 for the 3 genotypes)
library(stats)
modadd = glm(geno$casectrln ~rs743572_2 + factor(racegrp)+ factor(smokgp)+
factor(alcgp)+ factor(bmigp) + factor(ipssgp)
+ agebase, family="binomial")
drop1(mod,scale=0,test=c("Chisq"), x=NULL, k=2)
There is a great deal of missing data in this dataset for both snps and for
covariates--so I have been instructed not to simply drop all cases with
miss...