I am using Cuffcompare to identify novel and know transcripts from my cufflinks assemblies. In the paper from Nature protocols: http://www.nature.com/nprot/journal/v7/n3/full/nprot.2012.016.html there is a command in awk like this for: "printing a simple table for each assembly that lists how many transcripts in each assembly are complete matches to known transcripts".

$ for I in 'find . -name *.tmap'; do echo $i; awk 'NR  >  1 { s[$3] +  +  } END { \
                  for (j in s) { print j, s[j] }} ' $i; done;

Then, I get an error:

$ for I in 'find . -name *.tmap'; do echo $i; awk 'NR  >  1 { s[$3] +  +  } END {
for (j in s) { print j, s[j] }} ' $i; done;

find . -name *.tmap
awk: NR  >  1 { s[$3] +  +  } END {
awk:                        ^ syntax error

Unfortunately I don't know how to use awk and can't for my own know what can be causing this error. Maybe a different version of awk? Can someone help me?

$ awk --version
GNU Awk 3.1.7
Copyright (C) 1989, 1991-2009 Free Software Foundation.

Thank you for your support.