Supposons que j'ai variables aléatoires normales indépendantes
et . Comment pourrais-je caractériser la densité de si la distribution de chaque est tronquée à l'intérieur ? En d'autres termes, partir de distributions normales indépendantes , je échantillons qui ne se pas à moins de de chaque moyenne et je les additionne.
En ce moment, je fais cela avec le code R ci-dessous:
x_mu <- c(12, 18, 7)
x_sd <- c(1.5, 2, 0.8)
a <- x_mu - 2 * x_sd
b <- x_mu + 2 * x_sd
samples <- sapply(1:3, function(i) {
return(rtruncnorm(100000, a[i], b[i], x_mu[i], x_sd[i]))
})
y <- rowSums(samples)
Existe-t-il une méthode pour générer directement la densité de ?
Réponses:
Vous pouvez utiliser l'approximation par les méthodes de point de selle, pour la somme des normales tronquées. Je ne donnerai pas les détails maintenant, vous pouvez consulter ma réponse à la somme générale des distributions Gamma pour obtenir des conseils. Ce dont nous avons besoin, c'est de trouver la fonction de génération de moment pour une normale tronquée, ce qui est facile. Je vais le faire ici pour une normale standard tronquée à± 2 , qui a une densité
La fonction de génération de moment peut être calculée comme
la source
Je suis curieux de savoir pourquoi, mais oui, il existe un moyen simple de générer le pdf de cette somme de distributions:
la source