inst/doc/tips.R

In desctable: Produce Descriptive and Comparative Tables Easily

## ---- echo = F, message = F, warning = F--------------------------------------
library(desctable)

## -----------------------------------------------------------------------------
labels <- c(mpg   = "Miles/(US) gallon",
            cyl   = "Number of cylinders",
            disp  = "Displacement (cu.in.)",
            hp    = "Gross horsepower",
            drat  = "Rear axle ratio",
            wt    = "Weight (1000 lbs)",
            qsec  = "1/4 mile time",
            vs    = "Engine",
            am    = "Transmission",
            gear  = "Number of forward gears",
            CARBURATOR = "Number of carburetors")

mtcars %>%
  desc_table(.labels = labels) %>%
  desc_output("DT")

## ---- eval = F----------------------------------------------------------------
#  stats_auto <- function(data) {
#    data %>%
#      lapply(is.numeric) %>%
#      unlist() %>%
#      any -> numeric
#  
#    data %>%
#      lapply(is.factor) %>%
#      unlist() %>%
#      any() -> fact
#  
#    stats <- list("Min"  = min,
#                  "Q1"   = ~quantile(., .25),
#                  "Med"  = stats::median,
#                  "Mean" = mean,
#                  "Q3"   = ~quantile(., .75),
#                  "Max"  = max,
#                  "sd"   = stats::sd,
#                  "IQR"  = IQR)
#  
#    if (fact & numeric)
#      c(list("N" = length,
#             "%" = percent),
#        stats)
#    else if (fact & !numeric)
#      list("N" = length,
#           "%" = percent)
#    else if (!fact & numeric)
#      stats
#  }

## -----------------------------------------------------------------------------
stats = list(N = length,
             Mean = mean,
             SD = sd)

mtcars %>%
  desc_table(!!!stats) %>%
  desc_output("DT")


## -----------------------------------------------------------------------------
stats2 = list(N = length,
              mean,
              sd)

mtcars %>%
  desc_table(!!!stats2) %>%
  desc_output("DT")


## -----------------------------------------------------------------------------
default_stats <- function(data)
{
  list(N = length,
       mean,
       sd)
}

## ---- eval = F----------------------------------------------------------------
#  tests_auto <- function(var, grp) {
#    grp <- factor(grp)
#  
#    if (nlevels(grp) < 2)
#      ~no.test
#    else if (is.factor(var)) {
#      if (tryCatch(is.numeric(fisher.test(var ~ grp)$p.value), error = function(e) F))
#        ~fisher.test
#      else
#        ~chisq.test
#    } else if (nlevels(grp) == 2)
#      ~wilcox.test
#    else
#      ~kruskal.test
#  }

## -----------------------------------------------------------------------------
mtcars %>%
  group_by(am) %>%
  desc_table(mean, sd) %>%
  desc_tests(.default = ~t.test) %>%
  desc_output("DT")

## ---- warning = F-------------------------------------------------------------
mtcars %>%
  group_by(am) %>%
  desc_table(mean, sd, median, IQR) %>%
  desc_tests(.default = ~t.test, carb = ~wilcox.test) %>%
  desc_output("DT")

## -----------------------------------------------------------------------------
iris %>%
  group_by(Species) %>%
  desc_table(mean, sd) %>%
  desc_tests() %>%
  desc_output("DT", digits = 10)

## -----------------------------------------------------------------------------
iris %>%
  group_by(Species) %>%
  desc_table(mean, sd) %>%
  desc_output("DT", filter = "top")