# Examples

#### 2023-07-13

### Introduction

First, we apply labels on the well-known `mtcars` dataset:

    library(expss)
    data(mtcars)
    mtcars = apply_labels(mtcars,
                          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",
                          vs = c("V-engine" = 0,
                                 "Straight engine" = 1),
                          am = "Transmission",
                          am = c("Automatic" = 0,
                                 "Manual"=1),
                          gear = "Number of forward gears",
                          carb = "Number of carburetors"
    )

Table construction consists of at least of three functions chained with magrittr pipe operator: `%>%`. At first we need to specify variables for which statistics will be computed with `tab_cells`. Secondary, we calculate statistics with one of `tab_stat_*` functions. And last, we finalize table creation with `tab_pivot`: `dataset %>% tab_cells(variable) %>% tab_stat_cases() %>% tab_pivot()`. We can split our statistics by columns with `tab_cols` or by rows with `tab_rows`. After that we can sort table with `tab_sort_asc`, drop empty rows/columns with `drop_rc` and transpose with `tab_transpose`. Generally, resulting table is just a data.frame so we can use arbitrary operations on it. Statistic is always calculated on the last cells, column/row variables, weight, missing values and subgroup. To define new cell/column/row variables we can call appropriate function one more time. `tab_pivot` defines how we combine different statistics and where statistic labels will appear - inside/outside rows/columns.

### Simple column percent

    mtcars %>%
        tab_cells(cyl) %>%
        tab_cols(vs) %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        tab_caption("Simple column percent")

Simple column percent
| Engine
--- | ---
| V-engine | Straight engine
Number of cylinders
4 | 5.6 | 71.4
6 | 16.7 | 28.6
8 | 77.8 |
#Total cases | 18 | 14

### Split by columns and rows

    mtcars %>%
        tab_cells(cyl) %>%
        tab_cols(vs) %>%
        tab_rows(am) %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        tab_caption("Split by columns and rows")

Split by columns and rows
|  |  | Engine
--- | --- | --- | ---
|  |  |  | V-engine | Straight engine
Transmission
Automatic | Number of cylinders | 4 |  |  | 42.9
|  | 6 |  |  | 57.1
|  | 8 |  | 100.0 |
|  | #Total cases |  | 12 | 7
Manual | Number of cylinders | 4 |  | 16.7 | 100.0
|  | 6 |  | 50.0 |
|  | 8 |  | 33.3 |
|  | #Total cases |  | 6 | 7

### Multiple banners, table is sorted by total

    mtcars %>%
        tab_cells(cyl) %>%
        tab_cols(total(), vs, am) %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        tab_sort_desc() %>%
        tab_caption("Multiple banners, table is sorted by total")

Multiple banners, table is sorted by total
| #Total |  | Engine |  | Transmission
--- | --- | --- | --- | --- | ---
|  |  | V-engine | Straight engine |  | Automatic | Manual
Number of cylinders
8 | 43.8 |  | 77.8 |  |  | 63.2 | 15.4
4 | 34.4 |  | 5.6 | 71.4 |  | 15.8 | 61.5
6 | 21.9 |  | 16.7 | 28.6 |  | 21.1 | 23.1
#Total cases | 32 |  | 18 | 14 |  | 19 | 13

### Nested banners

    mtcars %>%
        tab_cells(cyl) %>%
        tab_cols(total(), vs %nest% am) %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        tab_caption("Nested banners")

Nested banners
| #Total |  | Engine
--- | --- | --- | ---
|  |  | V-engine |  | Straight engine
|  |  | Transmission |  | Transmission
|  |  | Automatic | Manual |  | Automatic | Manual
Number of cylinders
4 | 34.4 |  |  | 16.7 |  | 42.9 | 100
6 | 21.9 |  |  | 50.0 |  | 57.1 |
8 | 43.8 |  | 100 | 33.3 |  |  |
#Total cases | 32 |  | 12 | 6 |  | 7 | 7

### Multiple nested banners

    mtcars %>%
        tab_cells(carb) %>%
        tab_cols(total(), list(cyl, vs) %nest% am) %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        tab_caption("Multiple nested banners")

Multiple nested banners
| #Total |  | Number of cylinders |  | Engine
--- | --- | --- | --- | --- | ---
|  |  | 4 |  | 6 |  | 8 |  | V-engine |  | Straight engine
|  |  | Transmission |  | Transmission |  | Transmission |  | Transmission |  | Transmission
|  |  | Automatic | Manual |  | Automatic | Manual |  | Automatic | Manual |  | Automatic | Manual |  | Automatic | Manual
Number of carburetors
1 | 21.9 |  | 33.3 | 50 |  | 50 |  |  |  |  |  |  |  |  | 42.9 | 57.1
2 | 31.2 |  | 66.7 | 50 |  |  |  |  | 33.3 |  |  | 33.3 | 16.7 |  | 28.6 | 42.9
3 | 9.4 |  |  |  |  |  |  |  | 25.0 |  |  | 25.0 |  |  |  |
4 | 31.2 |  |  |  |  | 50 | 66.7 |  | 41.7 | 50 |  | 41.7 | 50.0 |  | 28.6 |
6 | 3.1 |  |  |  |  |  | 33.3 |  |  |  |  |  | 16.7 |  |  |
8 | 3.1 |  |  |  |  |  |  |  |  | 50 |  |  | 16.7 |  |  |
#Total cases | 32 |  | 3 | 8 |  | 4 | 3 |  | 12 | 2 |  | 12 | 6 |  | 7 | 7

### Multiple variable and multiple summary statistics

    mtcars %>%
        tab_cells(mpg, disp, hp, wt, qsec) %>%
        tab_cols(total(), am) %>%
        tab_stat_fun(Mean = w_mean, "Std. dev." = w_sd, "Valid N" = w_n) %>%
        tab_pivot() %>%
        tab_caption("Multiple variable and multiple summary statistics")

Multiple variable and multiple summary statistics
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic | Manual
Miles/(US) gallon
Mean | 20.1 |  | 17.1 | 24.4
Std. dev. | 6.0 |  | 3.8 | 6.2
Valid N | 32.0 |  | 19.0 | 13.0
Displacement (cu.in.)
Mean | 230.7 |  | 290.4 | 143.5
Std. dev. | 123.9 |  | 110.2 | 87.2
Valid N | 32.0 |  | 19.0 | 13.0
Gross horsepower
Mean | 146.7 |  | 160.3 | 126.8
Std. dev. | 68.6 |  | 53.9 | 84.1
Valid N | 32.0 |  | 19.0 | 13.0
Weight (1000 lbs)
Mean | 3.2 |  | 3.8 | 2.4
Std. dev. | 1.0 |  | 0.8 | 0.6
Valid N | 32.0 |  | 19.0 | 13.0
1/4 mile time
Mean | 17.8 |  | 18.2 | 17.4
Std. dev. | 1.8 |  | 1.8 | 1.8
Valid N | 32.0 |  | 19.0 | 13.0

### Multiple variable and multiple summary statistics - statistic lables in columns

    mtcars %>%
        tab_cells(mpg, disp, hp, wt, qsec) %>%
        tab_cols(total(), am) %>%
        tab_stat_fun(Mean = w_mean, "Valid N" = w_n, method = list) %>%
        tab_pivot() %>%
        tab_caption("Multiple variable and multiple summary statistics - statistic lables in columns")

Multiple variable and multiple summary statistics - statistic lables in columns
| #Total |  | Transmission
--- | --- | --- | ---
| Mean |  | Valid N |  | Automatic |  | Manual
|  |  |  |  | Mean | Valid N |  | Mean | Valid N
Miles/(US) gallon | 20.1 |  | 32 |  | 17.1 | 19 |  | 24.4 | 13
Displacement (cu.in.) | 230.7 |  | 32 |  | 290.4 | 19 |  | 143.5 | 13
Gross horsepower | 146.7 |  | 32 |  | 160.3 | 19 |  | 126.8 | 13
Weight (1000 lbs) | 3.2 |  | 32 |  | 3.8 | 19 |  | 2.4 | 13
1/4 mile time | 17.8 |  | 32 |  | 18.2 | 19 |  | 17.4 | 13

### Filter dataset and exclude empty columns

    mtcars %>%
        tab_subgroup(am == 0) %>%
        tab_cells(cyl) %>%
        tab_cols(total(), vs %nest% am) %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        drop_empty_columns() %>%
        tab_caption("Filter dataset and exclude empty columns")

Filter dataset and exclude empty columns
| #Total |  | Engine
--- | --- | --- | ---
|  |  | V-engine |  | Straight engine
|  |  | Transmission |  | Transmission
|  |  | Automatic |  | Automatic
Number of cylinders
4 | 15.8 |  |  |  | 42.9
6 | 21.1 |  |  |  | 57.1
8 | 63.2 |  | 100 |  |
#Total cases | 19 |  | 12 |  | 7

### Total at the top of the table

    mtcars %>%
        tab_cells(cyl) %>%
        tab_cols(total(), vs) %>%
        tab_rows(am) %>%
        tab_stat_cpct(total_row_position = "above",
                      total_label = c("number of cases", "row %"),
                      total_statistic = c("u_cases", "u_rpct")) %>%
        tab_pivot() %>%
        tab_caption("Total at the top of the table")

Total at the top of the table
|  |  | #Total |  | Engine
--- | --- | --- | --- | --- | ---
|  |  |  |  |  | V-engine | Straight engine
Transmission
Automatic | Number of cylinders | #number of cases |  | 19 |  | 12 | 7
|  | #row % |  | 100 |  | 63.2 | 36.8
|  | 4 |  | 15.8 |  |  | 42.9
|  | 6 |  | 21.1 |  |  | 57.1
|  | 8 |  | 63.2 |  | 100.0 |
Manual | Number of cylinders | #number of cases |  | 13 |  | 6 | 7
|  | #row % |  | 100 |  | 46.2 | 53.8
|  | 4 |  | 61.5 |  | 16.7 | 100.0
|  | 6 |  | 23.1 |  | 50.0 |
|  | 8 |  | 15.4 |  | 33.3 |

### Three different statistics in each cell - stat. labels in rows

    mtcars %>%
        tab_cells(am) %>%
        tab_cols(total(), vs) %>%
        tab_total_row_position("none") %>%
        tab_stat_cpct(label = "col %") %>%
        tab_stat_rpct(label = "row %") %>%
        tab_stat_tpct(label = "table %") %>%
        tab_pivot(stat_position = "inside_rows") %>%
        tab_caption("Three different statistics in each cell - stat. labels in rows")

Three different statistics in each cell - stat. labels in rows
|  |  | #Total |  | Engine
--- | --- | --- | --- | --- | ---
|  |  |  |  | V-engine | Straight engine
Transmission
Automatic | col % |  | 59.4 |  | 66.7 | 50.0
| row % |  | 100.0 |  | 63.2 | 36.8
| table % |  | 59.4 |  | 37.5 | 21.9
Manual | col % |  | 40.6 |  | 33.3 | 50.0
| row % |  | 100.0 |  | 46.2 | 53.8
| table % |  | 40.6 |  | 18.8 | 21.9

### Three different statistics in each cell - stat. labels in columns

    mtcars %>%
        tab_cells(am) %>%
        tab_cols(total(), vs) %>%
        tab_total_row_position("none") %>%
        tab_stat_cpct(label = "col %") %>%
        tab_stat_rpct(label = "row %") %>%
        tab_stat_tpct(label = "table %") %>%
        tab_pivot(stat_position = "inside_columns") %>%
        tab_caption("Three different statistics in each cell - stat. labels in columns")

Three different statistics in each cell - stat. labels in columns
| #Total |  | Engine
--- | --- | --- | ---
| col % |  | row % |  | table % |  | V-engine |  | Straight engine
|  |  |  |  |  |  | col % | row % | table % |  | col % | row % | table %
Transmission
Automatic | 59.4 |  | 100 |  | 59.4 |  | 66.7 | 63.2 | 37.5 |  | 50 | 36.8 | 21.9
Manual | 40.6 |  | 100 |  | 40.6 |  | 33.3 | 46.2 | 18.8 |  | 50 | 53.8 | 21.9

### Stacked statistics

    mtcars %>%
        tab_cells(cyl) %>%
        tab_cols(total(), am) %>%
        tab_stat_mean() %>%
        tab_stat_se() %>%
        tab_stat_valid_n() %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        tab_caption("Stacked statistics")

Stacked statistics
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic | Manual
Number of cylinders
Mean | 6.2 |  | 6.9 | 5.1
S. E. | 0.3 |  | 0.4 | 0.4
Valid N | 32.0 |  | 19.0 | 13.0
4 | 34.4 |  | 15.8 | 61.5
6 | 21.9 |  | 21.1 | 23.1
8 | 43.8 |  | 63.2 | 15.4
#Total cases | 32 |  | 19 | 13

### Stacked statistics with section headings

    mtcars %>%
        tab_cells(cyl) %>%
        tab_cols(total(), am) %>%
        tab_row_label("#Summary statistics") %>%
        tab_stat_mean() %>%
        tab_stat_se() %>%
        tab_stat_valid_n() %>%
        tab_row_label("#Column percent") %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        tab_caption("Stacked statistics with section headings")

Stacked statistics with section headings
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic | Manual
#Summary statistics
Number of cylinders
Mean | 6.2 |  | 6.9 | 5.1
S. E. | 0.3 |  | 0.4 | 0.4
Valid N | 32.0 |  | 19.0 | 13.0
#Column percent
Number of cylinders
4 | 34.4 |  | 15.8 | 61.5
6 | 21.9 |  | 21.1 | 23.1
8 | 43.8 |  | 63.2 | 15.4
#Total cases | 32 |  | 19 | 13

### Stacked statistics - different statistics for different variables

    mtcars %>%
        tab_cols(total(), am) %>%
        tab_cells(mpg, hp, qsec) %>%
        tab_stat_mean() %>%
        tab_cells(cyl, carb) %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        tab_caption("Stacked statistics - different statistics for different variables")

Stacked statistics - different statistics for different variables
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic | Manual
Miles/(US) gallon
Mean | 20.1 |  | 17.1 | 24.4
Gross horsepower
Mean | 146.7 |  | 160.3 | 126.8
1/4 mile time
Mean | 17.8 |  | 18.2 | 17.4
Number of cylinders
4 | 34.4 |  | 15.8 | 61.5
6 | 21.9 |  | 21.1 | 23.1
8 | 43.8 |  | 63.2 | 15.4
#Total cases | 32 |  | 19 | 13
Number of carburetors
1 | 21.9 |  | 15.8 | 30.8
2 | 31.2 |  | 31.6 | 30.8
3 | 9.4 |  | 15.8 |
4 | 31.2 |  | 36.8 | 23.1
6 | 3.1 |  |  | 7.7
8 | 3.1 |  |  | 7.7
#Total cases | 32 |  | 19 | 13

### Linear regression by groups

    mtcars %>%
        tab_cells(sheet(mpg, disp, hp, wt, qsec)) %>%
        tab_cols(total(), am) %>%
        tab_stat_fun_df(
            function(x){
                frm = reformulate(".", response = as.name(names(x)[1]))
                model = lm(frm, data = x)
                sheet('Coef.' = coef(model),
                      confint(model)
                )
            }
        ) %>%
        tab_pivot() %>%
        tab_caption("Linear regression by groups")

Linear regression by groups
| #Total |  | Transmission
--- | --- | --- | ---
| Coef. |  | 2.5 % |  | 97.5 % |  | Automatic |  | Manual
|  |  |  |  |  |  | Coef. | 2.5 % | 97.5 % |  | Coef. | 2.5 % | 97.5 %
(Intercept) | 27.3 |  | 9.6 |  | 45.1 |  | 21.8 | -1.9 | 45.5 |  | 13.3 | -21.9 | 48.4
`Displacement (cu.in.)` | 0.0 |  | 0.0 |  | 0.0 |  | 0.0 | 0.0 | 0.0 |  | 0.0 | -0.1 | 0.1
`Gross horsepower` | 0.0 |  | -0.1 |  | 0.0 |  | 0.0 | -0.1 | 0.0 |  | 0.0 | 0.0 | 0.1
`Weight (1000 lbs)` | -4.6 |  | -7.2 |  | -2.0 |  | -2.3 | -5.0 | 0.4 |  | -7.7 | -12.5 | -2.9
`1/4 mile time` | 0.5 |  | -0.4 |  | 1.5 |  | 0.4 | -0.7 | 1.6 |  | 1.6 | -0.2 | 3.4

### Subtotals

    mtcars %>%
        tab_cells(mpg) %>%
        tab_cols(total(), vs) %>%
        tab_rows(subtotal(cyl, 1:2, 3:4, "5 and more" = 5 %thru% hi)) %>%
        tab_stat_mean() %>%
        tab_pivot() %>%
        tab_caption("Subtotals in rows")

Subtotals in rows
|  |  | #Total |  | Engine
--- | --- | --- | --- | --- | ---
|  |  |  |  |  | V-engine | Straight engine
Number of cylinders
1 | Miles/(US) gallon | Mean |  |  |  |  |
2 | Miles/(US) gallon | Mean |  |  |  |  |
TOTAL 1/2 | Miles/(US) gallon | Mean |  |  |  |  |
3 | Miles/(US) gallon | Mean |  |  |  |  |
4 | Miles/(US) gallon | Mean |  | 26.7 |  | 26.0 | 26.7
TOTAL 3/4 | Miles/(US) gallon | Mean |  | 26.7 |  | 26.0 | 26.7
6 | Miles/(US) gallon | Mean |  | 19.7 |  | 20.6 | 19.1
8 | Miles/(US) gallon | Mean |  | 15.1 |  | 15.1 |
5 and more | Miles/(US) gallon | Mean |  | 16.6 |  | 16.1 | 19.1

### Subtotals at the bottom of the table

    mtcars %>%
        tab_cells(mpg, qsec) %>%
        tab_cols(total(), vs) %>%
        tab_rows(subtotal(cyl, 1:2, 3:4, "TOTAL 5 and more" = 5 %thru% hi, position = "bottom")) %>%
        tab_stat_mean() %>%
        tab_pivot() %>%
        tab_caption("Subtotals at the bottom of the table")

Subtotals at the bottom of the table
|  |  | #Total |  | Engine
--- | --- | --- | --- | --- | ---
|  |  |  |  |  | V-engine | Straight engine
Number of cylinders
1 | Miles/(US) gallon | Mean |  |  |  |  |
| 1/4 mile time | Mean |  |  |  |  |
2 | Miles/(US) gallon | Mean |  |  |  |  |
| 1/4 mile time | Mean |  |  |  |  |
3 | Miles/(US) gallon | Mean |  |  |  |  |
| 1/4 mile time | Mean |  |  |  |  |
4 | Miles/(US) gallon | Mean |  | 26.7 |  | 26.0 | 26.7
| 1/4 mile time | Mean |  | 19.1 |  | 16.7 | 19.4
6 | Miles/(US) gallon | Mean |  | 19.7 |  | 20.6 | 19.1
| 1/4 mile time | Mean |  | 18.0 |  | 16.3 | 19.2
8 | Miles/(US) gallon | Mean |  | 15.1 |  | 15.1 |
| 1/4 mile time | Mean |  | 16.8 |  | 16.8 |
TOTAL 1/2 | Miles/(US) gallon | Mean |  |  |  |  |
| 1/4 mile time | Mean |  |  |  |  |
TOTAL 3/4 | Miles/(US) gallon | Mean |  | 26.7 |  | 26.0 | 26.7
| 1/4 mile time | Mean |  | 19.1 |  | 16.7 | 19.4
TOTAL 5 and more | Miles/(US) gallon | Mean |  | 16.6 |  | 16.1 | 19.1
| 1/4 mile time | Mean |  | 17.2 |  | 16.7 | 19.2

### Nets

Net, contrary to `subtotal`, remove original categories.

    mtcars %>%
        tab_cells(mpg) %>%
        tab_cols(total(), vs) %>%
        tab_rows(net(cyl, 1:2, 3:4, "NET 5 and more" = 5 %thru% hi, prefix = "NET ")) %>%
        tab_stat_mean() %>%
        tab_pivot() %>%
        tab_caption("Nets in rows, custom prefix")

Nets in rows, custom prefix
|  |  | #Total |  | Engine
--- | --- | --- | --- | --- | ---
|  |  |  |  |  | V-engine | Straight engine
Number of cylinders
NET 1/2 | Miles/(US) gallon | Mean |  |  |  |  |
NET 3/4 | Miles/(US) gallon | Mean |  | 26.7 |  | 26.0 | 26.7
NET 5 and more | Miles/(US) gallon | Mean |  | 16.6 |  | 16.1 | 19.1

### Nets with complex grouping

    mtcars %>%
        tab_cells(net(mpg, "Low mpg" = less(mean(mpg)), "High mpg" = greater_or_equal(mean(mpg)))) %>%
        tab_cols(total(), am) %>%
        tab_stat_cases() %>%
        tab_pivot() %>%
        tab_caption("Nets with complex grouping")

Nets with complex grouping
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic | Manual
Miles/(US) gallon
Low mpg | 18 |  | 15 | 3
High mpg | 14 |  | 4 | 10
#Total cases | 32 |  | 19 | 13

### Significance testing on column percent

Letters marks cells which are significantly greater than cells in the appropriate columns. `-` and `+` marks values which are lower/greater than values in the first column. Significance testing on column percent should be applied on the result of `tab_stat_cpct` with total row.

    mtcars %>%
        tab_cells(cyl) %>%
        tab_cols(total(), vs) %>%
        tab_stat_cpct() %>%
        tab_pivot() %>%
        significance_cpct(compare_type = c("first_column", "subtable"), sig_level = 0.05) %>%
        tab_caption("Significance testing on column percent")

Significance testing on column percent
| #Total |  | Engine
--- | --- | --- | ---
|  |  | V-engine |  | Straight engine
|  |  | A |  | B
Number of cylinders
4 | 34.4 |  | 5.6 - |  | 71.4 + A
6 | 21.9 |  | 16.7 |  | 28.6
8 | 43.8 |  | 77.8 + |  |
#Total cases | 32 |  | 18 |  | 14

### Significance testing on means

Significance testing on means should be applied on the result of `tab_stat_mean_sd_n`.

    mtcars %>%
        tab_cells(mpg, disp, hp, wt, qsec) %>%
        tab_cols(total(), am) %>%
        tab_stat_mean_sd_n() %>%
        tab_pivot() %>%
        significance_means(compare_type = c("first_column", "subtable")) %>%
        tab_caption("Significance testing on means")

Significance testing on means
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Miles/(US) gallon
Mean | 20.1 |  | 17.1 - |  | 24.4 + A
Std. dev. | 6.0 |  | 3.8 |  | 6.2
Unw. valid N | 32.0 |  | 19.0 |  | 13.0
Displacement (cu.in.)
Mean | 230.7 |  | 290.4 B |  | 143.5 -
Std. dev. | 123.9 |  | 110.2 |  | 87.2
Unw. valid N | 32.0 |  | 19.0 |  | 13.0
Gross horsepower
Mean | 146.7 |  | 160.3 |  | 126.8
Std. dev. | 68.6 |  | 53.9 |  | 84.1
Unw. valid N | 32.0 |  | 19.0 |  | 13.0
Weight (1000 lbs)
Mean | 3.2 |  | 3.8 + B |  | 2.4 -
Std. dev. | 1.0 |  | 0.8 |  | 0.6
Unw. valid N | 32.0 |  | 19.0 |  | 13.0
1/4 mile time
Mean | 17.8 |  | 18.2 |  | 17.4
Std. dev. | 1.8 |  | 1.8 |  | 1.8
Unw. valid N | 32.0 |  | 19.0 |  | 13.0

### Multiple-response variables with weighting

Here we load data with multiple-responce questions. `mrset` means that we treat set of variables as multiple response varibale with category encoding. For dichotomy encoding use `mdset`.

    data(product_test)
    codeframe_likes = num_lab("
                              1 Liked everything
                              2 Disliked everything
                              3 Chocolate
                              4 Appearance
                              5 Taste
                              6 Stuffing
                              7 Nuts
                              8 Consistency
                              98 Other
                              99 Hard to answer
                              ")

    set.seed(1)
    product_test = product_test %>%
        let(
            # recode age by groups
            age_cat = recode(s2a, lo %thru% 25 ~ 1, lo %thru% hi ~ 2),
            wgt = runif(.N, 0.25, 4),
            wgt = wgt/sum(wgt)*.N
        )  %>%
        apply_labels(
            age_cat = "Age",
            age_cat = c("18 - 25" = 1, "26 - 35" = 2),
            a1_1 = "Likes. VSX123",
            b1_1 = "Likes. SDF456",
            a1_1 = codeframe_likes,
            b1_1 = codeframe_likes
        )

    product_test %>%
        tab_cells(mrset(a1_1 %to% a1_6), mrset(b1_1 %to% b1_6)) %>%
        tab_cols(total(), age_cat) %>%
        tab_weight(wgt) %>%
        tab_stat_cpct() %>%
        tab_sort_desc() %>%
        tab_pivot() %>%
        tab_caption("Multiple-response variables with weighting")

Multiple-response variables with weighting
| #Total |  | Age
--- | --- | --- | ---
|  |  | 18 - 25 | 26 - 35
Likes. VSX123
Nuts | 68.0 |  | 82.5 | 53.7
Taste | 39.3 |  | 48.4 | 30.3
Chocolate | 34.5 |  | 31.9 | 37.0
Appearance | 33.6 |  | 26.7 | 40.4
Stuffing | 31.5 |  | 23.6 | 39.3
Consistency | 12.9 |  | 7.7 | 18.0
Disliked everything | 0.7 |  |  | 1.4
Liked everything |  |  |  |
Other |  |  |  |
Hard to answer |  |  |  |
#Total cases | 150 |  | 70 | 80
Likes. SDF456
Nuts | 63.3 |  | 62.1 | 64.4
Taste | 31.6 |  | 35.2 | 28.0
Appearance | 30.3 |  | 28.8 | 31.8
Chocolate | 27.2 |  | 27.9 | 26.4
Stuffing | 21.5 |  | 14.3 | 28.7
Consistency | 8.2 |  | 1.7 | 14.7
Disliked everything | 2.9 |  | 3.3 | 2.5
Other | 1.0 |  |  | 1.9
Liked everything |  |  |  |
Hard to answer |  |  |  |
#Total cases | 150 |  | 70 | 80

### Side-by-side variables comparison

To make side-by-side comparison we use “|” to suppress variable labels and put these labels to the statistic labels. Statistics labels we place in columns with `tab_pivot`.

    product_test %>%
        tab_cols(total(), age_cat) %>%
        tab_weight(wgt) %>%
        # '|' is needed to prevent automatic labels creation from argument
        tab_cells("|" = unvr(mrset(a1_1 %to% a1_6))) %>%
        tab_stat_cpct(label = var_lab(a1_1)) %>%
        tab_cells("|" = unvr(mrset(b1_1 %to% b1_6))) %>%
        tab_stat_cpct(label = var_lab(b1_1)) %>%
        tab_pivot(stat_position = "inside_columns") %>%
        tab_caption("Side-by-side variables comparison")

Side-by-side variables comparison
| #Total |  | Age
--- | --- | --- | ---
| Likes. VSX123 |  | Likes. SDF456 |  | 18 - 25 |  | 26 - 35
|  |  |  |  | Likes. VSX123 | Likes. SDF456 |  | Likes. VSX123 | Likes. SDF456
Liked everything |  |  |  |  |  |  |  |  |
Disliked everything | 0.7 |  | 2.9 |  |  | 3.3 |  | 1.4 | 2.5
Chocolate | 34.5 |  | 27.2 |  | 31.9 | 27.9 |  | 37.0 | 26.4
Appearance | 33.6 |  | 30.3 |  | 26.7 | 28.8 |  | 40.4 | 31.8
Taste | 39.3 |  | 31.6 |  | 48.4 | 35.2 |  | 30.3 | 28.0
Stuffing | 31.5 |  | 21.5 |  | 23.6 | 14.3 |  | 39.3 | 28.7
Nuts | 68.0 |  | 63.3 |  | 82.5 | 62.1 |  | 53.7 | 64.4
Consistency | 12.9 |  | 8.2 |  | 7.7 | 1.7 |  | 18.0 | 14.7
Other |  |  | 1.0 |  |  |  |  |  | 1.9
Hard to answer |  |  |  |  |  |  |  |  |
#Total cases | 150 |  | 150 |  | 70 | 70 |  | 80 | 80

### Multiple tables in the loop with knitr

To make the task more practical we will create table with means for variables which have more than 6 unique elements. For other variables we will calculate column percent table. Note that you need to set `results='asis'` in the chunk options.

    # here we specify dataset and banner
    banner = mtcars %>%
        tab_cols(total(), am)

    for(each in colnames(mtcars)){
        # note ..$ which is used for indirect reference to variable
        # specify variable
        curr_table = banner %>%
            tab_cells(..$each)
        # calculate statistics
        if(length(unique(mtcars[[each]]))>6){
            curr_table = curr_table %>%
                tab_stat_mean_sd_n() %>%
                tab_pivot() %>%
                significance_means()
        } else {
            curr_table = curr_table %>%
                tab_stat_cpct() %>%
                tab_pivot() %>%
                significance_cpct()
        }
        # finalize table
        curr_table %>%
            tab_caption("Variable name: ", each) %>%
            htmlTable() %>%
            print()
    }

Variable name: mpg
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Miles/(US) gallon
Mean | 20.1 |  | 17.1 |  | 24.4 A
Std. dev. | 6.0 |  | 3.8 |  | 6.2
Unw. valid N | 32.0 |  | 19.0 |  | 13.0

Variable name: cyl
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Number of cylinders
4 | 34.4 |  | 15.8 |  | 61.5 A
6 | 21.9 |  | 21.1 |  | 23.1
8 | 43.8 |  | 63.2 B |  | 15.4
#Total cases | 32 |  | 19 |  | 13

Variable name: disp
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Displacement (cu.in.)
Mean | 230.7 |  | 290.4 B |  | 143.5
Std. dev. | 123.9 |  | 110.2 |  | 87.2
Unw. valid N | 32.0 |  | 19.0 |  | 13.0

Variable name: hp
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Gross horsepower
Mean | 146.7 |  | 160.3 |  | 126.8
Std. dev. | 68.6 |  | 53.9 |  | 84.1
Unw. valid N | 32.0 |  | 19.0 |  | 13.0

Variable name: drat
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Rear axle ratio
Mean | 3.6 |  | 3.3 |  | 4.0 A
Std. dev. | 0.5 |  | 0.4 |  | 0.4
Unw. valid N | 32.0 |  | 19.0 |  | 13.0

Variable name: wt
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Weight (1000 lbs)
Mean | 3.2 |  | 3.8 B |  | 2.4
Std. dev. | 1.0 |  | 0.8 |  | 0.6
Unw. valid N | 32.0 |  | 19.0 |  | 13.0

Variable name: qsec
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
1/4 mile time
Mean | 17.8 |  | 18.2 |  | 17.4
Std. dev. | 1.8 |  | 1.8 |  | 1.8
Unw. valid N | 32.0 |  | 19.0 |  | 13.0

Variable name: vs
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Engine
V-engine | 56.2 |  | 63.2 |  | 46.2
Straight engine | 43.8 |  | 36.8 |  | 53.8
#Total cases | 32 |  | 19 |  | 13

Variable name: am
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Transmission
Automatic | 59.4 |  | 100.0 |  |
Manual | 40.6 |  |  |  | 100.0
#Total cases | 32 |  | 19 |  | 13

Variable name: gear
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Number of forward gears
3 | 46.9 |  | 78.9 |  |
4 | 37.5 |  | 21.1 |  | 61.5 A
5 | 15.6 |  |  |  | 38.5
#Total cases | 32 |  | 19 |  | 13

Variable name: carb
| #Total |  | Transmission
--- | --- | --- | ---
|  |  | Automatic |  | Manual
|  |  | A |  | B
Number of carburetors
1 | 21.9 |  | 15.8 |  | 30.8
2 | 31.2 |  | 31.6 |  | 30.8
3 | 9.4 |  | 15.8 |  |
4 | 31.2 |  | 36.8 |  | 23.1
6 | 3.1 |  |  |  | 7.7
8 | 3.1 |  |  |  | 7.7
#Total cases | 32 |  | 19 |  | 13