10-13-2011
Hi tyler,
I used your instructions as follows:
C:\Perl>perl -F"\t" -lane 'BEGIN {@hdr = qw(chromo pos refalleles.SampleA refAllele.SampleA refCount.SampleA refFreq.SampleA altAllele.SampleA altCount.SampleA altFreq.SampleA alleles.SampleB refAllele.SampleB refCount.SampleB refFreq.SampleB altAllele.SampleB altCount.SampleB altFreq.SampleB alleles.SampleC refAllele.SampleC refCount.SampleC refFreq.SampleC altAllele.SampleC altCount.SampleC altFreq.SampleC alleles.SampleD refAllele.SampleD refCount.SampleD refFreq.SampleD altAllele.SampleD altCount.SampleD altFreq.SampleD)} $x{join"\t",@F[0..2]}.="\t".join"\t",@F[3..9]; { print join " ", @hdr; foreach $k (sort keys %x) {print "$k\t$x{$k}" if $k !~ /chromo/}}' samplea sampleb samplec sampled
But I get the following error message:
Can't find string terminator "'" anywhere before EOF at -e line 1.
C:\Perl>
10 More Discussions You Might Find Interesting
1. Shell Programming and Scripting
Hi ,
I want to read the data from 9 tables in oracle DB into 9 different files in the same connection instance (session). I am able to get data from one table to one file with below code :
X=`sqlplus -s user/pwd@DB <<eof
select col1 from table1;
EXIT;
eof`
echo $X>myfile
Can anyone... (2 Replies)
Discussion started by: net
2 Replies
2. Programming
I'm pretty new to the database world and I've run into a mental block of sorts. I've been unable to find the answer anywhere. Here's my problem: I have several tables and everything is as normalized as possible (as I've been lead to understand normalization.) Normalization has lead to some... (1 Reply)
Discussion started by: flakblas
1 Replies
3. Shell Programming and Scripting
I've hunted and hunted but nothing seems to apply to what I need. Any help will be much appreciated!
My input file looks like (Unix):
marker,allele1,allele2
RS1002244,1,1
RS1002244,1,3
RS1002244,3,3
RS1003719,2,2
RS1003719,2,4
RS1003719,4,4
Most markers are listed 3 times but a few... (2 Replies)
Discussion started by: Peggy White
2 Replies
4. Shell Programming and Scripting
Hi..
File 1:
1 aa rep
1 dd rep
1 kk rep
2 bb sad
2 ss sad
3 ee dam
File 2
1 apple fruit
2 mango tree
3 lilly flower
output:
1 aaple fruit aa,dd,kk rep (7 Replies)
Discussion started by: empyrean
7 Replies
5. Web Development
Say I have two tables like below..
status
HId sName dName StartTime EndTime
1 E E 9:10 10:10
2 E F 9:15 10:15
3 G H 9:17 10:00
logic
Id devName capacity free Line
1 E 123 34 1
2 E 345 ... (3 Replies)
Discussion started by: ilan
3 Replies
6. Shell Programming and Scripting
Hi all,
I have a complex (beyond my biological expertise) problem at hand.
I need to merge multiple files into 1 big matrix. Please help me with some code.
Inp1
Ang_0 chr1 98 T A
Ang_0 chr1 352 G A
Ang_0 chr1 425 C T
Ang_0 chr2 ... (1 Reply)
Discussion started by: newbie83
1 Replies
7. Shell Programming and Scripting
multiple files to load into different tables,
I have a script show below, but this script loads data from txt file into a table,
but i have multiple input files(xyzload.txt,xyz1load.txt,xyz2load.txt......) in the unix folder ,
can we load these files in diff tables (table 1, table2... (0 Replies)
Discussion started by: nani1984
0 Replies
8. Shell Programming and Scripting
multiple files to load into different tables,
I have a script show below, but this script loads data from txt file into a table,
but i have multiple input files(xyzload.txt,xyz1load.txt,xyz2load.txt......) in the unix folder ,
can we load these files in diff tables (table 1, table2... (1 Reply)
Discussion started by: nani1984
1 Replies
9. UNIX for Dummies Questions & Answers
Hi,
Please excuse me , i have searched unix forum, i am unable to find what i expect ,
my query is , i have 2 files of same structure and having 1 similar field/column , i need to merge 2 tables/files based on the one matched field/column (that is field 1),
file 1:... (5 Replies)
Discussion started by: karthikram
5 Replies
10. UNIX for Beginners Questions & Answers
Hello All,
just wanted to export multiple tables from oracle sql using unix shell script to csv file and the below code is exporting only the first table.
Can you please suggest why? or any better idea?
export FILE="/abc/autom/file/geo_JOB.csv"
Export= `sqlplus -s dev01/password@dEV3... (16 Replies)
Discussion started by: Hope
16 Replies
LEARN ABOUT DEBIAN
bio::graphics::glyph::allele_tower
Bio::Graphics::Glyph::allele_tower(3pm) User Contributed Perl Documentation Bio::Graphics::Glyph::allele_tower(3pm)
NAME
Bio::Graphics::Glyph::allele_tower - The "allele_tower" glyph
SYNOPSIS
See <Bio::Graphics::Panel> and <Bio::Graphics::Glyph>.
DESCRIPTION
This glyph draws a letter for each allele found at a SNP position, one above the other (i.e. in a column). For example:
A
G
See also http://www.hapmap.org/cgi-perl/gbrowse/gbrowse 'genotyped SNPs' for an example.
The common options are available (except height which is calculated based on the number of alleles). In addition, if you give the glyph
the minor allele frequency (MAF) and indicate which is the minor allele, the glyph will display these differences.
GETTING THE ALLELES
To specify the alleles, create an "Alleles" attribute for the feature. There should be two such attributes. For example, for a T/G
polymorphism, the GFF load file should look like:
Chr3 . SNP 12345 12345 . . . SNP ABC123; Alleles T ; Alleles G
Alternatively, you can pass an "alleles" callback to the appropriate section of the config file. This option should return the two alleles
separated by a slash:
alleles = sub {
my $snp = shift;
my @d = $snp->get_tag_values('AllelePair');
return join "/",@d;
}
OPTIONS
. Glyph Colour
. Different colour for alleles on the reverse strand
. Print out the complement for alleles on the reverse strand
. Major allele shown in bold
. Horizontal histogram to show allele frequency
GLYPH COLOR
The glyph color can be configured to be different if the feature is on the plus or minus strand. Use fgcolor to define the glyph color for
the plus strand and bgcolor for the minus strand. For example:
fgcolor = blue
bgcolor = red
For this option to work, you must also set ref_strand to return the strand of the feature:
ref_strand = sub {shift->strand}
REVERSE STRAND ALLELES
If the alleles on the negative strand need to be the complement of what is listed in the GFF files, (e.g. A/G becomes T/C), set the
complement option to have value 1
complement = 1
For this option to work, you must also set ref_strand to return the strand of the feature:
ref_strand = sub {shift->strand}
MAJOR/MINOR ALLELE
Use the 'minor_allele' option to return the minor allele for the SNP. If you use this option, the major allele will appear in bold type.
ALLELE FREQUENCY HISTOGRAMS
Use the 'maf' option to return the minor allele frequency for the SNP. If you use this option, a horizontal histogram will be drawn next
to the alleles, to indicate their relative frequencies. e.g.
A______
C__
Note: The 'label' option must be set to 1 (i.e. on) and the 'minor_allele' option must return a valid allele for this to work.
BUGS
Please report them.
SEE ALSO
Bio::Graphics::Panel, Bio::Graphics::Glyph, Bio::Graphics::Glyph::arrow, Bio::Graphics::Glyph::cds, Bio::Graphics::Glyph::crossbox,
Bio::Graphics::Glyph::diamond, Bio::Graphics::Glyph::dna, Bio::Graphics::Glyph::dot, Bio::Graphics::Glyph::ellipse,
Bio::Graphics::Glyph::extending_arrow, Bio::Graphics::Glyph::generic, Bio::Graphics::Glyph::graded_segments,
Bio::Graphics::Glyph::heterogeneous_segments, Bio::Graphics::Glyph::line, Bio::Graphics::Glyph::pinsertion, Bio::Graphics::Glyph::primers,
Bio::Graphics::Glyph::rndrect, Bio::Graphics::Glyph::segments, Bio::Graphics::Glyph::ruler_arrow, Bio::Graphics::Glyph::toomany,
Bio::Graphics::Glyph::transcript, Bio::Graphics::Glyph::transcript2, Bio::Graphics::Glyph::translation, Bio::Graphics::Glyph::allele_tower,
Bio::DB::GFF, Bio::SeqI, Bio::SeqFeatureI, Bio::Das, GD
AUTHOR
Fiona Cunningham <cunningh@cshl.edu> in Lincoln Stein's lab <steinl@cshl.edu>.
Copyright (c) 2003 Cold Spring Harbor Laboratory
This package and its accompanying libraries is free software; you can redistribute it and/or modify it under the terms of the GPL (either
version 1, or at your option, any later version) or the Artistic License 2.0. Refer to LICENSE for the full license text. In addition,
please see DISCLAIMER.txt for disclaimers of warranty.
perl v5.14.2 2012-02-20 Bio::Graphics::Glyph::allele_tower(3pm)