Insérer des équations

Packages utilisés

On utilise ici le package ‘equatiomatic’ qui va nous permettre de récupérer l’équation du modèle au format ‘latex’.

library(flextable)
library(equatiomatic)

Le flextable

On va générer un flextable depuis un modèle linéaire en utilisant la fonction as_flextable().

mod <- lm(mpg ~ cyl + disp, mtcars)
ft <- as_flextable(mod)
ft

	Estimate	Standard Error	t value	Pr(>\|t\|)
(Intercept)	34.661	2.547	13.609	0.0000	***
cyl	-1.587	0.712	-2.230	0.0337	*
disp	-0.021	0.010	-2.007	0.0542	.
Signif. codes: 0 <= '*' < 0.001 < '' < 0.01 < '*' < 0.05

Residual standard error: 3.055 on 29 degrees of freedom
Multiple R-squared: 0.7596, Adjusted R-squared: 0.743
F-statistic: 45.81 on 29 and 2 DF, p-value: 0.0000

Ajouter l’équation dans le tableau

La fonction equatiomatic::extract_eq() va être utilisée pour récupérer l’équation au format latex:

eq <- extract_eq(mod)
print(eq)
## $$
## \operatorname{mpg} = \alpha + \beta_{1}(\operatorname{cyl}) + \beta_{2}(\operatorname{disp}) + \epsilon
## $$

Il reste enfin à l’ajouter dans le tableau avec la fonction flextable::compose(). Il faut indiquer qu’il faut traiter la valeur comme une équation ‘latex’ en utilisant la fonction as_equation().

ft <- add_header_lines(ft, "", top = TRUE)

ft <- compose(
  x = ft, j = 1, i = 1, part = "header",
  value = as_paragraph(
    as_equation(eq, 
    width = 2, height = .5)
  )
)

ft <- align(ft, i = 1, part = "header", align = "center")
ft

$\operatorname{mpg} = \alpha + \beta_{1}(\operatorname{cyl}) + \beta_{2}(\operatorname{disp}) + \epsilon$
(Intercept)	34.661	2.547	13.609	0.0000	***
cyl	-1.587	0.712	-2.230	0.0337	*
disp	-0.021	0.010	-2.007	0.0542	.
Signif. codes: 0 <= '*' < 0.001 < '' < 0.01 < '*' < 0.05

Residual standard error: 3.055 on 29 degrees of freedom
Multiple R-squared: 0.7596, Adjusted R-squared: 0.743
F-statistic: 45.81 on 29 and 2 DF, p-value: 0.0000

Une autre illustration

eqs <- c(
  "(ax^2 + bx + c = 0)",
  "a \\ne 0",
  "x = {-b \\pm \\sqrt{b^2-4ac} \\over 2a}")
df <- data.frame(formula = eqs)
df
##                                   formula
## 1                     (ax^2 + bx + c = 0)
## 2                                a \\ne 0
## 3 x = {-b \\pm \\sqrt{b^2-4ac} \\over 2a}


ft <- flextable(df)
ft <- compose(
  x = ft, j = "formula",
  value = as_paragraph(as_equation(formula, width = 2, height = .5)))
ft <- align(ft, align = "center", part = "all")
ft <- width(ft, width = 2)
ft

formula
$(ax^2 + bx + c = 0)$
$a \ne 0$
$x = {-b \pm \sqrt{b^2-4ac} \over 2a}$