similar to: How to make a smooth ( linear ) CDF plot?

All, Thanks for all of your help & advice! I created a plot starting with qqnorm, if I remember right, but I had to add a LOT of extras. I was surprised that it doesn't include the probability (e.g. 50% at mean) and instead provides the standard deviations (e.g., 0 at mean) on the axis. Plus, there is no line through Q1 & Q3 to help give a sense of the normality of the graph.

On my machine (Mac OS 10.5.6, R 2.8.1) the following plot is drawn such that a substantial portion of the descender of the title is covered by the plot: library(ggplot2) qplot(x=1:10,y=1:10,main='p q j g') -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com Looking to arrange a meeting? Do so at: http://www.timetomeet.info/with/mike/ ~

Hi, I would like to overlap the cdf curve for observed and generated data  Here is my code: plot(cdf,main ="CDF of the sum for winter season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", xaxs="i",yaxs="i",col=c("black","red"), lty=c(1,1),ylab="Cumulative probability", xlim=c(0,800),lwd=1) lines(ecdf(datobs))

Hi all, I'm looking for an analytic method to obtain the PDF & CDF of the sampling distribution of a given correlation (rho) at a given sample size (N). I've attached code describing a monte carlo method of achieving this, and while it is relatively fast, an analytic solution would obviously be optimal. get.cors <- function(i, x, y, N){ end=i*N

Dear all, I have define a discrete distribution P(y_i=x_i)=p_i, which I want to plot a CDF plot. However, I can not find a function in R to draw it for me after searching R and R-archive. I only find the one for the sample CDF instead my theoretical one. I find stepfun can do it for me, however, I want to plot some different CDF with same support x in one plot. I can not manage how to do it with

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

My machine has 8GB memory. I had quit all other programs that might take a lot of memory when I try the script (before I post the first message in this thread). The cdf file is of only 741 MB. It is strange to me to see the error. On Mon, Dec 28, 2009 at 2:38 AM, Wolfgang Huber <whuber at embl.de> wrote: > Dear Peng Yu > > how big is the RAM of your computer? You could try with

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]]

I all, is there a way with R to perform an exact permutation test to replace the wilcoxon test to compare paired series and/or to perform pairwise multiple comparisons for related series after a Friedman test ? Thanks Gilles

Hi, I have a data set which is assumed to follow weibull distr'. How can I find of cdf for this data. For example, for normal data I used (package - lmomco) >cdfnor(15,parnor(lmom.ub(c(df$V1)))) Also, lmomco package does not have functions for finding cdf for some of the distributions like lognormal. Is there any other package, which can handle these distributions?

Hi all, Just a random shot in the dark. Does anyone have/know of a function for the CDF of a generalized error dist? -- Elliot Williams (ewilliams at ucsd.edu) Economics Department, UC San Diego -------------- next part -------------- An embedded message was scrubbed... From: Elliot Williams <ewilliams at ucsd.edu> Subject: [R] Generalized Error Distribution (Exponential Power) CDF?

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.

Hi, I've estimated a simple kernel density of a univariate variable with density(), but after I would like to find out the CDF at specific values. How can I do it? thanks for your help, with it I am very close to finish my first little bit more serious work in R, Viktor

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

Hi All, How can I get the expected value from a discrete cdf? Is there any R function that can do this? Thanks, Annie [[alternative HTML version deleted]]

Hello list, I wish to perform a randomization test on the F-statistics of a 2 way ANOVA but have not been able to find out how to do so - is there a package / function that can perform this that I am unaware of? FactorA has 6 levels (0,1,2,3,4,5) whereas FactorB has 3 (1,2,3). A sample: Resp. FactorA FactorB 2 0 2 3 1 2 1 2 2 0 3 2 0 4 2 0 5 2 4 0 1 6 1 1 1 2 1 0 3 1 1 4 1 0 5 1 2 0 2 3 1 2 1

Dear all, I have the cdf of the following power fuction distribution: F(y)=(y/350)^a               ,0<y<350, where " a " is some parameter with range a>0. I want to use it as the argument of the discretize function of the actuar package. So I think I need to define this function to R so that if I entered a=1, I get the following F(y)=(y/350) and if I entered a=4.5, I get the

Hi, I would like to use the following codes to plot the CDF for pareto distribution. Before doing this, I have plot the emperical one. x <- seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4477557,), col="red") Could anyone give me some advice whether the above codes are correct? Many thanks. -- View this message in context:

-- R Help List -- I am looking for a bivariate normal cdf routine in R. I have some fortran routines for this, which appear to be based on 15-point quadrature. Any guidance/suggestions on making these in loadable R-functions would be appreciated. Thanks, Dan =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Daniel A. Powers, Ph.D. Department of Sociology University of Texas at Austin

Dear Experts, I need to calculate the cdf of the standard normal distribution, i.e. H(x) = 1/sqrt(2*pi) integral(exp(-z^2/2) dz), where z is b/w -infi to infi. I know there should be a way to do it in R, but did not know to do it. I'd appreciate any help you could offer. Charlie Liu Graduate student intern at EPA/ECO