similar to: Multiple ecdf plots?

Hello Prof. Harrell and dear R-help mailing list, I wish to perform a non-parametric repeated measures anova. If what I read online is true, this could be achieved using a mixed Ordinal Regression model (a.k.a: Proportional Odds Model). I found two packages that seems relevant, but couldn't find any vignette on the subject: http://cran.r-project.org/web/packages/repolr/

I don't understand why nchar() gives different string lengths for vectors in a list than it does for vectors in a dataframe. Here's a snippet of code: > temp <- list(text=c("thug","jimbob","apple","thug"),numbers=1:4) > nchar(temp$text) [1] 4 6 5 4 > temp <- data.frame(temp) > nchar(temp$text) [1] 1 1 1 1 Could someone explain

Good day. We've got a question, as there's a need to mount NFS over RDMA. How we should define a pool on an already mounted NFS folder? The pool must be "shared" for all KVM nodes for live migration purposes. Thanks in advance.

Dear R-help: I am trying to plot the cumulative distribution function of a vector of around 33 million numeric observations. > plot.ecdf(myDataVector) R has been non-responsive for about an hour, and my guess is that it's probably not going to finish. Does anybody have a sense whether this a reasonable experience (and if so, is there a way to get the desired effect, or am I SOL)?

Hi All, following expression: x <- sort(rnorm(10)); e <- ecdf(x); d <- numericDeriv(e(x),"x"); makes d far from approximation of one dimensional pdf. What's wrong then here? Kind regards. --------------------------------------------------------------------------- Valery A.Khamenya Bioinformatics Department BioVisioN AG, Hannover

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

Dear R Gurus, Sorry if this has been asked before but I did not find it in the archives. I would like to add a horizontal display of mean and SEM on Hmisc ecdf plots done by group (ie variate~treatment). Has anyone written some code to do that ? Thanks and kind regards, Jean-Louis

Hi r-users,   I would like to compare the cdf between historical and predicted. My x.obs and x.pre are the frequency data in classes of 0-300. I tried: plot(ecdf(x.obs),ecdf(x.pre),type="l",col="red")   and it gives me: Error in plot.stepfun(x, ..., ylab = ylab, verticals = verticals, pch = pch) :   argument 4 matches multiple formal arguments   Thank you so much for any

Hi, all ecdf function (Empirical Cumulative Distribution Function) in "stats" package counts from smaller values to larger values. However, I want to draw it by counting from larger value to smaller values and I couldn't find options for this purpose. How can I draw ecdf or ecdf like graph by counting from larger values to smaller values. Thank you in advance. Hyunchul Kim

I'm replying to R-devel, the mailing list which should be used to discuss R feature enhancements. >>>>> "Norman" == Norman Warthmann <norman at warthmann.com> >>>>> on Wed, 18 Jan 2006 11:33:22 +0100 writes: Norman> .......... Norman> Is there a specific reason why in the ecdf-function Norman> the variable

Hi! I've noticed the ecdf() R code (R ver. 2.5.0) contains two call to sort: --- [R-code] --- ecdf <- function(x) x <- sort(x) n <- length(x) if (n < 1) stop("'x' must have 1 or more non-missing values") vals <- sort(unique(x)) rval <- approxfun(vals, cumsum(tabulate(match(x, vals)))/n, method

Dear R listers Could somebody please advise me how to draw two empirical cumulative distribution functions on the same set of axes? I know I should be using the ecdf command but I'm not sure what to add to it to force it to plot twice on the axes Thanks Jane ______________________________________________________ Check out all the latest outrageous email attachments on the Outrageous Email

Dear Prof. Brian Ripley, first of all thank you for your answer, I do appreciate how do you manage to keep successfully all your activities and answer posts in this forum! > An empirical CDF is a step function: it does not have a > derivative at the jump points, and has a zero > derivative everywhere else. of course! Let me add few words concerning my simple motivation. 1.

I am trying to use the ecdf function to find p-values (using a vector of numbers to represent my new distribution and a test specific t-statistic value). I am using : 1-ecdf(vector)(t-stat) vector<-c(5.386, 3.701717, 3.8289, 3.602, 4.469, 5.2087, 6.1613, 4.71181, 5.07716, 2.3517) ecdf(vector)(4.6604) [1] 0.5 R will only give me 1 significant digit after the decimal point and I am

Hi, I have a question for the group, perhaps someone can help me figure this out. I've already looked in the help files and they were no help to me. I have a vector of values and I am plotting an ecdf graph. 1. How can i draw a continuous line through the ecdf points? (lines and type for the plot with an ecdf object does not work) 2. Supposing I have this line drawn. I can add a

Hi Everyone, I want to generate some random numbers according to some empirical distribution. Therefore I am looking for the inverse of an empirical cumulative distribution function. I haven't found any in R. Can anyone give a pointer? Thanks, Geoffrey _______________________________________________________=0A= =0A= =0A= The information in this email or in any file attached

Hi all, I would like to plot part of the emperical CDF. Suppose the variable is x, I just need the part when x>1,therefore, I am using the following codes. tail <- x>1 plot(ecdf(x[tail]), do.points=FALSE, verticals=TRUE) The "x" value starts from 1, but the yaxs still begins from 0, not the corresponding value when "x" is 1. How can I make it match? Could anyone

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

I have a sample of observations: > yy [1] 0.00000000 2.39722222 4.35000000 -4.19722222 0.63611111 [6] 1.08055556 5.90555556 -1.87222222 2.13333333 -1.18055556 [11] 3.61666667 0.87777778 8.33888889 3.84166667 1.11111111 [16] -3.76111111 -11.67777778 -2.03055556 6.94444444 -11.76666667 [21] 4.81111111 -7.25833333 1.42222222 5.37222222 4.68055556

Hi. I'm trying to plot two ecdf's on the same graph using two different colors. I can plot using the same color, but it doesn't work when I change colors? Any suggestions? Thanks in advance for your help. x <- c(0.80, 0.83, 1.89, 1.04, 1.45, 1.38, 1.91, 1.64, 0.73, 1.46) y <- c(1.15, 0.88, 0.90, 0.74, 1.21) plot(ecdf(x)) # it works without col='blue', but doesn't