I have 1000+ protein sequences. I want to generate random sequences using a Markov model based on residue transitions found my sequences. I'm told Matlab will make a Markov chain based on multiple sequences, but I would like to use a free alternative to Matlab (python, ruby, R, etc). Can anyone provide me with a library or module?
1 answer
It seems Python programmers like writing Markov generators. I often see this topic pop up on Python blogs in the context of generating pseudo random text. A quick search shows a few hits:
I guess you would only need to change to tokenizer to split on letters rather than words.
•
written
by
Istvan Albert
For 1st and 2nd order chains you can use Sean Eddy's Squid lib. It's written in C. But, you could add some randomizations as a control. In this case you could use uShuffle to preserve high-level orderings. Then you'll see why biology is very much local . . .
All links, but the last one, are dead.
pretty sad state of the internet - you can't trust anything to be around
Log in to answer this question.
-
Building a hidden markov model from synthetic protein designs
written by v.berriosfarias •
-
Trouble with Biopython residue depth
written by Anand KR •
-
How to extract subset of protein structure (PDB format) file based on a subsequence of protein
written by gundalav
-
Biopython: Retrieving atomic coordinates of a model from a pdb file
written by souparnoa91 •
-
multiple sequence alignment using HMM and simulated annealing
written by Rojhin •
-
How to run multiple alignment and SNP-calling of WGS data in .gb and .fasta using Python or Ruby/J…
written by fashiondesignrussian •
-
How Do I Generate A -Bfile For Meme?
written by a1ultima
-
Simulation Protein Sequences
written by User 1933 •
-
Why Doesn'T The "Markov Property" Disqualify Markov Chains For Sequence Analysis?
written by Ib84 •
-
Visualize Effect Of Snp Variant On Protein Structure
written by Hanif Khalak
