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

From |
"Martin Weiss" <martin.weiss1@gmx.de> |

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

Subject |
st: AW: How to write a fast code to obtain a specific calculation |

Date |
Wed, 18 Feb 2009 13:56:45 +0100 |

<> If it is all about speed, ************* qui summ WT, mean ************* is the way to go. Contrary to its name, it does leave a few things in its wake, among them -r(sum)- HTH Martin -----Ursprüngliche Nachricht----- Von: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] Im Auftrag von Tiago V. Pereira Gesendet: Mittwoch, 18. Februar 2009 00:53 An: statalist@hsphsun2.harvard.edu Betreff: st: How to write a fast code to obtain a specific calculation Dear statalisters, I usually need to compute a value X that satifies a specific condition. A typical example is shown below. But this is only a single example, and ordinarily I have similar, but more complex aims. Hence, I would like to know your expert opinion on how one can make codes these kind of codes faster. Perhaps only Mata helps? In this example, I want to compute the value of tau^2 that imposes the lower limit of my confidence interval to be 0. . . . scalar lower_limit = 999999999 while lower_limit>0 { local tau2 = `tau2'+0.0001 cap drop T_i V_i W_i WT gene T_i= ln(_ES) gene V_i = _selogES^2 gene W_i = 1/(V_i+`tau2') gene WT = W_i*T_i qui summ WT scalar sum_WT = r(sum) qui sum W_i scalar sum_W_i = r(sum) scalar summary_random = sum_WT/sum_W_i scalar lower_limit = summary_random-(1.96/sqrt(sum_W_i)) } Cheers! Tiago * * 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/

**References**:**st: How to write a fast code to obtain a specific calculation***From:*"Tiago V. Pereira" <tiago.pereira@incor.usp.br>

- Prev by Date:
**st: How to write a fast code to obtain a specific calculation** - Next by Date:
**st: multivariate normal regression with ml, linear form** - Previous by thread:
**st: How to write a fast code to obtain a specific calculation** - Next by thread:
**st: RE: How to write a fast code to obtain a specific calculation** - Index(es):

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