similar to: ROC curve using epicalc (after logistic regression) (re-sent)

Dear R-help list, I'm attempting to use the ROC routine from the epicalc package after performing a logistic regression analysis. My code is included after the sessionInfo() result. The datafile (GasketMelt1.csv) is attached. I updated both R and the epicalc packages and tried again before sending this request. sessionInfo result: R version 2.9.1 (2009-06-26) i386-pc-mingw32 locale:

Hello list, I am trying to avoid "Rifying" some of my SAS code to generate ROC plots, and the logistic.display() and lroc() functions in the epicalc package do what I want. However, I must generate my logistic model with bigglm because I have 1) limited hardware, 2) ~2.5 million rows, and 4 categorical and 2 continuous independent variables. When I attempt to invoke epicalc's

Dear R-friends, ? I am using the epicalc package and the manual by V. Chongsuvivatwong "Analysis of epidemiological data using R and Epicalc" to get the hang of some basic epidemiological analyses.??? ? After running all the analyses of chapter 7, one is supposed to wrap it up by saving the data writing: ? ? > save(.data, file = "Chapter7.Rdata") ? ...?after writing the

Hi everyone, I'm working with the Epical library, specicatly using the power test in proportions. I think this test is not working like in the book: Biostatistical Analysis (4th Edition): Jerrold H. Zar In the example 23.25. (I attach this Pic) It's not the same answer. Using the follow command don't give the same answer. library(epicalc) power.for.2p(0.75, 0.50, 50, 45, alpha =

I know how to compute the ROC curve and the empirical AUC from the logistic regression after fitting the model. But here is my question, how can I compute the standard error for the AUC estimator resulting form logistic regression? The variance should be more complicated than AUC based on known test results. Does anybody know a reference on this problem? [[alternative HTML version deleted]]

Hi everyone, I'm studying the manual name: Analysis of Epidemiological Data Using R and Epicalc, maked by: Virasakdi Chongsuvivatwong and Edward McNeil. And I can't find the data base that they use in some examples, this are the names: Chapter7.Rdata,Chapter8.Rdata,Chapter9.Rdata Somebody can tell me how can I have this files? Thk! Jos? O__ ---- Jos? Bustos M. c/ /'_ ---

Hi, I'm trying to set up a prediction software, now i testing the performance of my method, so i need to calculate a ROC curve, specially auc, cut-off, sens and spec, i just looking at ROCH package, but it's a mass for me, i'm not a math guy and I'm getting lost Could any of you recommend me an easy-to-use package to do this task? i just have a list of positive/negative samples

Hi, I use the epicalc package which provides the function aggregate.numeric. Unfortunately aggregate.numeric produces warnings when aggregate is used by functions not under my control on a numeric value. If I don't load epicalc, aggregate.default is used instead by these functions and does not produce any warning. However I need epicalc. So to get around this, what I would do is firstly

On 10 October 2009 at 05:25, lmramba at gmail.com wrote: | Dear all, | I installed my R-2.9.2 on my ubuntu version 9.04 successfully using= the | command |=20 | sudo apt-get install r-base-dev |=20 | The problem is that I cannot install any package. Huh? This has _nothing_ to do with the error below. =20 | See my details below: |=20 |=20 |=20 | > install.packages("epicalc") |

Buenas noches Javier y José, Estoy en contra de usar attach(), asi que propongo la siguiente alternativa con with(): # paquete require(epicalc) # los argumentos en ... pasan de epicalc:::cc # ver ?cc para mas informacion foo <- function(var1, var2, var3, ...){ or1 <- cc(var1, var2, ...) or2 <- cc(var1, var3, ...) list(or1 = or1, or2 = or2) } # datos x <-

Dear R Users; I am an R user who has recently bought a new laptop;Toshiba Satellite U405 running on both Windows Vista and Ubuntu. I have problems on the wondows vista when installing packages. The argument 'lib' is missing. How do i solve this problem? Illustration: > install.packages("epicalc") Warning in install.packages("epicalc") : argument 'lib'

Dear all, I installed my R-2.9.2 on my ubuntu version 9.04 successfully using the command sudo apt-get install r-base-dev The problem is that I cannot install any package. See my details below: > install.packages("epicalc") Warning in install.packages("epicalc") : argument 'lib' is missing: using '/home/lmramba/R/i486-pc-linux-gnu-library/2.9'

Dear R-help, I am attempting to plot data using standard R plot utilities. The data was retrieved from FRED (St. Louis Federal Reserve) using the package quantmod. My question is NOT about quantmod. While I retrieve data using quantmod, I am not using its charting utility. I have been having success using the standard R "plot" utilities to this point with this type of data.

Estimados Estoy entrenando hacer funciones que respondan a comandos, en esta caso en la salida gráfica se observa que dice : Exposure=var3 y outcome=var 1 quisiéramos que se reflejan los nombres de la base de datos : var1=estado, var2=cake, var3=chocolate Espero haberme explicado adecuadamente Adjunto tabla con datos #################################### #Comando que llama

Dear all, I get an error shown below when i run the command "use" found in package epicalc on R2.7.1 and R2.7.0 but dont get the error when i run it on R 2.6.2 What is likely to be the problem? library(epicalc) Loading required package: foreign Loading required package: survival Loading required package: splines use("malvac.dta") Error in fromchar(x) : character string

Hola José A. Entiendo que cargas lso datos del fichero con algo parecido a esto.... misdatos=read.table("sanda.csv",sep=",",header=TRUE) luego, si no quieres referirte a las variables así misdatos$IMC tendrás que hacer un attach(misdatos) por lo que ya funcionarían cosas como summary(IMC) [en vez de sumemry(misdatos$IMC) en tu caso... require(epicalc)

José, La función survfit espera como argumento una fórmula. Al pie del help de la función survfit dice: "Versiones anteriores del código también permiten la especificación de una sola curva al omitir el lado derecho de la fórmula, es decir, ~ 1. El manejo de este caso requiere algunas manipulaciones poco comunes y frágiles, este caso ya no es compatible." (sepan disculpar la

En el libro EPICALC (pagina 229-230) en el que está el siguiente script, todo nos funciona bien, pero cuando vamos a life table, ya allí no avanza, lo señalamos en el script, por favor quizá se nos haya ido algún detalle, pero fuimos siguiéndolo por el libro paso a paso y no no hemos percatado Todos los de el paquete survival de la ayuda del R funcionan perfectamente

Hi, I need to build ROC curve in R, can you please provide data steps / code or guide me through it. Thanks and Regards Rithesh M Mohan [[alternative HTML version deleted]]

Hola, El problema está con el nombre de las columnas y las filas del fichero que importas. Mira este detalle tras hacer la importación: > names(mort) [1] "X" "X15.24..79." "X25.34..79." [4] "X35.44..79." "X45.54..79." "X55.64..79." [7] "X65.74..79." "X75.84..79."