Question: Position Weight Matrix in R or Python
2
gravatar for esantosh10
5.4 years ago by
esantosh1030
Germany
esantosh1030 wrote:

Hello All,

I would like to create a position weight matrix in R or Python for some peptide sequences . I found that there are packages in R to create a position weight matrix for gene sequences(ATGC) like Biostrings but was difficult to find them for proteins/peptides. Any suggestions?

peptide pwm matrix • 8.4k views
ADD COMMENTlink modified 4.2 years ago by kamil_cygan20 • written 5.4 years ago by esantosh1030
2
gravatar for madk00k
5.4 years ago by
madk00k340
Heidelberg
madk00k340 wrote:

Biopython allows to create Position Weight Matrices using motifs module:

http://biopython.org/DIST/docs/tutorial/Tutorial.html#sec244

As far as I know, there are know constrains for the alphabet, so it could polypeptide.

 

ADD COMMENTlink modified 5.4 years ago • written 5.4 years ago by madk00k340
2
gravatar for kamil_cygan
4.2 years ago by
kamil_cygan20
United States
kamil_cygan20 wrote:

Actually, you can create PSSM for amino acids in a similar manner as for DNA/RNA sequences.

library(Biostrings)
sequences <- c("ERLRELETKIGSLQEDLC","CERLRELETKIGSLQEDL","ERLRECLETKIGSLQEDL")
string_set <- AAStringSet(sequences)
PSSM <- consensusMatrix(string_set, as.prob = TRUE)
ADD COMMENTlink modified 4.2 years ago • written 4.2 years ago by kamil_cygan20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 870 users visited in the last hour