[R] Quantiles of sums of independent discrete random variables

Greetings, Minimal reproducible example as requested by the technical expert Jeff Newmiller: library(bayesmeta) # density of $(1/10)*\sum_{j=1}{10}N(j,0.01$ # (convex sum of normal distributions) # f <- Vectorize(function(s) sum(vapply(1:10,    FUN = function(j) dnorm(s,mean=j,sd=0.01)/10, FUN.V

[R] Quantiles of sums of independent discrete random variables

Greetings, I have the following  Problem: Given k (=10) discrete independent random variables X_i with n_i (= 5 to 20) values each,compute quantiles of the distribution of the sum X = X_1+...+X_k. Here X has n=n_1 x n_2 ... n_k distinct values which is too large to list them all together with t