Search Results

cat pattern.txt
chr1 1000 2000 chr1 10000 11000
chr1 1000 2000 chr1 11000 12000
chr1 1000 2000 chr1 12000 13000
cat DB_file.txt
chr1 1000 ...

yes, it would. Try to pass the file name correctly and handle it in code

try also:
awk '
NR>1{l[$1]=$1; c[$1]++;
for (i=2; i<=NF; i++) a[$1,i]+=$i;
}
END {
for (g in l) {
printf g " ";
for (i=2; i<=NF; i++) printf ("%.1f ", (a[g,i]/c[g]));
...

l[$1]=$1 is a useless value; l[$1] alone defines the key (no value).
Or store the length, i.e. allow an individual length for each gene type:
awk '
NR>1 {
L[$1]=NF; c[$1]++
for (i=2;...

This is by reading same file twice, you can also process this in END block

$ cat file
Gene Sample_1 Sample_2 Sample_3
gene_A 2 4 5
gene_B 1 2 3
gene_A 0 5 7
gene_B 4 5 6
gene_A 11 12 13...

awk '{if(!($1 $2 in a) && !($2 $1 in a)){print};{a[$1 $2]++}}'

Try: awk -F, '$3<=-0.5 || $3>=0.5' file



--
BTW your input does not correspond to your output. That should probably be:
A,B,0.6
C,D,-0.7

awk '$1 != $2 && $1=$1' OFS=, myFile

Output order will be same as that of edge_list.txt
Try :
$ awk 'FNR==NR{A[$1];next}($1 in A || $2 in A)' input.txt edge_list.txt
A B
J H
C A
G H
H A
K G

Hi, this is a fairly frequent question. If the file is sorted try:
awk 'p!=$1{if(s)print s; p=$1; s=$0; next} {s=s OFS $2} END{print s}' file

May try this...for given input this will work

$ cat file
LOC_Os01g01870 GO:0006139
LOC_Os01g01870 GO:0009058
LOC_Os01g02570 GO:0006464
LOC_Os01g02570 GO:0009987
LOC_Os01g02570 ...

$ grep -f nodes.txt edge.txt
LOC_Os11g37970 LOC_Os04g47140
LOC_Os03g19480 LOC_Os02g47060
LOC_Os03g49630 LOC_Os03g19480
LOC_Os04g50860 LOC_Os11g37970
Although you don't get the header...

Using awk:
awk 'NR==FNR{A[$1];next}$1 in A || $2 in A' node.txt edge.txt

Using awk:
awk '$3 >= 0.02' file

Another way with sort:
$ sed -n '2,$p' file | sort -n -k1,1 -k3,3 | sort -ur -k1,1
sed to suppress first line (column name).
Regards.

Put this into "script.awk":function abs(x){
return ((x < 0.0) ? -x : x)
}
NR>1{
if (max[$1]<abs($3)){
max[$1]=abs($3)
a[$1]=$2
b[$1]=$3
}
}
END{
for (i in a){
...

Using awk:
awk '{if(A[$1]<$2||!(A[$1])) A[$1]=$2}END{for(k in A) print k,A[k]}' file