Difference between revisions of "Dr. Rasmus Nielsen Laboratory"
From Christoph's Personal Wiki
(→Overview of research: Added more categories to "MrBayes") |
(format) |
||
Line 9: | Line 9: | ||
==Overview of research== | ==Overview of research== | ||
− | * Taxonomy (from GenBank files) | + | * '''Taxonomy''' (from GenBank files) |
− | * ClustalW (align top 50; ignore ''E''-value; remove redundancy) | + | * '''ClustalW''' (align top 50; ignore ''E''-value; remove redundancy) |
− | ** clustalw all.fsa | + | ** <tt>clustalw all.fsa</tt> |
− | ** clustalx all.aln | + | ** <tt>clustalx all.aln</tt> |
− | ** njplot all.ph | + | ** <tt>njplot all.ph</tt> |
− | * MrBayes | + | * '''MrBayes''' |
** [[Bayesian inference]] | ** [[Bayesian inference]] | ||
** [[Phylogenetics]] | ** [[Phylogenetics]] | ||
Line 20: | Line 20: | ||
** [[Markov chain Monte Carlo]] | ** [[Markov chain Monte Carlo]] | ||
** [[NEXUS file format]] | ** [[NEXUS file format]] | ||
− | * Find probability of each group | + | * '''Find probability of each group''' |
Nucleotide model: GTR (option: +gamma) | Nucleotide model: GTR (option: +gamma) |
Revision as of 23:32, 26 December 2005
Dr. Rasmus Nielsen Laboratory is where I am currently doing research (December 2005-present). It is located at the Centre for Bioinformatics, København Universitet, Denmark.
Contents
Professor Information
- Rasmus Nielsen (Ole Roemer Fellow)
- Office: 314
- Phone: +45 3532 1279
- E-mail: rasmus@binf.ku.dk
Overview of research
- Taxonomy (from GenBank files)
- ClustalW (align top 50; ignore E-value; remove redundancy)
- clustalw all.fsa
- clustalx all.aln
- njplot all.ph
- MrBayes
- Find probability of each group
Nucleotide model: GTR (option: +gamma) Constrain all phylogenetic groups to be monophyletic Make option: strict clock trees (uniform) (prset brlenspr=clock:uniform)
1,000,000 updates (e.g. cycles) discard first 50,000 as burn-in sample a total of 10,000 trees (say)
Then process output to get probabilities of each possible phylogenetic assignment of query sequence.
In MrBayes - if you run the program with the constraints of monophyletic groups - you are forcing the query sequence _not_ to be part of these groups. So that won't work. Instead, run it without the constraints and simply check how often the query sequence is member of a partition that only contains one particular phylogenetic group as members (at any phylognetic level, e.g. species, genus, family, order, etc). Remember, MrBayes will output the probabililities of specific groups (or partitions) directly. So you don't have to do anything with trees yourself.
External links
Taxonomy
- PLANTS Database — by the USDA
- Wikispecies
ClustalW
MrBayes
- MrBayes: Bayesian Inference of Phylogeny — Official Homepage
- MrBayes Wiki — ('Christoph')