-- %TEACHINGWEB%.TaiNg - 20 Nov 2011
Information Content Homework
To implement:
info.pl FILE
I stored a lot of the information in a matrix. if you want to see the matrix, unsupress print in line 191-198
The matrix printed out will be in the form:
SEQUENCE1
SEQUENCE2
SEQUENCE3
Number of A
Number of U
Number of C
Number of G
Number of Gap
Probability of A
Probability of U
" " of C
" " of G
" " of Gap
Entropy
To see the matrix, unsupress line 191-198
To see the joint distribution: Unsupress line 149-1521
To see mutual information, unsupress 170-173
AnalysisOfData
Extra Credit
This 3D graph plots the Column i and Column j number as well as the mutual information shared between them. This is done by MATLAB surface fitting tool. As you can see, the spikes are where the mutual information is highest between the two columns. The peaks are where pairing is most likely to occur because they share the most information in most of the sequences. .
This 2D plot plots Column number versus the entropy of that column. Entropy is lowest when the that area of the sequence is the most conserved. As one can see, the most conserved area is probably position 61-71. Important structures of the tRNA may be encoded there.

Copyright © 2008-2013 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki?
Send feedback