similar to: plotting least-squares regression against x-axis

I've got a data frame describing comments on an electronic journal, wherein each row is a unique comment, like so: commentID author articleID 1 1 smith 2 2 2 jones 3 3 3 andrews 2 4 4 jones 1 5 5 johnson 3 6 6 smith 2 I want know the number of unique authors per article. I can get a

Dear R helpers,   I have two queries.   (1) If the dataset contains some variables having all the entries = 0 and while analysing I want to delete those pericular columns, how do acheive this. i.e.   dataset1   sr_no    var1      var2       var3        var4        var5 1           5           0          3              1           0 2           3           0          2              9           0

I tried to use BATCH command to run an input file (a simple text file) which contained a few R commands. I first saw the BATCH help page using      > ?BATCH The help page mentioned the Usage as: Batch Execution of RR CMD BATCH [options] infile [outfile] The Detailed section mentioned the following: Use R CMD BATCH --help to be reminded of the usage. My input file is called batchchk, which

I need to detect outliers in a large data set which is highly right-skewed. I plan to use medcouple-based outlier detection. Is there any support for medcouple-based outlier detection in R? Are there any other routines in R to perform outlier detection in highly right-skewed data? Manuj Sharma See the Web's breaking stories, chosen by people like you. Check out Yahoo! Buzz.

I need to run a set of R commands in batch mode. Can somebody please provide pointers to information on how to write R scripts to execute commands in batch mode? Thanks, Manuj Try the new Yahoo! India Homepage. Click here. http://in.yahoo.com/trynew [[alternative HTML version deleted]]

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

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

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

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

Dear all I am trying to fit Variable Coefficients Models on Unbalanced Panel Data. I managed to fit such models on balanced panel data (the example from the "plm" vignette), but I failed to do so on my real, unbalanced panel data. I can reproduce the error on a modified example from the vignette: > require(plm) > data("Hedonic") > Hed <- pvcm(mv ~ crim + zn + indus

Hello to everybody, I'd like to submit a problem I'm dealing with, and I can't get an answer to by myself. I have to test if my data come from a specific probability distribution, of which I know the analytic form both of the p.d.f. and the c.d.f. Namely, it is the hypoexponential distribution, sum of two exponentials with different parameters. Is there any way I can manage the

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

The need comes from the PISA data. (http://www.pisa.oecd.org) In the data there are many cases and each of them carries a numeric variable that signifies it's weight. In SPSS the command would be "WEIGHT BY" In simpler words here is an R sample ( What is get VS what i want to get ) > data.recieved <- data.frame( + kindergarten_attendance = factor(c(2,1,1,1), labels =

Hi, I've got two vectors with ratings from two coders, like this: x<-c("red", "yellow", "blue", "red") #coder number 1 y<-c("red", "blue", "blue", "red") #coder number 2 I want to find Cohen's Kappa using the wkappa function in the psych package. The only example in the docs is using a matrix, which

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

> 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

On Fri, Jun 26, 2009 at 10:14 AM, Brent Jones<brent at servuhome.net> wrote: > On Thu, Jun 25, 2009 at 12:00 AM, James Lever<j at jamver.id.au> wrote: >> >> On 25/06/2009, at 4:38 PM, John Ryan wrote: >> >>> Can I ask the same question - does anyone know when the 113 build will >>> show up on pkg.opensolaris.org/dev ? >> >> On