[R] SSPIR problem

2006-04-29 Thread Fernando Saldanha
I am having a problem with the package SSPIR. The code below illustrates it. I keep getting the message: Error in y - f : non-conformable arrays. I tried to tweak the code below in many different ways, for example, substituting rbind for cbind, and sometimes I get a different error message, but I

[R] Where do I find Cohen´s kappa???

2006-04-29 Thread jones christian
Hello, I´m looking for a way to measure the goodness of fit of my model with Cohen´s Kappa (scaling between 0 and 1). The kappa function does not give the results I´m looking for. Heres the code: z-glm(x~y,binomial) kappa(z, exact = T) Does anyone know more? many thanks Christian --

Re: [R] Where do I find Cohen´s kappa???

2006-04-29 Thread ronggui
kappa(base) is used to Estimate the Condition Number,so it is NOT what you want. And several packages have function related to your question: cohen.kappa(concord)kappa reliability coefficient for nominal data kappa2(irr) Cohen's Kappa and weighted Kappa for two raters

[R] question on residuals df in weighted linear regression

2006-04-29 Thread Antonio, Fabio Di Narzo
Hi all. I have a doubt with weighted linear regression. I've noted that supplying integer weights gives different residuals df than giving 'double' weights. Here an example: x - 1:5 y - -1 + 2*x + rnorm(length(x))*0.1 y - c(y, x + rnorm(length(x))*0.1) dat - data.frame(x=rep(x,2), y = y)

Re: [R] How to get a grid behind a boxplot

2006-04-29 Thread Camarda, Carlo Giovanni
Hi, likely it's not the most elegant way, but it's straight-forward and the first which popped in my mind: a - rnorm(100) # either boxplot(a, col=white, border=white) my.grid - seq(-4, 4, 0.5) abline(h=my.grid, col=grey, lty=2) boxplot(a, col=red, add = TRUE) # or boxplot(a, col=white,

Re: [R] How to get a grid behind a boxplot

2006-04-29 Thread Neuro LeSuperHéros
Simply add the same boxplot over your first graph using par(new=TRUE) x -1:10 y -1:5 boxplot(x, y,col=blue) grid(nx=NA, ny=NULL) #grid over boxplot par(new=TRUE) boxplot(x, y,col=blue)#grid behind boxplot Neuro From: Toby Popenfoose [EMAIL PROTECTED] Reply-To: Toby Popenfoose [EMAIL

[R] Resolving conflicting function definitions {was Error in rm.outlier}

2006-04-29 Thread Martin Maechler
Marc == Marc Schwartz (via MN) [EMAIL PROTECTED] on Fri, 28 Apr 2006 14:38:53 -0500 writes: Marc Sachin, I don't have a definitive thought, but some Marc possibilities might be a conflict somewhere in your Marc environment with a local function or with one in the Marc

[R] new computer for running R, Coore Duo or Athlon X2? + Operating system choice

2006-04-29 Thread Jason Liao
Dear R experts here, I am considering a new PC mainly to run R simulation. I will buy a double core machine as I regularly run several R sessions at a time. I am choosing between Core Duo and Athlon X2. I would like your input. I also need to hook up R with a Fortran 90 subroutine of about 1000

Re: [R] Resolving conflicting function definitions {was Error in rm.outlier}

2006-04-29 Thread Marc Schwartz
On Sat, 2006-04-29 at 16:21 +0200, Martin Maechler wrote: Marc == Marc Schwartz (via MN) [EMAIL PROTECTED] on Fri, 28 Apr 2006 14:38:53 -0500 writes: Marc Sachin, I don't have a definitive thought, but some Marc possibilities might be a conflict somewhere in your Marc

Re: [R] help with box-tidwell

2006-04-29 Thread John Fox
Dear Xuemei Bian, This is, I guess, the box.tidwell() function in the car package. To know exactly what happened, it would help to see the data, but I suspect that you have negative values among the predictors. The log of a negative number is not defined, producing NaN (which means not a

[R] predicted possibilities for classification tree??

2006-04-29 Thread jones christian
Hello! I would like to use a classification tree as a species modelling method. The dependent variable is binary (0,1),the independent variables are numeric. So I have to surround the dependent variable with factor(). That tells tree() that I want a logistic, or classification, tree rather

[R] How to rep a matrix by row?

2006-04-29 Thread Jiantao Shi
Hi, i have a dataframe like this, SourceTreatDrugReplicate control0A1 control10A2 control30A3 10A1 And i want to rep this dataframe 3 times by row,the resulting matrix as follow, SourceTreatDrugReplicate control

Re: [R] How to rep a matrix by row?

2006-04-29 Thread Gabor Grothendieck
Try this where DF is your data frame: DF[rep(seq(nrow(DF)), each = 3), ] On 4/29/06, Jiantao Shi [EMAIL PROTECTED] wrote: Hi, i have a dataframe like this, SourceTreatDrugReplicate control0A1 control10A2 control30A3 10A1

[R] How to get predicted probabilities out of a class.tree

2006-04-29 Thread jones christian
Hello! I would like to use a classification tree as a species modelling method. The dependent variable is binary (0,1),the independent variables are numeric. So I have to surround the dependent variable with factor(). That tells tree() that I want a logistic, or classification, tree rather than

[R] splitting and saving a large dataframe

2006-04-29 Thread Dimitri Szerman
Hi, I searched for this in the mailing list, but found no results. I have a large dataframe ( dim(mydata)= 1297059 16, object.size(mydata= 145280576) ) , and I want to perform some calculations which can be done by a factor's levels, say, mydata$myfactor. So what I want is to split this

Re: [R] R-News submission

2006-04-29 Thread John Fox
Dear Roman, I've attached the latex file for the copy-edited version of your manuscript -- correcting English grammar and style, and making the paper conform more closely to the R News style. Can you check this over to confirm that I've not altered the meaning of your text and to approve the

[R] Making R talk to Win/OpenBUGS in Linux (again)

2006-04-29 Thread Paul Johnson
I'm back! I've just learned that, on a fully updated Fedora Core Linux5 sytem, the working solution to access Winbugs under wine via the R package rbugs no longer works. Here was my last post on this topic (with the formerly working solution) from January.

Re: [R] Reshaping genetic data from long to wide

2006-04-29 Thread Farrel Buchinsky
Gabor Grothendieck ggrothendieck at gmail.com writes: http://news.gmane.org/gmane.comp.lang.r.general or one of these: http://dir.gmane.org/gmane.comp.lang.r.general Yes but when I hit Post this article it send something to gMane (I think) but not to R-help@stat.math.ethz.ch

[R] Writing responses to the R-Help list

2006-04-29 Thread Farrel Buchinsky
A while back Gabor Grothendieck suggested that I try http://news.gmane.org/gmane.comp.lang.r.general. This was after I asked how to easily reply to posts on the listserve. Ideally I would like the functionality that I find in Microsoft Outlook Express newsreader for usenet groups or what I find in

Re: [R] Resolving conflicting function definitions {was Error in rm.outlier}

2006-04-29 Thread Richard M. Heiberger
I agree that conflicts() belongs in the posting guide. I devoted a small section in my book to that very topic. Here is the lead paragraph which you are welcome to use Inexplicable Error Messages In general, weird and inexplicable errors mean that there are masked function names. That's the

Re: [R] Reorganizing rows and columns

2006-04-29 Thread kewley
That did the job. thanks, akn Suppose your data frames are A and B: AB - merge(A, B, c(Doc, Query), all=TRUE) AB[order(AB$Query, AB$Doc),] gets the answer you are asking for. (Not sure why you want to sort it to use a correlation test, as those are indifferent to ordering.) On Sun, 23

Re: [R] Reshaping genetic data from long to wide

2006-04-29 Thread Gabor Grothendieck
You might want to check the double check the list archives https://www.stat.math.ethz.ch/pipermail/r-help/ to see if your posts got through or not just in case its just some problem in displaying your own posts. On 4/29/06, Farrel Buchinsky [EMAIL PROTECTED] wrote: Gabor Grothendieck