similar to: trend in incidence rate

Dear sir, I am Shibu John from Thrombosis Research Institute India. It is a multidisciplinary organisation concerned with the interrelated problems of thrombosis and atherosclerosis. I was searching for Cochran armitage trend test program in R. Then I had seen your R coding for C-A trend test. I tried that in the R software. But I can?t run the program due the [Error: could not find function

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)

Hi, further down from the problem I asked about yesterday I encounter the following problem (please appologize for the lengthy mail - I have not found a brief example making the issue clear): - using R packages making use of namespaces & - S4 classes when using --with-package-versions in the install call: Pkgs install ok and the first pkg 'base' can be loaded in R

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

Hallo! Does anybody know how to calculate a sample size estimation for proportions with continuity correction? ? I only found EpiR which seems to calculate without continuity correction: library(epiR) epi.studysize(treat = .65, control = .50, n = NA, sigma = NA, power = 0.80, r = 1, conf.level = 0.95, sided.test = 2, method = "cohort") $n [1] 340 Thanks! Karl

Hello R developers, installing R 2.0.0 from source on a Linux Debian system (gcc 3.2.2) - which worked just fine - I encountered the following problem when trying to update/install packages with the option --with-package-versions: example: RUnit 0.4.0 (also on CRAN) R20 CMD INSTALL -l /mnt/local/R/R-2.0.x-libs-EpiR --with-package-versions RUnit_0.4.0.tar.gz or after untaring R20 CMD INSTALL

A new package 'epiR' is available on CRAN. Package description as follows: Package: epiR Version: 0.9-3 Date: 2008-03-24 Title: Functions for analysing epidemiological data Author: Mark Stevenson <M.Stevenson at massey.ac.nz> with contributions from Telmo Nunes, Javier Sanchez, and Ron Thornton. Maintainer: Mark Stevenson <M.Stevenson at massey.ac.nz> Description: A package

A new package 'epiR' is available on CRAN. Package description as follows: Package: epiR Version: 0.9-3 Date: 2008-03-24 Title: Functions for analysing epidemiological data Author: Mark Stevenson <M.Stevenson at massey.ac.nz> with contributions from Telmo Nunes, Javier Sanchez, and Ron Thornton. Maintainer: Mark Stevenson <M.Stevenson at massey.ac.nz> Description: A package

Anyone know if coin can run a permutation test based on a (user-defined) statistic other than the mean difference? The function independence_test does the permutation t-test via difference in means. I'm wondering if it's possible to use independence_test to run a permutation test for some other statistic than the difference in means. For example, I'd like to run a permutation test

Hi, I wonder whether a package with namespace & dynamic library can be installed with the version number attached, ie. with the argument --with-package-versions. Is this currently possible? Using R 1.9.1 on Debian 3.0 I encounter a problem when trying to load a package installed with R91 CMD INSTALL --with-package-versions -l /mnt/local/R/R-1.9.x-libs-EpiR RIO in R started with R91

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

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 ideas very very welcome thank

Hello, Would you tell my how to extract a result from a test - it's justified because I need to run this test many times. Here is an example from authors' test: > library("coin") > lungtumor <- data.frame(dose = rep(c(0, 1, 2), c(40, 50, 48)), tumor = c(rep(c(0, 1), c(38, 2)), rep(c(0, 1), c(43, 7)), rep(c(0, 1), c(33, 15)))) > ca.test<-independence_test(tumor ~

Hi, I was wondering if there was some sort of package or function that calculated asymmetric confidence intervals for small proportions. I thought of both the epicalc and epitools package, but I am hoping to find something where you can just plug in a standard error and point estimate and it will output the upper and lower CI bounds. Thanks! Sarah

Hi, since a couple of days ago I see failures when using R CMD INSTALL -l lib path-to-src/pkg.tar.gz This occurs in automated library updating via a shell script and is caused by the immediate preceding call to leave an 00LOCK folder even though the package installation of that call succeeded. This occurs at random (no obvious/systematic correlation to package name or or such) and only for

I've recently set up a new server running RH 6.1 and Samba 2.0.5a. All our older servers are running 1.9.18p10, and were left that way on the "If it ain't broke" principle. On the new server, I've noticed that DOS wildcards don't work correctly for some (but not all) old 16-bit DOS programs running under Win95/98. The symptom is that the program only

HI R-users, I am trying to repeat an example from Rayner and Best "A contingency table approach to nonparametric testing (Chapter 7, Ice cream example). In their book they calculate Durbin's statistic, D1, a dispersion statistics, D2, and a residual. P-values for each statistic is calculated from a chi-square distribution and also Monte Carlo p-values. I have found similar p-values

Hi everybody, I'm trying to create a numerical data frame on which to perform PRCC. So far I have created a data frame that consists of function/vector output that displays in numerical form, but when I try and run PRCC (from epiR package) I get the following error message: "Error in solve.default(C) : Lapack routine dgesv: system is exactly singular" It appears this is because

Hi, I am new to bioconductor, trying to install KEGGSOAP package, but got warnings() when installing and error message when trying to load the package, can anyone suggest what went wrong? many thanks John > source("http://bioconductor.org/biocLite.R") Bioconductor version 2.11 (BiocInstaller 1.8.3), ?biocLite for help > biocLite("KEGGSOAP") BioC_mirror:

Hi, I recently updated to R 3.0.1. I'm running linux ubuntu 12.04. I realized that I have to update all the installed packages so I run > update.packages(checkBuilt=TRUE) as described here http://www.r-bloggers.com/r-3-0-0-is-released-whats-new-and-how-to-upgrade/ . The first thing it did was telling me that it will not update several packages. When it was finished I used the warnings