[Bioperl-l] trouble with pairwise_kaks.PLS on Cygwin/XP platform
Caleb F. Davis
cdavis at bcm.tmc.edu
Fri Sep 1 00:00:59 UTC 2006
Hi folks,
Thank you for all your hard work with Bioperl. My name is Caleb Davis,
a first year grad student at Baylor College of Medicine. I'm trying to
use bioperl on a WinXP/cygwin configuration. I've got it working great
for everything up to this point. I can read and write files, use the
wrappers around bl2seq, etc. However, I am now getting starting using
PAML, and I am stuck. I could not find anyone with my problem on the
mailing list so I thought I'd ask you.
I'm using pairwise_kaks.PLS code and sample data to test my machine
configuration first. The code fails in various ways, and I can't figure
out why. I'm guessing it has something to do with the tmpdir
environmental variable and the differences between windows and unix
paths, but that's a stab in the dark really.
Will you point me in the right direction? Thanks,
--Caleb
Included below are:
********
output of the command:
pairwise_kaks.PLS -i ~/src/bioperl/bioperl-live/t/data/worm_fam_2785.cdna
********
my environmental variables set in .bashrc
********
//////////////////////////////////////////////
******************************************
CLUSTAL W (1.83) Multiple Sequence Alignments
Sequence format is Pearson
Sequence 1: CBG10100 363 aa
Sequence 2: F22B7.13 525 aa
Sequence 3: C38C10.4 525 aa
Start of Pairwise alignments
Aligning...
Sequences (1:2) Aligned. Score: 9
Sequences (1:3) Aligned. Score: 8
Sequences (2:3) Aligned. Score: 96
Guide tree file created:
[C:/cygwin/home/cdavis/tmp/qsGfDVyxqJ/wez3HD46tR.dnd]
Start of Multiple Alignment
There are 2 groups
Aligning...
Group 1: Sequences: 2 Score:11215
Group 2: Delayed
Sequence:1 Score:2745
Alignment Score 3187
GCG-Alignment file created
[c:/cygwin/home/cdavis/tmp/qsgfdvyxqj/ec6g4yeljx]
Can't remove directory C:/cygwin/home/cdavis/tmp/Yqsvs4wNAX: Directory
not empty at /usr/lib/perl5/site_perl/5.8/Bio/Root/IO.pm line 811
Use of uninitialized value in pattern match (m//) at
/usr/lib/perl5/site_perl/5.8/Bio/Tools/Phylo/PAML.pm line 615, <GEN2>
line 105.
Can't call method "get_MLmatrix" on an undefined value at
./pairwise_kaks.PLS line 177, <GEN2> line 105.
******************************************
#BioPerl settings
export BIOPERL_INDEX_TYPE=SDBM_File
export BIOPERL_INDEX=~/bioperlindex
export TMPDIR=C:/cygwin/home/cdavis/tmp/
#Blast settings
export PATH=~/src/ncbi/blast-2.2.14/bin:${PATH}
export BLASTDIR=~/src/ncbi/blast-2.2.14
export BLASTDB=~/src/ncbi/blast-2.2.14/data
#other program settings
export CLUSTALDIR=~/bioperl-runnables/clustalw1.83
export PAMLDIR=~/bioperl-runnables/paml3.15/bin
#FTP settings
export FTP_PASV=1
********************************************
More information about the Bioperl-l
mailing list