search for: epi

Hello, I was wondering who to contact to see about somehow interconnecting or integrating the programs R and Epi Info. I have been involved in the development of Epi Info which is a software program developed by the Centers for Disease Control and Prevention (CDC) in Atlanta (www.cdc.gov/epiinfo<http://www.cdc.gov/epiinfo>). Epi Info allows for creating data entry systems and performs certain epidemio...

?? Hola. Estoy intentando ordenar el gráfico tal que las columnas se apilen en el siguiente orden: debajo, los EPI; los SemiEPI en el medio y los NoEpi arriba de todo. Estos son mis datos: > tabla.barras Rangos Clase Asistentes 1 De 45d EPI 405 2 De 1a EPI 812 3 De 2a EPI 1639 4 De 3a EPI 2202 5 De 4a EPI 1994 6 De 5a EPI 255 7 D...

Hola, Mira la versión de ggplot2 que estás usando por tenerla actualizada a la última... A mi me funciona introduciendo algún cambio salvo la última línea cuando cambias el tema... #----------- Lines <- " id Rangos Clase Asistentes 1 45d EPI 405 2 1a EPI 812 3 2a EPI 1639 4 3a EPI 2202 5 4a EPI 1994 6 5a EPI 255 7 45d SemiEPI 614 8 1a SemiEPI 881 9 2a SemiEPI 2174 10 3a SemiEPI 3254 11 4a SemiEPI 2934 12 5a SemiEPI...

The name conflicts in Epi and ROC packages (2 'ROC' functions are the problem) cause the following code to work once, but not twice: library(MASS); data(cats); x = cats[,2] y = ifelse(cats[,1]=='F',0,1) library(Epi); ROC(x,y,grid=0)$AUC library(ROC); AUC(rocdemo.sca(y, x, dxrule.sca)) What i...

Dear all, I would like to calculate the optimal cut off (threshold) of a test using the Epi package. Here I am presenting some data based on the output of two tests. I am interested in identifying the optimal cut off an its 95% CI. Running the ROC() function with the Epi package I obtain a nice picture that returns what I interpret as the optimal cut off with lr.eta=0.431. would be this...

This is a really simple question, I'm sure,but I can't make EpiR work! I keep getting the following: > epi.2by2(47, 263483, 282, 935028, method="cohort.time", conf.level=0.95) Error in epi.2by2(47, 263483, 282, 935028, method = "cohort.time", conf.level = 0.95) : unused argument(s) (935028) and I really don't know why!. Any ide...

Hello: Does anybody knows if Epi Info, from CDC (www.cdc.gov/epiinfo) runs on Wine? I did try AppDB first but there are no entries for Epi Info. Epi Info is a data base for epidemiology and public-heath developed by Centers for Disease Control and Prevention (USA) and heavilly used in the world. It is based on Access 97 or 2000...

...r linux. Here are the commands I have been using: postscript("crocheinterp.ps", family="Times", height=7, width=4, horizontal=FALSE, paper="letter") par(mfcol=c(2,2)) par(mfg=c(1,1,2,2)) par(mai=c(0.5, 0.5, 0.5, 0.5)) plot(croche.day.epi.grid$xy, type="n", xlim=c(0,100), ylim=c(0,200), xlab="", ylab="") image(croche.day.epi.grid$x, croche.day.epi.grid$y, matrix(croche.day.epi.grid$krige$zhat, length(croche.day.epi.grid$x), length(croche.day.epi.grid$y)), add=T, col=grays)...

...matrix of the coefficients. I know how to get corr from survreg, say, summary(f,corr=T), but how to get covariance matrix of the coefficient ? can I put them in a matrix and output ? thanks below is what I did <mailto:r-help@stat.math.ethz.ch> n=100; beta=1.0 z=runif(n) u=runif(n) epi=log(u/(1-u)) d=rep(1,n) t=z*beta+epi f<-survreg(Surv(t,d)~z+epi, dist="logistic") summary(f, corr=T) below is the result > summary(f, corr=T) ....... ........ Correlation of Coefficients: (Intercept) z epi z -0.882 epi 0.16...

There is now an Epi-package on CRAN. It is intended for epidemiological analysis in R. It has its own homepage, http://www.pubhealth.ku.dk/~bxc/Epi The package has been used at the course "Statistical practise in Epidemiology with R", see http://www.pubhealth.ku.dk/~bxc/SPE. A mailing list for people usi...

There is now an Epi-package on CRAN. It is intended for epidemiological analysis in R. It has its own homepage, http://www.pubhealth.ku.dk/~bxc/Epi The package has been used at the course "Statistical practise in Epidemiology with R", see http://www.pubhealth.ku.dk/~bxc/SPE. A mailing list for people usi...

Dear all, I would like to calculate the optimal cut off (threshold) of a test using the Epi package. Here I am presenting some data based on the output of two tests. I am interested in identifying the optimal cut off and its 95% CI. Running the ROC() function with the Epi package I obtain a nice picture that returns what I interpret as the optimal cut off with lr.eta=0.491. would be thi...

The stat.table function in the Epi package won't do standard deviations. It didn't seem that it would be difficult to add an "sd" function to the stat.table function. Following the example for the mean, I set up a similar function for the sd (and included it as an options) but it just won't work. (I tried sendi...

...with 3 genes, there are 3 comparisons, 4 genes=6 comparisons) in a new data.frame or matrix so that I can then test each pair with a chi-square test in a loop. Below is some sample data of the form I will be using. > lets<-c("ab","ac","bc","bd") > epi<-data.frame(cbind("M1"= c(sample(lets,10, replace=TRUE)),"M2"=c(sample(lets,10,replace=TRUE)), "M3"=c(sample(lets,10, replace=TRUE)))) > print(epi) M1 M2 M3 1 ac bc bd 2 ac ac bd 3 bd bd bd 4 ab ac bd 5 ac bc bd 6 bd bd bc 7 ab ac ab 8 bc bd ab 9 bd a...

A new major upgrade of the Epi package for Epidemiological data analysis has been put on CRAN, it is now at version 1.0.7. It contains an entirely new way of representing follow-up data on multiple timescales and multiple states. See the function Lexis(). Plus a lot of other useful stuff for epidemiological analysis. See more...

Hi all, I have a question and I do not know If I am at the right place to ask this question. But is there someone that has the formula of CKD-Epi in code in R? I have tried a lot of loops, but none of the approaches give me the right answer. Is there someone who has this formula coded? Thank you! [[alternative HTML version deleted]]

Hi all, I want to know how to perform the test Breslow-Day test for homogeneity of odds ratios (OR) stratified for svytable. This test is obtained with the following code: epi.2by2 (dat = daty, method = "case.control" conf.level = 0.95, units = 100, homogeneity = "breslow.day", verbose = TRUE) where "daty" is the object type table svytable consider it, but when I run the code does not throw the homogeneity test. Thanks....

...haven" (and possibly others) no longer access files on shared network drives. The problem appears to be in the normalizePath() function. The file can be read from a local drive or by functions that don't call normalizePath(). The error thrown is Error: path[1]="\\Hzndhhsvf2/data/OCPH/EPI/BHSDM/Group/17.xls": The system cannot find the file specified Here's my session: library(readxl) library(XLConnect) # attempting to read file from network drive df1 <- read_excel("//Hzndhhsvf2/data/OCPH/EPI/BHSDM/Group/17.xls") # pathname is fully qualified, but error thro...

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 library(Epi) ROC(form=y~x,plot="ROC") library(rms) rcorr.cens(x,y) lrm(y~x) As you can see, ROC() gave an AUC 0.782, while both rcorr...

...>> problem appears to be in the normalizePath() >> function. The file can be read from a local drive or by >> functions that don't call normalizePath(). The error >> thrown is >> >> Error: >> path[1]="\\Hzndhhsvf2/data/OCPH/EPI/BHSDM/Group/17.xls": >> The system cannot find the file specified >> >> Here's my session: >> >> library(readxl) library(XLConnect) >> >> # attempting to read file from network drive df1 <- >> read_excel(...