similar to: Question abou pROC package

Dear R users, pROC is a package to compare, visualize, and smooth receiver operating characteristic (ROC) curves. The package provides the following features: * Partial or full area under the curve (AUC) computation * Comparison of two ROC curves (curves and AUC) * Calculating and plotting confidence intervals * Smoothing of the ROC curve * Coordinates extraction ('coords' function).

Dear R users, pROC is a package to compare, visualize, and smooth receiver operating characteristic (ROC) curves. The package provides the following features: * Partial or full area under the curve (AUC) computation * Comparison of two ROC curves (curves and AUC) * Calculating and plotting confidence intervals * Smoothing of the ROC curve * Coordinates extraction ('coords' function).

Packages pROC and ROCR both calculate/approximate the Area Under (Receiver Operator) Curve. However the results are different. I am computing a new variable as a predictor for a label. The new variable is a (non-linear) function of a set of input values, and I'm checking how different parameter settings contribute to prediction. All my settings are predictive, but some are better. The AUC i

Dear useRs: I use pROC package to compute the bootstrap C.I. of AUC. The command was as follows: roc1<-roc(all$D,all$pre,ci=TRUE,boot.n=200) However, the result was: Area under the curve: 0.5903 95% CI: 0.479-0.7016 (DeLong) Why the C.I. was computed by the Delong Method? Yao Zhu Department of Urology Fudan University Shanghai Cancer Center Shanghai, China [[alternative HTML version

Dear R-philes, I am plotting ROC curves for several cross-validation runs of a classifier (using the function below). In addition to the average AUC, I am interested in obtaining a confidence interval for the average AUC. Is there a straightforward way to do this via the ROCR package? plot_roc_curve <- function(roc.dat, plt.title) { #print(str(vowel.ROC)) pred <-

Hi, I am using logistic regression to classify a binary psychometric data. using glm() and then predict.glm() i got the predicted odds ratio of the testing data. Next i am going to plot ROC curve for the analysis of my study. Now what i will do: 1. first select a cut-off (say 0.4) and classify the output of predict.glm() into {0,1} segment and then use it to draw ROC curve using ROCR package

I might be missing something but I thought that AUC was a measure for comparing ROC curves, so there is nothing else needed to "compare" them. The larger AUC is the higher correlation of 2 variables compared. No other measures or calculations are needed. Jarek Tuszynski -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On

Dear R-philes, I am having some trouble averaging across multiple runs of a classifier in an ROC Curve. I am using the ROCR package and the plot() method. First, I initialize a list with two elements where each element is a list of predictions and labels: vowel.ROC <- list(predictions=list(), labels=list()) For every run of the classifier, I append the scores and labels to their

Hi Abou, Try this: library(plotrix) curve(rescale(dnorm(x ,mean=mean(Lizard.tail.lengths),sd=sd(Lizard.tail.lengths)), c(0,6)),add=TRUE, col=2, lwd = 2) Jim On Mon, Sep 25, 2017 at 9:35 AM, AbouEl-Makarim Aboueissa <abouelmakarim1962 at gmail.com> wrote: > Dear All: > > One more thing. > > I want to add the normal curve to the histogram. Is there away to stretch > the

Dear all, I try to compare the performances of several parameters to diagnose lameness in dogs. I have several ROC curves from the same dataset. I plotted the ROC curves and calculated AUC with the ROCR package. I would like to compare the AUC. I used the following program I found on R-help archives : From: Bernardo Rangel Tura Date: Thu 16 Dec 2004 - 07:30:37 EST

Hello I am trying to make a graph of 10 different lines built each from 4 different segments and to add a darker line that will represent the average of all graphs - all in the same plot.Actually each line is a ROC plot The code I'm using for plotting one line is as follows: logit.roc.plot <- function(r, title="ROC curve") { old.par <- par(no.readonly = TRUE);

Hi, I have a dataset (see attached) with 2 variables "Y" is binary, "x" is a continuous variable. I want to calculate area under the curve (AUC) for the ROC curve, but I got different AUC values using ROC() from Epi package vs. rcorr.cens() from rms package: test<-read.table("test.txt",sep='\t',header=T,row.names=NULL) y<-test$y x<-test$x

Dearh all, I have following question: a script (using pROC functions) that works when run in Rgui, failes when run through rscript. This is the script: library(pROC) hits <- c("T", "D", "T", "D", "T", "D", "T", "D", "T", "D", "T", "D", "T", "D",

*Dear Bert, * Thank you very much for your feedback and the useful link https://rseek.org/ and https://www.r-bloggers.com/calculating-auc-the-area-under-a-roc-curve/. Actually, I want to know different performance between Stata and R, in multilevel logistic regression. For this purposes, I replicate ".do" file use Stata in

Hello I am trying to make a graph of 10 different lines built each from 4 different segments and to add a darker line that will represent the average of all graphs - all in the same plot.Actually each line is a ROC plot The code I'm using for plotting one line is as follows: logit.roc.plot <- function(r, title="ROC curve") { old.par <- par(no.readonly = TRUE);

Hello, I would like to use the ROCR package to draw ROC curves and compute AUC values. However, in the specific context of my application, the true positive rates and false positive rates are already provided by some upstream method. Of course, I can draw a ROC plot with the following command : plot(x=FPrate, y=TPrate, "o", xlab="false positive rate", ylab="true

Hello, I am comparing two ROC curves with bootstraping. However, some runs return "p-value = NA," and I have no clue why. Is this anyhow related to like sample size or no sufficient numbers of bootstraping? I used the default value (i.e. boot.n=2000), and the number of observations are quite big since I am comparing maps (e.g., the largest has more than 9 million observations).

Hi Brian, Your underlying dataset for the ROC curve only has 4 unique values for specificity, even though there are 23 elements in the vector, hence the step function nature of the first plot. The default smoothing in the smooth() function is "binormal". You might try one of the other smoothing options to see the result and whether they make visual sense. In the absence of smoothing,

Hi everybody, I am having a problem building a ROC curve with my data using the ROCR package. I have 10 lists of proteins such as attached (proteinlist.xls). each of the lists was calculated with a different p-value. The goal is to find the optimal p-value for the highest number of true positives as well as lowaest number of false positives. As far as I understood the explanations from the

I'm using randomSurvivalForest to predict survival from a rather small sample. As it's not enough data to have training and validation set, I rely on the "Estimate of error rate" computed by the randomForest. If I understand the method correctly, it repeatedly partitions the data into varying training/validation sets during the learning steps, which also produces the