similar to: [Plea to the R Gods] Theoretical and Empirical CDFs

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

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

> x <- runif(100) > y <- runif(100) > ks.test(x,y) Two-sample Kolmogorov-Smirnov test data: x and y D = 0.11, p-value = 0.5806 alternative hypothesis: two-sided ok I expected that, but: > ks.test(runif(100), "runif") One-sample Kolmogorov-Smirnov test data: runif(100) D = 0.9106, p-value < 2.2e-16 alternative hypothesis: two-sided How

I have a small LAN in my basement, with my PC running RH5 Linux and a Win95 PC connected by Samba. Everything works fine, but when I smbmount a CD (burned with xcdroast) from the Linux box (CD is in the Win95 machine) I do not get the proper filenames. All files are uppercase, and the TRANS.TBL files appear to be ignored. I have tried playing with name mangling, but it didn't help. Is this a

Hi, I have to generate 10 cdfs in a graph. I need to compare the cdf's nature by plotting ten cdfs in a graph. Thus, I need multiple plots in a graph. I would appreciate if you could give some solution to the problem asap. Thanking you, Sincerely, Ajay. -- Ajay Singh Research Scientist, SOM, IIT-Bombay, Powai, MUMBAI-400076, MH (INDIA).

Hello I just downloaded a great program called Transcribe (to slow down music so you can learn it) available at http://www.seventhstring.demon.co.uk/ and it worked smoothly and without problem through Wine on my Gentoo 1.4 laptop. The problem is that the program wants to open .cda files to access the CD tracks and I don't know about a way to mount audio CDs in Linux. I could of course rip the

Hi, I want to plot the residuals of a least-squares regression. plot(lm(y~x), which=1) does this, but it plots the y-axis of my data on the x-axis of the residuals plot. That is, it plots the residual for each y-value in the data. Can I instead use the x-axis of my data as the x-axis of the residuals plot, showing the residual for a given x? Thanks! Jason Priem University of North

Hi, some time ago I tried to boot an isolinux boot loader from another one to be able to boot e.g. a Knoppix CD with its original boot loader and menus from a "main" isolinux on a multiboot DVD. This works for isolinux <= 3.07 now! However, newer isolinux boot loaders have some problems and the isolinux (e.g. of the Knoppix CD) loaded from the multiboot DVD's main isolinux

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 want to generate random samples from truncated normal say Normal(0,1)Indicator((0,1),(2,4)). It has more than one intervals. In the library msm, it seems to me that the 'lower' and 'upper' arguments can only be a number. I tried rtnorm(1,mean=0,sd=1, lower=c(0,2),upper=c(1,4)) and it didn't work. Can you tell me how I can do truncated normal at more than one intervals?

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

Hello r-help, I notice that r includes distribution, and quantile values for the Chi-Square, F, t, and normal functions. Could someone give me some good reference(s) for how these probabilities and quantiles are calculated. I have already calculated the normal probability values using a Taylor expansion (http://www.sisweb.com/math/tables.htm). I realize these calculations must be very

Dear all Does anybody know of updated Custom CDFs with binaries for MacOS. I'm trying to do affy-analysis using R-studio for Mac. Thanks Petter

Hi everyone, I have pushed out Syslinux 3.81 to the usual places. Again, special thanks to my employer, Intel, for letting me spend time working on this stuff ;) -hpa Syslinux 3.81 is primarily a bug fix release. The main new features are halt-on-idle support and some features to isohybrid, including the ability to boot from an image written to a partition. Changes in 3.81: *

Hi, I have just started using R. Now I have the following problem: I want to create an Empirical Cumulative Distribution Function and I only came so far: F10 <- ecdf(x) plot(F10, verticals= TRUE, do.p = TRUE, lwd=3) x=c(1.6,1.8,2.4,2.7,2.9,3.3,3.4,3.4,4,5.2) Now I'd like to use arguments such as xlabs and main but I don't know how to integrate them. I hope someone can help me, I am

I am trying to do a MLE fit of the weibull to some data, which I attach. fitweibull<-function() { rt<-scan("r/rt/data2/triam1.dat") rt<-sort(rt) plot(rt,ppoints(rt)) a<-9 b<-.27 fn<-function(p) -sum( log(dweibull(rt,p[1],p[2])) ) cat("starting -log like=",fn(c(a,b)),"\n") out<-nlm(fn,p=c(a,b), hessian=TRUE)

3.80 is a bug fix and internal restructuring release. In particular, the entire shuffler subsystem has been rewritten largely from scratch, as has the Multiboot module. Changes in 3.80: * New shuffler mechanism and API. * Rewritten mboot.c32 module. * The syslinux_boot_linux() function has been simplified. * Don't hang trying to boot a "menu

Dear Masters, I have a question to submit consider the following script m<-4.95 obs<-rpois(36,m) # i generate 36 realization from a poisson(m) hist(obs,freq=F) curve(dpois(x,m),add=T,col="red") #i wish to overlay on the histogram the theorical poisson density function errors are returned saing the x vector doesn't contain integers.... really bizarre i can't give

Hi, I would like to know, if there's any R function, which allows calculation of probability values (0,1) from multivariate normal densities. I would be grateful for any output. Cheers, MG

I am having consistent difficulty with write access and installing applications. The errors show up as write protect errors or not enough disk space errors. 'winefile' shows the disks as being NTFS. Is this correct or is it a winefile anomaly? Are these conditions related? What is the correct file system type? How can this be managed? 'winecfg' does not seem to allow