Move all output formatting from analyses() to table_out()

[gerasy1987]

Currently treat_only, stars and round_digits are arguments for usefulr::analyses() which does not make sense, since the function is mainly for estimation. Consider moving this functionality and all output formatting to usefulr::table_fun()

This most likely will mean that the following chunk should be moved to usefulr::table_fun()

out <-
    dplyr::mutate(estout,
                  printout =
                    ifelse(is.nan(estimate), "-- [--]",
                           paste0(fround(estimate, digits = round_digits),
                                  ifelse(stars, add_stars(p.value), ""),
                                  " [", fround(std.error, digits = round_digits), "]")),
                  estimate =
                    round(estimate, digits = 3), std.error = round(std.error,
                                                                   digits = round_digits),
                  p.value = round(p.value, digits = round_digits))

out <- dplyr::select(.data = out,
                                  term, estimate, std.error,
                                  printout, p.value)