Re: st: Ologit storing p-values up to 10 decimal places

[Nick Cox <njcoxstata@gmail.com>]

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 <debs_stata@yahoo.com> 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/