Tag: GLM
Ejecuto un glm() utilizando errores estándar robustos. Para una comparación de modelos posterior, calculo la diferencia entre dos modelos de regresión (coeficientes y errores estándar). Para ese cálculo, utilizo la función summary(). Sin embargo, la función summary de los modelos muestra errores estándar diferentes a los que obtengo de coeftest(). . . . Read more
Estoy intentando optimizar manualmente un modelo de regresión binomial negativa utilizando el paquete optim en R para predecir una variable de conteo y utilizando una matriz de factores X utilizando el siguiente código: “`{r} <h1>generando datos ficticios</h1> n <- 1000 X <- matrix(NA, ncol = 5, nrow = n) X[,1] . . . Read more
Quiero encontrar un valor umbral para poder clasificar qué observación se considera éxito o no. Pero estoy confundido acerca de cómo usar el paquete para regresión logística. Por ejemplo, mi modelo glm es: library(vcd) library(ggplot2) library(dplyr) library(MLmetrics) library(pROC) d=read.delim(“http://dnett.github.io/S510/Disease.txt”) d$disease=factor(d$disease) d$ses=factor(d$ses) d$sector=factor(d$sector) model=glm(disease~age+sector, family=binomial(link=logit), data=d) Donde éxito se define como . . . Read more
He realizado una regresión logística para un resultado (un tipo de efecto secundario – si los pacientes lo tienen o no). La fórmula y los resultados de este modelo se muestran a continuación: modelo <- glm(side<em>effect</em>G1 ~ edad + bmi + tipoCirugía1 + tipoCirugía2 + tipoCirugía3 + cvd + rt_áxila, . . . Read more
¿Alguien sabe cómo extraer las estimaciones de parámetros de un término aleatorio cuando se utiliza la sintaxis (1 | …) en un modelo glmer (incluyendo se, t ratio y p-value)? Solo puedo acceder a la varianza promedio y la desviación estándar con la función summary. Algunos antecedentes: Utilicé términos aleatorios . . . Read more