Statalist


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

Re: st: Re: Monte Carlo Simulation, help!


From   STATA Dndee <[email protected]>
To   [email protected]
Subject   Re: st: Re: Monte Carlo Simulation, help!
Date   Sun, 18 Jan 2009 19:03:05 +0000 (GMT)

Thank you Martin!
I did that : "list" or looking at data editor shows 10,000 empty rows!

Kindest regards



--- On Sun, 18/1/09, Martin Weiss <[email protected]> wrote:

> From: Martin Weiss <[email protected]>
> Subject: Re: st: Re: Monte Carlo Simulation, help!
> To: [email protected]
> Date: Sunday, 18 January, 2009, 6:52 PM
> Type -list- in Stata and you see your ouput (10.000 rows, so
> you may prefer 
> a -summarize-)
> 
> HTH
> Martin
> _______________________
> ----- Original Message ----- 
> From: "STATA Dndee"
> <[email protected]>
> To: <[email protected]>
> Sent: Sunday, January 18, 2009 7:46 PM
> Subject: Re: st: Re: Monte Carlo Simulation, help!
> 
> 
> > Many thanks Martin and Carlo,
> > You have already put a lot of time for this thread.
> >
> > I tried:
> >
> > program define lnsim, rclass
> >        version 9.2
> >        syntax [, obs(integer 1) mu(real 0) sigma(real
> 1) ]
> >        drop _all
> >        set obs `obs'
> >        tempvar z
> >        gen `z' = exp(`mu' +
> `sigma'*invnorm(uniform()))
> >        summarize `z'
> >        return scalar mean = r(mean)
> >        return scalar Var  = r(Var)
> >    end
> >
> >
> > and then as in the stata hekp file example:
> >
> > simulate mean=r(mean) var=r(Var), reps(10000): lnsim,
> obs(50) mu(-3)
> >        sigma(7)
> >
> >
> > for each observation or for all of them no error this
> time but I got just 
> > rows like this(from 50 to 10,000):
> >
> > Simulations (10000)
> > ----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
> > .xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx   
> 50
> >
> > without any data output.
> >
> > Anyway, thank you again for your time.
> > Kindest regards
> >
> >
> >
> >
> > --- On Sun, 18/1/09, Martin Weiss
> <[email protected]> wrote:
> >
> >> From: Martin Weiss <[email protected]>
> >> Subject: Re: st: Re: Monte Carlo Simulation, help!
> >> To: [email protected]
> >> Date: Sunday, 18 January, 2009, 5:17 PM
> >> This thread is getting weirder by the minute:
> StataDndee
> >> seems to think that
> >> Stata will read each line of this dataset, and
> enter his
> >> information into
> >> the -simulate- command. This is not going to
> happen. If you
> >> want to enter
> >> these values one after another and run -simulate-,
> look at
> >> -help foreach-...
> >>
> >> HTH
> >> Martin
> >> _______________________
> >> ----- Original Message ----- 
> >> From: "STATA Dndee"
> >> <[email protected]>
> >> To: <[email protected]>
> >> Sent: Sunday, January 18, 2009 5:47 PM
> >> Subject: Re: st: Re: Monte Carlo Simulation, help!
> >>
> >>
> >> > Thanks dear Martin
> >> > I tried that and got the same error. Then
> arranged the
> >> data in data editor
> >> > like this
> >> >
> >> >     +-----------------------+
> >> >     | obs   mean        var |
> >> >     |-----------------------|
> >> >  1. |   1     71   .0007049 |
> >> >  2. |   2    259    .000269 |
> >> >  3. |   3    121     .00066 |
> >> >  4. |   4    262   .0003177 |
> >> >  5. |   5    192      .0005 |
> >> >     |-----------------------|
> >> >  6. |   6    120   .0011873 |
> >> >     +-----------------------+
> >> > and then this command with vluaes of the mean
> obs from
> >> the last table:
> >> >
> >> > simulate mean=r(mean) var=r(Var),
> reps(10000): lnsim,
> >> obs(6) mu(713
> >> > ) sigma(0.000102877)
> >> >
> >> > And got the same error
> >> > loc not found
> >> > an error occurred when simulate executed
> lnsim
> >> > r(111);
> >> >
> >> > And that as before everything in data editor
> was gone.
> >> there should be
> >> > something with "lnsim" in command.
> >> >
> >> > Thank you so much for your help.
> >> >
> >> >
> >> > --- On Sun, 18/1/09, Martin Weiss
> >> <[email protected]> wrote:
> >> >
> >> >> From: Martin Weiss
> <[email protected]>
> >> >> Subject: Re: st: Re: Monte Carlo
> Simulation, help!
> >> >> To: [email protected]
> >> >> Date: Sunday, 18 January, 2009, 3:55 PM
> >> >> "loc not found" sounds weird,
> the real
> >> trouble is
> >> >> that you put a "%" sign in
> >> >> the "sigma" specification. Just
> put
> >> 0.027
> >> >> there...
> >> >>
> >> >> HTH
> >> >> Martin
> >> >> _______________________
> >> >> ----- Original Message ----- 
> >> >> From: "STATA Dndee"
> >> >> <[email protected]>
> >> >> To:
> <[email protected]>
> >> >> Sent: Sunday, January 18, 2009 4:50 PM
> >> >> Subject: Re: st: Re: Monte Carlo
> Simulation, help!
> >> >>
> >> >>
> >> >> > Hi
> >> >> >
> >> >> > I have tried the following commande:
> >> >> >
> >> >> > simulate mean=r(mean) var=r(Var),
> >> reps(10000): lnsim,
> >> >> obs(100)
> >> >> >
> >> >> > from STATA help and for Inst. F
> >> >> >
> >> >> > my command was:
> >> >> > simulate mean=r(mean) var=r(Var),
> >> reps(10000): lnsim,
> >> >> obs(273) mu(71)
> >> >> > sigma(2.7%)
> >> >> >
> >> >> > and I got the following error:
> >> >> >
> >> >> > loc not found
> >> >> > an error occurred when simulate
> executed
> >> lnsim
> >> >> > r(111);
> >> >> > When it resolved I should run the
> analysis
> >> for all
> >> >> INSs and the mean?
> >> >> >
> >> >> > the SE in the floowing table was SD.
> >> >> >
> >> >> >
> >> >> > I have copied the  table again.
> >> >> >
> >> >> > Many thanks for any comment
> >> >> > Bests
> >> >> >
> >> >> >
> >> >> >
> >> >>
> >>
> +-------------------------------------------------------+
> >> >> >     | inst   cases   total   outcome
>     lci
> >>    uci
> >> >>   sd |
> >> >> >
> >> >>
> >>
> |-------------------------------------------------------|
> >> >> >  1. |    f      71     273    26.01%
>   20.8%
> >>  31.2%
> >> >> 2.7% |
> >> >> >  2. |    e     259     809    32.01%
>   28.8%
> >>  35.2%
> >> >> 1.6% |
> >> >> >  3. |    c     121     345    35.07%
>   30.0%
> >>  40.1%
> >> >> 2.6% |
> >> >> >  4. |    d     262     726    36.09%
>   32.6%
> >>  39.6%
> >> >> 1.8% |
> >> >> >  5. |    b     192     480    40.00%
>   35.6%
> >>  44.4%
> >> >> 2.2% |
> >> >> >
> >> >>
> >>
> |-------------------------------------------------------|
> >> >> >  6. |    a     120     204    58.82%
>   52.1%
> >>  65.6%
> >> >> 3.4% |
> >> >> >  7. | mean     713    2153    33.12%
>   31.1%
> >>  35.1%
> >> >> 1.0% |
> >> >> >
> >> >> >
> >> >> >
> >> >> >
> >> >> >
> >> >> > --- On Sun, 18/1/09, STATA Dndee
> >> >> <[email protected]> wrote:
> >> >> >
> >> >> >> From: STATA Dndee
> >> <[email protected]>
> >> >> >> Subject: Re: st: Re: Monte Carlo
> >> Simulation, help!
> >> >> >> To:
> [email protected]
> >> >> >> Date: Sunday, 18 January, 2009,
> 1:20 PM
> >> >> >> Thanks for your kind and quick
> reply
> >> Martin
> >> >> >>
> >> >> >> I tried it before but
> couldn't
> >> understand the
> >> >> commands
> >> >> >> in Help file. I arranged the
> data with
> >> the CIs.
> >> >> >>
> >> >> >>
> >> >> >>
> >> >>
> >>
> +-------------------------------------------------------+
> >> >> >>      | inst   cases   total  
> outcome
> >> lci
> >> >> uci
> >> >> >> se |
> >> >> >>
> >> >> >>
> >> >>
> >>
> |-------------------------------------------------------|
> >> >> >>   1. |    f      71     273   
> 26.01%
> >> 20.8%
> >> >> 31.2%
> >> >> >> 2.7% |
> >> >> >>   2. |    e     259     809   
> 32.01%
> >> 28.8%
> >> >> 35.2%
> >> >> >> 1.6% |
> >> >> >>   3. |    c     121     345   
> 35.07%
> >> 30.0%
> >> >> 40.1%
> >> >> >> 2.6% |
> >> >> >>   4. |    d     262     726   
> 36.09%
> >> 32.6%
> >> >> 39.6%
> >> >> >> 1.8% |
> >> >> >>   5. |    b     192     480   
> 40.00%
> >> 35.6%
> >> >> 44.4%
> >> >> >> 2.2% |
> >> >> >>      |
> >> >> >>   6. |    a     120     204   
> 58.82%
> >> 52.1%
> >> >> 65.6%
> >> >> >> 3.4% |
> >> >> >>   7. | mean     713    2153   
> 33.12%
> >> 31.1%
> >> >> 35.1%
> >> >> >> 1.0% |
> >> >> >>
> >> >> >>
> >> >>
> >>
> +-------------------------------------------------------+
> >> >> >>
> >> >> >> Regards
> >> >> >>
> >> >> >>
> >> >> >>
> >> >> >> --- On Sun, 18/1/09, Martin
> Weiss
> >> >> >> <[email protected]>
> wrote:
> >> >> >>
> >> >> >> > From: Martin Weiss
> >> >> <[email protected]>
> >> >> >> > Subject: st: Re: Monte
> Carlo
> >> Simulation,
> >> >> help!
> >> >> >> > To:
> [email protected]
> >> >> >> > Date: Sunday, 18 January,
> 2009,
> >> 11:49 AM
> >> >> >> > A very warm welcome to you
> in the
> >> Stata
> >> >> community,
> >> >> >> good to
> >> >> >> > see you :-) Try -help
> simulate- for
> >> a start!
> >> >> >> >
> >> >> >> > HTH
> >> >> >> > Martin
> >> >> >> > _______________________
> >> >> >> > ----- Original Message
> ----- From:
> >> >> "STATA
> >> >> >> Dndee"
> >> >> >> >
> <[email protected]>
> >> >> >> > To:
> >> <[email protected]>
> >> >> >> > Sent: Sunday, January 18,
> 2009 12:46
> >> PM
> >> >> >> > Subject: st: Monte Carlo
> Simulation,
> >> help!
> >> >> >> >
> >> >> >> >
> >> >> >> > > Hi
> >> >> >> > > I want to rank
> institutions for
> >> and
> >> >> outcome given
> >> >> >> in
> >> >> >> > outcome% based on Monte
> Carlo
> >> simulation and
> >> >> produce a
> >> >> >> table
> >> >> >> > like table 2 in
> >> >> >> >
> >> >>
> http://plus.maths.org/issue46/risk/index.html.
> >> The
> >> >> >> method
> >> >> >> > is also described by
> Spiegelhalter
> >> in BMJ,
> >> >> 1998,
> >> >> >> > 316(1701-1705).
> >> >> >> > > I am a beginner in
> STATA. Could
> >> anyone
> >> >> help me in
> >> >> >> > performing the analysis in
> STATA? I
> >> need the
> >> >> commands
> >> >> >> and
> >> >> >> > number of iterations. 
> Number of
> >> institutions
> >> >> is more
> >> >> >> than
> >> >> >> > in the sample.
> >> >> >> > >
> >> >> >> > > Any help will be
> really
> >> appreciated.
> >> >> >> > > Kind regards
> >> >> >> > >
> >> >> >> > > Inst. cases total
> outcome%
> >> >> >> > > f 71 273 26.01%
> >> >> >> > > e 259 809 32.01%
> >> >> >> > > c 121 345 35.07%
> >> >> >> > > d 262 726 36.09%
> >> >> >> > > mean 713 2153 33.12%
> >> >> >> > > b 192 480 40.00%
> >> >> >> > > a 120 204 58.82%
> >> >> >> > >
> >> >> >> > >
> >> >> >> > >
> >> >> >> > >
> >> >> >> > > *
> >> >> >> > > *   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/
> >> >> >
> >> >> >
> >> >> >
> >> >> >
> >> >> > *
> >> >> > *   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/
> >> >
> >>
> >> *
> >> *   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/


      

*
*   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–2024 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index