1 Verosimiglianza profilo.

Vengono caricati i packages necessari per realizzare questo documento

1.1 Verosimiglianza totale

Verosimiglianza rispetto a \(\mu,\sigma^2\) per un campione di ampiezza 20 estratto da una normale standard, con osservazioni \(x_i, \ (i=1,2,..,20)\):

 [1] -1.2071  0.2774  1.0844 -2.3457  0.4291  0.5061 -0.5747 -0.5466 -0.5645 -0.8900 -0.4772 -0.9984 -0.7763  0.0645
[15]  0.9595 -0.1103 -0.5110 -0.9112 -0.8372  2.4158

You must enable Javascript to view this page properly.

1.2 Verosimiglianza profilo

Verosimiglianza profilo: per ciascun valore fissato \(\mu\) stimiamo \(\sigma^2\) da un campione di ampiezza 20 mediante il classico stimatore di massima verosimiglianza, ma in funzione di un particolare valore di \(\mu\), e per questo espresso come \(\hat{\sigma}^2(\mu)\):

\[ \hat{\sigma}^2(\mu)=\frac{\sum_{i=1}^{n} (x_i-\mu)^2}{n}\] ad esempio per il valore di \(\mu = -0.51\) la verosimiglianza è massima per \[ \hat{\sigma}^2(\mu)=\frac{\sum_{i=1}^{n} (x_i-(-0.51))^2}{n} = 1.05 \]

You must enable Javascript to view this page properly.

Nel grafico corrisponde a fare una sezione in corrispondenza del valore \(\mu = -0.51\) e calcolare il valore di \(\hat{\sigma}^2(\mu)\) per la quale è massima questa curva ottenuta come intersezione fra il piano con \(\mu\) costante e la verosimiglianza totale (in colore blu)

1.3 Rappresentazione dell’intera verosimiglianza profilo

A questo punto calcoliamo \(\hat{\sigma}^2(\mu)\) per tutti valori di \(\mu\) e rappresentiamo la sulla verosimiglianza la curva \(L(\hat{\sigma}^2(\mu),\mu )\), che è ora solo funzione di , e comprende i punti sulla superficie di \(L(\cdot)\), in corrispondenza di ciascun valore di \(\mu\) e dei corrispondenti ottimi \(\hat{\sigma}^2(\mu)\) (in rosso nel grafico)

glX 
  4