add standard error calculation, but it seems to not give the same results as paml, possibly because paml uses the plain kappa and omega and branch lengths instead of the log transformations of these values