Search Archives

Hi R users, I want to do a data reshape from long to wide, I thought it was easy using tidyverse spread function, but it did not work well. Can you help me? Thank you, Ding test1 data frame is...

Hi Bert, Yes, your code works very well. Thank you very much!! Ding ________________________________ From: Bert Gunter [bgunter.4567 at gmail.com] Sent: Wednesday, March 18, 2020 1:02 PM To: Yuan Chun Ding Cc: r-help mailing list Subject: Re...

When you do x1g$fdr, you adjusted for a total of 15568 tests, so FDR is high for those first four entries. When you do p.adjust(pval,method="BH"), you assumed there were only a total of 4 multiple...

Hi Jeff, Thank you for your quick response!! I made this following summary file, however, I want to add the count number sequentially into the original file, because I have several more columns to explain or annotate the first two...

Hi Thierry, the values in the example data frame are fake numbers, my original data frame has hundreds of row and values are in wide range, not min or max of two variables, also the number 23 is also different...

Hi Rui, Thank you so much!! You code works well and I am looking into the pivot_wider function. Yuan Ding -----Original Message----- From: Rui Barradas [mailto:ruipbarradas at sapo.pt] Sent: Thursday, January 7, 2021 12:19 PM To...

Hi Bert, I am very sorry that I made a mistake for the second count value. The original data are much more complicated, I just made a toy table. test <- data.frame(group1=c("g1", "g1", "g1", "g2", "g2", "g2...

Hi R users, I want to create a new variable, Ravg, in data frame tem2 based on values of two other variables m1 and m2. the condition: if m1 = 23 and m2 =23 then Ravg =23; else if m1 != 23...

Hi Jeff, Your code generated the count what I wanted. As I emailed to Bert, yes, I made a mistake when generating the toy table, the second count value should be 2. Correct file is test <- data.frame(group1=c...

Hi Jim, Thank you for the message! Yes, you are right, I contacted Dr. Xu Zongli, he said he removed the ComBat.mc function from the new version. Ding From: Jim Lemon [mailto:drjimlemon at gmail.com] Sent: Wednesday, November...

Hi Bert, Sorry that I was in a hurry going home yesterday afternoon and just posted my question and hoped to get some advice. Here is what I got yesterday before going home. --------------------------------------------------------------- setwd("C:/Awork/VNTR/GETXdata/GTEx_genotypes...

You are right. I also just thought about that, no intercept is not applicable to my case. Ding From: Bert Gunter <bgunter.4567 at gmail.com> Sent: Saturday, August 10, 2024 1:06 PM To: Yuan Chun Ding <ycding at...

Hi Denes, thank you very much for your mesasge!! I would like to learn from you. We generated preliminary data for a new project, so performing explorative analysis. since only sample size is relative small, I probably will both Cox...

Hi Richard, Nice to know you had similar experience. Yes, your understanding is right. all correlations are negative after removing double-zero rows. It is consistent with a heatmap we generated. 1 is for a cancer patient with a specific...

Hi Jim, Yes, you are right. I sorted the tem4$Var1 first, then find rising peaks in Freq variable from left to right. I guess I probably need to define the minimal rising and drop on both side of a...

Hi Bert, No, this Is not home work related. Original data have 87352 rows. I used the standard reshape function and got warning message. So I reformatted the wide format to meet my research purpose. mut2 <-mut[,c("Tumor_Sample...

Hi Avi, Thank you for your patient with me!! Your code below generated what I want. test %>% group_by(group1,group2) %>% mutate(count2 = ifelse(row_number()==1, TRUE, FALSE)) %>% ungroup() your code reminded me of the SAS first. And last...

Hi Abby, Thank you so much for your effort!! I really appreciate your help!! I modified your code a little to get both maxv and corresponding alleles in var1. For now, I do not care too much for identifying small...

Hi Bert and Rui, Thank you very much! I had thought for every condition of matrix ==0 , rnrom(n=1, 1, 0.1) will randomly generate a number with mean 1 and sd 0.1. Yuan From: Bert Gunter [mailto...

Hi Rui, Thank you very much! Yes, I verified using real data, it worked correctly as expected after adding tidyr:: to the pivot_longer function and dplyr:: to the group_by and summarize Function. I did not know how to...