[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

From |
"Nick Winter" <nwinter@policystudies.com> |

To |
<statalist@hsphsun2.harvard.edu> |

Subject |
st: RE: help with tables |

Date |
Mon, 19 Aug 2002 11:52:36 -0400 |

> -----Original Message----- > From: Neil Shephard [mailto:mdeasnds@fs1.ser.man.ac.uk] > Sent: Monday, August 19, 2002 9:51 AM > To: statalist@hsphsun2.harvard.edu > Subject: st: help with tables > > > Hi all, > > I was wondering if anyone knows of an .ado file that might solve the > problem I have. > > My dataset consists of LOD scores from whole genome screens, > and I have carried out analysis on various stratifications of the > whole data set. > > I would like to produce a table that gives the maximum LOD score > on each chromosome for each replicate, with the position of that > LOD score. I can achieve the first part of this using > > -table chr strat, c(max lod)- > > which produces > > ------------------------------------------ > Chrom | strat > | 0 1 2 > -------------+------------------------------- > Chr01 | 1.021457 .744813 1.516707 > Chr02 | .484565 .526793 .633966 > Chr03 | .582743 1.079223 .228413 > Chr04 | .515382 1.269104 .107555 > Chr05 | .290357 .15492 .935334 > Chr06 | 3.610746 6.724064 .210437 > Chr07 | .658136 .859743 .266631 > Chr08 | .460963 .786725 .294109 > Chr09 | .47856 .520341 .121645 > Chr10 | .673608 .62035 1.421104 > Chr11 | .100738 .555981 .933838 > . > . > > > however I can't work out how to include the pos(ition) that > corresponds to each of the maximum LOD scores. None of the > clist arguments seem to allow for this. > > I've tried a few different -findit- searches but nothing > obvious stuck > out, I'd be grateful for any suggestions. > > Thanks in advance > > Neil There may be a fancier way, but one way is to first flag the relevant observation(s), then table only those: egen double ismax = max(lod) , by(chr strat) replace ismax = lod==ismax table chr strat if ismax==1, c(mean lod mean pos n pos) The first line simply generates a variable with the maximum lod within chr and strat; then the next line converts this to a 1/0 variable that flags the maximum. The final line tables only the maximum observations, listing the mean lod and mean pos. Note the inclusion also of the number of maximum observations; this will tip you off if there is more than one that takes on the maximum value, in which case you might want to examine more closely. Nick Winter * * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

- Prev by Date:
**RE: st: Generating a logistically distributed random variable.** - Next by Date:
**st: First announcement and call for papers, NASUG'2003** - Previous by thread:
**st: help with tables** - Next by thread:
**st: First announcement and call for papers, NASUG'2003** - Index(es):

© Copyright 1996–2017 StataCorp LLC | Terms of use | Privacy | Contact us | What's new | Site index |