Dato una variabile casuale con funzione di densità di probabilità f (x), come calcolare il valore atteso di questa variabile casuale in R?Dato una variabile casuale con funzione di densità di probabilità f (x), come calcolare il valore atteso di questa variabile casuale in R?
Se si vuole calcolare il valore atteso, basta calcolare:
E (X) = integrale di xf (x) dx su tutto il dominio di X.
usandoL'integrazione può facilmente essere fatto la funzione integra().
dici che stai avendo una normale funzione di densità (è possibile definire facilmente il proprio funzione di densità):
f <- function(x){
1/sqrt(2*pi)*exp((-1/2)*x^2)
}
si calcola il valore atteso semplicemente:
f2 <- function(x){x*f(x)}
integrate(f2,-Inf,Inf)
Prestare attenzione, a volte devi usare Vectorize() per la tua funzione. Questo è necessario per integrarsi al lavoro. Per maggiori informazioni, consultare le pagine della guida di integrate() e Vectorize().
Aiuta a sapere che l'aspettativa E è l'integrale di x*f(x) dx per x in (-inf, inf)?
Perché è necessario 'Vectorize'? 'f2' è già una funzione vettoriale. O mi sono perso qualcosa? – Marek
In questo caso, lo è. Volevo solo segnalare l'errore evidente quando si utilizza l'integrazione. L'ho visto molte volte prima. Correggerò –