search for: epicalc

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....

...il, and for any inconvenience. I have also sent the question to the package maintainer, as suggested in the posting guide. Regards, Cliff ---------- Forwarded message ---------- From: Clifford Long <gnolffilc at gmail.com> Date: Sun, Jul 26, 2009 at 8:46 PM Subject: Fwd: ROC curve using epicalc (after logistic regression) To: cvirasak at medicine.psu.ac.th Dear Virasakdi Chongsuvivatwong, After sending the message below to the R-help mailing list, it occurred to me that I probably should also have sent a copy to you, per R posting guidance. I would be interested in any thoughts or sug...

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-...

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, fi...

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 <- read.csv("~/Downloads/OR.csv") head(x) # resultados SIN g...

...ibrary, 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 = 0.05) with this command the Power = 0.6402. But in the Book we find Power=0.72. ?what it's the correct answer?...?Zar or Chongsuvivatwong? Thank guys! O__ ---- Jos? Bustos M. c/ /'_ --- Master Apllied Stat Program (*) \(*) -- University...

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 "lroc" function on the bigglm model R returns the err...

...e de datos : var1=estado, var2=cake, var3=chocolate Espero haberme explicado adecuadamente Adjunto tabla con datos #################################### #Comando que llama a una función rm(list=ls()) #setwd("D:/DEMO_new/demo_scripts/OR/") #setwd("D:/Public/Documents/R/EPICALC/funciones/OR/") data= mydata<-read.csv("OR.csv",header=TRUE, sep=",", dec=".") use(data) attach(data) var1=estado var2=cake var3=chocolate library(epicalc) source("function_or.r") odratios(data,var1,var2,var3) #función odratios <...

...rs; 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' is missing: using 'C:\Users\lazurus\Documents/R/win-library/2.7' --- Please select a CRAN mirror for use in this session --- trying URL 'http://cran.za.r-project.org/bin/windows/contrib/2.7/epicalc_2.7.1...

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' --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done Error in read.dcf(file = tmpf) : Line starting &...

Dear All, I installed R 2.7.0 and tried to call a dataset i had ealier own called on R2.6.2 but i keep on getting an error: use("maltreat.dta") Error in fromchar(x) : character string is not in a standard unambiguous format Tried doing the same with R2.7.1 but i get the same error. However if i call the same on R 2.6.2, there is no error: use("maltreat.dta") > des()

...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") | Warning in install.packages("epicalc") : | argument 'lib' is missing: using | '/home/lmramba/R/i486-pc-linux-gnu-library/2.9' This is suspicuous.=20 If you used the Ubuntu package, it would use /usr/local/lib/R/site-libr= ary by default. So if you have a loca...

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...

#Imagino que hay un camino más corto para resolver mi duda #Adjunto unos datos en csv library (epicalc) #Gráficos exploratorios #este comando no se ejecuta porque no ve las variables summ(IMC, by=escolaridad ) ##?porqué no ve las variables? #tuve que hacer esto que me parece que se podría hacer de manera más sencilla IMC <- (sanda[,c(1)]) escolaridad <- (sanda[,c(2)])...

Hi all: My data has 3 variables: age(3levels : <30y=1 30-50y=2, >50y=3) gender(Male=0, Female=1) CD4 cell count(raw lab measurement) y(1:death 0:alive) I perform logistic regression to find out the factors that influence y. result<-glm(y ~ factor(age) + factor(gender) + CD4,family = binomial) >From the result,I can get OR(Odds Ratio) of gender via exp(Estimate of Female,

Estimados Cuando existia epicalc, había una manera muy fácil de determinar la media de una variable (en esta caso Gain) por grupos, en este caso (Diet). ?Como se puede hacer ahora? Diet Gain 1 d1 270 2 d1 300 3 d1 280 4 d1 280 5 d1 270 6 d2 290 7 d2 250 8 d2 280 9 d2 290 10 d2 280 1...

Estimados ¿Existe en R alguna manera de que las tablas de estadística descriptiva como las de codebook de epicalc, salgan en una tabla como las de SPSS ? Saludos José -- Nunca digas nunca, di mejor: gracias, permiso, disculpe. Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La perso...

Hi, I am trying to calculate confidence intervals using ci.numeric from epicalc package. If I generate a normal set of data and find the 99% and 95% CI, they seem too narrow to me. Am I doing something wrong?? The IQR goes from -0.62 to 0.62, so I thought the CI limits should be more extreme than these values. x<- rnorm(200,0,1) ci.numeric(x=mean(x),n=200,sds=sd(x),alpha=...

...> > survfit(surv.marr ~ 1) > > Repito no se si tiene sentido, quizás tú lo encuentres o alguien más pueda > opinar sobre este uso en este caso. > > Daniel Merino > > > 2013/10/4 Dr. José A. Betancourt B. <jbetancourt en iscmc.cmw.sld.cu> > >> 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 percata...

...argo cuando mandamos a leer desde la tabla csv original de los autores no hace el análisis porque algo falta y no nos percatamos de que es. Adjunto tabla original Saludos cordiales #ESTO ES LO QUE VIENE EN LA AYUDA Y TRABAJA BIEN rm(list = ls()) library(FactoMineR) data(mortality) library(epicalc) des(mortality) #Agregue esto para comparar las dos datas res<-MFA(mortality,group=c(9,9),type=c("f","f"), # REALIZA EL ANÁLISIS BIEN name.group=c("1979","2006")) plot(res,choix="freq",invisible="ind",habillage="group...