On Wed, Aug 26, 2009 at 5:39 AM, Jane
Maguire<[email protected]> wrote:
> 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.
What does "rubbish" mean?
> I am not a biostatistician- am a Phd student in medicine and Public Health.
You'd benefit from getting your head round the basics of statistics in
the long run. A good book is Intuitive Biostatistics by Harvey
Motulsky. Not too heavy on the maths, but explains things very well.
> 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
As per above, what does "rubbish" mean? If you could show the output
of running these two commands it would be more informative for those
on the list who might then be able to comment on what the problem is.
Neil
--
"The combination of some data and an aching desire for an answer does
not ensure that a reasonable answer can be extracted from a given body
of data." ~ John Tukey (1986), "Sunset salvo". The American
Statistician 40(1).
Email - [email protected]
Website - http://slack.ser.man.ac.uk/
Photos - http://www.flickr.com/photos/slackline/
*
* 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/
