similar to: multiple plots in a graph

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

Hi I am looking for a structural equation modeling package in R which can be used for categorical data. Is anyone aware of the existence of such a package? Would appreciate any help on this. Thank you Upasna -- --------------------------------------------------------------------- Upasna Sharma Research Scholar Shailesh J. Mehta School of Management, Indian Institute of Technology, Bombay

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

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 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 all, I calculate the test statistic for the KS test outside R, and wish to use R only to calculate the corresponding p-value. Is there a way for doing this? (as far as I see, ks.test() requires raw data as input). Alternatively, is there a way to provide the ks.test() the two CDFs (two samples test) rather than the (x, y) data vectors? Thanks in advance, Rani

Hi. I have the following two ecdf plots in one graph: plot( ecdf(....), do.points=FALSE, verticals=TRUE, main=paste("Ecdf of distances ",DIM,sep=""), col="red" ); lines( ecdf(....), do.points=FALSE, verticals=TRUE ); How do I change the color of the resulting graph? Adding col="red" to either plot or lines results in an error

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

I have a problem where I need to calculate the corresponding cohort percentile ranks for each of several variables. Essentially, what I need is a function that will calculate the distribution-free percentiles from each variable's data vector, returning a corresponding vector of percentiles: e.g.: percentile.my.data<-/function/(my.data) I tried to make ecdf() perform this task but

Hi I have used the hetcor function on dataset which consists of about 40 variables. Since the output is quite large entire output is not visible on the R screen. Is there a way, I can retrieve the output that is not visible? Thank you Regards Upasna -- --------------------------------------------------------------------- Upasna Sharma Research Scholar Shailesh J. Mehta School of Management,

Hi I need to know which package in R calculates the Goodman Kruskal's tau statistic for nominal data. Also is there any implementation for multiple classification analysis (Andrews at al 1973) in R? Any information on this would be greatly appreciated. Thank you Upasna -- --------------------------------------------------------------------- Upasna Sharma Research Scholar Shailesh J. Mehta

Dear All, I have to get the value of say 90th percentile of precipitation time series.. The series is of daily precipitation value of 96 years, I have to to get 90the percentile value of daily precipitation each year. If you know the R code or command for this please let me know. I would appreciate your early response. Thanking you, Sincerely, Ajay.

Dear List, I have some discrete data and want to calculate the percentiles and the percentile ranks for each of the unique scores. I can calculate the percentiles with quantile(). I know that "ecdf" can be used to calculate the empirical cumulative distribution. However, I don't know how to exact the cumulative probabilities for each unique element. The requirement is similar

Heeeelp Given this data x <- c(1,5,100,300,250,200,550,900,1000) > quantile(x) 0% 25% 50% 75% 100% 1 100 250 550 1000 When I run the quantile, I can only know the value of the nth percentile I want to know what's the percentile position of each items in the list Sample 1 = 100% on the list has 1 or more 5 = more than x% on the list has 5 or more 100 = more than x% on the

Given a dataset x, the ecdf is ecdf(x). Then I can use ecdf(x)(y) to find the percentile of y. Given the ecdf is there a way to determine what is the value of y that is the boundary of let's say 95 percentile? In other words, is there a function I can call on the ecdf like: fomeFunc( ecdf( x ), 0.95 ) Which will return the highest value of y, for which ecdf( y ) < 0.95? The only solution

Greetings, I have two questions that I could not answer from the documentation. A - ecdf and confidence intervals : Is there a (simple) way to generate confidence intervals (95%) for a ecdf? B - cross-validation of rpart trees : a colleague is using S to generate decision tree and mentioned to me the use of cross-validation. Is this function enabled in R ? if so, how should one proceed to

Dear all, I am familiar with obtaining the value corresponding to a chosen probability via the quantile function. Now I am facing the opposite problem I have a value an want to know it's corresponding percentile in the distribution. So is there a function for this as well? Thank you for your support in advance, Felix

Hi, can anyone help me to figure out how to compute the percentile of an individual observation with respect to a reference distribution. What I mean is. Let's assume I have a vector consisting of 10 numbers {3,5,8,1,9,5,4,3,5.5,7} and I want figure out what percentile the number 4.9 corresponds to. I failed to find any reference to such a function, although I would assume this must

Dear All, we have: t<-seq(0,24,1) a<-10*exp(-0.05*t) b<-10*exp(-0.07*t) c<-10*exp(-0.1*t) d<-10*exp(-0.03*t) z<-data.frame(a,b,c,d) res<-t(apply(z, 1, quantile, probs=c(0.3))) my goal is to do a 'reverse" of the function here that produces "res" on a data frame, ie: to get the answer 0.3 back for the percentile location when I have

x <- c(6.6493705109108, 7.1348436721241, 8.76886994525624, 6.12907548096037, 6.88379118678109, 7.17841879427688, 7.90737237492867, 7.1207373264833, 7.82949407630692, 6.90411547316105) plot(ecdf(x), log="x") It does the plot fine, but complains: Warning message: In xy.coords(x, y, xlabel, ylabel, log) : 1 x value <= 0 omitted from logarithmic