similar to: Special Functions calculations

Forwarded Message: > To: Steve Arthur <sarthur at crick.protogene.com> > From: Thomas Lumley <thomas at biostat.washington.edu> > Subject: Re: [R] couldn't find function "..." > Date: Wed, 5 Jul 2000 16:16:25 -0700 (PDT) > ----- > On Wed, 5 Jul 2000, Steve Arthur wrote: > > > Hi, > > > > R 1.0.1, on Linux(mandrake/RedHat).

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

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

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") >

Hi, This may belong more to r-develop, but general discussion may be useful (for the how many-th time ?) seq(2,5,-2) seq(5,2,2) both result in Error in seq.default(2, 5, -2) : wrong sign in 'by' argument But often, if not always, mathematicians and programmers want a behaviour e.g. in for loops, where this statement results in an empty statement, that is for (ii in seq(2,5,-2))

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

try: library(foreign) read.ssd("J:\\QM\\Reports\\Sarthur\\SAS_Application\\SAS_Data_Sets","use") instead, hth, Merry Christmas, Bernhard -----Original Message----- From: Stephen Arthur [mailto:sarthur67 at yahoo.com] Sent: 20 December 2002 16:55 To: ripley at stats.ox.ac.uk; rossini at blindglobe.net Cc: r-help at stat.math.ethz.ch; stvjc at channing.harvard.edu Subject:

On Sun, Jul 16, 2017 at 06:02:15PM +0100, Pete Biggs wrote: > > > > The physicists and mathematicians who count there need high durations. > > Yes. I too run HPC clusters and I have had uptimes of over 1000 days - > clusters that are turned on when they are delivered and turned off when > they are obsolete. It is crucial for long running calculations that you > have a

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

Halo Valeri, let's think about what a hpc cluster is for. Second, one should always ask the question where security is to be applied,then one can come to the following decision: - The firewall is placed in front of the cluster. - After you have found a safe base for this, you freeze it. - We have a rsync of centos and epel on the head node. >From here, we can always reinstall a node

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)

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 all, I am not sure if this feature exists in the R-console command line prompt: In Matlab, if I want to enter a command which is similar to what I have entered before, I can enter a few prefix, then press "->", the previous command that matches with this prefix will then appear on this command line, and it saves a lot of our time. For example: > abline(lm(new~old)) >

Rer's, Get ready for this intellectual fest! Linux i686 2.2.14-1mdklinus R-1.0.1 When I try to q() from my workspace, I get the following responses: ############################### > q() Save workspace image? [y/n/c]: y Error: can't save data -- unable to open ./.RData then: > save.image() Error in save(list, file, ascii, oldstyle) : unable to open file and finally:

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

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