Search Results

How come your ch1 400 line is "flanked" by the ch2 567 line? If the chromosome doesn't need to match, why then wouldn't ch3 456 flank the 400 line?

---------- Post updated at 21:30 ----------...

How about awk '
FNR==NR {if ($3 > 0.7) EXP[$1] = EXP[$1] "," $2
next
}
{split ($1 EXP[$1], U, ",")
split ($2 EXP[$2], V,...

Can't tell - even though scores below 0.7 won't be collected into awk arrays, 345 E6 records may be too many. And, I'm afraid you can't avoid three loops if you want to permute three columns.

Try awk '
NR==FNR {a[$1]
next
}
FNR==1
{split ($2, T, "/")}
T[1] in a
' file1 file2
col1 col2 col3
1 a/e aa
2 b/e ...

Why don't you tell us WHAT's going wrong in lieu of leaving it up to us to guess?

For the data file, NR can't be 1. User FNR instead.

---------- Post updated at 22:01 ---------- Previous update...

Like so:awk ' {LN[$1]; HD[$2]; MX[$1,$2]++}
END { printf "%10s", ""; for (i in HD) printf "%10s", i; print "";
for (j in LN) {printf "%10s",j; ...

Brilliant! I seem to have been a bit worn out last night...

Could I suggest these changes in red for m calculation:

for (g in GN) {printf "%-40s", g
for (m in M) {
...

Try awk 'FNR==1 {FILE++}

FILE==1 {CHR[$1]=$2
POS[$1]=$3
next
}

FILE==2 ...

You're not printing anything in the first awk so there's nothing to pipe to the second. Even if there were, the second awk just reads file1 but not stdin from the pipe (you can use "-" to supply...

If you save the following awk script in a file named Serials.awk:
BEGIN { FS = OFS = ","
}
FNR == NR {
n[$1,++nc[$1]] = $2
next
}
FILENAME == "Serials.csv" {
nextfile
}
FNR == 1 {...

Close. To make it really read n columns, tryawk 'NR==1 {for (i=4; i<=NF; i++) READ[i]=$i
COLS=NF
next
}
...

There's a well tested solution for converting single column data to tables/matrices published several times in these fora; some upfront conditioning of your data and we can use that ideally:awk ...

Could this help you ?
script name - Forum1.pl
#!/usr/bin/perl

my %hashTree;
my @array;

while (<>) {
chomp;
@array=split;
$hashTree{$array[0]}{$array[3]}=[$array[1] ,$array[2]];
}...

Is you input file (other than the header line) always in sorted order with the Group (field 1) being the primary sort key and the Serial Number (field 4) being the secondary sort key (as in your...

This is probably no better than what you already have, but what the heck!

#!/bin/bash

# sort -k 1,1 -k 4,4 -o infile impute.data
infil=infile
outfil=outfile
wrkfil=/tmp/impute
touch...

Another one, perhaps even more clumsy, but works with any awk (but the Solaris /usr/bin/awk -> oawk)
awk '
BEGIN {OFS="\t"}
function prlast(){
f=""
if (T1=="") {
t=substr(T2,1,1)
...

Well, this again is far from elegant, but tryawk 'function prlast( ) {
for (y in Y)
if (T[y]) {sub (T[y], "", C[LAST])
...

Strange, complex conditions ... there might be a more elegant solution, but tryawk 'NR==1 {print $0, "similar to"}
NR==FNR {if ($2 !~/\*/) {L[$1,$2]=$3
...