I made the change in the code and attached the output, but it looks like the value is column 1 is off and that the values are shifted by a column. Thank you .
Code:
#!/bin/perl
use strict;
# Accept the input and output files as parameters
my $input_file = $ARGV[0];
my $output_file = $ARGV[1];
# Set the header columns to be added to the left
# and to the right of the header in the input file
my @left = (
"Index",
"Chromosome Position",
"Gene",
"Inheritance",
"RNA Accession",
"Chr",
"Coverage",
"Score",
"A(#F,#R)",
"C(#F,#R)",
"G(#F,#R)",
"T(#F,#R)",
"Ins(#F,#R)",
"Del(#F,#R)",
"SNP db_xref",
"Mutation Call",
"Mutant Allele Frequency",
"Amino Acid Change"
);
my @right = (
"HP",
"SPLICE",
"Pseudogene",
"Classification",
"HGMD",
"Disease",
"Sanger",
"References"
);
# open the input file, read the header line and sandwich it
# between @left and @right arrays
my $final_header;
open (FH, "<", $input_file) or die "Can't open $input_file: $!";
chomp(my $hdr=<FH>);
$final_header = sprintf("%s\t%s\t%s\n", join("\t", @left), $hdr, join("\t",@right));
# final header is set, print it to the output file
open (OF, ">", $output_file) or die "Can't open $output_file: $!";
print OF "$final_header";
# close (FH) or die "Can't close $output_file: $!";
my @colsleft = map "Null",(0..$#left);
my @colsright = map "Null",(0..$#right);
while(<FH>) { # puts row of input file into $_
chomp;
print OF join("\t",$.,@colsleft,$_,@colsright),"\n"; # row data is set, print it to the output file
}
Greetings.
I'm having a time of it with this file. I'm trying to do a script that will take two command line inputs, string1 and string2 and use sed to change the text over files in the current directory. This is what I have so far. It appears to work a little, it does create the... (3 Replies)
I can not for the life of me figure out how to iterate over this array
{
'name1' => {
'a' => .
'b' =>
},
'name2' => {
'a' => .
'b' =>
}
}
I want a for loop to iterate through the first element (name1 and name2 in my example) but I can't figure it out. Help... (2 Replies)
Hi,
I have a directory call Test, which contains files "a", b", "c", etc. I want to go through all of the files inside Test and remove any empty file. How would I do that with shell csh?
So far I got...
#!/bin/csh
if (($#argv == 0) || ($#argv > 1)) then
echo "no argument or too... (2 Replies)
Hi,
I am writing a comparator script, which comapre two txt files(column by column)
below are the precondition of this comparator
1)columns of file are not seperated
Ex.
file1.txt
8888812341181892
1243548895685687
8945896789897789
1111111111111111
file2.txt
9578956789567897... (2 Replies)
Hi
i need to find x in the following equation such that it satisfies this condition:
y/x-ln(x)-1.24=0
how can i write a script to iterate to give random x to satisfy this equation.
y is different each time too. any help with awk/shell script will be awesome!
thanks (1 Reply)
I'm trying to send the file list as parameter to another job and execute it.
But the loop doesn't work, the inner job is running only once and not twice as expected
for filelist in $(ls -rt *.txt | tail -2)
do
echo $filelist
export filelist
cmd="$Program -config $configfile -autoexec... (11 Replies)
Have two 3 files which has list of servers,users and location and base url which is common on every server
A = server1 server2 server3
B = user1 user2 user3
C = dom1 dom2 dom3
baseurl=/opt/SP/
and what i have to achieve is below via ssh from REMOTE SERVER
for it's first iteration it... (7 Replies)
I would like to iterate over `dirs`in a script, but the script will never show more than one (current) folder
#! /bin/bash
for i in `dirs`
do
echo ${i}
done
echo ++++++++++++++++++
for i in $( dirs -p )
do
echo ${i}
done
echo ------------------
dirscontent=`dirs`
echo... (5 Replies)
I have a job that produces a file of barcodes that gets added to every time the job runs
I want to check the list to see if the barcode is already in the list and report it out if it is. (3 Replies)
Discussion started by: worky
3 Replies
LEARN ABOUT SUSE
pdl::reduce
Reduce(3) User Contributed Perl Documentation Reduce(3)NAME
PDL::Reduce -- a "reduce" function for PDL
DESCRIPTION
Many languages have a "reduce" function used to reduce the rank of an N-D array by one. It works by applying a selected operation along a
specified dimension. This module implements such a function for PDL by providing a simplified interface to the existing projection
functions (e.g. "sumover", "maximum", "average", etc).
SYNOPSIS
use PDL::Reduce;
$a = sequence 5,5;
# reduce by adding all
# elements along 2nd dimension
$b = $a->reduce('add',1);
@ops = $a->canreduce; # return a list of all allowed operations
FUNCTIONS
reduce
reduce dimension of piddle by one by applying an operation along the specified dimension
$a = sequence 5,5;
# reduce by adding all
# elements along 2nd dimension
$b = $a->reduce('add',1);
$b = $a->reduce('plus',1);
$b = $a->reduce('+',1); # three ways to do the same thing
[ As an aside: if you are familiar with threading you will see that this is actually the same as
$b = $a->mv(1,0)->sumover
]
NOTE: You should quote the name of the operation (1st arg) that you want "reduce" to perform. This is important since some of the names are
identical to the names of the actual PDL functions which might be imported into your namespace. And you definitely want a string as
argument, not a function invocation! For example, this will probably fail:
$b = $a->reduce(avg,1); # gives an error from invocation of 'avg'
Rather use
$b = $a->reduce('avg',1);
"reduce" provides a simple and unified interface to the projection functions and makes people coming from other data/array languages
hopefully feel more at home.
$result = $pdl->reduce($operation [,@dims]);
"reduce" applies the named operation along the specified dimension(s) reducing the input piddle dimension by as many dimensions as supplied
as arguments. If the dimension(s) argument is omitted the operation is applied along the first dimension. To get a list of valid operations
see canreduce.
NOTE - new power user feature: you can now supply a code reference as operation to reduce with.
# reduce by summing over dims 0 and 2
$result = $pdl->reduce(&sumover, 0, 2);
It is your responsibility to ensure that this is indeed a PDL projection operation that turns vectors into scalars! You have been warned.
canreduce
return list of valid named "reduce" operations Some common operations can be accessed using a number of names, e.g. '+', "add" and "plus"
all sum the elements along the chosen dimension.
@ops = PDL->canreduce;
This list is useful if you want to make sure which operations can be used with "reduce".
AUTHOR
Copyright (C) 2000 Christian Soeller (c.soeller@auckland.ac.nz). All rights reserved. There is no warranty. You are allowed to redistribute
this software / documentation under certain conditions. For details, see the file COPYING in the PDL distribution. If this file is
separated from the PDL distribution, the copyright notice should be included in the file.
perl v5.12.1 2009-10-17 Reduce(3)
04-15-2015
.