Dear all
Data is case control, snps are coded 0,1,2
One gene, two snps. Haplotype for hap00 is quite rare and when interaction term (cox2use) is added the output is rubbish. I am not a biostatistician- am a Phd student in medicine and Public Health. However, I think the problem is the numbers are too small...unless I have made a coding error? any advice?
This command works- without interaction,
xi: haplologit stroke i.cox2use i.aspirin i.nsaid i.gender i.agegroup2 i.bmigroup2 i.Diabetes4 i.HT i.AF2 ///
(i.Smoker) if ~(inlist(kozak,.) & inlist(thrmet,.)), snpvars(kozak thrmet) inher(d) ///
riskhap1("11") ///
happrefix(_hap) or noemtable noemshow
But this one doesn't,
set more off
xi: haplologit stroke i.cox2use i.aspirin i.nsaid i.gender i.agegroup2 i.bmigroup2 i.Diabetes4 i.HT i.AF2 ///
(i.Smoker) if ~(inlist(kozak,.) & inlist(thrmet,.)), snpvars(kozak thrmet) ///
riskhap1("11", inter(cox2use)) ///
inher(d) happrefix(_hap) or
Many Thanks
Jane
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
