Coleção modelos
Regressão linear
Ajusta um modelo linear (lm) pelas colunas ou por uma fórmula digitada.
O que o bloco faz
models/lm Ajusta uma regressão por mínimos quadrados e produz models/fit. A saída é models/fit.
Antes de continuar
Termos que esta página usa, em poucas palavras. As definições completas estão no glossário.
- Graus de liberdade
- Quantas informações independentes sobram para estimar uma variabilidade depois de estimados os outros parâmetros. Definem qual distribuição t, F ou qui-quadrado se usa no teste; com poucos graus de liberdade a inferência fica mais incerta.
- Resíduo
- A diferença entre o valor observado e o valor que o modelo ajustado prevê para aquela observação. É nos resíduos que se conferem os pressupostos do modelo: eles devem parecer ruído, sem padrão.
- Modelo linear
- Um modelo em que a resposta é a soma de efeitos das variáveis explicativas, cada um multiplicado por um coeficiente, mais um erro aleatório. Regressão e análise de variância são casos dele; “linear” se refere aos coeficientes, não à forma da curva.
- Fator
- Uma variável explicativa categórica cujo efeito se quer estudar, como adubo, cultivar ou dose. Num experimento fatorial, dois ou mais fatores são combinados.
- Interação
- Há interação quando o efeito de um fator depende do nível de outro: o adubo aumenta a produção numa cultivar e não em outra. Com interação, olhar o efeito médio de cada fator sozinho pode enganar.
Quando usar
Use quando a resposta contínua e os preditores puderem ser descritos por uma combinação linear; use a fórmula para interações, curvatura ou transformações.
Pressupostos
A relação entre a resposta e os preditores é a da fórmula (linear nos parâmetros, com as interações e curvaturas escritas).
Verificar:Diagnóstico dos resíduos
Se falhar: Acrescente termos (
poly(x, 2),a * b) ou transforme; teste o termo a mais com omodels/compare.As observações são independentes.
Se falhar: Medidas repetidas ou agrupadas (mesmo animal, parcela, local) pedem o
models/lmer.Os erros são normais.
Verificar:Normalidade dos resíduosDiagnóstico dos resíduos
Se falhar: Transforme a resposta, ou use o
models/glmcom a família da resposta.A variância do erro é constante (não cresce com o ajustado nem muda entre grupos).
Verificar:Breusch-PaganLeveneDiagnóstico dos resíduos
Se falhar: Transforme a resposta (log) ou use o
models/glmcom família gama.Sem colinearidade forte entre os preditores: com ela os coeficientes ficam instáveis e os erros padrão, grandes.
Verificar:Coeficientes
Configuração
Use resposta e preditores para efeitos principais; preencha Fórmula para interações, transformações ou codificação explícita de fatores. Fórmula prevalece sobre as colunas.
Exemplo
library(trama)
reg <- tr_registry()
tr_use("trama.models", registry = reg)
tr_flow(reg) |>
tr_add("dados", "models/example", dataset = "mtcars") |>
tr_add("ajuste", "models/lm", formula = "mpg ~ wt + hp", from = "dados")mtcars oferece consumo mpg e preditores wt e hp.
Como interpretar
R² e R² ajustado resumem variação explicada; o teste F global avalia preditores em conjunto. Coeficientes dependem das unidades e dos contrastes dos fatores.
Implementação
stats::lm()versão 4.6.0Mínimos quadrados; texto e lógico viram fator (contraste de tratamento, o padrão do R), número fica número.
No código O bloco chama trama.models::tr_models_lm(), em trama.models/R/ajustar.R:59.
A conta é feita nesta chamada:
stats::lm(f, data = p$dados)
Outras funções de pacotes que o bloco usa (6)
stats::as.formula(texto)stats::as.formula(paste(.tr_models_bt(resposta), "~", rhs))stats::complete.cases(d[, vars, drop = FALSE])stats::df.residual(ajuste)stats::nobs(ajuste)stats::coef(ajuste)
Referências
Livro-texto
- Rencher, A. C.; Schaalje, G. B. (2008). Linear Models in Statistics. 2. ed. Hoboken: Wiley.
- Montgomery, D. C. (2017). Design and Analysis of Experiments. 9. ed. Hoboken: Wiley.