similar to: Cochran-Armitage-trend-test

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

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

Hi, my name is Giacomo. I would like to know how to create a Trend variable in a regression using R. Thank you for your help. My best regards, Giacomo --------------------------------- [[alternative HTML version deleted]]

Dear R-Users, is there any statement to fit a orthogonal regression in R environment? Many thanks in advance. Best regards, Vito Diventare costruttori di soluzioni Became solutions' constructors "The business of the statistician is to catalyze the scientific learning process." George E. P. Box "Statistical thinking will one day be as necessary for efficient

When using fitdistr() with the exponential, log-normal and beta distributions, you get the relevent rate, mean, standard deviation, shape1 and shape2 but you get a number bellow those that are in () and I was wandering what exactly those numbers represent and how they relate to the data. Many thanks Mark Miller

Dear all, I'm looking for some material on data mining with R. I have something from Luis Torgo but I'd like to see something else. If anybody could help me I'll be thankful Adri??n

Hello, anyone got an idea on how to use stl() so that the remainder eventually becomes white noise? i used stl repeatedly but there is autocorrelation in the remainder that i can't get rid of. os: linux suse9.3 ------------------------------------------------ Sebastian Leuzinger Institute of Botany, University of Basel Sch??nbeinstr. 6 CH-4056 Basel ph 0041 (0) 61 2673511 fax 0041 (0)

hi all i have a another stats question. i would like to solve the following question: y(i)=a+b*x(i)+e(i) i.e. estimate a and b (they should be fixed) but i dont want to specify the standard density to the straight line. this can be done using kernel regression. the fitted line is however fitted locally. does anyone have a reference that will help me with my problem. i am still new to

Hello to everybody, I'd like to submit a problem I'm dealing with, and I can't get an answer to by myself. I have to test if my data come from a specific probability distribution, of which I know the analytic form both of the p.d.f. and the c.d.f. Namely, it is the hypoexponential distribution, sum of two exponentials with different parameters. Is there any way I can manage the

Does p-value on Box.test(data,lag=l) returns probability, that H0: cor(1)=cor(2)=..=cor(l)=0 holds? Thanks. [[alternative HTML version deleted]]

To the R-help list, In the documentation of the prop.trend.test function in the stats package, no bibliography has been provided which would allow one to find out the theoretical basis of that function and/or details of its implementation. May I suggest that some bibliography be included, as it generally happens with other statistical functions. I currently use R 2.8.0 running on Windows XP.

Hello, I'm trying to calculate a chi-squared test to see if my data are different from the theoretical distribution or not: chisq.test(rbind(c(79, 52, 69, 71, 82, 87, 95, 74, 55, 78, 49, 60),c(80,80,80, 80, 80, 80, 80, 80, 80, 80, 80, 80))) Pearson's Chi-squared test data: rbind(c(79, 52, 69, 71, 82, 87, 95, 74, 55, 78, 49, 60), c(80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80,

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

hi, rgev function gives me random deviates and I have a data set which I am fitting to an EVD,IS there a way I can plot both observed and ideal evd on the same plot thankyou Rangesh

Hi there, I'm a newbie, plesae bear with me. I have a dataset with about 10000 ~ 30000 data points. Would like fit to both Gamma and Normal distribution to see which one fits better. How do I do this in R? Or I could do a normality test of the data, if it's normal, I then will do a normal fit, otherwise, a gamma fit. But again, I don't know how to do this either. Please help! David

Hello, I am trying to simulate 10 relicates of 100-tables. Each table is a 2 x 3 and 80% pf the tables are true nulls and 20% are non-nulls. The nulls follow the Hardy Weinberg distribution (ratio) 1:2:1. I have the code below but the p-values are not what I am expecting. I want to use the Cochran Armitage trend test to get the p-values. num.reps=10 num.vars=1000 pi0 = 80 num.subjects = 100

I haven't been an R lister for a bit, but I hope to enlist someone's help here. I think this is a simple question, so I hope the answer is not much trouble. Can you please respond directly to this email address in addition to the list (if responding to the list is warranted)? I have a histogram and I want to see if the data fit a Poisson distribution. How do I do this? It is

Hi, I have been trying to use the function cochran.test from the Outliers package to test for homogeneity of variance. This works well except when I use transformed data. Would anyone have an idea why it doesn't work and how I could do the cochran test on transformed data? Thanks, Stephanie >library(outliers) > set.seed(1234) > x=rnorm(100) >

Dear all, I am trying to look for a built in function that performs the cochran Q test. that is, cochranq.test(X) where X is a contingency table (maybe a matrix or data.frame). The output will naturally be the test statisitcs, p-value, etc. A quick search on Google gives me the cochran.test in the 'outlier' package, but I had a look at the description of the test and it doesn't look

Hello, I'm new in this R world and I don't know much about statistics, but now I have to analize some data and I've got some first queries yet: I have 5 sets of area mesures and each set has 5 repetitions. My first step is to check data looking for outliers. I've used the outliers package. I have to use the cochran test and the grubbs test in case I find any outlier. The problem