similar to: how to plots two pairwise data sets into a same graph

Dear All, I am desperately trying to find any R package that fits a mixture survival models also know as a cure models. This are survival model where the survival function is improper which also means that a fraction of subjects are expected not to expreience the event. A Huge literature has been developed for thes type of models but I couldn't find any R package that fits this type of

Hi! I''m experiencing such error while running simple model rspec just to make sure everything is working. $ script/spec spec/models/site_user_spec.rb /var/lib/gems/1.8/gems/activesupport-1.3.1/lib/active_support/dependencies.rb:100:in `const_missing'': uninitialized constant Spec::Rails (NameError) from

I am familiar with pairwise t-tests, corrections for multiple testing, etc. however I have a problem whose answer I have not found after extensive R-help archive and Google searching. What I have done in the past: I have N items which are measured, exposed to a condition, and then measured again. I wish to know if the condition changes the items so I can perform a t-test. Better yet I can

Dear all, I am looking for an equivalent of the pairwise.t.test() function using a complex survey design. Any suggestions on this or how to account for weights using the pairwise.t.test itself? Thanks, V [[alternative HTML version deleted]]

An embedded and charset-unspecified text was scrubbed... Name: ?????? ?? ????????. URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130412/7bbdecd1/attachment.pl>

Dear list I am comparing the results of several different experimental setups. With kruskal.test() I can test if there is any difference at all in any of them, if I understand it correctly. But now, when there is a difference, how do I generate a (half-) table of pairwise comparisons, using e.g. wilcox.test(), to find the ones where the difference actually occurs. I guess I don't have to

Does anybody know if there is some sort of "pairwise" option for handling missing cases in lm and computing the relevant statistics? I would be much obliged if anyone could help... Regards Alan Simpson Roberts Research Group

Hi, I have a number of genes (columns) for which I want to examine pairwise associations of genotypes (each row is an individual)...For example (see data below), I would like to compare M1 to M2, M2 to M3, and M1 to M3 (i.e. does ac from M1 tend to be found with bc from M2 more often than expected.) Down stream I will be performing chi square tests for each pair. But I am looking for a way to

I want to create a MDS plot with pairwise Fst values derived from a population genetics project. My Fst values are in a tab-delimited file (lower triangle only) that I view with Excel. When I use the cmdscale command I get the message: Error in cmdscale(x) : Distances must be result of dist or a square matrix In addition: Warning messages: 1: "^" not meaningful for factors in:

Hello R listing, I did two-way anova on lm. Further question the investigator interested in is what two treatments are different? I am looking for a command which could do pairwise comparison for every treatment. Could anyone help me out? Thanks a bunch Kevin

How can I compute pairwise p-values in Kaplan-meier plots for three or more groups? bin.1<-cut(score,c(-1000,-1,1,1000),c("low","intermediate","high")) I use "km.coxph.plot" currently which reports one p-value. Thanks very much. [[alternative HTML version deleted]]

Dear R gurus, to compute the correlation matrix of "n" variables with "n_obs" observations each, possibly including NA, I use cor(M, use="pairwise.obs") where m is a "n" x "nobs" matrix. Now I want to know the number of observations actually used in this computation, namely for each pair of columns in M, say pair (i,j), I want to compute sum(

There are two line 216 and 218 Three development stages 5 WEEK (5W), 7W, 9W. Three tissue: Ca, Co, Pa each with 2 biological replicate. With two biological replicate. I want to do differential gene expression analysis using DESeq2 so I tried these codes after reading about DESeq2: ,my aim is to do the pairwise comparison. how to make colData and design formula. library("DESeq2")

Hi, I'm wondering if i can use the function "TukeyHSD" to perform the all pairwise comparisons of a "aov()" model with one factor (e.g., GROUP) and one continuous covariate (e.g., AGE). I did for example: library(multcomp) data('litter', package = 'multcomp') litter.aov <- aov(weight ~ gesttime + dose, data = litter) TukeyHSD(litter.aov, which =

Dear R Core, pairwise.wilcox.test does not handle "paired = TRUE" correctly; e.g. set.seed(13) x <- rnorm(20) g <- c(rep(1, 10), rep(2, 10)) wilcox.test(x ~ g)$p.value # 0.075 pairwise.wilcox.test(x, g)$p.value # 0.075, o.k wilcox.test(x ~ g, paired = TRUE)$p.value # 0.105 pairwise.wilcox.test(x, g, paired = TRUE)$p.value # 0.075, wrong The line wilcox.test(xi, xj,

Hi, I have just encountered a strange behaviour from 'cor' with regards to the treatment of NAs when calculating Spearman correlations. I guess it is a subtle bug. If I understand the help page correctly, the two modes 'complete.obs' and 'pairwise.complete.obs' specify how to deal with correlation coefficients when calculating a correlation _matrix_. When calculating

Hello, this is probably more a statistical question than an R-specific problem, but I'll risk it. I've fitted a Cox Proportional hazard model with one factor Treatment (seven levels) as a predictor variable. The general Null hypothesis (all groups show the same survival behaviour) is clearly rejected. Now, is there any (statistically sensible) way of doing pairwise comparisons and/or

Hello, I have a linear model with different slopes for different treatment groups. I need to pairwise compare the different slope estimates for the different treatment groups. Is there a package that does pairwise comparisons of slope coefficients, making the appropriate adjustments in the P values? Thanks, John. -- ========================================== John Janmaat Department of

R-List, I'm trying to compute pairwise distances among pairs of observations, which each pair containing data from 2 groups. There are more than 100000 unique pairs. I have programmed a distance function that has three parameters, a vector of covariates from the ith observation in Group 1, a vector of covarites from the jth observation in Group 2, and a weighting matrix. I have used

hi - sorry for a newbie question, but I've tried to go through the documentation and couldn't find anything that would address my specific need. in a nutshell, I have a txt file containing a data matrix with 10 columns of data. I would like to generate pairwise plots of the data, i.e. 1 column against the other 9. Since this will produce 50 plots, I'd like to do it using a