Welcome

Who Am I,

I am a biologist who, from the onset, favored command line prompts over pipettes. My doctoral training is nearly complete, and I wish to part unto you knowledge I learned along my journey.


What to expect,

 In these pages you will find tutorials for
  • writing programs for genetic analysis
  • explaining common file formats for biologists
  • helpful time-saving tips
  • fun projects using 23andMe/Ancestry data

I work in a UNIX environment. Hence, most of the tutorials should apply for GNU/Linux / MacOS users. For Windows users, I would recommend a virtual machine running GNU/Linux.

I write primarily in Perl and Python, but will post items on C/C++.

What's the deal with all this Magus talk?

Magus is the Latin word for "Magician", derived from the Greek word "Magos". The Greeks associated Zoroaster (the founder of the Persian faith) with magic, alchemy, and astrology. Subsequently, anyone who practiced in such mysterious arts became known as Magi.

Programming is magic. It's an art that dazzles the minds of those not versed in it. And for those who do practice the art of programming realize that truly magnificent code is undeniably magic.

Feel the power of the command line, let your verbal instructions translate into electrical current. Wield the power of an intelligence far superior than yours in calculation. Accept the magic of programming and may you always conjure benevolent spirits with your Bio Sigils.

Comments

Post a Comment

Popular posts from this blog

Summoning Singletons with plink

Singletons are variants private to individuals, a personal mutation. Sometimes these variants can carry disease risk when passed to offspring. Especially when the risk is additive and multiple disease variants are transmitted to an affected child. Summoning singletons is easy with plink . Here's a basic guideline of how to do it. 1. Convert the VCF into a plink bfile 2. Run plink --freq counts That's it. I included a test case you can download from GitHub $ git clone --recursive https://github.com/biosigil/plink_singletons.git All following commands are assuming you installed plink and changed directory to plink_singletons/ The VCF is located in input/1kgp.vcf.gz and consists of three trios (child and parents). These data are a sampling from the 1000 Genome Project .  Convert the VCF $ plink --vcf input/1kgp.chr21.vcf.gz --make-bed --out 1kgp plink will now make three files, 1kgp.bed, 1kgp.bim, 1kgp.fam The fam file is not correct, the relation inform...
Read more

Custom Usage Message for Python

Spellcrafting in Python is approachable and enjoyable. The library argparse facilitates the parsing of arguments for code destined for wide distribution. Below is an example of basic argparse usage #!/usr/bin/env python3 """ in this example get input from the user and manlipulate it """ from argparse import RawTextHelpFormatter import argparse,sys def get_args (): parser = argparse.ArgumentParser(formatter_class = RawTextHelpFormatter) arg_grp = parser.add_argument_group( 'arguments' ) # add an argument to take an int # set the default to None to check for errors, required arg_grp.add_argument( '-n' ,type = int,default = None) # add a required argument for a string # set the required option to True # this argument has two flags, -s --string arg_grp.add_argument( '-s' , '--string' ,type = str,required = True) # this command loads the args # object with user ...
Read more

Push to GitHub with the Terminal

Push to GitHub with the Terminal GitHub is an online version control repository for code. It facilitates sharing and improvement of code. GitHub uses the Git platform to update code from a terminal. However, setting up a server to push changes to GitHub can be quite a pain. Below is a guide on one way to set up such a server. 1. Generate a deploy key Deploy keys allow a server to push changes to a GitHub repository. Open a terminal window and type the following command. $ ssh-keygen -t rsa -C "deploy key for repo" You will be given a prompt to change the file name. Change it to reflect the name of your repository. For this example, the name will be "repo". Generating public/private rsa key pair. Enter file in which to save the key (/home/usr/.ssh/id_rsa): /home/usr/.ssh/repo_rsa # change the file name to match your repo Enter passphrase (empty for no passphrase): # hit enter for no passphrase Enter same passphrase again: # hit enter for no pa...
Read more