Re: [R] ROCR for combination of markers

2012-04-22 Thread suo
Hi R users, I am sort of new on those. In the following logistic regression function "tigol( 0.3 + 3*x + 5*z + 7*w ) ", why the intercept was set 0.3, and the coefficients were set 3 ,5 and 7? Thanks. -- View this message in context: http://r.789695.n4.nabble.com/ROCR-for-combination-of-marke

Re: [R] ROCR for combination of markers

2012-04-22 Thread suo
Hi Eik or other who might help: I got this error: "Error in roc.formula(form = y1 ~ x + z, plot = "ROC") : Invalid formula: exactly 1 predictor is required in a formula of type response~predictor. " when I ran "out=ROC( form = y1 ~ x + z, plot="ROC") " from your code. How to fix it? Thanks. -

Re: [R] ROCR for combination of markers

2011-04-28 Thread Rasanga Ruwanthi
Thanks Eik. As you said both packages give the same result, except labelling of the x axis. --- On Thu, 28/4/11, Eik Vettorazzi wrote: From: Eik Vettorazzi Subject: Re: [R] ROCR for combination of markers To: "Rasanga Ruwanthi" Cc: "R Help" Date: Thursday, 28 Apr

Re: [R] ROCR for combination of markers

2011-04-28 Thread Eik Vettorazzi
t for (1-sensitivity) vs specificity? > > Thanks > Rasanga > > > > --- On *Thu, 28/4/11, Eik Vettorazzi > //* wrote: > > > From: Eik Vettorazzi > Subject: Re: [R] ROCR for combination of markers > To: "Rasanga Ruwanthi" > Cc: &

Re: [R] ROCR for combination of markers

2011-04-28 Thread Eik Vettorazzi
ol( 0.3 + 3*x + 5*z + 7*w ) ) > out=ROC( form = y1 ~ x + z, plot="ROC",MI=FALSE) > > But this function does not produce SE or CI of the AUC or any other > statistics. Any suggestion to get these? > > Thanks again > Rasanga > > > --- On *Thu, 28/4/11, E

Re: [R] ROCR for combination of markers

2011-04-28 Thread Rasanga Ruwanthi
tigol( 0.3 + 3*x + 5*z + 7*w ) ) out=ROC( form = y1 ~ x + z, plot="ROC",MI=FALSE)   But this function does not produce SE or CI of the AUC or any other statistics. Any suggestion to get these?   Thanks again Rasanga --- On Thu, 28/4/11, Eik Vettorazzi wrote: From: Eik Vettorazzi

Re: [R] ROCR for combination of markers

2011-04-28 Thread Frank Harrell
Be careful not to use cutpoints at any stage. improveProb does not require this. ROC analysis on the other hand is dangerous in the sense that it tempts one to find cutpoints, which are not replicable or consistent with decision theory. Frank Eik Vettorazzi wrote: > > Hi Rasanga, > you may ha

Re: [R] ROCR for combination of markers

2011-04-28 Thread Eik Vettorazzi
... and additionally, 'ROC' from the Epi package does the second step all in one. Am 28.04.2011 13:01, schrieb Eik Vettorazzi: > Hi Rasanga, > you may have a look at the 'improveProb' function from the Hmisc > package. There you can compare the increase in prognostic power for > several combinatio

Re: [R] ROCR for combination of markers

2011-04-28 Thread Eik Vettorazzi
Hi Rasanga, you may have a look at the 'improveProb' function from the Hmisc package. There you can compare the increase in prognostic power for several combinations of markers. You can create a ROC curve for a combination of markers by using the predicted risks eg. from a logistic regression model

[R] ROCR for combination of markers

2011-04-27 Thread Rasanga Ruwanthi
Dear list   I have 5 markers that can be used to detect an infection in combination. Could you please advise me how to use functions in ROCR/ other package to produce the ROC curve for a combination of markers?   I have used the following to get ROC statistics for each marker. pred <- prediction(