Statalist


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

R: st: a possible - foreach - question?


From   "Carlo Lazzaro" <carlo.lazzaro@tiscalinet.it>
To   <statalist@hsphsun2.harvard.edu>
Subject   R: st: a possible - foreach - question?
Date   Tue, 14 Oct 2008 17:18:06 +0200

Dear Austin,
many thanks for your helpful hints. I have just modified your syntax (too
advanced for me) in order to approach my ultimate goal, that is calculating,
for each one of the 10,000 iterations, how many times ICER_round*Delta_Q was
> Delta_C (please, see below). I was not absolutely clear on this point in
my previous messages as it was a subsequent thought of mine.
Once again I would like to express my grateful thankfulness  to Austin, Ben,
Johannes, and Martin (in alphabetical order) for their support and for
having increased my poor skills about Stata programming. This experience
makes me think that the famous Bayesian saying "today's posterior will be
tomorrow's prior" applies very well when it comes to knowledge improvement
that can be obtained from more experienced Statalisters.

Thanks a lot again and Kind Regards,
Carlo
  

g double d=.
(10000 missing values generated)

. edit
- preserve

. qui forv i=1/`=_N' {
.  if  ICER[`i']<. {
. g double tmp=  ICER_round[`i']* Delta_Q> Delta_C
.   su tmp, meanonly
. replace d=r(mean) in `i'
.   drop tmp
.   }
.   }
-----Messaggio originale-----
Da: owner-statalist@hsphsun2.harvard.edu
[mailto:owner-statalist@hsphsun2.harvard.edu] Per conto di Austin Nichols
Inviato: martedì 14 ottobre 2008 15.31
A: statalist@hsphsun2.harvard.edu
Oggetto: Re: st: a possible - foreach - question?

Carlo Lazzaro <carlo.lazzaro@tiscalinet.it>
Not sure I've understood, but if so, here is one method of calculation.

use bp* tcresult using http://www.stata-press.com/data/r9/nhanes2, clear
ren bpsystol a
ren bpdiast b
ren tcresult c
g double d=.
qui forv i=1/`=_N' {
 if a[`i']<. {
  g double tmp=a[`i']*b-c
  su tmp, meanonly
  replace d=r(mean) in `i'
  drop tmp
  }
 }

On Tue, Oct 14, 2008 at 4:50 AM, Carlo Lazzaro
<carlo.lazzaro@tiscalinet.it> wrote:
> Dear Statalisters,
> I have three variables (a; b; C with 10,000 observations each) in long
> format.
> For each value of a I would like to obtain the following one:
>
> (a1 * b1) - c1
> .
> .
> .
> .
> (a1*b10,000)-c10,000;
>
> Repeat the above reported would-be 10,000 times (i.e., from a1 to a10,000)
>
> Create a new variable d in which the mean of each one of the 10,000
> iterations are stored and can be displayed.
>
> I suppose there's way of doing this with Stata 9.2/SE, but I cannot figure
> out the proper syntax.
*
*   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–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index