Computational biology and genomics
Can the R bioconductor package “WGCNA” find artefactually created modules? Firstly some (subpar) code to generate an artefactual list of genes with co-expression patterns (modules): Running the code: ./generate_random_module.pl 10 1000 20 > 10_sample_1000_list_20_module.tsv Patterns: 1 0 0 0 0 1 0 0 0 1 0 0 0 0 0 0 1 0 1 1…
Is there any nucleotide bias with the -40 region of RefSeqs? Taking all hg19 RefSeqs that mapped to assembled chromosomes (36,004) and extracting the nucleotide sequences 40 bp upstream of the RefSeq gene model, I generated a sequence logo. No obvious TATA box enrichment, which was expected since only 10-20% of genes in eukaryotes have…
Update 15th May 2018: I recommend using the pheatmap package for creating heatmaps. Heatmaps are great for visualising large tables of data; they are definitely popular in many transcriptome papers. Here are the basic commands for making your own heatmap: The colour schemes from the function display.brewer.all() Returning the values used for the heatmap Taken…
Create your own dissimilarity matrix and save it in a file, say test.matrix column1 column2 column3 0 1 3 1 0 2 3 2 0
Foreach item in the array, start a fork. $pid returns 0 if it is the child process. This way you can spawn 4 child processes from one parent. This page explains it all and the code shown below is an adaptation of code shown from the link. I just added code to show the current…