Per dettagli sul metodo utilizzato dall'Istituto Superiore di Sanità per la stima dell'efficacia vaccinale, vedere GLM Poisson.
I dati sono tratti dall'ultimo aggiornamento disponibile all'indirizzo EpiCentro: Ultimi Aggiornamenti
Gli aggiornamenti ISS fornisco media e intervalli di confidenza al 95% dell'efficacia vaccinale stimata.
Per ottenere una distribuzione Beta dell'efficacia da questi valori non possiamo utilizzare direttamente gli intervalli di confidenza. Si può però sfruttare l'equivalenza asintotica di una distribuzione Beta ad una distribuzione Normale per $\alpha$ e $\beta$ "sufficientemente grandi" ovvero
\begin{equation} \tag{1} \mathbf{Beta}(\alpha, \beta) \sim_{(\alpha, \beta) \rightarrow (\infty, \infty)} \mathbf{Norm}(\mu, \sigma) \end{equation}Dato quindi che
\begin{equation} \tag{2} \int_{\mu-2\sigma}^{\mu+2\sigma} \mathbf{Norm}(\mu, \sigma) \approx 0.95 \end{equation}possiamo approssimare la deviazione standard della distribuzione Beta come
\begin{equation} \label{sigma} \tag{3} \sigma = \frac{1}{2} \frac{|E_{lo} - E_{up}|}{2} \end{equation}dove $E_{lo}$ e $E_{up}$ sono gli estremi dell'intervallo di confidenza dell'efficacia al 95%.
Sappiamo che la media di una distribuzione Beta è pari a
\begin{equation} \tag{4} \mu = \frac{\alpha}{\alpha + \beta} \end{equation}da cui possiamo facilmente ricavare $\beta$ in funzione di $\alpha$ e $\mu$
\begin{equation} \label{beta} \tag{5} \beta = \frac{\alpha(1-\mu)}{\mu} \end{equation}Sappiamo che la varianza di una distribuzione Beta è pari a
\begin{equation} \tag{6} \sigma^2 = \frac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)} \end{equation}da cui, sostituendo il valore di $\beta$ della (\ref{beta}) e con semplici passaggi algebrici, otteniamo
\begin{equation} \tag{7} \sigma^2 = \frac{\mu^2(1-\mu)}{\alpha+\mu} \end{equation}e possiamo così ricavare $\alpha$ in funzione di $\mu$ e $\sigma$
\begin{equation} \tag{8} \alpha = \frac{\mu^2(1-\mu)}{\sigma^2} - \mu \end{equation}e, infine, $\beta$ in funzione di $\alpha$ e $\mu$ dalla (\ref{beta}).
La distribuzione Beta dell'efficacia $E$ sarà pertanto
\begin{equation} \tag{9} E \sim \mathbf{Beta}\left( \alpha= \frac{\mu^2(1-\mu)}{\sigma^2} - \mu , \beta = \frac{\alpha(1-\mu)}{\mu} \right) \end{equation}dove $\mu$ è la media e $\sigma$ è la deviazione standard, approssimata dalla (\ref{sigma}) con gli estremi dell'intervallo di confidenza al 95%, dei valori di efficacia vaccinale stimata forniti nell'aggiornamento ISS.
Exported from Italia/Vaccini_Efficacia_ISS.ipynb
committed by maxdevblock on Sun May 22 12:01:14 2022 revision 109, dc77c872