<>
*************
clear*
input Month s1 u1 s2 u2
s3 u3
1 .0198312 .0016526 .0189071 .00157559 .0386583
.00322152
2 .01853398 .0015445 .01767033 .00147253 .03612954
.00301079
3 .01732161 .00144347 .01651445 .0013762 .03376619
.00281385
end
compress
list, noobs
egen float monthmeanui =/*
*/ rowmean(u1 u2 u3)
forv i=1/3{
egen totals`i'=total(s`i')
}
list Month monthmeanui/*
*/ totals1 totals2 totals3,/*
*/ noobs
*************
HTH
Martin
-----Ursprüngliche Nachricht-----
Von: [email protected]
[mailto:[email protected]] Im Auftrag von Carlo Lazzaro
Gesendet: Montag, 22. Juni 2009 15:00
An: [email protected]
Betreff: st: R: RE: rowmean within a loop
Dear Martin and Nick,
thanks a lot for your kind reply and sorry for not being clear in my
previous message.
I am doing a probabilistic sensitivity analysis (PSA) with 10,000
simulations for each variable.
I would like to obtain 10,000 s`i' simulations (that is, 10,000 s`i',
expressing the monthly probability of surviving for people affected by a
given disease and followed-up from 50 to 100 years of age);
Then I multiply this 10,000 probabilities for the surviving expressed in a
monthly base (LYS), in order to get 10,000 LYS weighted for s`i'(ie,
s`i'*LYS) (where LYS =.0833333).
My problem is: I would like to collect, for each month, the rowmean value of
u`i', calculated across all the 10,000 simulations.
Assuming for sake of simplicity that my PSA stops after 3 simulations, I
would figure out a file.dta like the following one:
Month s1 u1 s2 u2
s3 u3
1 .0198312 .0016526 .0189071 .00157559
.0386583 .00322152
2 .01853398 .0015445 .01767033 .00147253
.03612954 .00301079
3 .01732161 .00144347 .01651445 .0013762
.03376619 .00281385
By the way, for another research purpose, I am interested in storing in
another variable, the total() per column of each u`i'.
Thanks a lot for your kindness and for your time.
Kind Regards,
Carlo
-----Messaggio originale-----
Da: [email protected]
[mailto:[email protected]] Per conto di Nick Cox
Inviato: lunedì 22 giugno 2009 13.03
A: [email protected]
Oggetto: st: RE: rowmean within a loop
The intent of the code you give is difficult to understand. It commits
you to generating 20000 new variables, s1-s10000 and u1-u10000, which I
doubt is what you really want.
egen v`i' = rowmean(u`i')
within a loop is just going to put the rowmean of one variable (that
variable, identically) in another variable, which is pointless.
Do start again from the beginning:
1. Explain your data.
2. Explain your problem.
Nick
[email protected]
Carlo Lazzaro
after typing
qui forval i =1/10000 {
qui g double s`i' = Surv_Pos_IFI[`i']*(1-Monthly_Pr_Death_AML_MDS_Prob)
in 1
replace s`i' = s`i'[_n-1]*(1-Monthly_Pr_Death_AML_MDS_Prob) if s`i'==.
qui g double u`i' = s`i'*LYS
}
I would like to obtain a variable where I can store the rowmeans for all
the
u`i'.
I tried unsuccessfully to use egen v`i'=rowmean(u`i'), but this way I
(obviously, I would suppose) get v`i' for each u`i', and this not help
the
achievement of my research purpose.
*
* 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/
