[R] Drawing from an empirical distribution

roger koenker rkoenker at uiuc.edu
Wed Jan 7 00:25:10 CET 2009


Nit-picking about syntax does seem needed, mea culpa,  I
intended something more like:

	Qn <- ecqf(x)
	Qn(runif(whatever))

On Jan 6, 2009, at 5:06 PM, roger koenker wrote:

> Sure, but it would be more 'fun' to modify ecdf() slightly to produce
> an ecqf()  function -- essentially reversing the arguments to  
> approxfun()--
> and then use
>
> 	ecqf(runif(whatever))
>
> no nit-picking about efficiency, please.
>
>
> url:    www.econ.uiuc.edu/~roger            Roger Koenker
> email    rkoenker at uiuc.edu            Department of Economics
> vox:     217-333-4558                University of Illinois
> fax:       217-244-6678                Champaign, IL 61820
>
>
>
> On Jan 6, 2009, at 4:42 PM, Antonio, Fabio Di Narzo wrote:
>
>> If the ecdf is 'ecdf(x)', do just:
>>> sample(x, size=whatever, replace=TRUE)
>>
>> HTH,
>> Antonio.
>>
>> 2009/1/6 culpritNr1 <ig2ar-saf1 at yahoo.co.uk>:
>>>
>>> Hi All,
>>>
>>> Does anybody know if there is a simple way to draw numbers from an  
>>> empirical
>>> distribution?
>>>
>>> I know that I can plot the empirical cumulative distribution  
>>> function this
>>> easy:
>>> plot(ecdf(x))
>>>
>>> Now I want to pick a number between 0 and 1 and go back to domain  
>>> of x.
>>> Sounds simple to me.
>>>
>>> Any suggestion?
>>>
>>> Thank you,
>>>
>>> Your culprit
>>> (everybody needs a culprit)
>>>
>>>
>>> --
>>> View this message in context: http://www.nabble.com/Drawing-from-an-empirical-distribution-tp21320810p21320810.html
>>> Sent from the R help mailing list archive at Nabble.com.
>>>
>>> ______________________________________________
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>
>>
>>
>> -- 
>> Antonio, Fabio Di Narzo
>> Ph.D. student at
>> Department of Statistical Sciences
>> University of Bologna, Italy
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.




More information about the R-help mailing list