Error en el Modelo de Efectos Fijos

Question

Quiero hacer un modelo de efectos fijos en un conjunto de datos ficticio que tiene 3 columnas que representan la gravedad de los pacientes y su tasa de supervivencia en diferentes hospitales.

La TASA DE SUPERVIVENCIA es la variable dependiente y la GRAVEDAD como mi variable independiente con la columna HOSPITAL como covariable con efecto fijo. Quiero solo un modelo de efectos fijos y no efectos aleatorios ya que estoy tratando de entender la diferencia entre efectos fijos y ols.

Intenté usar la función lme y estoy obteniendo este error

model_fixed1<-lme(SURVIVALRATE~SEVERITY,datos=data_hos)
Error en getGroups.data.frame(dataMix, groups) : 
  fórmula no válida para grupos

También intenté la función lmer

model_fixed2<-lmer(SURVIVALRATE ~ SEVEREITY,datos=data_hos)
Error: No se especificaron términos de efectos aleatorios en la fórmula

¿Está mal mi sintaxis? ¿cómo lo hago? También, si es posible, guíame sobre cómo hacer efectos aleatorios a partir de ahí.

Answer 1

No existe tal cosa como un modelo mixto lineal con solo efectos fijos. Eso es solo regresión lineal/mínimos cuadrados ordinarios, y se ajusta utilizando la función lm:

modelo_fijo1<-lm(TASADESOBREVIVENCIA~GRAVEDAD,data=datos_hospital)