similar to: Help me in Cochran armitage trend test Coding

Dear R-enthusiasts, I am trying to do a Cochran-Armitage test for trend in R. After consulting google I found Torsten Hothorn's remark that the 'coin' library could be used. 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)),

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 am searching for the Cochran-Armitage-trend-test. Is it included in an R-package? Thank you! --

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

Dear R helpers, I am trying to understand how to use the independence_test function in the coin package. I think I suffer from a misunderstanding about what the package does. Either that or I do not understand how to use it properly. Specifically, I cannot understand if I can test independence of arbitrary statistics. Take the following example: set.seed(10) d <- data.frame(y = c(rnorm(10,

I partly answered my question since independence_test() function in coin package apparently do Cochran-Armitage trend test just like Eric Lecoutre's function tabletrend() - slightly modified here: > independence_test(pheno ~ geno, data = dat2, teststat = "quad", scores = list(geno = c(0, 1, 2))) Asymptotic General Independence Test data: pheno by groups 1 < 2

Dear List, I'm interested in performing a permutation test on the interaction between a binary treatment indicator and a covariate (either continuous or categorical). I'm interested in the p-value of the interaction effect from a permutation test, and I'm using the coin package for that purpose. As I haven't seen any examples like this in the package documentation (or anywhere

Are there internal methods available for R2.6 (I'm using the mac os x gui) for extracting (or subsetting) S4 objects? Using the independence_test() function implemented in the COIN package, I can't seem to select out p-values upon completion of each iteration of a loop. Sorry if my search was incomplete beforehand. -- View this message in context:

Hi all, Since upgrading to perl 5.8.5-36.el4_5.2 I have had no end of problems, starting with needing compress::zlib which I got from rpmforge, then I had to reinstall scalar::util to get mailscanner working but now I can't install or upgrade any perl modules. Has anyone else had problems or can someone tell me what I've done wrong please? regards tom

Hi Centos Users I would like to allow outgoing traffic to Google Analytics servers (destination port 80). I wish to do a iptables rule. How to whitelist all Google Analytics servers? cheers Simon -- XMPP: sjolle at swissjabber.org

Is it possible to implement the Jonckheere-Terpstra test for ordered alternatives using the coin package: Conditional Inference Procedures in a Permutation Test Framework? I found jonckheere.test{clinfun}, but it uses a normal approximation when ties are present in the data. To make this concrete, I've include a small dataset. Thanks. --Dale Hollander and Wolfe, 1999 Table 6.6, pg 205

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

Dear friend , when i try to change my spooler directory to on of network drive ,iam getting this error "Server spool directory location is invalid".please help me for this problem. Thanks & Regards Shibu Varghese U Bahrain Specialist Hospital Information Technology P O Box 10588 Manama, Bahrain Tel: +973 812028 Fax: +973 812012 _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _

Hello everyone, I run a 64bit system here and install flash plugins for firefox using nspluginwrapper. I get the movies but there is no sound. Flashplayer is: flash-plugin-9.0.48.0-release.i386.rpm nspluginwrappers are: nspluginwrapper-0.9.91.5-1.x86_64.rpm nspluginwrapper-i386-0.9.91.5-1.x86_64.rpm Any help in tracking this down would be great. Thnx. -- Regards Robert Smile...

Hello, I've got a CentOS 5 box. I want to install the uucp rpm to give access to the cu command. I do a yum -y install uucp and yum is hanging on the transaction test. I get the file downloaded fine, but on transaction test yum freezes and i have to kill it's pid to stop it. Any ideas? I'm looking for a program for serial terminals. Thanks. Dave.

I'm interested in doing centralized logging & analysis of logs from my CentOS boxes. I messed around with syslog-ng and it seems like it's a better syslog. But I noticed most of the "usual suspects" of third-party repo maintainers (Dag, Axel, etc.) don't include it. Karanbir has an el4 version, but all my boxes are CentOS 5. I was going to go ahead and roll my own

I am using "apply" to run exact permutation t-tests by columns using the coin package. For example: library(coin) dat <- matrix(rnorm(7*35),7,35) fun <- function(x) { pvalue( independence_test(x~f, data=data.frame(x, f=factor(c(rep("a",4),rep("b",3)))), distribution = "exact")) } p.vals <- apply(dat, 2, fun) Some small-scale

Dear R-list members, This is in fact a question about statistics, not directly about the R software. The coin package, for conditional inference through permutation methods, has as it main function the function independence_test. One of its parameters is teststat, about which the package documentation says: teststat: "a character, the type of test statistic to be applied: either a

Dear R Users, I am hoping you can help me. I have received code from a colleague who uses Matlab. I need to translate it into R. I am wondering if there is a randomization t-test (from non-parametric statistics) function already defined in R. (In Matlab the function is randtest.m.) ************************************************************************ ********** QUESTION: Is anyone

I would like to compare the incidence rates of three groups. They are supposed to have different risks so I would like to test whether there is a increasing trend in the incidence rates. Does R or any packages provide a trend test for incidence rates? I checked epiR and epitools. It seems they do not have this function. Thank you for the help. -- View this message in context: