similar to: Error in plotmath

Hello all, I would like to plot the emperical CDF, normal CDF and pareto CDF in the same graph and I amusing the following codes. "z" is a vector and I just need the part when z between 1.6 and 3. plot(ecdf(z), do.points=FALSE, verticals=TRUE, xlim=c(1.6,3),ylim=c(1-sum(z>1.6)/length(z), 1)) x <- seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4373,-0.2398), col="red") y

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

Dear UseRs, maybe is a silly question: how can I get Empirical CDF values from an object created with ecdf()?? Using print I obtain: Empirical CDF Call: ecdf(t) x[1:57] = 4.1, 4.4, 4.5, ..., 491.3, 671.27 Thanks in advance. Regards, Vito Diventare costruttori di soluzioni Became solutions' constructors "The business of the statistician is to catalyze the scientific

Hi all, is there any package in R capable of smooth approximation of CDF basing on given sample? (Thus, I am not speaking about ecdf) In particular, I expect very much that the approximation should subject to the property: f(x0)<=f(x1) for x0<x1, where x0 and x1 belong to range of the sample given. Polynomial approximation could be OK for me as well. P.S.

Good Afternoon, I'm trying to create a cdf plot, with the following code. It works well, but I have little doubt, if you can help solve. When I create the plot, like the graph line would still not appear with point #cdf x<-table(Dataset$Apcode) View(s) hist(s) *plot(ecdf(x))* x<-1 37607 2 26625 3 5856 4 25992 5 30585 6 16064 7 9850 .. ... .. 186 52 -- View this message in

Hi, I'm bit confused about ecdf (read the help files but still not sure about this). I have an analytical expression for the pdf, but want to get the empirical cdf. How do I use this analytical expression with ecdf? If this helps make it concrete, the pdf is: f(u) = \sum_{t = 1}^T 1/n_t \sum_{i = 1}^{n_t} 1/w K((u - u_{it})/w) where K = kernel density estimator, w = weights, and u_{it} =

Dear List: I am trying to find the area between two ECDFs. I am examining the gap in performance between two groups, males and females on a student achievement test in math, which is a continuous metric. I start by creating a subset of the dataframe male<-subset(datafile, female="Male") female<-subset(datafile, female="Female") I then plot the two CDFs via

Dear R Users, I am trying to suppress the information printed by the ecdf function during an assignment. Various alternatives have failed me so far: > a=summary(ecdf(rnorm(100)))["1st Qu."] Empirical CDF: 100 unique values with summary > invisible(a=summary(ecdf(rnorm(100)))["1st Qu."]) Empirical CDF: 100 unique values with summary > (function()

Does R do cumulative frequency distribution plots? -- Tommy E. Cathey, Senior Scientific Application Consultant High Performance Computing & Scientific Visualization SAIC, Supporting the EPA Research Triangle Park, NC 919-541-1500 EMail: cathey.tommy at epa.gov My e-mail does not reflect the opinion of SAIC or the EPA. Federal Contact - John B. Smith 919-541-1087 - smith.johnb at epa.gov

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") >

Dear list, This is a copy of a mail sent to Max Kuhn, original author and maintainer of the odfWeave package, which seems not to have received it. It reports a problem that seems to be very implementation specific (reproductible on three Debian testing amd64 machine, does *not* happen on two i686 Debian testing systems, does *not* happen on an Ubuntu 11.06 amd64 machine) and therefore not

Hi r-users, I would like to overlap 2 ecdf plots.  I tried this below and it gives me two plots of ecdf but just both just in black. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs)) lines(ecdf(gam_sum_gen)) Then I try to add colors etc and also the legend but fail. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i")

There seems to be 2 functions call ecdf... http://lib.stat.cmu.edu/S/Harrell/help/Hmisc/html/ecdf.html http://127.0.0.1:11885/library/stats/html/ecdf.html How do I get the one ecdf {Hmisc} to run instead of the ecdf {stats} A pointer in the right direction would be greatly appreciated. Tried to instal Hmisc but got this message, so I assume I have it > utils:::menuInstallPkgs() Warning:

Dear Group, I am using the ecdf function as follows: cawa.cdp <- ecdf(cawaocc$LEFF80) summary(cawa.cdp) Empirical CDF: 223 unique values with summary Min. 1st Qu. Median Mean 3rd Qu. Max. 0.07918 1.35700 1.68600 1.61000 1.91200 2.70000 I can see by the summary that the y-value for the 3rd quartile is 1.912. How can I obtain the x-value for a specified y-value (e.g., 0.8)?

Dear R users, I'm using Ecdf (Hmisc library) to plot four cdf in a same graphic. In this graphic I also plot the 0.99 quantile for these cdf. I successfully plot cdfs using different types of line to distinguish them, but I can't determine the type of lines showing 0.99 quantile. Is there a way to assign different line types for quantile lines in Ecdf plot? Best regards, -- Mateus da

Good day everyone, I want to assess the error when fitting a Gram-Charlier CDF to some data 'ws', that is, I want to calculate: Err = |ecdf(ws) - GCh_ser(ws)| The problem is, I cannot get the F(x) values from the ecdf. 'Summary(ecdf())' returns some of the x-axis values, but how do you get the F(x) values? Thank you for any help you can provide. Regards, Augusto

Hi,   I would like to distinguish my plots using colors but I got error message.  How do I correct that?     plot(ecdf(z), main ="CDF for observed and simulated weighted sum",type="l",lwd=2,col="blue", xlab="Weighted sum (mm)", ylab="Cumulative Percent", xlim=c(0,15), xaxs ='i', yaxs ='i',ylim=c(0,1)) par(new=TRUE)

Hello all, I have two columns of numbers. I would like to do the following: (1) Plot both cdfs, F1 and F2 on the same graph. (2) Find smoothed approximations of F1 and F2 lets call them F1hat and F2hat (3) Find values for F1hat when we substitue a value of x in it. (4) Find the corresponding densities of the cdfs. Any ideas? -- Thanks, Jim. [[alternative HTML version deleted]]

Dear all, I would like for a given vector to calculate the pdf and the cdf and then plot these two graphs. Is there in R any toolkit that can help me do that fast? Regards Alex