Displaying 1 result from an estimated 1 matches for "nconvolutions".
Did you mean:
convolutions
2024 Jan 24
0
Quantiles of sums of independent discrete random variables
...#
f <- Vectorize(function(s) sum(vapply(1:10,
?? FUN = function(j) dnorm(s,mean=j,sd=0.01)/10, FUN.VALUE=0
)))
g <- function(s) dnorm(s,mean=0,sd=0.01)
cat("\n\n")
for(i in 1:5){
? ?
? cat("Doing convolution ",i,"\n")
? g <- convolve(g,f)$density
}
cat("\nConvolutions finished, plotting density.")
s <- seq(0,100,length.out=1024)
matplot(s,g(s),type="l")
??
?
Michael Meyer
[[alternative HTML version deleted]]