similar to: Bonferroni correction for multiple correlation tests

Hi there, I have subscribed to R-help but am not sure how to view or post questions? I think this is the right way. I am planning on doing a multivariate regression investigating the relationship between depression (a continuous variable) and social support variables (mostly continuous, some categorical) among older people. I have a number of demographic and health-related variables that I am

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

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 am performing a sequential bonferroni adjustment on the results of an ANCOVA but the equation I have for calculating p-values from the t-values is not working. I can't seem to find it anywhere else. This is the code I have now: 2*(1-pt(t,df)) where t=t-value and df=degrees of freedom from the ANCOVA. Is there an error in the code? maybe another way to find the p-values? Help! -- View this

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

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"

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

I ungraded my Linux Operating system from Ubuntu 16.04 LTS to 8.04 LTS Samba will not start after an upgrade to 18.04. My samba is running as an NT style member server server with smbd nmbd and winbindd all running. None of them will start. It worked this morning. This morning it was running Version 4.3.11-Ubuntu This afternoon after the upgrade it's Version 4.7.6-Ubuntu The logs in

Dear readers, S-Plus have a rts function (Regular Time Series), which is used like so: fveks<-read.csv('http://louise.hoffman.googlepages.com/veks.csv',header=TRUE,sep=',') attach(fveks) acf(ts.intersect(rts(HC.f),rts(Ta.f),rts(GR.f),rts(W.f))) Warning the csv file is 750kB. Can the same be done in R? Lots of love =) Louise

Please can someone explain why I seem to get these contradictory results? choose(5,2) [1] 10 gamma(6)/(gamma(3)*gamma(4)) [1] 10 gamma(6)/(gamma(3)*gamma(4)) == choose(5,2) [1] TRUE # all's well so far. # now look what happens: gamma(21)/(gamma(6)*gamma(16)) == choose(20,5) [1] FALSE # check individual terms: gamma(21)/(gamma(6)*gamma(16)) [1] 15504 choose(20,5) [1] 15504 # so they are the

Hi! I have been working several years with R but it's my first public question. I hope I'll be clear :) . This question is related to obtaining letter-based representation of non-parametric pairwise comparisons. I have a dataframe with this structure (but with quite more rows and cols): A B C factor 1 2 2 one 2 1 2 one 2 2 3 two 2 3 2 two 1 4 2 three 9 8 1 three I have no normality,

Hi, I am performing GEV analysis on temperature/precipitation data and want to use the ks.boot function but I am unsure of how to implement it using the pgev distribution. For example: ks.test(data,pgev,shape,location,scale) ks.boot(data,(distribution parameters?),alternative="pgev",nboots=1000) Any advice? Apologies in advance if I have used the wrong email address. Regards, Louise

Sorry for speaking too soon. It is missing from 3.2 which is the version I have been using. If it is already implemented than it is very good news. Thanks. -----Original Message----- From: Eli Bendersky [mailto:eliben at google.com] Sent: Monday, November 26, 2012 9:25 AM To: Manny Ko Cc: Albert Graef; Benjamin Kramer; LLVM Developers Mailing List Subject: Re: [LLVMdev] Removing old JIT

I have two samba servers running Ubuntu 10.04 Samba Version 3.4.7 One server acts as domain controller and stores user ids in a .tdb Somehow I've ended up with a duplicate user name. On the Domain Controller # pdbedit -w -L|grep debbie debbie:1005:XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX:84DEC6FE3B018B0FB977EDDF5009742C:[U ]:LCT-4D4B086F: On the other Server running winbind I get #

Dear list fellows, I want to study time series and use R to analyse time series of fishing data from several species (landings and cpue) investigating the correlation between them and with environmental factors (water temperature, wind, etc.). Searching at Amazon I found three books with examples in R: Time Series Analysis: With Applications in R by Jonathan D. Cryer and Jonathan D. Cryer

Thanks Andrew for the update. Does it support NotifyFunctionEmitted? Cheers. -----Original Message----- From: Kaylor, Andrew [mailto:andrew.kaylor at intel.com] Sent: Monday, November 26, 2012 1:06 PM To: Manny Ko Cc: LLVM Developers Mailing List Subject: RE: [LLVMdev] Removing old JIT CodeEmitters for ARM and PPC The current 3.2 code branch has events being broadcast from the MCJIT engine

Hi, I want to fit some psychophysical data with cumulative gaussians. There is quite a convenient toolbox for matlab called 'psignifit' (formerly known as 'psychofit'). It allows the lower bound of the sigmoid to vary slightly from zero, aswell as the upper bound to vary from one. with these two free parameters, the fitted function is less sensitive to noisy data and outliers.

I see. If I am writing my own listener can I still get back the (llvm::Function, void* Code, size_t Size)? Ciao. -----Original Message----- From: Kaylor, Andrew [mailto:andrew.kaylor at intel.com] Sent: Monday, November 26, 2012 1:24 PM To: Manny Ko Cc: LLVM Developers Mailing List Subject: RE: [LLVMdev] Removing old JIT CodeEmitters for ARM and PPC No. It adds a new event,