similar to: creating 3-way tables for mantelhaen.test

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 listers, I am trying to compute the exact conditional test given strata margins of a 2 by 2 by K array using the mantelhaen.test function to get a common odds ratio estimate. The estimate for the test on the following data is 0, which in my opinion dosen't make any sense. x <- array(c(53, 6098, 1006, 4521, 63, 8070, 1163, 6137), dim=c(2,2,2)) x , , 1 [,1] [,2] [1,] 53 1006

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 R users Having both a faster CPU and more memory will boost computing power. I was wondering if only adding more memory (1GB -> 2GB) will significantly reduce R computation time? Taka, _________________________________________________________________ Get FREE company branded e-mail accounts and business Web site from Microsoft Office Live

Dear R-users I have four simple linear models, which are all in the form of a*X+b The estimated parameter vectors are a <- c(1,2,3,4) b <- c(4,3,2,1) My goal is to draw a plot where x-axis is X (range from -100 to 100) and y-axis is the sum of all four linear models X <- seq(-100,100, length=10000) plot(X, sum of the four linear functions) I started with a function for summing

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

Dear R-users I need to read some text files produced by some other software. I used readLines (with n = -1 ) command and then tried to find some numbers I liked to extract or some line numbers I like to know. The problem is that there is no empty line at the end of the text files. R gave me this warning below; In addition: Warning message: incomplete final line found by readLines on

Dear R-users >prob <- c(0.5,0.4,0.3,0.1,0.0) >cal <- prob * log(prob,base=2) >cal [1] -0.5000000 -0.5287712 -0.5210897 -0.3321928 NaN Is there any way to change NaN to zero ? I did come up with this by applying Ripley's relpy to my previous question cal <-prob*log(pmax(prob,0.00000001),base=2) Any suggestion ? Thank you Taka

Dear R-users, I am looking for a way to scan a directory and then make a string vector consisting of the file names in the directory. For example, under c:\temp\ there are 4 files a.txt, b.txt, c.txt, and d.txt I would like to have a string vector c("a.txt","b.txt","c.txt","d.txt") How do I do that? Thanks Taka,

Hi R users Which version of S does the current R version use? Venables and Ripley's book (2000) says R is version 3 of S language. Is there any change since 2000 ? I searched R-help archive for using C or C++ with R. Although there are some postings, I am looking for up-to-date answers. Which one ( C and C++) works better with R? What complier is the best suited for R ( gcc or visual

Hi R users Here is what I got with help from Petr Pikal (Thanks Petr Pikal). I modified Petr Pikal's code to a little to meet my purpose. I created a function to generate a matrix generate.matrix<-function(n.variable) { mat<-matrix(0,n.variable,(n.variable/2)/5+1) #matrix of zeroes dd<-dim(mat) # actual dimensions mat[1:(dd[1]/2),1]<-1 #put 1 in first half of first column

Hi R-users I have a data set. There are 10 products and the numbers of people who ranked the products. The format of the data set is productID rank1 rank2 rank3 rank4 rank5 rank6 rank7 rank8 rank9 rank10 ------------------------------------------------------------------------------------------------------- 1 10 2 3 3 6 4 2 5

Hi, I am analyzing a data set with greater than 1000 independent cases (collected in an unrestricted manner), where each case has 3 variables associated with it: one, a factor variable with 0/1 levels (called XX), another factor variable with 8 levels (X) and a third response variable with two levels (Y: 0/1). I am trying to see if X1 has an effect on the relationship between X2 and the

Hi R user I need to read in some values from a computer program output. I can't change the output format because the developer of the program doesn't allow to change the format of output. There are two formats. First one looks like this if I have 10 variables, ------------------------------------------------------------------------------------------------------ [ 1]

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

Dear R-users I have 5 dependent variables (y1 to y5) and one independent variable (x) and 3 conditioning variables (m, n, and 0). Each of the conditioning variables has 2 levels. I created 2*4 panel plots. xyplot(y1+y2+y3+y4+y5 ~ x | m*n*o,layout = c(4,2)) I would like to reorder the 8 panels. I tried to use index.cond (e.g., index.cond = list(c(1,3,2,4,5,7,6,8)) but it didn't work out.

Hi R users I like to extract lower diagonal elements of a matrix in such a way like, data[1,2], data[1,3], ...., data[5,6] are extracted from a matrix called 'data' This short script below is what I have written so far. ########################################## data <- matrix(rnorm(36,0,1),nrow=6) temp<-c() for (i in 1:(nrow(data)-1)) { for (j in (i+1):nrow(data)) {

Hi R users I like to generate some strings (a character vector) in a special way like If i have 5 variables "002.001", "003.001", "003.002", "004.001", "004.002", "004.003", "005.001", "005.002", "005.003", "005.004" so the created string vector's elements are "002.001",

Hi R users This looks a simple question Is there any difference between between rnorm(1000,0,1) and running rnorm(500,0,1) twice in terms of outcome ? TM