diff --git a/10-fatorial-duplo-1-adicional.Rmd b/10-fatorial-duplo-1-adicional.Rmd
index 94b7db62461408eef6fa35b94ddde711b044e7b4..e34b4bb1c9ba9120ba9275860f90f840130d978f 100644
--- a/10-fatorial-duplo-1-adicional.Rmd
+++ b/10-fatorial-duplo-1-adicional.Rmd
@@ -192,7 +192,8 @@ feita.
 ```{r}
 # Linhas distintas da matriz do modelo com as colunas de efeitos
 # estimáveis.
-unname(unique(model.matrix(m0)[, which(!is.na(coef(m0)))]))
+s <- names(coef(m0))[!is.na(coef(m0))]
+unname(unique(model.matrix(m0)[, s]))
 ```
 
 O 2º coeficiente corresponde ao efeito do fator A. Ele possui 2 níveis