Utiliser awk pour calculer la somme de plusieurs colonnes

J'ai une liste avec les 8 colonnes, dont les 6 premiers sont les mêmes.

6   99999715    99999771    NM_001013399    0   -   23  0.0714286
6   99999715    99999771    NM_001013399    0   -   24  0.0178571
6   99999715    99999771    NM_001013399    0   -   25  0.1250000

J'ai besoin de calculer la moyenne de la colonne 7 et column8, ainsi que 7 $*$8, et obtenir le format comme:

6   99999715    99999771    NM_001013399    0   -   ave($7) ave($8) ave($7*$8)

Comment dois-je faire?
thx

OriginalL'auteur LookIntoEast | 2012-05-05

  1. 7

    N'ai pas essayé, mais il doit être juste:

    {sum7+=$7; sum8+=$8; mul+=$7*$8} END {print sum7/NR,sum8/NR,mul/NR}

    En réponse à la demande populaire, je vais ajouter le printf.

    {sum7+=$7; sum8+=$8; mul+=$7*$8}
END {printf "%s %4i %10.7f %10.7f\n", substr($0,0,49),sum7/NR,sum8/NR,mul/NR}
    C'est vraiment une réponse utile pour montrer certains awk calculs, mais ne donne pas le format de sortie demandée.

    OriginalL'auteur stark

  2. 6
    awk '
{
  if( common == "" ) { 
    fn=1   # field number
    cn=1   # column number
    tmp=$0
    f=0 
    while( match( tmp, / *|$/) && f<=NF ) 
    {  f+=1
       cnA[fn]=cn           # column number of start of field fn
       cnZ[fn]=cn+RSTART-1  # column number of   end of field fn
       ++fn
       cn+=RSTART+RLENGTH-1
       tmp=substr( tmp, RSTART+RLENGTH )
    }
    common = substr($0,1,cnA[7]-1)
    dlim78 = substr($0,cnZ[7], cnZ[7]-cnA[7])
  }  
  print $0
  (f7+=$7)
  (f8+=$8)
}
END {
  p7=".0" # decimal places ($7)
  p8=".7" # decimal places ($8)
  pP=".7" # decimal places ($7*$8) 
  printf( "%s%"p7"f%s%"p8"f%s%"pP"f\n" ,
           common, f7/NR, dlim78, f8/NR, dlim78,f7*f8/NR )
}
' <<'EOF'
6   99999715    99999771    NM_001013399    0   -   23  0.0714286
6   99999715    99999771    NM_001013399    0   -   25  0.1250000
EOF

    De sortie: 

    6   99999715    99999771    NM_001013399    0   -   23  0.0714286
6   99999715    99999771    NM_001013399    0   -   25  0.1250000
6   99999715    99999771    NM_001013399    0   -   24  0.0982143  4.7142864
    Pour whowever il était qui ont marqué cette réponse vers le bas. J'aimerais leur montrer aucun d'autre moyen de maintian le champ de séparation, à l'aide d'un script awk... Qui semble être un point essentiel de la question.
    ... et c'est encore commenté. Beau travail 🙂

    OriginalL'auteur Peter.O