Penguins Tutorial: Easy to Use, Easy to Read

This tutorial shows the full easytable workflow using palmerpenguins.

1) Fit a sequence of models

library(easytable)
library(palmerpenguins)

m1 <- lm(body_mass_g ~ flipper_length_mm, data = penguins)
m2 <- lm(body_mass_g ~ flipper_length_mm + species, data = penguins)
m3 <- lm(body_mass_g ~ flipper_length_mm + species + island, data = penguins)

2) Create a readable baseline table

easytable(
  m1, m2, m3,
  model.names = c("Baseline", "With Species", "Full Model")
)

Coefficient

Baseline

With Species

Full Model

(Intercept)

-5780.83 ***
(305.81)

-4031.48 ***
(584.15)

-4047.52 ***
(585.52)

flipper_length_mm

49.69 ***
(1.52)

40.71 ***
(3.07)

41.09 ***
(3.09)

species:Chinstrap

-206.51 ***
(57.73)

-205.59 ***
(70.35)

species:Gentoo

266.81 ***
(95.26)

199.81 *
(109.71)

island:Dream

-59.75
(75.74)

island:Torgersen

-101.92
(77.65)

N

342

342

342

R sq.

0.76

0.78

0.78

Adj. R sq.

0.76

0.78

0.78

Significance: ***p < .01; **p < .05; *p < .1

By default, each coefficient cell is two lines:

  1. Estimate + significance stars
  2. Standard error in parentheses

3) Highlight significant coefficients

easytable(
  m1, m2, m3,
  model.names = c("Baseline", "With Species", "Full Model"),
  highlight = TRUE
)

Coefficient

Baseline

With Species

Full Model

(Intercept)

-5780.83 ***
(305.81)

-4031.48 ***
(584.15)

-4047.52 ***
(585.52)

flipper_length_mm

49.69 ***
(1.52)

40.71 ***
(3.07)

41.09 ***
(3.09)

species:Chinstrap

-206.51 ***
(57.73)

-205.59 ***
(70.35)

species:Gentoo

266.81 ***
(95.26)

199.81 *
(109.71)

island:Dream

-59.75
(75.74)

island:Torgersen

-101.92
(77.65)

N

342

342

342

R sq.

0.76

0.78

0.78

Adj. R sq.

0.76

0.78

0.78

Significance: ***p < .01; **p < .05; *p < .1

4) Collapse controls into indicator rows

easytable(
  m1, m2, m3,
  model.names = c("Baseline", "With Species", "Full Model"),
  control.var = "island",
  highlight = TRUE
)

Coefficient

Baseline

With Species

Full Model

(Intercept)

-5780.83 ***
(305.81)

-4031.48 ***
(584.15)

-4047.52 ***
(585.52)

flipper_length_mm

49.69 ***
(1.52)

40.71 ***
(3.07)

41.09 ***
(3.09)

species:Chinstrap

-206.51 ***
(57.73)

-205.59 ***
(70.35)

species:Gentoo

266.81 ***
(95.26)

199.81 *
(109.71)

island

Y

N

342

342

342

R sq.

0.76

0.78

0.78

Adj. R sq.

0.76

0.78

0.78

Significance: ***p < .01; **p < .05; *p < .1

control.var is useful when models include many factor levels or fixed effects.

5) LaTeX output for PDF workflows

easytable(
  m1, m2, m3,
  model.names = c("Baseline", "With Species", "Full Model"),
  output = "latex",
  control.var = "island",
  highlight = TRUE
)

6) Robust standard errors (optional dependency path)

if (requireNamespace("lmtest", quietly = TRUE) &&
    requireNamespace("sandwich", quietly = TRUE)) {
  easytable(
    m1, m2, m3,
    model.names = c("Baseline", "With Species", "Full Model"),
    robust.se = TRUE
  )
}

Coefficient

Baseline

With Species

Full Model

(Intercept)

-5780.83 ***
(292.56)

-4031.48 ***
(542.24)

-4047.52 ***
(536.56)

flipper_length_mm

49.69 ***
(1.44)

40.71 ***
(2.86)

41.09 ***
(2.84)

species:Chinstrap

-206.51 ***
(52.21)

-205.59 ***
(66.68)

species:Gentoo

266.81 ***
(95.08)

199.81 *
(112.03)

island:Dream

-59.75
(82.45)

island:Torgersen

-101.92
(83.82)

N

342

342

342

R sq.

0.76

0.78

0.78

Adj. R sq.

0.76

0.78

0.78

Significance: ***p < .01; **p < .05; *p < .1

Note: Robust Standard Errors

7) Marginal effects (optional dependency path)

if (requireNamespace("margins", quietly = TRUE)) {
  easytable(
    m2, m3,
    model.names = c("With Species", "Full Model"),
    margins = TRUE
  )
}

Coefficient

With Species

Full Model

flipper_length_mm

40.71 ***
(3.07)

41.09 ***
(3.09)

species:Chinstrap

-206.51 ***
(57.73)

-205.59 ***
(70.35)

species:Gentoo

266.81 ***
(95.26)

199.81 *
(109.71)

island:Dream

-59.75
(75.74)

island:Torgersen

-101.92
(77.65)

N

342

342

R sq.

0.78

0.78

Adj. R sq.

0.78

0.78

Significance: ***p < .01; **p < .05; *p < .1

Note: Average Marginal Effects (AME)

8) Export outputs

easytable(
  m1, m2, m3,
  model.names = c("Baseline", "With Species", "Full Model"),
  highlight = TRUE,
  export.word = "penguins_table.docx",
  export.csv = "penguins_table.csv"
)

Design Notes

easytable keeps these display invariants:

  1. Two-line coefficient cells.
  2. Zebra only in coefficient rows.
  3. One divider between coefficient rows and model-stat rows.
  4. Control indicators in the model-stat block.

These defaults are intentional so tables remain legible in long workflows.