03-24-2011
@Bartus11
THANK YOU VERY MUCH like always. I was expeting part of the previous code to repeat to get to the value of genotype. I need to better understand the loops. Just one question .... why did you use a hash and not just a simple array to hold the values??
Cheers
... hv a nice day
10 More Discussions You Might Find Interesting
1. Shell Programming and Scripting
Hello,
I have a bash shell script and I use awk to print certain columns of one file and direct the output to another file. If I do a less or cat on the file it looks correct, but if I email the file and open it with Outlook the lines outputted by awk are concatenated.
Here is my awk line:... (6 Replies)
Discussion started by: xadamz23
6 Replies
2. Shell Programming and Scripting
Hi,
I have 2 files
file1
********
01-05-09|java.xls|
02-05-08|c.txt|
08-01-09|perl.txt|
01-01-09|oracle.txt|
********
file2
********
01-02-09|windows.xls|
02-05-08|c.txt|
01-05-09|java.xls|
08-02-09|perl.txt|
01-01-09|oracle.txt|
******** (8 Replies)
Discussion started by: shekhar_v4
8 Replies
3. UNIX for Advanced & Expert Users
Hi All,
I have a very huge file (4GB) which has duplicate lines. I want to delete duplicate lines leaving unique lines. Sort, uniq, awk '!x++' are not working as its running out of buffer space.
I dont know if this works : I want to read each line of the File in a For Loop, and want to... (16 Replies)
Discussion started by: krishnix
16 Replies
4. Shell Programming and Scripting
Hi friends,
I have multiple files. For now, let's say I have two of the following style
cat 1.txt
cat 2.txt
output.txt
Please note that my files are not sorted and in the output file I need another extra column that says the file from which it is coming. I have more than 100... (19 Replies)
Discussion started by: jacobs.smith
19 Replies
5. UNIX for Dummies Questions & Answers
hi
i have used comm -13 <(sort 1.txt) <(sort 2.txt) option to get the unique lines that are present in file 2 but not in file 1. but some how i am getting the entire file 2. i would expect few but not all uncommon lines fro my dat. is there anything wrong with the way i used the command?
my... (1 Reply)
Discussion started by: anurupa777
1 Replies
6. Shell Programming and Scripting
hi
my problem is little complicated one. i have 2 files which appear like this
file 1
abbsss:aa:22:34:as akl abc 1234
mkilll:as:ss:23:qs asc abc 0987
mlopii:cd:wq:24:as asd abc 7866
file2
lkoaa:as:24:32:sa alk abc 3245
lkmo:as:34:43:qs qsa abc 0987
kloia:ds:45:56:sa acq abc 7805
i... (5 Replies)
Discussion started by: anurupa777
5 Replies
7. Shell Programming and Scripting
Hello everyone,
Maybe somebody could help me with an awk script.
I have this input (field separator is comma ","):
547894982,M|N|J,U|Q|P,98,101,0,1,1
234900027,M|N|J,U|Q|P,98,101,0,1,1
234900023,M|N|J,U|Q|P,98,54,3,1,1
234900028,M|H|J,S|Q|P,98,101,0,1,1
234900030,M|N|J,U|F|P,98,101,0,1,1... (2 Replies)
Discussion started by: Ophiuchus
2 Replies
8. Shell Programming and Scripting
file 1
Sun Mar 17 00:01:33 2013 submit , Name="1234"
Sun Mar 17 00:01:33 2013 submit , Name="1344"
Sun Mar 17 00:01:33 2013 submit , Name="1124"
..
..
..
..
Sun Mar 17 00:01:33 2013 submit , Name="8901"
file 2
Sun Mar 17 00:02:47 2013 1234 execute SUCCEEDED
Sun Mar 17... (24 Replies)
Discussion started by: aravindj80
24 Replies
9. UNIX for Dummies Questions & Answers
I would like to print unique lines without sort or unique. Unfortunately the server I am working on does not have sort or unique. I have not been able to contact the administrator of the server to ask him to add it for several weeks. (7 Replies)
Discussion started by: cokedude
7 Replies
10. UNIX for Beginners Questions & Answers
I have a directory of files, I can show the number of lines in each file and order them from lowest to highest with:
wc -l *|sort
15263 Image.txt
16401 reference.txt
40459 richtexteditor.txt
How can I also print the number of unique lines in each file?
15263 1401 Image.txt
16401... (15 Replies)
Discussion started by: spacegoose
15 Replies
LEARN ABOUT DEBIAN
bio::popgen::genotype
Bio::PopGen::Genotype(3pm) User Contributed Perl Documentation Bio::PopGen::Genotype(3pm)
NAME
Bio::PopGen::Genotype - An implementation of GenotypeI which is just an allele container
SYNOPSIS
use Bio::PopGen::Genotype;
my $genotype = Bio::PopGen::Genotype->new(-marker_name => $name,
-individual_id => $indid,
-alleles => @alleles);
DESCRIPTION
This object will contain alleles for a given marker for a given individual.
The class variable BlankAlleles (accessible through $Bio::PopGen::Genotype::BlankAlleles = 'somepattern') can be set to a regexp pattern
for identifying blank alleles which should no be counted (they are effectively missing data). By default it set to match white space, '-',
'N' or 'n', and '?' as blank alleles which are skipped.
FEEDBACK
Mailing Lists
User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to the
Bioperl mailing list. Your participation is much appreciated.
bioperl-l@bioperl.org - General discussion
http://bioperl.org/wiki/Mailing_lists - About the mailing lists
Support
Please direct usage questions or support issues to the mailing list:
bioperl-l@bioperl.org
rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address
it. Please include a thorough description of the problem with code and data examples if at all possible.
Reporting Bugs
Report bugs to the Bioperl bug tracking system to help us keep track of the bugs and their resolution. Bug reports can be submitted via the
web:
https://redmine.open-bio.org/projects/bioperl/
AUTHOR - Jason Stajich
Email jason-at-bioperl.org
CONTRIBUTORS
Matthew Hahn, matthew.hahn-at-duke.edu
APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
new
Title : new
Usage : my $obj = Bio::PopGen::Genotype->new();
Function: Builds a new Bio::PopGen::Genotype object
Returns : an instance of Bio::PopGen::Genotype
Args : -marker_name => string representing name of the marker
-individual_id => string representing individual id (optional)
-alleles => arrayref with each item in the array being an allele
marker_name
Title : marker_name
Usage : my $name = $genotype->marker_name();
Function: Get the marker name for a genotype result
Returns : string
Args : [optional] marker name value to store
individual_id
Title : individual_id
Usage : my $indid = $genotype->individual_id();
Function: Gets the individual id associated with a genotype
This is effectively a back reference since we will typically
associate a genotype with an individual with an
individual HAS-A genotype relationship.
Returns : unique id string for an individual
Args : none
get_Alleles
Title : get_Alleles
Usage : my @alleles = $genotype->get_Alleles();
Function: Get the alleles for a given marker and individual
Returns : array of alleles (strings in this implementation)
Args : $showblank - boolean flag to indicate return ALL alleles not
skipping the coded EMPTY alleles
Note : Uses the class variable $BlankAlleles to test if alleles
should be skipped or not.
add_Allele
Title : add_Allele
Usage : $genotype->add_Allele(@alleles);
Function: Add alleles to the genotype, at this point there is no
verification to insure that haploid individuals only have 1
allele or that diploids only have 2 - we assume that is
done by the user creating these objects
Returns : count of the number of alleles in genotype
Args : Array of alleles to store
reset_Alleles
Title : reset_Alleles
Usage : $genotype->reset_Alleles;
Function: Resets the stored alleles so the list is empty
Returns : None
Args : None
perl v5.14.2 2012-03-02 Bio::PopGen::Genotype(3pm)