+
We use the Poisson(5) distribution as an importance sampler. Note constant exp(-5) needs to be accounted for. #generate a sample of n from the Poisson(5) n=100000 x=rpois(n,5) est=0 s2=0 for(i in 1:n) { est=est+sin(x[i]**2) s2=s2+(sin(x[i]**2))**2 } est=exp(5)*est/n s2=exp(5)*exp(5)*s2/n-est**2 cat("estimate +- 3 standard errors","\n") cat(n,est-3*sqrt(s2/n),est,est+3*sqrt(s2/n),"\n") 10000 -42.25446 -39.57449 -36.89451 1e+05 -39.7936 -38.94073 -38.08785
April 10 1-3pm
Join Zoom Meeting
https://utoronto.zoom.us/j/89664260184
Meeting ID: 896 6426 0184
Passcode: 654227
April 16 1-3pm
Join Zoom Meeting
https://utoronto.zoom.us/j/81802845727
Meeting ID: 818 0284 5727
Passcode: 756900
April 17 1-3pm.
Join Zoom Meeting
https://utoronto.zoom.us/j/83135124574
Meeting ID: 831 3512 4574
Passcode: 171127
First class is Monday, January 8.
The following topics will be covered. The material on Monte Carlo will be interspersed throughout the course.