Research | Teaching | Fall11
 Biowiki > Fall11 > JacobVoganHomework8a

 Search all webs web groupthis web topic names Advanced search...

 Topics BioE131/231 BioE241 %SESSION_IF_AUTHENTICATED% %GET_SESSION_VARIABLE{AUTHUSER}% %SESSION_ELSE% Login Register %SESSION_ENDIF% Changes Help
-- JacobVogan - 24 Nov 2011

Homework #8a - Bayesian Analysis of DNA sequence Origin

• ORFprob.pl: ORF Bayesian Perl script - outputs posterior probability on sequence origin

ORFprob.pl Features:

• Reads FASTA formatted sequence file as input
• Computes the probability distribution, P(x), for each nucleic acid "x" in sequence "S".
• Outputs the posterior probability P(G=1|S) and the log ratio of P(G=1|S) / P(G=0|S)

Examples:

• perl ORFprob.pl file.fasta        # Outputs the posterior probability and log ratio of posterior probabilities
• perl ORFprob.pl file.fasta -h        # Prints program information
• perl ORFprob.pl file.fasta -move [#]        # Option to shift the reading frame (e.g. -move 1 will shift reading frame by 1 nucleotide)
• perl ORFprob.pl file.fasta -stop [X]        # Option to enter unique STOP codon (e.g. TGG)

Computation:

Background on sequence origin probabilities P(G=0) and P(G=1) can be found here.

This program first computes the posterior probability P(G=0|S) from the probability P(S|G=0). The prior probability P(S|G=0) is computed using the probability P(G=0) would not generate a STOP codon for the given sequence codon-length and reading frame (this also depends on P(x=T), P(x=A), and P(x=G)). From this, P(G=0|s) is computed by:

The posterior probability P(G=1|S) is then:

Terminal Example of ORFprob.pl:

For test.fasta file with sequences:

I Attachment History Action Size Date Who Comment
txt ORFprob.pl.txt r4 r3 r2 r1 manage 7.0 K 2011-11-30 - 02:04 JacobVogan ORF Bayesian Perl script - outputs posterior probability on sequence origin
 Actions:
 Fall11.JacobVoganHomework8ar27 - 2011-11-30 - 02:04:46 - JacobVogan Biowiki content is in the public domain. Comments on this site? JacobVoganHomework8a">Send feedback