TODO: Write a project description
TODO: Describe the installation process
| ##################################################### | |
| ## (A) Add fingerprint folding argument to desc2fp ## | |
| ########################################## | |
| ## Intersect and Venn Diagram Functions ## | |
| ########################################## | |
| ## Author: Thomas Girke | |
| ## Last update: March 24, 2012 | |
| ## Utilities: | |
| ## (1) Venn Intersects | |
| ## Computation of Venn intersects among 2-20 or more sample sets using the typical |
| ##################################################### | |
| ## (A) Add fingerprint folding argument to desc2fp ## | |
| ########################################## | |
| ## Intersect and Venn Diagram Functions ## | |
| ########################################## | |
| ## Author: Thomas Girke | |
| ## Last update: March 24, 2012 | |
| ## Utilities: | |
| ## (1) Venn Intersects | |
| ## Computation of Venn intersects among 2-20 or more sample sets using the typical |
| my_frame[!duplicated(my_frame[,2]),] # Removes rows with duplicated values in selected column. | |
| my_frame[my_frame$y2 > my_frame$y3,] # Prints all rows of data frame where values of col1 > col2. Comparison operators are: == (equal), != (not equal), >= (greater than or equal), etc. Logical operators are & (and), | (or) and ! (not). | |
| x <- 0.5:10; x[x<1.0] <- -1/x[x<1.0] # Replaces all values in vector or data frame that are below 1 with their reciprocal value. | |
| x <-data.frame(month=month.abb[1:12], AB=LETTERS[1:2], no1=1:48, no2=1:24); x[x$month == "Apr" & (x$no1 == x$no2 | x$no1 > x$no2),] # Prints all records of frame 'x' that contain 'Apr' AND have equal values in columns 'no1' and 'no2' OR have greater values in column 'no1'. | |
| x[x[,1] %in% c("Jun", "Aug"),] # Retrieves rows with column matches specified in a query vector. | |
| x[c(grep("\\d{2}", as.character(x$no1), perl = TRUE)),] # Possibility to print out all rows of a data frame where a regular expression matches (here all double digit values in col 'no1'). | |
| x[c(g |
| ############################# | |
| ## Source R code from Gist ## | |
| ############################# | |
| ## Display code snippet as raw by pressing Raw button on top right corner of Gist | |
| ## Save URL in address bar of browser | |
| # library(RCurl) | |
| ## (a) Source specific commit of code | |
| # source(textConnection(getURL("https://gist.githubusercontent.com/tgirke/c26daac0f647f1732a58/raw/521292eca2d4daeab8911df2552f0c55b0e58e6b/Source_Gist_in_R"))) | |
| ## (b) Source latest version of code (delete commit ID in URL after .../raw/...) | |
| # source(textConnection(getURL("https://gist.githubusercontent.com/tgirke/c26daac0f647f1732a58/raw/Source_Gist_in_R"))) |
| ############################ | |
| ## appendCounter Function ## | |
| ############################ | |
| ## Author: Thomas Girke | |
| ## Last update: 04-Oct-15 | |
| ## Function to append occurrence counter to entries in character | |
| ## vector and return the results as named vector where the | |
| ## original data are in the same order in the data slot | |
| ## and the counting result in the name slot. |
| ##################################### | |
| ## Quality Trimming of FASTQ Reads ## | |
| ##################################### | |
| ## Author: Thomas Girke | |
| ## Last update: May 30, 2016 | |
| ## Usage of below function combined with preprocessReads form systemPipeR: | |
| # qcTrim <- "qualityTrimming(fq, phred_cutoff=20, cutoff_occurrences=1, N_cutoff=1, minreadlength=100)" | |
| # preprocessReads(args=args, Fct=qcTrim, batchsize=100000, overwrite=TRUE, compress=TRUE) | |
| ## Arguments: |
| ######################################## | |
| ## Import of DrugBank Annotation Data ## | |
| ######################################## | |
| ## Function to import DrugBank XML to data.frame | |
| ## Last step gives error. To debug, the following function may help. | |
| ## Note, this functions needs some major speed improvements. Ideally, | |
| ## it should be replaced with a standard XML import method. | |
| ## (1) Download | |
| ## - download DrugBank XML (https://www.drugbank.ca/releases/latest) | |
| ## - name uncompressed file 'drugbank.xml' |
| x <- 1:12 |