{total+=$1; totalp2+=($1)*($1)} END{if (NR==0) {avg=0; std_err=0} else { avg=total/NR; std_dev=sqrt((totalp2-total*total/NR)/(NR-1)); std_err=std_dev/sqrt(NR) } printf("\t%s\t%s\n",avg,std_err) }