El problema de la estimación inversa

Supongamos que tenemos unos niños de los que sabemos las edades x_i y las alturas y_i. Supongamos además que podemos estimar las segundas en función de las primeras con un modelo lineal clásico

y_i \sim N(a_0 + a_1 x_1, \sigma).

Este modelo nos permite, dada una edad, estimar la altura y los correspondientes intervalos de confianza. Pero, dada una altura, ¿qué nos dice de la edad? Este es el problema conocido como de la estimación inversa.

En este caso concreto es relativamente sencillo. Pero, ¿qué si y_i depende de más de una variable predictora? ¿O si la regresión no es lineal? ¿O si…?

Estos problemas y algunas maneras de afrontarlos se discuten en investr: An R Package for Inverse Estimation. Artículo que, cuyo nombre bien indica, viene acompañado de un paquete de R, investr. Que me hubiese resultado muy útil en más de una ocasión.

Nota: Desde hace unas semanas vengo prefieriendo representar el modelo lineal como arriba en lugar de la manera más tradicional, y_i = a_0 + a_1 x_1 + \epsilon_i. Cuando tenga más claro el motivo, os lo cuento.