similar to: calculating p-values from t-values for a Bonferroni adjustment

The Partek package (www.partek.com) allows only two selections for Multiple Test Correction: Bonferroni and Dunn-Sidak. Can anyone suggest why Partek implemented Dunn-Sidak and not the other methods that R has? Is there any particular advantage to the Dunn-Sidak method? R knows about these methods (in R 2.1.1): > p.adjust.methods [1] "holm" "hochberg" "hommel"

Hi! I need to run a wilcoxon (Mann-whitly, in fact) test with bonferroni correction, as I am running 10 consecutive wilcoxon test not independent, and I know that bonferroni will partially correct for this problem, but I have no idea how to do it with R, I have been looking in the archive but couldn't understand how to do it. The format I am using at the moment is r4_o <- [1]

Hi, I'm having no luck figuring out how to find Bonferroni simultaneous confidence intervals to obtain a family of estimates in R. Does anyone know how to do this? Thank you!

Hola a todos/as estoy analizando un diseño factorial con una variable intrasujeto y una inter-sujeto. Ambas resultan significativas en el anova y quiero hacer un análisis por pares, pero aún no sé como puedo hacerlo en R. Me interesan ambos Bonferroni y Tukey. ¿Alguien podría indicarme? Gracias [[alternative HTML version deleted]]

I append below a suggested update for p.adjust(). 1. A new method "yh" for control of FDR is included which is valid for any dependency structure. Reference is Benjamini, Y., and Yekutieli, D. (2001). The control of the false discovery rate in multiple testing under dependency. Annals of Statistics 29, 1165-1188. 2. I've re-named the "fdr" method to "bh" but

Hello- I have run an ANOVA on 4 treatments with unequal sample sizes (n=9,7,10 and 10). I want to determine where my sig. differences are between treatments using a Bonferroni test, and have run the code: pairwise.t.test(Wk16, Treatment, p.adf="bonf") I receive an error message stating that my arguments are of unequal length: Error in tapply(x, g, mean, na.rm = TRUE) :

Please can someone advise me how I can adjust correlations using bonferroni's correction? I am doing manny correlation tests as part of an investigation of the validity/reliability of a psychometric measure. Help would be so appreciated! Cheers, Louise

Hello. We have a question concerning the nonparametric analysis of a dataset, which resulted in rejection of the null hypothesis (Kruskal-Wallis-test = H-test). In order to find out which sample means actually are statistically different, we want to do multiple comparisons with the Wilcoxon rank sum test (= U-test); the p-level should be corrected according to Bonferroni-Holm. Thus we decided to

Hi! as my subject says I am struggling with the different of a two-way ANOVA and a (two-way) ANCOVA. I found the following examples from this webpage: http://www.statmethods.net/stats/anova.html # One Way Anova (Completely Randomized Design) fit <- aov(y ~ A, data=mydataframe) # Randomized Block Design (B is the blocking factor) fit <- aov(y ~ A + B, data=mydataframe) # Two Way

I've run an ancova, edadysexo is a factor with 3 levels,and log(lcc) is the covariate (continous variable) I get this results > ancova<-aov(log(peso)~edadysexo*log(lcc)) > summary(ancova) Df Sum Sq Mean Sq F value Pr(>F) edadysexo 2 31.859 15.9294 803.9843 <2e-16 *** log(lcc) 1 11.389 11.3887 574.8081 <2e-16 ***

I'm running the following ANCOVA model in R: X=a*b*d + d where d is the covariate. I want to produce adjusted means for X for all possible treatment combinations. a has 6 levels, b has 2 levels, c has 2 levels, so this is a 6x2x2 = 24 possible treatment combinations. When I run: effect("d", Ancova.Model) I am given only 20 adjusted means (also, I'm not exactly sure of how the

Hello R gurus, I am a beginner with R. I am doing an ANCOVA analysis using 'aov,' and need some help understanding how 'aov' works. I have a dataset (taken from http://faculty.vassar.edu/lowry/ch17pt2.html) looking at hypnotic induction. The variable 'X' is a measure of how susceptible the subject is to being hypnotized, the variable 'Y' is how well the

I have some data with two categorises plus/minus (p53) and a particular time (Time) and the outcome is a continuous vairable (Result). I set up a maximum model. ancova <- lm(Result~Time*p53) > summary(ancova) .. Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.05919 0.55646 0.106 0.916 Time -0.02134 0.01785 -1.195 0.241 p53plus

Dear all, I am trying to run Wilcoxon Rank-Sum Test with Bonferroni's correction. I have two lists: l0, l1: mapply(function(x,y)wilcox.test(x,y)$p.value, l0, l1) How do I run Bonferroni's correction on mapply? Any help is much apperciated. Thanks, -Raj

Dear R-help, I have a question about pairwise.t.test and adjustment for multiple comparisons for paired data points. I have the following data: n=c("x", "x", "x", "x", "x", "x", "x", "x", "x", "x", "y", "y", "y", "y", "y", "y",

Hi R-Users, I am working with sm.ancova (in the package sm) and I have two problems with the graph, which is automatically generated when sm.ancova() is run. 1-Besides of the fitted lines, the observed data appeared automatically in the graph. I prefer that only fitted lines appear. I check the sm.options, but I could not find the way that the observed data do not appear in the graph. 2-I

Hi, I was trying to figure out how to do post-hoc tests for Two Way ANOVAs and found the following 2 approaches: a. Do pairwise t-tests (bonferroni corrected) if one finds significance with the ANOVA. Link- http://rtutorialseries.blogspot.com/2011/01/r-tutorial-series-two-way-anova-with.html b. Do TukeyHSD on an aov model Link-

Hi, i am trying to run an ANCOVA and a bootstrapped ANCOVA analysis on a specific data set. I am using the ancova and ancboot functions as in the following code: setwd("C:/Users/User/Desktop/Rdatabilingualstudy2012") bilingualismdata<-read.spss("bilingualdataforconferences2012.sav", use.value.labels = TRUE, to.data.frame = TRUE)

How to Ancova in R? I know this has got to be an FAQ, because I see it asked in the lists, but I haven't seen an answer to it. I see the R-sm has the ancova thing happening, but I kind of doubt that what I'm trying to do is "smoothing"... -- Pete Hurd phurd at uts.cc.utexas.edu http://www.zo.utexas.edu/research/phurd Section of Integrative Biology, University of Texas, Austin

Dear list members: I searched the R-help for packages to calculate power for an ANCOVA problem I have. I have found power.t.test, power.anova.test. But it seems that I can not found one with ANCOVA. I have two datasetsets with variables: univariate response(one data with continous response and one with 0 and 1), treatment(two levels), covariates(x1,x2,x3). I would appreciate your help. Tony