similar to: Asymmetric CIs

mantelhaen.test() gives the exact conditional p-value (for independence) and confidence intervals (CIs)for the common odds ratio for a stratified 2x2 table. The epitools package by Tomas Aragon (available via CRAN) contains functions which use fisher.test() to calculate mid-p exact p-values and CIs for the CMLE odds ratio for a single 2x2 table. The mid-p p-value for independence for a stratified

At work after updating to R 2.5 I get an error using epitab from package epitools, when at home (R 2.4) I get no error. Could someone help me? Thanks Pietro Bulian Servizio di Onco-Ematologia Clinico-Sperimentale I.R.C.C.S. Centro di Riferimento Oncologico Via Franco Gallini 2 33081 Aviano (PN) - Italy phone: +39 0434 659 412 fax: +39 0434 659 409 e-mail: pbulian at cro.it (at work)

R-experts: A quick question, please. >From a lab exp, I got 12 positives out of 50. To get 90% CI for this , I think binom.test might be the one to be used. Is there a better way or function to calculate this? > binom.test(x=12, n=50, p=12/50, conf.level = 0.90) Exact binomial test data: 12 and 50 number of successes = 12, number of trials = 50, p-value = 1 alternative

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

The concrete problem is that I am refereeing a paper where a confidence interval is presented for the risk ratio and I do not find it credible. I show below my attempts to do this in R. The example is slightly changed from the authors'. I can obtain a confidence interval for the odds ratio from fisher.test of course === fisher.test example === > outcome <- matrix(c(500, 0, 500, 8),

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

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

Hello, I have data which I have arcsin transformed to analyse. I want to plot my data with error bars however as my data is back-transformed my standard errors are uneven. Is there a simple way to draw these asymmetric error bars in R? Thanks for your help. [[alternative HTML version deleted]]

> Dear expeRts, > > Currently, my colors are as follows: > mycol <- > c("blue1","blue2","blue3","blue4","black","yellow4","yellow3","yellow2","y > ellow1") > heatmap(snp, Rowv=NA, Colv=NA, col=mycol) > > However, I would like to have the following colors: > bright blue ->

# is there a bug in the calculation of the odds ratio in fisher.test? # Nicholas Horton, nhorton at smith.edu Fri Jan 22 08:29:07 EST 2010 x1 = c(rep(0, 244), rep(1, 209)) x2 = c(rep(0, 177), rep(1, 67), rep(0, 169), rep(1, 40)) or1 = sum(x1==1&x2==1)*sum(x1==0&x2==0)/ (sum(x1==1&x2==0)*sum(x1==0&x2==1)) library(epitools) or2 = oddsratio.wald(x1, x2)$measure[2,1] or3 =

On Fri, Dec 19, 2014 at 2:40 PM, Brian Mathis <brian.mathis+centos at betteradmin.com> wrote: > > > GPG is really what you want to be using for this. OpenSSL is a general > toolkit that provide a lot of good functions, but you need to cobble some > things together yourself. GPG is meant to handle all of the other parts of > dealing with files. > > I will expand on

We are running tinc (v. 1.0.26) in switch mode with bmx6 (another mesh protocol) running on top of the tap0 interface on about 25 devices. The asymmetric routing of UDP packets is causing my firewall and I presume others to drop some of the packets, since there are no outbound SYN packets originating from the device running tinc. Is there any way to mitigate this issue besides enabling tcponly

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 =

Hello all, I've been using binconf (package Hmisc) at a range of alpha values and noticed that using the 'Wilson' method when alpha is larger (i.e. narrow CI), results in the upper value being smaller than the lower value. The 'exact' and 'asymptotic' methods give results in the realm I'd expect. But the help file suggests: "Following Agresti and Coull, the

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

Hello CentOS list, I have a requirement that I need to use encryption technology to encrypt very large tar file on a daily basis. The tar file is over 250G size and those are data backup. Every night the server generated a 250G data backup and it?s tar?ed into one tarball file. I want to encrypt this big tarball file. So far I have tried two technologies with no success. 1) generating RSA 2048

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