Distribuzione di Kumaraswamy

Template:F In teoria della probabilità la distribuzione di Kumaraswamy è una distribuzione di probabilità continua, definita sull'intervallo [0,1] e dipendente da due parametri. È simile alla variabile casuale beta, ma è più semplice da usare grazie alle semplici espressioni chiuse della funzione di densità di probabilità e della frequenza cumulata. Porta il nome di Poondi Kumaraswamy che la descrisse per primo.

Le funzioni di densità di probabilità delle variabile casuale di Kumaraswamy per alcuni valori dei parametri.

Confronto tra le variabili casuali beta e di Kumaraswamy per una scelta dei parametri.

Confronto tra le variabili casuali beta e di Kumaraswamy per una scelta dei parametri.

Caratteristiche

La funzione di densità di probabilità è definita da

$f(x|a;b)=ab{x}^{a-1}(1-x^a{)}^{b-1}$, dove a e b sono i due parametri e $x\in [0,1]$

si ottiene così che la cumulata è

$F(x|a;b)=[1-(1-x^a{)}^b]$

e il valore atteso diventa

$\mu =\frac{b\mathrm{\Gamma }(1+1/a)\mathrm{\Gamma }(b)}{\mathrm{\Gamma }(1+1/a+b)}$

mentre la mediana è

${(1-{\left(\frac{1}{2}\right)}^{1/b})}^{1/a}$

e la moda

${\left(\frac{a-1}{ab-1}\right)}^{1/a}.$

I momenti di ordine n sono calcolabili con

$m_n=\frac{b\mathrm{\Gamma }(1+n/a)\mathrm{\Gamma }(b)}{\mathrm{\Gamma }(1+b+n/a)}=bB(1+n/a,b).$

dove $\mathrm{\Gamma }$ e $B$ sono rispettivamente la funzione gamma e la funzione beta di Eulero.

Relazione con altre distribuzioni

• Se $X\sim \text{Kumaraswamy}(1,1)$ allora $X\sim U(0,1)$
• Se $X\sim U(0,1)$ (distribuzione continua uniforme) allora ${(1-{(1-X)}^{\frac{1}{b}})}^{\frac{1}{a}}\sim \text{Kumaraswamy}(a,b)$
• Se $X\sim \text{Beta}(1,b)$ (variabile casuale beta) allora $X\sim \text{Kumaraswamy}(1,b)$
• Se $X\sim \text{Beta}(a,1)$ (variabile casuale beta) allora $X\sim \text{Kumaraswamy}(a,1)$
• Se $X\sim \text{Kumaraswamy}(a,1)$ allora $(1-X)\sim \text{Kumaraswamy}(1,a)$
• Se $X\sim \text{Kumaraswamy}(1,a)$ allora $(1-X)\sim \text{Kumaraswamy}(a,1)$
• Se $X\sim \text{Kumaraswamy}(a,1)$ allora $-ln(X)\sim \text{Exponential}(a)$
• Se $X\sim \text{Kumaraswamy}(1,b)$ allora $-ln(1-X)\sim \text{Exponential}(b)$
• Se $X\sim \text{Kumaraswamy}(a,b)$ allora $X\sim \text{GB1}(a,1,1,b)$, la distribuzione beta generalizzata di primo ordine.

Implementazioni in software

In R tramite il pacchetto extraDistr sono disponibili le seguenti funzioni

dkumar(x, a = 1, b = 1, log = FALSE)
pkumar(q, a = 1, b = 1, lower.tail = TRUE, log.p = FALSE)
qkumar(p, a = 1, b = 1, lower.tail = TRUE, log.p = FALSE)
rkumar(n, a = 1, b = 1)

rispettivamente funzione di densità, di probabilità, dei quantili e generatore di numeri casuali.

Bibliografia

• Template:Cita pubblicazione
• Template:Cita pubblicazione
• Template:Cita pubblicazione
• Template:Cita pubblicazione