seqLogo package: sequence logo graph with background frequency?
0
0
Entering edit mode
Tim Smith ★ 1.1k
@tim-smith-1532
Last seen 7.9 years ago
Hi, I had a question about the package/function seqLogo. I had a set of sequences for which I want to plot the sequence Logo. However, the background frequencies are not 1/4. I have another matrix with the background frequencies. How can I incorporate these in the seqLogo too? For example, my sequence frequency matrix is: > freqmat         1       2       3       4       5       6       7       8 9      10 A 0.20000 0.20000 0.06667 0.13333 0.06667 0.33333 0.26667 0.13333 0.20000 0.40000 G 0.13333 0.06667 0.13333 0.33333 0.26667 0.20000 0.20000 0.20000 0.06667 0.33333 C 0.53333 0.66667 0.46667 0.40000 0.46667 0.40000 0.46667 0.53333 0.66667 0.06667 T 0.13333 0.06667 0.33333 0.13333 0.20000 0.06667 0.06667 0.13333 0.06667 0.20000 and my background frequency matrix is: > background_freq_mat         1       2       3       4       5       6       7       8 9      10 A 0.22182 0.22320 0.22473 0.22234 0.22258 0.22234 0.22053 0.22102 0.22335 0.22099 G 0.27851 0.28159 0.27729 0.28225 0.28118 0.28125 0.28052 0.28357 0.27801 0.28031 C 0.28063 0.27739 0.28000 0.27885 0.27989 0.27691 0.28134 0.27925 0.27756 0.28011 T 0.21904 0.21782 0.21797 0.21655 0.21635 0.21950 0.21761 0.21615 0.22109 0.21859 I know that I can get the seqLogo using an assumed background frequency of 1/4 for each base by: > seqLogo(freqmat) But how do I incorporate the background frequencies too? If I cannot do it with seqLogo, is there any other package available that might do the trick? thanks!! [[alternative HTML version deleted]]
seqLogo seqLogo • 949 views
ADD COMMENT

Login before adding your answer.

Traffic: 225 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6