[Bioperl-l] seqboot and random seeds
june tantoolvesm
junet9876 at yahoo.com
Mon Jun 14 00:33:10 EDT 2004
Hi,
I was just wondering if someone could please tell me
how random seeds are handled in the SeqBoot module.
For example, is there a way to find out what random
number is being used, is it the same each time?
Thank you.
Best regards,
June
--- Joe Felsenstein <joe at gs.washington.edu> wrote:
> Date: Sun, 13 Jun 2004 08:38:09 -0700
> From: Joe Felsenstein <joe at gs.washington.edu>
> To: june tantoolvesm <junet9876 at yahoo.com>
> Subject: Re: bootstrapping related question
>
>
> June Tantoolvesm --
>
> > I have automated a few of your programs into a
> > pipeline using bioperl
> > (seqboot->protdist->neighbor->consense). My
> problem
> > is, the final trees I get from when I run the
> programs
> > manually and from when I use bioperl, are not the
> > same. I tried several times.
> >
> > Is it normal to not get the same tree everytime we
> run
> > the pipeline of programs
> > (seqboot->protdist->neighbor->consense). (But
> actually
> > I do get the same tree when I run the programs
> > manually a few times, so I'm confused.) Also, does
> > this have anything to do with the random seed
> number
> > we input for seqboot?
>
> You should get the same tree each time when you run
> the
> Bioperl script, provided that it properly renames
> files,
> and provided that the random number seeds used are
> the
> same each time. There is no reason to expect them
> to
> be different.
>
> Check on the issue of random number seeds, and
> whether the
> output files have been renamed to input files
> properly.
>
> J.F.
> ----
> Joe Felsenstein joe at gs.washington.edu
> Department of Genome Sciences and Department of
> Biology,
> University of Washington, Box 357730, Seattle, WA
> 98195-7730 USA
>
__________________________________
Do you Yahoo!?
Friends. Fun. Try the all-new Yahoo! Messenger.
http://messenger.yahoo.com/
More information about the Bioperl-l
mailing list