Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: AW: Formatting tables with estpost & esttab


From   Eduardo Nunez <[email protected]>
To   [email protected]
Subject   Re: st: AW: Formatting tables with estpost & esttab
Date   Mon, 5 Oct 2009 11:10:19 -0400

Hi Martin,

Stata/MP 11.0 for Windows (64-bit x86-64)
Born 26 Aug 2009
Copyright (C) 1985-2009

Total physical memory:     4191812 KB
Available physical memory: 2227852 KB

Single-user 2-core Stata perpetual license:
       Serial number:  XXXXXXXXXXXX
         Licensed to:  Eduardo Nunez
                       Personal license


I tried to adapt the code for means-SD and ttest to medians IQR and
median unsuscefully. Below is the code. I get the p-value for each
variable but not the descriptive statistics.
What am i doing wrong?

Here is the code:

**********************
use "C:\Users\Eduardo\Desktop\HF dataset - August 24,
2009\Longitudinal_lymphocyte_analysis.dta", clear
tempname hdle
capt erase info.dta
postfile `hdle' str10 var Median_0 IQR25_0 IQR75_0 /*
*/  Median_1 IQR25_1 IQR75_1 pval Overallmedian OverallIQR25 OverallIQR75 /*
*/ using info

foreach var of varlist los edad_b fc_ingreso tas_ingreso tad hb leucos
neutrofilos linfocitos monocitos cr acurico colesterol ldl hdl tga na{
       sum `var' if !sexo
       loc Median_0=r(p50)
       loc IQR25_0=r(p25)
	   loc IQR75_0=r(p75)
       sum `var' if sexo
       loc Median_1=r(p50)
       loc IQR25_1=r(p25)
	   loc IQR75_1=r(p75)
       median `var', by(sexo) exact
       local pval=r(p_exact)
       sum `var'
       post `hdle' ("`var'") (`Median_0') (`IQR25_0') (`IQR75_0') /*
       */ (`Median_1') (`IQR25_1') (`IQR75_1') (`pval') /*
       */ (r(p50)) (r(p25)) (r(p75))
}

postclose `hdle'

use info, clear
l, abbrev(12) noobs sep(0)
*********************************************************
Thisw is the result I got:
+----------------------------------------------------------------------------------------------------------------------------------+
  |        var   Median_0   IQR25_0   IQR75_0   Median_1   IQR25_1
IQR75_1       pval   Overallmed~n   OverallIQR25   OverallIQR75 |
  |----------------------------------------------------------------------------------------------------------------------------------|
  |        los          .         .         .          .         .
    .   .4098835              .              .              . |
  |     edad_b          .         .         .          .         .
    .   .0000138              .              .              . |
  | fc_ingreso          .         .         .          .         .
    .   .6778672              .              .              . |
  | tas_ingres          .         .         .          .         .
    .    .000103              .              .              . |
  |        tad          .         .         .          .         .
    .   .4437288              .              .              . |
  |         hb          .         .         .          .         .
    .   7.24e-11              .              .              . |
  |     leucos          .         .         .          .         .
    .   .0460077              .              .              . |
  | neutrofilo          .         .         .          .         .
    .     .00087              .              .              . |
  | linfocitos          .         .         .          .         .
    .    .833713              .              .              . |
  |  monocitos          .         .         .          .         .
    .   1.11e-06              .              .              . |
  |         cr          .         .         .          .         .
    .   1.36e-08              .              .              . |
  |    acurico          .         .         .          .         .
    .   .0039871              .              .              . |
  | colesterol          .         .         .          .         .
    .   7.56e-08              .              .              . |
  |        ldl          .         .         .          .         .
    .   .0002156              .              .              . |
  |        hdl          .         .         .          .         .
    .   .0047893              .              .              . |
  |        tga          .         .         .          .         .
    .   1.35e-09              .              .              . |
  |         na          .         .         .          .         .
    .   .6747179              .              .              . |
  +----------------------------------------------------------------------------------------------------------------------------------+


One more thing...
The code you sent me for discrete variables:
***
sysuse auto, clear

estpost tab rep for, chi2
esttab ., main(b) aux(colpct)  /*
*/ nostar unstack noobs nonum nonote  /*
*/ sca("chi2 Chi2 Value" "p corr. p-value")  /*
*/ sfmt(%4.3fc) wide compress

***
How do I loop over several variables? I used tab2, firstonly but it
seems estpost doesn't recognize it.



Note: In biomedical sciences these are the traditional way to express
table 1 for the manuscript where you show how baseline characteristics
are distributed over your exposure variable.
I imagine a lot of statalister that work in biomedical science would
take advantage of this automatic process to save time.

I am very grateful of your help and excuse my ignorance in Stata programing.

Best wishes,

Eduardo





2009/10/5 Martin Weiss <[email protected]>:
>
> <>
>
> Not sure how you would get the upper panel into esttab, but you can sure get
> it into a Stata file which you could export in any way you see fit...
>
>
> *************
> sysuse auto, clear
>
> tempname hdle
>
> capt erase info.dta
>
> postfile `hdle' str10 var dommean domsd /*
> */ formean forsd pval overallmean overallsd /*
> */ using info
>
>
> foreach var of varlist price mpg{
>        sum `var' if !for
>        loc dommean=r(mean)
>        loc domsd=r(sd)
>        sum `var' if for
>        loc formean=r(mean)
>        loc forsd=r(sd)
>        ttest `var', by(foreign)
>        local pval=r(p)
>        sum `var'
>        post `hdle' ("`var'") (`dommean') (`domsd') /*
>        */ (`formean') (`forsd') (`pval') /*
>        */ (r(mean)) (r(sd))
> }
>
> postclose `hdle'
>
> use info, clear
> l, abbrev(12) noobs sep(0)
> *************
>
>
>
> HTH
> Martin
>
> -----Ursprüngliche Nachricht-----
> Von: [email protected]
> [mailto:[email protected]] Im Auftrag von Eduardo Nunez
> Gesendet: Montag, 5. Oktober 2009 01:18
> An: [email protected]
> Betreff: Re: st: AW: Formatting tables with estpost & esttab
>
> Hi Martin,
>
> It did work pretty good for discrete variables. For continuous
> variables, I would like a code that provide me side by side the
> meanpuls-minusSD for the two categories of Foreign and a third column
> with the p-value from ttest.
>
> Thank you so much in advance,
>
> Eduardo
>
>
>
> On Sun, Oct 4, 2009 at 5:36 PM, Martin Weiss <[email protected]> wrote:
>>
>> <>
>>
>>
>>
>> You can -estpost- all these results in a quite straightforward manner, and
>> then -esttab- them with the matrix function of this command.  For
> instance,
>> the lower panel of your table would be
>>
>> ***
>> sysuse auto, clear
>>
>> estpost tab rep for, chi2
>> esttab ., main(b) aux(colpct)  /*
>> */ nostar unstack noobs nonum nonote  /*
>> */ sca("chi2 Chi2 Value" "p corr. p-value")  /*
>> */ sfmt(%4.3fc) wide compress
>>
>> ***
>>
>>
>>
>> HTH
>> Martin
>>
>>
>> -----Ursprüngliche Nachricht-----
>> Von: [email protected]
>> [mailto:[email protected]] Im Auftrag von Eduardo Nunez
>> Gesendet: Sonntag, 4. Oktober 2009 22:33
>> An: [email protected]
>> Betreff: st: Formatting tables with estpost & esttab
>>
>> Dear statalister,
>>
>> I have been trying to get the right code in estpost and esttab
>> (without success) to format tables as the example below (using auto
>> dataset):
>>
>>                         Domestic (N=52)                Foreign (N=22)
>> P-value         Total (N=74)
>>
>> Price                       6179.3ą3189.0               6070.1ą2221.0
>> 0.015           6146.0ą2912.4
>>
>> Mileage (mpg)                  19.54ą4.753              25.29ą6.310
>> 0.006           21.29ą5.866
>>
>> Repair Record 1978, n(%)
>> *0.000
>> 1                                2 (4.17)
>> 0 (0)
>>        2 (2.9)
>> 2                               8 (16.67)
>> 0 (0)
>>        8 (11.59)
>> 3                               27 (56.25)
>> 3 (14.29)                               30 (43.48)
>> 4                               9 (18.75)
>> 9 (42.86)                               18 (26.09)
>> 5                               2 (4.17)
>> 9 (42.86)                               11 (15.94)
>>
>> *Omnibus P-value
>>
>>
>>
>>
>> The meanąSD  [or median (IQR) when appropriate] for each variable and
>> for the total population with the p-value corresponding (either a
>> t-test or ranksum test);
>> For discrete variables, the n (%)...count(col percentage) and the
>> p-value corresponding to chi2.
>>
>> I would really appreciate any help on this issue.
>>
>> Best wishes,
>>
>> Eduardo
>>
>> *
>> *   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/
>>
>>
>> *
>> *   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/
>>
>
> *
> *   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/
>
>
> *
> *   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/
>

*
*   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/



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