03-17-2016
I do not understand what output you are trying to produce.
For each input gene line do you want:
- one output line with the gene number and 10 averages from all 13 samples,
- thirteen output lines with the gene number and the 10 averages from one sample on each line, or
- one output line with the gene number and 130 averages where each set of 10 averages comes from one sample?
Can you show us the exact output you're hoping to produce from the data provided in your sample input for genes 1 and 2?
Was the data for gene 3 in your sample truncated, or will some inputs have missing fields that should be treated as zero values when computing the averages?
This User Gave Thanks to Don Cragun For This Post:
10 More Discussions You Might Find Interesting
1. Shell Programming and Scripting
Hi All,
I need the modification for the below mentioned code (found in one more post https://www.unix.com/shell-programming-scripting/27161-script-generate-average-values.html) to find the average values for all the columns(but for a specific rows) and print the averages side by side.
I have... (4 Replies)
Discussion started by: ks_reddy
4 Replies
2. Shell Programming and Scripting
I have a lot of input files that have the following form:
Sample Cq Sample Cq Sample Cq Sample Cq Sample Cq
1WBIN 23.45 1WBIN 23.45 1CVSIN 23.96 1CVSIN 23.14 S1 31.37
1WBIN 23.53 1WBIN 23.53 1CVSIN 23.81 1CVSIN 23.24 S1 31.49
1WBIN 24.55 1WBIN 24.55 1CVSIN 23.86 1CVSIN 23.24 S1 31.74 ... (3 Replies)
Discussion started by: isildur1234
3 Replies
3. Shell Programming and Scripting
Hi,
I have a file which looks like this:
FID IID MISS_PHENO N_MISS N_GENO F_MISS
12AB43131 12AB43131 N 17774 906341 0.01961
65HJ87451 65HJ87451 N 10149 906341 0.0112
43JJ21345 43JJ21345 N 2826 906341 0.003118I would... (11 Replies)
Discussion started by: kayakj
11 Replies
4. Shell Programming and Scripting
Hi to all,
I have two files. File1 has no header, two columns:
sample1 A
sample2 B
sample3 B
sample4 C
sample5 A
sample6 D
sample7 D
File2 has a header, except for the first 3 columns (chr,start,end). "sample1" is the header for the 4th ,5th ,6th columns, "sample2" is the header... (4 Replies)
Discussion started by: aec
4 Replies
5. Shell Programming and Scripting
Hello,
I have some tab delimited text files with a three header rows. The headers look like, (sorry the tabs look so messy).
index group Name input input input input input input input input input input input... (9 Replies)
Discussion started by: LMHmedchem
9 Replies
6. Shell Programming and Scripting
Hi experts,
I want to group by average, for multiple columns starting column $7 until NF,
group by ($1-$5), please help
For just 7th column, I can do
awk '
NR>1{
arr += $7
count += 1
}
END{
for (a in arr) {
print a, arr/count
... (10 Replies)
Discussion started by: ritakadm
10 Replies
7. UNIX for Beginners Questions & Answers
I have this code below that only prints out certain columns from the first two rows (doesn't affect rows 3 and beyond). How can I do the same on a partial header pattern “G_TP” instead of having to know specific column numbers (e.g. 374-479)? I've tried many other commands within this pipe with no... (4 Replies)
Discussion started by: aachave1
4 Replies
8. Shell Programming and Scripting
Hello,
I have to fish out some specific columns from a file based on the header value. I have the list of columns I need in a different file. I thought I could read in the list of headers I need,
# file with header names of required columns in required order
headers_file=$2
# read contents... (11 Replies)
Discussion started by: LMHmedchem
11 Replies
9. Shell Programming and Scripting
hello, I have three files in the following order
==> File1 <==
1 20977000 20977000 A C 1.00 0,15 15 45
1 115829313 115829313 G A 0.500 6,7 13 99
==> File2 <==
1 20977000 20977000 A C 1.00 0,13 13 39
1 115829313 ... (5 Replies)
Discussion started by: nans
5 Replies
10. UNIX for Beginners Questions & Answers
I have files that have the following columns
chr pos ref alt sample 1 sample 2 sample 3
chr2 179644035 G A 1,107 0,1 58,67
chr7 151945167 G T 142,101 100,200 500,700
chr13 31789169 CTT CT,C 6,37,8 0,0,0 15,46,89
chr22 ... (3 Replies)
Discussion started by: nans
3 Replies
LEARN ABOUT DEBIAN
bio::seqfeature::gene::genestructurei
Bio::SeqFeature::Gene::GeneStructureI(3pm) User Contributed Perl Documentation Bio::SeqFeature::Gene::GeneStructureI(3pm)
NAME
Bio::SeqFeature::Gene::GeneStructureI - A feature representing an arbitrarily
complex structure of a gene
SYNOPSIS
#documentation needed
DESCRIPTION
A feature representing a gene structure.
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 one
of the Bioperl mailing lists. 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 the bugs and their resolution. Bug reports can be submitted via the
web:
https://redmine.open-bio.org/projects/bioperl/
AUTHOR - Hilmar Lapp
Email hlapp@gmx.net
APPENDIX
The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _
transcripts
Title : transcripts()
Usage : @transcripts = $gene->transcripts();
Function: Get the transcript features/sites of this gene structure.
See Bio::SeqFeature::Gene::TranscriptI for properties of the
returned objects.
Returns : An array of Bio::SeqFeature::Gene::TranscriptI implementing objects
representing the promoter regions or sites.
Args :
promoters
Title : promoters()
Usage : @prom_sites = $gene->promoters();
Function: Get the promoter features/sites of this gene structure.
Note that OO-modeling of regulatory elements is not stable yet.
This means that this method might change or even disappear in a
future release. Be aware of this if you use it.
Returns : An array of Bio::SeqFeatureI implementing objects representing the
promoter regions or sites.
Args :
exons
Title : exons()
Usage : @exons = $gene->exons();
@inital = $gene->exons('Initial');
Function: Get all exon features or all exons of specified type of this gene
structure.
Refer to the documentation of the class that produced this gene
structure object for information about the possible types.
See Bio::SeqFeature::Gene::ExonI for properties of the
returned objects.
Returns : An array of Bio::SeqFeature::Gene::ExonI implementing objects
representing the exon regions.
Args : An optional string specifying the type of the exon.
introns
Title : introns()
Usage : @introns = $gene->introns();
Function: Get all introns of this gene structure.
Returns : An array of Bio::SeqFeatureI implementing objects representing the
introns.
Args :
poly_A_sites
Title : poly_A_sites()
Usage : @polyAsites = $gene->poly_A_sites();
Function: Get the poly-adenylation features/sites of this gene structure.
Returns : An array of Bio::SeqFeatureI implementing objects representing the
poly-adenylation regions or sites.
Args :
utrs
Title : utrs()
Usage : @utr_sites = $gene->utrs();
Function: Get the UTR features/sites of this gene structure.
See Bio::SeqFeature::Gene::ExonI for properties of the
returned objects.
Returns : An array of Bio::SeqFeature::Gene::ExonI implementing objects
representing the UTR regions or sites.
Args :
perl v5.14.2 2012-03-02 Bio::SeqFeature::Gene::GeneStructureI(3pm)