Sigma becomes matrix when number of network is 1
moved extractTail and computeSummaryTable to function.r
rewrite log as logarithm