similar to: Trouble using Ecdf () from the Hmisc library

Hi, (please Cc me) In xyplot (), type = "l" (or one that includes "l", *el*) is (generally) meaningful only when the 'x' variable is sorted. In practice, one either sorts the data frame before hand or writes a tiny panel function which sorts the supplied x and then calls the default panel.xyplot(). Trouble arises when there is a conditional variable as well as a

Full_Name: P Kapat Version: 2.1.1 (2005-06-20) OS: GNU/Linux 2.6.8-2-386, Debian testing Submission from: (NULL) (65.24.56.41) Relevant Bugs Ids : 7987, 7989 H is a 100x100 singular but symmetric matrix (a matrix defining the neighbourhood structure for a spatial data) available from.. http://www.stat.ohio-state.edu/~pkapat/miscl/bugMatrix.RData > load(file="bugMatrix.RData") >

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 again, this is a follow-up of my previous post : I've noticed that > attr(first.list[[1]]$size.ecdf,"call") gives ecdf(test.moms[, "m.pxs"]) and that there is no "test.moms" field in my records (a local variable named test.moms was used inside a function, but it was stored under another name in the list, outside the function) May this be the answer (and

Hi all, I'm working with data to which I'm fitting three-parameter weibull distributions (shape, scale & shift). The data are of low sample sizes (between 10 and 80 observations), so I'm reluctant to check my fits using chi-square (also, I'd like to avoid bin choice issues). I'd use the Kolmogorov-Smirnov test, but of course this is invalid when the distribution

Hi, i want to plot empirical cumulative density functions for two variables in one plot. For better visualizing the differences in the two cumulative curves I'd like to log-scale the axis. So far I found 3 possible functions to plot ecdf: 1) ecdf() from the package 'stats'. I don't know how to successfully set the log.scales? Combining two plots is not a problem: plot(ecdf(x1))

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

Hi, Newbie here. I read the R for Beginners but i still don't get this. I have the following data (this is just an example) in a CSV file: courseid numstudents 101 209 141 13 246 140 263 8 321 10 361 10 364 28 365 25 366 23 367

Dear R users, I am trying to minimize the distance between my data points and theoretical gamma distribution over shape and scale parameters. the function "mde" from actuar package does it for empirical distribution function and theoretical gamma distribution. However, I would like to minimize the distance by using only the data between 0.1 and 0.9 quantiles. I cannot use ecdf in this

Hi, I have two variables ranging both from 0 to 1 (n=500 each). Now I am interested in plotting them both in one plot (using ggplot2). So far I used ecdf() (from an example I found with google) to get values for the cumulatice distribution function which gives a relative curve. I also want to do the same plot but using absolute cumulative values instead of relative. Can that be done with ecdf or

Hi All, I have codes as follows to get the ecdf plots: > day.hos2<-c(6,4,6,6,4,6,5,4,7,5,6,6,8,6,17,9,8,4,6,3,5,8,7,12,5,10,6,4,6 ,13,7,6,6,25,4,9,96,6,6,6,6,9,4,5,5,4,10,5,7,6) > day.hos3<-c(5,6,7,6,4,5,6,6,6,6,19,7,5,9,8,8,7,5,6,20,40,5,8,7,7,5,6,13, 11,9,4,6,9,16,6,7,6) > f<-ecdf(day.hos2) > plot(f,col.p='red',col.h='red') >

Hi all, given this example #start a<-c(0,70,50,100,70,650,1300,6900,1780,4930,1120,700,190,940, 760,100,300,36270,5610,249680,1760,4040,164890,17230,75140,1870,22380,5890,2430) length(a) b<-c(0,0,10,30,50,440,1000,140,70,90,60,60,20,90,180,30,90, 3220,490,20790,290,740,5350,940,3910,0,640,850,260) length(b) out<-ks.test(log10(a+1),log10(b+1)) # max distance D

Hi, I am using ecdf-function and want to use the ecdf()-data-points for nls() as data-parameter. nls() expects 'list' or 'environment' as a data-type, knots(ecdf(mydata)) gives me 'numeric'. What should I do now? Thanks in advance - Jochen Here is the code: ################################################# # --- Fit --- # Gumbel-Dist-Function, cumulative,

I have a csv file that has approximately 50k rows. In the first value of each row, a file name is listed, and there are 162 different file names. At the end of each row, there is a number value. What I would like to be able to do is for the 162 different files (or we could call them categories), is compute the ecdf for the values within that category. Then plot the ecdf for each file on the

Hello! I want to plot a P-P plot. So I've implemented this function: ppplot <- function(x,dist,...) { pdf <- get(paste("p",dist,sep=""),mode="function"); x <- sort(x); plot( pdf(x,...), ecdf(x)(x)); } I have two questions: 1. Is it right to draw as reference line the following: xx <- pdf(x,...); yy <- ecdf(x)(x); l <- lm(

Dear all, I would like to calculate the complementary cumulative distribution function. As it is known, the ccdf is the 1-ecdf(X)==1-F(x). (From ?ecdf help is shown that ecdf returns a function which is the F(x) I would like to use that function inside my function and after call curve (which accepts function as input). My code looks like that onVector<-seq(1:100)

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

Dear R users, Quick question on the use of ecdf for producing a cumulative distribution function. After having produced a cumulative distribution using ecdf, I would like to extract the exact x-value for a given percentile. For example I can ascertain the 1st and 3rd quartile using the summary function: cawa.cdp <- ecdf(cawaocc$LEFF80) plot(cawa.cdp, ylab="Fn(x)")

Hi all R users I would like to know how acess the coordinates of the ECDF object. I look for the example, in this part: ###################### print(ls.Fn12 <- ls(env= environment(Fn12))) ###################### but I do not know to extract the Y coordinate and put it in other variable. My objective is to make a plot and identify the points with labels. ############# Example by

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