search for: smokgp

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...