A más gripe, ¿menos mortalidad? En determinados submundos frecuentistas, sí

Estos días he tenido que adaptar y ejecutar con datos españoles una serie de modelos para medir la virulencia de diversos subtipos de gripe. Y todo bien, salvo que para uno de ellos y determinados grupos de edad… a mayor prevalencia, menor mortalidad. ¡Estupendo!

Todo sucede porque un coeficiente que debería haber sido necesariamente positivo fue estimado como negativo (además, significativamente).

Y el coeficiente tenía el signo cambiado (¡error de tipo S!) debido a una serie de problemas sobradamente conocidos:

  • Alta correlación entre los proxies de prevalencia de gripe y otras causas de mortalidad, como el frío.
  • Alta correlación entre los proxies de prevalencia de distintos subtipos de gripe.
  • Una ratio señal/ruido no particularmente favorable.

Todo lo cual pide a gritos un modelo más informativo. Que, como poco, fuerce la restricción >0 en aquellos coeficientes que no puedan ser de otra manera.