similar to: Games-Howell, Gabriel, Hochberg

Dear R users, I am conducting multiple comparisons among 12 groups (after a significant F-test) that are heteroscedastic (as judged by a significant Levene's test). It seems from the literature that the Games-Howell post-hoc test is the most appropriate for these data - but, I can't seem to locate a function on the R-site or CRAN or the help lists that implements this test in any R

Hello, I am trying to use the p.adjust function for multiple testing. here is what i have 9997 201674_s_at 0.327547396 9998 221013_s_at 0.834211067 9999 221685_s_at 0.185099475 I import them from excel have have the gene symbol as well as the pvalue here is the issue > pa<-p.adjust(pt,method="BH") Error in p[nna] : object is not

Hello, I use t.test for normal distributed and wilcox.test for non-normal distributed samples. It is easy to write a function for t.test that calculates the effect size, because all parts of the formula are available from the t.test result: r = sqrt(t*t / (t*t + df)) However, for Wilcoxon tests, the formula for effect sizes is: r = Z / sqrt(N) I wonder how I can calculate the Z-score in R for

How are we all documenting complex dialplan? Is there something similar to Doxygen? I've got around 20 config files covering around 60 contexts and 40 variables. Of course, I've maintained a basic list of the major stuff, and documented the code throughout, but it's grown to the stage where it needs to be better documented, have a proper flowchart etc. Talking of flowcharts, I see

Hi, thanks - that looks really good! I was about to embark on some non-visual stuff using Ragic, but this looks great. Is there a binary anywhere, or any instructions to compile? I've never compiled C# code before, and although a quick google suggests it shouldn't be too hard, I might need to know a few things like what version of .net it should be compiled with. The readme just points

I am posting this to both R and BioC communities because I believe there is a lot of confusion on this topic in both communities (having searched the mail archives of both) and I am hoping that someone will have information that can be shared with both communities. I have seen countless questions on the BioC list regarding limma (Bioconductor) and its calculation of FDR. Some of them involved

I am posting this to both R and BioC communities because I believe there is a lot of confusion on this topic in both communities (having searched the mail archives of both) and I am hoping that someone will have information that can be shared with both communities. I have seen countless questions on the BioC list regarding limma (Bioconductor) and its calculation of FDR. Some of them involved

ALCON I'm trying to figure out how to rename groups in a data frame after groups by selected variabels. I am using the dplyr library to group my data by 3 variables as follows # group by lat (StoreX)/long (StoreY) priceStore <- LapTopSales[,c(4,5,15,16)] priceStore <- priceStore[complete.cases(priceStore), ] # keep only non NA records priceStore_Grps <- priceStore %>%

Hello, See if this is it: priceStore_Grps$StoreID <- paste("Store", seq_len(nrow(priceStore_Grps)), sep = "_") Hope this helps, Rui Barradas On 5/26/2018 2:03 PM, Jeff Reichman wrote: > ALCON > > > > I'm trying to figure out how to rename groups in a data frame after groups > by selected variabels. I am using the dplyr library to group my

I have a large dataset and I would like to make some kind of flowchart from this dataset. The idea is to show rowcounts from data subsets: data1 rows= 100000 | / \ / \ males females, rows=500000 rows=500000 /\ / \

Hi, I produce a series of diagrams with R in order to include them in my documents (LaTeX). However, there is a white border around the diagrams. For some that do not have anything written at the very bottom, the white border is relatively large. The rather big space between figure and caption at the final document looks not nice. It would be best not to have any white border. I played with

Hi, I have a data frame named "questions" that I use to get a subset and then a table: failcondonly <- subset(questions, errorreason=="condition") failcondtab <- table(failcondonly$type, failcondonly$qid, exclude=c("ORGA", "skipped")) The failcondtab looks like this: 6 11 12 13 14 15 17 26 30 31 39 41 gave up 0 1 1 1 1 0 0 0 0 8 0

Hello, Sorry, but I think my first answer is wrong. You probably want something along the lines of sp <- split(priceStore_Grps, priceStore_Grps$StorePC) res <- lapply(seq_along(sp), function(i){ sp[[i]]$StoreID <- paste("Store", i, sep = "_") sp[[i]] }) res <- do.call(rbind, res) row.names(res) <- NULL Hope this helps, Rui Barradas On 5/26/2018

Hello, I have a table (1) of the form q1 q3 q4 q8 q9 A 5 2 0 1 3 B 2 0 2 4 4 I have another table (2): q1 q2 q3 q4 q5 q6 q7 q8 q9 C 10 7 4 2 6 9 3 1 2 I would like to divide the numbers in table (1) by the number of the appropriate column in table (2): q1 q3 q4 q8 q9 A 5/10 2/4 0/2 1/1 3/2 B 2/10 0/4 2/2 4/1 4/2

Hello. Hope you can help me further. I can vaguely remember that once I saw a reference how to preconfigure debconf in puppet manifests. So you can define answers in debconf inside of manifests. No idea where anymore. Anybody knows? greetinx Christo

Hello dear R help group, I wish to create a flowchart of my R code (specifically for profiling the connection of different functions). Maybe something similar to this: http://www.fatesoft.com/s2f/ But for R, and that it will be free (open source would be best). After some searching, I found this package: http://cran.r-project.org/web/packages/diagram/ But couldn't find a way to this

Why doesn''t the yum provider support version, I tried added ensure -> ''ver'' and it barfs up (on clients) with: <snip> err: //stan.gbuild.org/any-host/java/Package[j2sdk]/ensure: change from 1.4.2_12-fcs to 1.4.2_13-fcs failed: Package provider yum does not support specifying versions at /etc/puppet/manifests/classes/ java.pp:13 </snip>

Hello everybody, I use the sweetpotato database included in R package: data(sweetpotato) This dataset contains two variables: yield(continous variable) and virus(factor variable). Due to Levene test is significant I cannot assume homogeneity of variances and I apply Welch test in R instead of one-way ANOVA followed by Tukey posthoc. Nevertheless, the problems come from when I apply posthoc

Hello, I have been running multiple exec scripts that could run in parallel. Does puppet run as much as possible in parallel? If it does not, is there a way to specify that actions may be run in parallel? Also, is puppet a multiple or single thread application? Thank you, -- Rob -- ____________________________________________________________________________________ Sucker-punch spam

Dear List, We are planning a genotyping study to be analyzed using false discovery rates (FDRs) (See Storey and Tibshirani PNAS 2003; 100:9440-5). I am interested in learning if there is any consensus as to how many features (ie. how many P values) need to be studied before reasonably reliable FDRs can be derived. Does anyone know of a citation where this is discussed? Bill Dupont William D.