Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: st: Ologit storing p-values up to 10 decimal places
From
Nick Cox <[email protected]>
To
[email protected]
Subject
Re: st: Ologit storing p-values up to 10 decimal places
Date
Mon, 8 Aug 2011 19:19:16 +0100
Your code is very confused, mixing up macros and variables and
assuming that values exist which you never assigned. However, with the
big caveat that no P-value can be trusted out to the tenth decimal
place, the technique in this may be a start at what you seek.
webuse fullauto
foreach var of var length-gratio {
capture ologit rep78 `var'
di "`var' {col 20}" %12.10f e(p)
}
For more manipulation, look at -postfile-.
If you are inspecting that many P-values, also consider building in a
penalty for chance hits.
On Mon, Aug 8, 2011 at 6:59 PM, Debs Majumdar <[email protected]> wrote:
> I am trying to write a code which will go through around 500 ordinal logistic regressions and the output is the id for each of the rs # and the corresponding p-values, where the p-values are sorted in ascending order. i want the p-values reported up to 10 decimal places. This is my code so far. I may be wrong here. How do I make Stata report p-values up to 10 decimal places? Do I need to specify e(b) under display if I want the beta coefficients printed out for each rs #?
>
>
> ~~~~~~~~
> local list "age sex c1 c2 c3"
> local pvalue = e(p)
>
> foreach `x’ of varlist rs1-rs500 {
>
> quietly ologit outcome `x' `list'
>
> sort pvalue
> display `x' pvalue
> }
>
> ~~~~~~~~
>
> I am using Stata 11.1.
>
> I want the output to show the following:
>
> ~~~~
>
> rs5 5e-7
> rs493 4e-6
> ...
> ...
> rs257 0.33
>
> ~~~
*
* 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/