similar to: Sort in ecdf

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

Full_Name: Martin Frith Version: R-2.0.0 OS: linux-gnu Submission from: (NULL) (134.160.83.73) I have large vectors containing 100,000 to 20,000,000 numbers. However, they only contain a few hundred *distinct* numbers (e.g. positive integers < 200). When I do ecdf(v), it either runs out of memory, or it succeeds, but when I plot the ecdf with postscript, the output is unnecessarily bloated

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

I would like you consider that the function ecdf could be extended in the following way to handle weights when computing Empirical distribution Functions. There exist particular cases that supports this kind of extension, see for example: Rao, C. R., 1997. Statistic and True. Putting chance to work. World Scientific Publishing. Cox, D. R., 1969. Some Sampling Problems in Technology. New

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

Dear all, I have run into a problem when running some code implemented in the Bioconductor panp-package (applied to my own expression data), whereby gene expression values of known true negative probesets (x) are interpolated onto present/absent p-values (y) between 0 and 1 using the *approxfun - function*{stats}; when I have used R version 2.8, everything had worked fine, however, after updating

Hi It is usually not recommended but if you insist maybe library(plotrix) ?twoord.plot twoord.plot(lx=D[,1],ly=D[,2], rx=D[,1], ry=D[,3]) or plot.yy(x=D[,1],yright=D[,3], yleft=D[,2]) which allows only one x axis (see below). Cheers Petr plot.yy <- function (x, yright, yleft, yleftlim = NULL, yrightlim = NULL, xlab = NULL, yylab = list(NA, NA), pch = c(1, 2), col = c(1,2), linky

Thank you very much peter. It worked out nicely. I have additional question. How can I get Y-axis on log-scale? Thank you very much in Advance, Eliza UoS PP ________________________________ From: PIKAL Petr <petr.pikal at precheza.cz> Sent: 23 November 2017 16:22:39 To: Eliza Botto; r-help at r-project.org Subject: RE: adding percentage secondary y-axis Hi It is usually not

#See David Williams' book "Weighing the odds", p286 y <- c(1.21, 0.51, 0.14, 1.62, -0.8, 0.72, -1.71, 0.84, 0.02, -0.12) ybar <- mean(y) ylength <- length(y) ybarv <- rep(ybar, ylength) x <- 1:ylength plot(x,y,asp=1,xlab="position",ylab="ybar",type="n",ylim=c(-1,1)) segments(x[1], ybar, x[ylength], ybar) segments(x,ybarv,x,y)

Dear R-devel, The following line crashes R > approx(1, 1, 0, method='const', rule=2, f=0, yleft=NULL, ties='ordered')$y Process R:2 exited abnormally with code 5 at Tue Jul 21 14:18:09 2009 > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 9.1 year

I am working on a project where we start with start with 2 long, equal-length vectors, massage them in various ways, and end up with a function mapping one interval to another. I'll call that function "f1." The last step in R is to generate f1 as the value of the approxfun function. I would like to put f1 into a package, but without having the package redo the creation of

Hi people, I am new in this list and could not find a FAQ for it in particular, furthermore I could not find my question answered in the official R FAQ or docs. I have simply something like this: > f<-approxfun(data[,1],data[,2]) and f is: > f function (v) .C("R_approx", as.double(x), as.double(y), as.integer(n), xout = as.double(v), as.integer(length(v)),

Readers, According to the help '?approxfun', the function can be used to obtain the interpolated values. The following test was tried: > testinterpolation<-read.csv('test.csv',header=FALSE) > testinterpolation V1 V2 1 10 2 2 20 NA 3 30 5 4 40 7 5 50 NA 6 60 NA 7 70 2 8 80 6 9 90 9 10 100 NA >

Petr Pikal said: > I am not really a R specialist but for this task I use function: and he pasted his code into the email. I reindented the code, and wrote a fragment to experiment with it. Here it is: --------------------------------------------------------------------------- plot.yy <- function(x, yright, yleft, yleftlim=NULL, yrightlim = NULL,

Dear useRs, I have this dataset (D) with three columns. > dput(D) structure(c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 2.990484802, 3.005018792, 3.019552781, 3.03408677, 3.048620759, 3.063154749, 3.077688738, 3.092222727, 3.106756717, 3.121290706, 3.135824695, 3.150358684, 3.164892674, 3.179426663, 3.193960652, 3.208494642, 3.223028631, 3.23756262,

How do I plot using multiple(2) y-axes? I have two series that use the same x-data, but have very different scales. Appreciate any feedback, Per Bak

hi all i have another probably simple question. I have three variables say x, y and z. x and y are quite large and z is relative small. how can one plot the three variables on the same graph with two separate axis? (one for x and y and the other for z) e.g. x<-c(101,110,150,167,120) y<-c(120,135,175,95,200) z<-c(0.001, 0.15, 0.6, 0.8, 1) regards Allan

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