similar to: question on cdf compare in R

Dear R-helper, I am working with ROCR of Tobias Sing et. al. to compare the performances of logistic and nnet models on a binary response. I had the performance plots, but I have problem finding out other performance statistics (eg. MSE/ASE, AUC). Any help on this? Thanks Ilham [[alternative HTML version deleted]]

Dear list, I want to combine several plots in one graph. I did this: plot(a1); plot(a2, add=TRUE); ...plot(a5, add=TRUE) The problem is the more plot we put, the more complex the graph. Is there any way to label each line; or other way just to make sure I know which one which? Thank you for the help, Ilham [[alternative HTML version deleted]]

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 everybody, I tried to analyze a custom Affymetrix 3'-biased Array. So I wanted to make a cdf package. (My CDF file size is 1.12Go). I tried several methods but the same error occured Method 1 > #Set the working directory > setwd("D:/Analyse R/Cel files") > #library to create cdf env > library("makecdfenv") >#Create cdf environment >pkgpath

Hi R help, I run my data in nnet with skip layer, factor response (with 0 & 1 values) and explicitly put softmax=T to compare the result of the default nnet with no softmax specification. I assume this should give me the same result. I got the result the default one, but not the softmax version and I got the error message that I did not quite understand. test6.nn.skipT.softm.Yfac <-

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

> Almost any method of fitting a density estimate would work on > integrating (numerically) the result. it is a nice idea concerning the monotony property, which will be obtained automatically, but I am going to use results of approximation analytically > In particular, look at package polspline, where > p(old)logspline does the integration for you. thank you, I am going to

Hi, I have a 324X15 matrix (No of years vs. heavy precipitation days) and I want to calculate the cdf at 5 different data points for each row. I tried by the following codes but it's not working. heavyprec <- read.csv (file="heavyprecdays_CSV.csv",header=TRUE,sep=",") a <- heavyprec pdf <- density (a, bw="SJ", kern= "gaussian") f <-

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

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, Ia m working in discretized data. Here my data: x <- c(2,1,3, 5), and I want to make (0,1) data based on the length of each component in x. So the new data should like: y = (0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1). I spent too much time with "seq", "rep". Still didn't get it. Any help? Thanks Ilham [[alternative HTML version deleted]]

This question might seem silly, because I felt that it MUST be in the mailing list archives or help files somewhere, but I simply couldn't find it. I want to make some simple CDF (cumulative distribution function) plots to check whether distributions are Gaussian / normal. But in order to check how "normal" the distribution is, I really need the y-axis to be Gaussian as well

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

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, 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:

Dear R-list, I want to combine several plots in one graph. I did this: plot(a1); plot(a2, add=TRUE); ...plot(a5, add=TRUE) The problem is the more plot we put, the more complex the graph. Is there any way to label each line; or other way just to make sure I know which one which? Thanks, Ilham [[alternative HTML version deleted]]

Dear all, I am calculating the bivariate skew normal cdf in "sn" package using "pmsn" function. Although it is quite convenient ( thanks to prof. Azzalini) but it seems to be slow. For example, it takes about 1 minute in calculation of 100k of such cdf values. I am thinking to write a c++ code for this although not very familiar with it. Any other idea?    Thanks in advance,

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.