similar to: p-adjust using Benjamn and Hochberg

Hello, I read a book about statistics in psychology. The authors use SPSS. They talk about post hoc tests after ANOVA finds significant effects: - Gabriel's procedure (for equal or slightly different sample sizes) - Hochberg's GT2 (for different sample sizes) - Games-Howell procedure (for populations with unequal variances) I could not find them in R. Do they not exist in R

Hi, I have two questions want to ask. 1. If I have a matrix like this, and I want to figure out the rows whose value in the 3rd column are less than 0.05. How can I do it with R. hsa-let-7a--MBTD1 0.528239197 2.41E-05 hsa-let-7a--APOBEC1 0.507869409 5.51E-05 hsa-let-7a--PAPOLA 0.470451884 0.000221774 hsa-let-7a--NF2 0.469280186 0.000231065 hsa-let-7a--SLC17A5

Hi all, I'm just getting started with R and I would appreciate some help. I'm having trouble creating a boxplot with whiskers at the 95th and 5th percentiles instead of at 1.5 * IQR. I have read the relevant documentation, and checked existing mails on this topic. I found a small modification that should work : https://stat.ethz.ch/pipermail/r-help/2001-November/016817.html and tried to

> On May 18, 2016, at 2:38 PM, Rowland penny <rpenny at samba.org> wrote: > > On 18/05/16 19:44, Smith, Jarrod A wrote: >> We just upgraded to 4.3.9 (from 4.1.x) and are experiencing a few issues/differences around permissions on files written from Windows clients authenticated from winbind/AD. One specific issue that we have is directories with permissions like: >>

Dear all, I constructed this function called my.boxplot.stats by replacing fivnum() with quantile() in function boxplot.stats(). So I can try different quantile methods in bwplot(). The problem is I couldn't pass different values to the "type" argument to my.boxplot.stats, which in turn is an argument in bwplot(). Now I just have to manually change the "type" value in

Dear list, I have a problem with freq from prettyR. Please have a look at my syntax with a litte example: library(prettyR) #Version 1 test.df<-data.frame(q1=sample(1:4,8,TRUE), gender=sample(c("f","m"),8,TRUE)) test.df freq(test.df) #No error message #Version 2 test.df<-data.frame(gender=sample(c("f","m"),8,TRUE), q1=sample(1:4,8,TRUE)) test.df

We just upgraded to 4.3.9 (from 4.1.x) and are experiencing a few issues/differences around permissions on files written from Windows clients authenticated from winbind/AD. One specific issue that we have is directories with permissions like: drwxrws---+ 9 myapp 9997 2048 May 16 17:38 . It's owned by user "myapp" and GID 9997 and as you can see we have the SGID bit set on this

Hello, I know this has been covered on here before, but as a complete novice, I need a little more guidance. I would like to produce boxplots with the whiskers extending to the 10 and 90th percentiles. I found this code: myboxplot.stats <- function (x, coef = NULL, do.conf = TRUE, do.out = TRUE) { nna <- !is.na(x) n <- sum(nna) stats <- quantile(x, c(.1,.25,.5,.75,.9), na.rm

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- Is anyone working on multiple comparisons of means or post-hoc tests (ie: Tukey, Bonferroni) for R? I saw in the winter 98' archives of the R mailing lists that these tests had not been implemented yet, I was just wondering if I could look foward to having them. ;-) I also looked through the contributed packages and didn't see anything that offered such tests. I guess I could check

Hi All.. Am working on continous backup of mailboxes using rsync(for e.g by running rsync every 2 min) Things gets more complicated when users create Subfolders in INBOX , SENT , etc.. If anyone among you did that plz guide & advise. Thanks CoolAtt _________________________________________________________________ Your E-mail and More On-the-Go. Get Windows Live Hotmail Free.

Dear R community I sent a message out a while ago asking for help with multiple comparison tests for ANOVA's, but haven't had any response yet. I'm sending a final desperate plea. If I can't get this done in R I'm going to have to redo a whole lot of stuff in a commercial package, which I'm REALLY not keen to do! My problem is how to implement these tests in R. Below

Matt, Here's a Bonferroni-corrected multiple one-sample t-test that I wrote some years ago. It took a while to get it into R, as na.omit doesn't seem to handle vectors and I had to write a quick kludge (na.remove). Another more general point was that I discovered that the help page for t.test gives the name "parameters" for the degrees of freedom, as in S. However, the name

Hi R folks, Can anyone tell me how to read in a pipe ("|") delimited text file? I've tried the following: read.delim("c:/junk/junk.txt",sep="|", skip=7, check.names=FALSE,quote = "", header=F) The file looks something like the following: RD|I|04|013|9997|68103|5|7|017|830|20000221|00:00|12.6||6|||||||||||||

https://bugzilla.samba.org/show_bug.cgi?id=9997 Summary: File listed only whose attributes have changed for single '-v' option Product: rsync Version: 3.0.9 Platform: All OS/Version: All Status: NEW Severity: normal Priority: P5 Component: core AssignedTo: wayned at

Hi everyone, I think I have found a minor issue with the R function "boxplot.stats". But before I make such a rash comment, I'd like to check my facts by fixing what I think is the problem. However, why I try to do this, R does not behave as I expect. Can you tell me what I'm doing wrong? If I use the command: debug(boxplot.stats) I am allowed to step through the code as

L.S. In the current version of ?p.adjust.Rd, one needs to scroll down to the examples section to find confirmation of one's guess that "fdr" is an alias of "BH". Please find a patch in attachment which mentions this explicitly. Best, Tobias -------------- next part -------------- A non-text attachment was scrubbed... Name: p.adjust.Rd.patch Type: text/x-patch Size: 633

Dear all, It was a pleasure to meet you at Iowa State University. Two days ago I submitted two experimental packages to CRAN (hope it will be there soon): rindex: quick indexing of large objects (currently only character, see ?index) regtest: some first support for automated regression testing (heavily used in \dontshow{} section of ?index) With rindex you can for example i <-

Full_Name: Ludo Pagie Version: 2.8.1 OS: linux Submission from: (NULL) (194.171.7.39) p.adjust in stats seems to have a bug in handling n>length(p) for (at least) the methods 'holm' and 'hochberg'. For method 'holm' the relevant code: i <- 1:n o <- order(p) ro <- order(o) pmin(1, cummax((n - i + 1) * p[o]))[ro] where p is the

Dear R Helpers, I have run into a problem trying to perform a number of actions on a set of quantmod data objects through a loop and I am hoping that this is an easy problem for someone else as opposed to very difficult for me. The example task is to get the first three objects of the quarterly balance sheet for a number of companies from the getFinancials object and put them together into a