Regressione nonlineare

Approssimazione di un set di osservazioni tramite polinomi di diversi gradi

In statistica la regressione non lineare è un metodo di stima di una curva interpolante un modello della forma:

$\ Y=f(X;\vartheta )+\varepsilon$

su un insieme di osservazioni (eventualmente multi-dimensionali), concernenti le variabili $\ X$, $\ Y$.

Metodi di stima

Diversamente da quanto accade nel caso della regressione lineare, non esiste un metodo generale per determinare i valori dei parametri che garantiscono la migliore interpolazione dei dati. A tal fine, si ricorre a classi di algoritmi numerici di ottimizzazione, che a partire da valori iniziali, scelti a caso o tramite un'analisi preliminare, giungono a punti ritenuti ottimali. Si potrebbero avere dei massimi locali della bontà del fitting, in contrasto ancora con il caso della regressione lineare, in cui il massimo è di natura globale.

Linearizzazione

Diversi modelli nonlineari possono essere linearizzati (cioè trasformati in modelli lineari, riducendo così l'onerosità numerica del problema di stima). Un esempio è dato dai modelli:

$\ Y=\alpha \exp \left\{\beta X\right\}$

e

$\ Y=\alpha X_{1}^{\beta }X_{2}^{\gamma }$

Il primo emerge naturalmente in una varietà di contesti, come soluzione di equazioni differenziali ordinarie; il secondo è tipico dell'ambito dell'economia e dell'econometria, come modello per la funzione di produzione. In entrambi i casi è possibile linearizzare i modelli applicando una trasformazione logaritmica.

I logaritmi nella regressione

Caso	Regressione	Interpretazione di $\beta _{1}$
lin-log	$Y_{i}=\beta _{0}+\beta _{1}\ln(X_{i})+u_{i}$	Una variazione percentuale dell'1% in X determina una variazione pari a 0,01 $\beta _{1}$ in Y.
log-lin	$\ln(Y_{i})=\beta _{0}+\beta _{1}(X_{i})+u_{i}$	Una variazione di un'unità in X (ΔX = 1) determina una variazione pari al 100 $\beta _{1}$ % in Y.
log-log	$\ln(Y_{i})=\beta _{0}+\beta _{1}\ln(X_{i})+u_{i}$	Una variazione pari all'1% in X determina una variazione pari al $\beta _{1}$ % in Y (elasticità).

Ulteriori metodi

Modelli di maggiore complessità, quali ad esempio quelli caratterizzati da equazioni trascendenti come $\ Y=A_{1}e^{B_{1}X}+A_{2}e^{B_{2}X}$, sono stimati tramite algoritmi più sofisticati. Diversi software matematici contengono librerie di ottimizzazione: Gauss, GNU Octave, MATLAB, Mathematica; sono inoltre ampiamente disponibili librerie di ottimizzazione per linguaggi avanzati quali C++ o Fortran.

Chiarimenti sull'ambito di applicazione del metodo

Sovente si considera - erroneamente - che l'impiego del metodo dei minimi quadrati per stimare i parametri $\ \alpha$, $\ \beta$, $\ \gamma$ in un modello del tipo:

$\ Y_{i}=\alpha X_{i}^{2}+\beta X_{i}+\gamma +\varepsilon _{i}$

costituisca un caso di regressione nonlineare. In realtà, l'aggettivo (non-)lineare è riferito ai parametri, non alla(e) variabile(i) dipendente(i), per cui il modello sopra è stimato tramite i minimi quadrati ordinari come un modello di regressione lineare; si veda al riguardo il relativo articolo.

Voci correlate

• Interpolazione
• Interpolazione polinomiale
• Interpolazione spline
• Regressione lineare
• Algoritmo di Gauss-Newton

Controllo di autorità	Thesaurus BNCF 57284 · GND (DE) 4251077-6