[Bioperl-l] Problem with Bio::PopGen::Individual

Sheri Simmons sheris at eps.berkeley.edu
Fri Jul 13 18:42:32 UTC 2007


Hi,
I have a collection of sequencing reads aligned with a consensus sequence that 
I input into a Bio::PopGen::Population object in order to calculate allele 
frequencies. The consensus sequence is included to force clustalw to give a 
better alignment. However,  I need to remove the consensus sequence before 
calculating allele frequencies in the individual reads. I'm having trouble 
with this part of it. I get the following error message:

"Can't locate object method "person_id" via package "Bio::PopGen::Individual" 		
at /usr/share/perl5/Bio/PopGen/Population.pm line 260, <GEN0> line 49."

Here is the code snippet producing the error. $pop is a 
Bio::PopGen::Population object.

	my @consensus = "gene_consensus";
	$pop->remove_Individuals(@consensus);

I also tried:
	my @consensus = $pop->get_Individuals(-unique_id => "gene_consensus"); 
	$pop->remove_Individuals(@consensus);

which produced the same error. Can anyone send me in the right direction? I 
suspect this is a simple problem.

Sheri

-- 
Sheri Simmons
Department of Earth and Planetary Sciences
University of California, Berkeley
Berkeley, CA 94720-4767



More information about the Bioperl-l mailing list