similar to: how to generate a random data from a empirical distribition

On 7/27/2010 6:00 AM, r-help-request at r-project.org wrote: > Date: Mon, 26 Jul 2010 11:36:29 -0700 (PDT) > From: xin wei<xinwei at stat.psu.edu> > To:r-help at r-project.org > Subject: [R] how to generate a random data from a empirical > distribition > Message-ID:<1280169389379-2302716.post at n4.nabble.com> > Content-Type: text/plain; charset=us-ascii > >

http://r.789695.n4.nabble.com/file/n3567636/ecdfs.jpg ecdfs.jpg http://r.789695.n4.nabble.com/file/n3567636/ecdf_curve.gif ecdf_curve.gif Hello, I have generated a plot of two empirical CDFs (attachment 1). As a result, they are stepwise when plotted. The following code was used: > plot(ecdf(mut), do.points=FALSE, verticals=TRUE, xlim=range(mut, non), > col="red") >

Hi folks, I know that density function will give a estimated density for a give dataset. Now from that I want to have a percentage estimation for a certain range. For examle: > y = density(c(-20,rep(0,98),20)) > plot(y, xlim=c(-4,4)) Now if I want to know the percentage of data lying in (-20,2). Basically it should be the area of the curve from -20 to 2. Anybody knows a simple

Hello, I wanted to know if there is a simple way of getting the inverse cdf for a KDE estimate of a density (using the ks or KernSmooth packages) in R ? The method I'm using now is to perform a numerical integration of the pdf to get the cdf and then doing a search for the desired probablity value, which is highly inefficient and very slow. Thanks, -fj [[alternative HTML version deleted]]

Forum, I'm a grad student in Civil Eng, took some Stats classes that required students learn R, and I have since taken to R and use it for as much as I can. Back in my lab/office, many of my fellow grad students still use proprietary software at the behest of advisers who are familiar with the recommended software (Statistica, @Risk (Excel Add-on), etc). I have spent a lot of time learning

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:

Dear R-list, I have made a simple kernel density estimation by x <- c(2,1,3,2,3,0,4,5,10,11,12,11,10) kde <- density(x,n=100) Now I would like to know the estimated probability that a new observation falls into the interval 0<x<3. How can I integrate over the corresponding interval? In several R-packages for kernel density estimation I did not found a corresponding function. I

Hi, I have just started using R. Now I have the following problem: I want to create an Empirical Cumulative Distribution Function and I only came so far: F10 <- ecdf(x) plot(F10, verticals= TRUE, do.p = TRUE, lwd=3) x=c(1.6,1.8,2.4,2.7,2.9,3.3,3.4,3.4,4,5.2) Now I'd like to use arguments such as xlabs and main but I don't know how to integrate them. I hope someone can help me, I am

Hello,   Could someone help me please and to tell how to get the probability from empirical distribution (not parametric) for each data value (R function). For example, for normal distribution there is such a function like:   “pnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)”   I need the same function only for the empirical distribution (which does not correspond to any typical

Hello,   Could someone help me please and to tell how to get the probability from empirical distribution (not parametric) for each data value (R function). For example, for normal distribution there is such a function like:   “pnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)”   I need the same function only for the empirical distribution (which does not correspond to any typical

I'm trying to produce empirical bayes estimates based on the lognormal model in disease mapping Is there a way this can be done in R? thanks Oarabile

Is it possible to get the empirical (sandwich) S.E. estimates for the fixed effects in lme () (thus allowing possibly correlated errors within the group)? In SAS you can get it by the 'empirical' option to PROC MIXED. Cheers, Michael -- Na (Michael) Li, Ph.D. Division of Biostatistics A443 Mayo Building, MMC 303 School of Public Health 420 Delaware

Hi, are there any function to plot the empirical probability distribution function? I just don't want to reinvent the wheel... Best wishes, Ott -- Ott Toomet PhD Student Dept. of Economics ?rhus University Building 322 Universitetsparken 8000 ?rhus C Denmark otoomet (a) econ au dk ph: (+45) 89 42 20 27 ------------------------------------------- (o_ (*_ (O_

Dear Sir/Madam, I apologize for any cross-posting. I got a simple question, which I thought the R list may help me to find an answer. Suppose we have Y_1, Y_2, ., Y_n ~ Poisson (Lambda_i) and Lambda_i ~Gamma(alpha_i, beta_i). Empirical Bayes Estimator for hyper-parameters of the gamma distr, i.e. (alpha_t, beta_t) are needed. y=c(12,5,17,14) n=4 What about a Hierarchal B ayes

Hi, I was trying to draw an empirical distribution function with uniform confidence bands. So I tried to find a way to calculate values of the Kolmogorov-Smirnov Distribution but failed. I guess it must be hidden somewhere (since the ks-test is implemented), but I was unable to find it. Is there any way to do this? Thanks Leif Boysen

Hi, I would like to implement the empirical copula in R, does anyone know if it is included in a package? I know it is not in the "Copula" package. This one only includes a gof-test based on the empirical copula process. Thanks for your help! Gregor -- View this message in context: http://www.nabble.com/Empirical-copula-in-R-tf4018319.html#a11412335 Sent from the R help mailing list

I've got the graphic with this: a=c(120,40,75,85,55,75,55,90,90,55,155) plot(ecdf(a)) How can I draw lines to show the lower quatile and upper quartil? I only know there is the argument lines(), but I don't know how to use it -- View this message in context: http://r.789695.n4.nabble.com/How-can-I-draw-a-line-in-empirical-distribution-function-tp3090675p3090675.html Sent from the R

Hi I have created the following plot over the empirical returns.. What I now want to do is to overlay a curve/line with the normal density as a comparison of the two. Does anyone know how to do this? (NB the last two lines are the problem, and are wrong, I know). Thank you in advance! Rikke http://r.789695.n4.nabble.com/file/n3768783/S%26P_500_spot_and_return_2010.csv

Dear R Users, Does R have an "inverse empirical cumulative distribution" function, something one can use to invert ecdf ? Thanks in advance, Tolga Generally, this communication is for informational purposes only and it is not intended as an offer or solicitation for the purchase or sale of any financial instrument or as an official confirmation of any transaction. In the event you

Dear list, I would like to plot the empirical cumulative distribution of the time needed by a treatment to attain a certain goal. A number of experiments is run with a strict time limit. In some experiments the goal is attained before the time limit, in other experiments time expires before the goal is attained. The situation is very similar to survivial analysis with censored data. I tryed