[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 |
AW: st: Bootstrapping new observations to add to an existing dataset |

Date |
Mon, 22 Jun 2009 15:12:11 +0200 |

<> Note that you can -inspect- your data and the -return- list contains the value " r(N_unique)" which does give you a hint as to the kind of var you are facing. If it is 2, it`s a binary one and so on... ************* clear* set obs 200 input x y 0 1 0 2 1 3 0 6 1 8 end gen z = rnormal() foreach var of varlist _all{ ins `var' ret li } ************* HTH Martin -----Ursprüngliche Nachricht----- Von: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] Im Auftrag von Davide Cantoni Gesendet: Montag, 22. Juni 2009 14:59 An: statalist@hsphsun2.harvard.edu Betreff: Re: st: Bootstrapping new observations to add to an existing dataset Thank you, gentlemen. To clarify: yes, In a sense I want more of the same observations. I just don't want to create new observations by filling up the dataset with -invnorm(uniform())- because I want to preserve the properties of the variables: some are binary, some are integers, and some are just any kind of number. The values of these variables are irrelevant for my regressions, their underlying distribution is not: I want binary variables to remain binary etc. Since there are more than 300 variables in the dataset, I do not want to do this one-by-one, finding out whether var4 consists only of integers between 1 and 10, var231 is binary etc. So one way to go ist just take more of the same observations. E.g., obs 201 is the same as obs 133, obs 202 is the same as obs 78 and so on. Or, and that's the other thing I was thinking of, I could create obs 201 by drawing a new value from the existing distribution (given by obs1-obs200) of var1, then drawing a new value from the distribution of var2 and so on... Davide 2009/6/22 Martin Weiss <martin.weiss1@gmx.de>: > > <> > > Davide said that he wanted to keep "the same (unknown to me) data generating > process". Every advice so far has assumed that this means that he just wants > "more of the same observations". If that was the case, he could also use a > random frequency weight for every observation which would reduce the size of > his dataset. > > Davide could clarify whether he merely wants to duplicate observations > randomly or whether he really wants "new" observations... > > > > HTH > Martin > > > -----Ursprüngliche Nachricht----- > Von: owner-statalist@hsphsun2.harvard.edu > [mailto:owner-statalist@hsphsun2.harvard.edu] Im Auftrag von Austin Nichols > Gesendet: Montag, 22. Juni 2009 14:33 > An: statalist@hsphsun2.harvard.edu > Betreff: Re: st: Bootstrapping new observations to add to an existing > dataset > > Davide Cantoni <davide.cantoni@gmail.com> : > You don't say how many more obs you want--let's assume you want about > 100 times as many: > > expand 100 > > will do it, or > > g u=round(uniform()*200) > expand u > > for a random-sized sample about 100 times as big with the same DGP. > You could also > > loc n=_N*100 > g u=round(uniform()*1000) > expand u > drop u > g u=uniform() > sort u > drop if _n>`n' > > for a sample 100 times as big but with random numbers of replications > of each obs. > > On Sun, Jun 21, 2009 at 11:58 PM, Davide Cantoni > <davide.cantoni@gmail.com> wrote: >> >> Hello, I am stuck while thinking about this issue and I would >> appreciate your suggestions. I have a dataset which I use for >> simulation purposes, to test whether my do-files run correctly. The >> issue is that this dataset is too short for many applications, as it >> has only 200 observations. >> >> What I want to do is expand this dataset to include more observations, >> but keeping the same (unknown to me) data generating process that >> created the first 200 observations. So I was thinking to proceed in a >> bootstrapping manner, by drawing the values for each one of the >> variables (var1, var2 etc etc) for the new observations from the >> empirical distributions of var1, var2,... in the first 200 >> observations. Yet, I have no idea on how to implement this. I'm >> grateful for any idea. Thanks for your interest, >> >> Davide > * > * 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/

**References**:**st: Bootstrapping new observations to add to an existing dataset***From:*Davide Cantoni <davide.cantoni@gmail.com>

**Re: st: Bootstrapping new observations to add to an existing dataset***From:*Austin Nichols <austinnichols@gmail.com>

**Re: st: Bootstrapping new observations to add to an existing dataset***From:*Davide Cantoni <davide.cantoni@gmail.com>

- Prev by Date:
**Re: st: Bootstrapping new observations to add to an existing dataset** - Next by Date:
**st: AW: dropping one unit at a time for test** - Previous by thread:
**Re: st: Bootstrapping new observations to add to an existing dataset** - Next by thread:
**R: st: Bootstrapping new observations to add to an existing dataset** - Index(es):

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