This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to automate a script using for?

Hey guys, All good?

I use this simple script to generate roc curves:

 fit1=glm(Class ~ ACTB, data = mydata, family='binomial')

 pred1 = predict(fit1)

roc1=roc(data$Class ~ pred1)

plot(roc1, col='black', print.auc=TRUE,

     print.auc.cex = 0.8,

      print.auc.y=0.1,

     print.auc.x=0.2)

My input table is as follows:

Column 1 contains the name of the samples, column 2 contains the class of the sample, and the other columns with the intensity of the proteins that I want to make the ROC curve.

I would like to do everything automatically, because the way I do it now, I need to do it one by one. Could someone help me, please?

Thanks in advance!!

roc curve

Please make the data reproducible, add example data via dput().

0 answers

No answers yet.

Log in to answer this question.