hello,
How can I find position specific probability matrix from a multi fasta file in R.
I know BioStrings in R as consensusMatrix. But what I would like to do is read all the sequences from a single file and generate a PSPM matrix as output.
Any suggestions? How to achieve this?
•
written
by
k.nirmalraman
1 answer
Given that they are aligned, all you need to do is loop through the strings, and at each position add up the number of each base.
e.g.
12345678
GATAGACC
GTTAGACG
GAAAGACG
-
--| 1 2 3 4 5 6 7 8
--+----------------
A | 0 2 1 3 0 3 0 0
T | 0 1 2 0 0 0 0 0
C | 0 0 0 0 0 0 3 1
G | 3 0 0 0 3 0 0 2
You could also add a laplacian correction (http://cis.poly.edu/~mleung/FRE7851/f07/naiveBayesianClassifier.pdf) by adding a constant to each cell. Then you may normalize each column if you wish.
Log in to answer this question.
-
Extracting the aligned sequence from a BLAST output in a multi-fasta file
written by Ming •
-
Create PSI-BLAST PSSM from position-specific weight matrix (probabilities)
written by M •
-
BAC annotation strategy
written by abhijit.synl
-
PPM, PSSM, and PWM
written by bharata1803
-
RSAT problem with matrix conversion
written by chemcehn
-
Position Frequency Matrix to Position Weight Matrix, whose colsums need to be 1.
written by Zhilong Jia
-
Position Weight Matrix in R or Python
written by esantosh10 •
-
miRNA isomiR analysis using Biostrings
written by Palgrave •
-
Fasta Length
written by k.nirmalraman
-
How To Create Genbank File(.Gbk/.Gb) Format From Multi-Fasta Format File?
written by Kits •
