Search Archives

Hi R users, I want to find multiple modes (10, 8, 149) for the following vector. freq =1,2,5,5 10,4,4,8,1,1,8,8,2,4,3,1,2,1,1 138 149 14,1...

Hi R Users, I tried to use sql in R; test1 works well, library(sqldf) test1 <- sqldf("select * from r where V1 like '%GCCATGTCAGCACACTACC%TGAAACCTTTAACTATTT%'") test2 <- sqldf("select * from r where V1 like markerinfo$flank1[1] ") but I want to...

Hi R users, Does any one know a R package for genetic analysis of VNTR (variable number of tandem repeats) genotype data? For example, performing Hardy Weinberg equilibrium test, calculating heterozygosity and polymorphism information content etc. I am sorry if...

Thank you for pointing me there. Ding -----Original Message----- From: peter dalgaard [mailto:pdalgd at gmail.com] Sent: Thursday, March 05, 2020 2:21 PM To: Yuan Chun Ding Cc: r-help at r-project.org Subject: Re: [R] how...

Thank you, Sarah, I modified your suggestion a little and it works. Ding -----Original Message----- From: Sarah Goslee [mailto:sarah.goslee at gmail.com] Sent: Sunday, June 28, 2020 2:01 PM To: Yuan Chun Ding Cc: r-help at...

Hi R users, in a Shiny test code, choices = list(1, 2, 3), it generates three components in the choice list, my code works fine and generated tables and plots I want. now I want to modify to choices = list...

Dear R users, I am running the following code below, the gem751be.rpkm is a dataframe with dim of 751 samples by 35164 variables, 73 phenotypic variables in the furst to 73rd column and 35091 genomic variables or genes in...

Hi R users, We like to use shiny program to share data view between collaborators in different universities. When running shiny, I got the following errors, so I just found that I am allowed to upload up to five applications...

Hi R users, I generated 2972 list files in R, each list includes four data frame files , file names for those list file are VNTR13576, VNTR14689, etc. the second data frame in each list has the same 11 column names...

Hi R users, In this test file, test <- data.frame(group1=c("g1", "g1", "g1", "g2", "g2", "g2", "g2", "g2", "g2"), group2=c("k1", "a2", "a2", "c5", "n6", "n6", "n6", "m10","m10"), count= c( 1, 1, 2, 1, 2, 2...

Dear R user, I want to reshape a long data frame to wide format, I made the following example files. Can you help me? Thank you, Yuan Chun Ding sample <-c("xr" , "xr" , "fh" , "fh" , "fh" , "uy" , "uy" , "uy" , "uy...

Hi R users, I am analyzing miRNA sequence data for a special network analysis, I need to replace zero values in a matrix as random numbers with mean of 1 and standard deviation of 0.1. er.miRCounts[er.miRCounts...

Hi R users, The following code worked well to summarize four data groups in a dataframe for three variables (t_depth, t_alt_count, t_alt_ratio), 12 columns of summary, see attached. However, after running another 2000 lines of...

Dear R user, The following qqunit.plot function generated correct qq plot, however, I want to make axis label (1 ,2 ,.....8) have larger size for publication. I tried to add cex.axis =1.2 code following the pch =20...

Hi R users, I thought the follow R code should work, but I got error, Can you fix my code? Thank you, Ding outlier_tcga_MAD3 <- outlier_tcga %>% filter(n_two >0) %>% mutate(freqMAD3_gain2ratio = N_MAD3_gain2/n_two...

Hi Rolf, Yes, choices = as.list(seq(1,5211,1)) is what I want. Thank you so much, Ding ________________________________________ From: Rolf Turner [r.turner at auckland.ac.nz] Sent: Wednesday, January 1, 2020 3:29 PM To: Yuan Chun Ding...

Hi R users, I have a folder (called genotype) with 652 files; the file names are GTEX-1A3MV.out, GTEX-1A3MX.out, GTEX-1B8SF.out, etc; in each file, only one column of data without a header as below 201...

Hi Rui, Thanks a lot, i got this error, I have library(tidyverse). Ding Error in pivot_wider(., id_cols = "vntr1", names_from = "group", names_prefix = "a", : could not find function "pivot_wider" ________________________________________ From: Rui Barradas [ruipbarradas at sapo.pt...

Hi All, I am very sorry to bother you. This week, I updated my R from 3.5.3? to 4.0.3, reinstalled all required R packages and reran the following R code; then got the following error: could...

Dear R users, I generated a boxplots in combination of dotplots using the R code below for the attached test file. boxplot(value~score, data = test, outpch = NA, xlab="",ylab="",xaxt='n', cex.lab=1.2, cex.axis=1.2...