similar to: exact conditional mantelhaen.test estimate is 0 ?!

Hi R users I have serveral binary variables (e.g., X1, X2, X3, X4, X5, X,6, and X7) and one continuous variable (e.g., Y1). I combined these variables using data.frame() mydata <- data.frame(X1,X2,X3,X4,X5,X6,X7,Y1) after that, I sorted this data.frame rank.by.Y1<-order(mydata[,8]) sorted.mydata<-mydata[rank.by.Y1,] after that, I replaced Y1's values with values ranging from 1

Full_Name: Chien-yu Peng Version: 2.0.1 OS: Windows XP Professional Submission from: (NULL) (140.109.72.181) Dear all: Although I don't know you, I am thankful for your help. When I use the function mantelhaen.test for R x C x K (R, C > 2) table, the output is not the same as SAS's. I don't know that the result consist with one of SAS's. But it works correctly for 2

Hello, I was using the mantelhaen.test function (2x2 J tables for conditional independence testing). I noticed that the option for the continuity correction (correction=T or correction=F) sometimes made a difference while sometimes it did give the same results regardless of correction=T or correction=F. Does the continuity correction apply only in certain cases or in certain strata with some

Hi, I am new to these functions. I'm wondering if there is anyway to save the entire results (all attributes of the result object) from the chisq.test or mantelhaen.test functions? For example, from chisq.test function, you will have statistic, parameter, p.value, expected, etc. in the result list. How can I save all of them in one shot to, says, a text file or csv file? Thank you. - adschai

Hi, My apologies for the naive question! I have three overlapping sets and I want to find the probability of finding a larger/greater intersection for 'A intersect B intersect C' (in the example below, I want to find the probability of finding more than 135 elements that are common in sets A, B & C). For a two set problem, I guess I would do a Fisher or chi-square test. Here is what

Dear All, I'd like to perform the generalized Cochran-Mantel-Haenszel-Test (as described in Agresti (1990), Categorical Data Analysis) for my nominal data. My problem is that I have more than two groups. In fact I think I'd need an 5-dimensional array for the response variable, the control variable and three group variables. Could you please tell me if this is possible in R - and if

This is an OpenPGP/MIME signed message (RFC 2440 and 3156) --------------enig5A700CC15DEABD01445A1B70 Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Hi, I have modified the code originally posted to include capability for arbitrary weightings for the rows and columns (which apply in the nominal-ordinal and the ordinal-ordinal cases). The code is pasted

Hello, I have some satellite tag time-at-depth (TAD) frequency data that I would like some help with. The data was transmitted via satellite as percent time spent in each of 7 depth bins (0m, 0-1m, 1-10m, 10-50m etc.), binned over 6-hour intervals. I categorized each row of data corresponding to a date and time into summer vs. winter, and day vs. night, and then summed and averaged the given

Hi all: For stratified count data,how to perform regression analysis? My data: age case oc count 1 1 1 21 1 1 2 26 1 2 1 17 1 2 2 59 2 1 1 18 2 1 2 88 2 2 1 7 2 2 2 95 age: 1:<40y 2:>40y case: 1:patient 2:health oc: 1:use drug 2:not use drug My purpose: Anaysis whether case and

In SAS, for a two-way (or 3-way, stratified) table, the CMH option in SAS PROC FREQ gives 3 tests that take ordinality of the factors into account, for both variables, just the column variable or neither. Is there an equivalent in R? The mantelhaen.test in stats gives something quite different (a test of conditional independence for *nominal* factors in a 3-way table). e.g. I'd like to

Hi all, I was quite succesfully working with the RODB package on a 32 bit linux box to connect to a MSSQL Server via the freeTSL driver. After changing to a 64 bit environment I ran into some segmentation faults using function sqlUpdate on large database operations (the actual seg faults occured in the call to the C function ODBCUpdate). I just briefly looked into the code and found some extensive

How i can perform a Pearson's Chi-squared Test in this data set: | Outcome -----------------+-----------+----------------------------------+ Treatment | Sex | None |Some | Marked | Total -----------------+------------+--------+--------+-------------+ Active | Female | 6 | 5 | 16 | 27

tl;dr is the R bug tracker down or am I being an idiot? Help please ... ---- I decided I would follow up on https://stat.ethz.ch/pipermail/r-devel/2018-January/075410.html (reporting/suggesting a patch for an issue in stats::mantelhaen.test() with large data sets) Reading the instructions at https://www.r-project.org/bugs.html suggests that if I'm sure I have found something worthy of

Hi! I would like to conduct a Cochran`s Q Test in R, but have not found any suitable function. My first idea was: J <- as.table(matrix(c(6,16,2,4),ncol=2, dimnames = list("C" = c("Favorable","Unfavorable"),"Drug A Favorable"=c("B Favorable","B Unfavorable")))) L <- as.table(matrix(c(2,4,6,6),ncol=2, dimnames =

Hi everybody, I am trying to write a simple progress display based on a tcltk toplevel. My first approach was to use the progressBar widget from the BWidget library but since this is not available on every system (missing on at least almost all windows systems, I guess...) I wanted to have a backup there. So my second strategy was to use a simple toplevel with a label and update the tclvariable

{again a message that was sent to owner-r-help (which is me, currently) why on earth ???!??!? reply to R-help or the original sender Brett Presnell; } I'm teaching a categorical data analysis course this term, and a minor "problem" has resurfaced that I have often thought about before. This applies equally to Splus I suppose, but my undergrads aren't using Splus. It

Dear listers, After activating the name space for my bioconductor package (prada) I successfully ran R CMD check. However when loading the package in R and running the examples the imported function brewer.pal from package RColorBrewer is not found. I can directly call brewer.pal from the RColorBrewer name space typing RColorBrewer::brewer.pal, but it is not imported into my prada name space. When

Dear listers, After activating the name space for my bioconductor package (prada) I successfully ran R CMD check. However when loading the package in R and running the examples the imported function brewer.pal from package RColorBrewer is not found. I can directly call brewer.pal from the RColorBrewer name space typing RColorBrewer::brewer.pal, but it is not imported into my prada name space. When

Dear all, I have been trying to find a function to calculate the Breslow-Day test for homogeneity of the odds ratio in R. I know the test can be preformed in SAS but i was wondering if anyone could help me to perform this in r. In addition i have the fullrefman file to search for functions in the basic R packages, does anyone have any suggestions of an efficient way of searching for

Just a general question concerning the woolf test (package vcd), when we have stratified data (2x2 tables) and when the p.value of the woolf-test is below 0.05 then we assume that there is a heterogeneity and a common odds ratio cannot be computed? Does this mean that we have to try to add more stratification variables (stratify more) to make the woolf-test p.value insignificant? Also in the