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 ?