Let g:R→R be an unknown function, twice continuously differentiable with ∣g′′(x)∣⩽M for all x∈R. For some x0∈R, we know the value g(x0) and we wish to estimate its derivative g′(x0). To do so, we have access to a pseudo-random number generator that gives U1∗,…,UN∗ i.i.d. uniform over [0,1], and a machine that takes input x1,…,xN∈R and returns g(xi)+εi, where the εi are i.i.d. N(0,σ2).
(a) Explain how this setup allows us to generate N independent Xi=x0+hZi, where the Zi take value 1 or −1 with probability 1/2, for any h>0.
(b) We denote by Yi the output g(Xi)+εi. Show that for some independent ξi∈R
Yi−g(x0)=hZig′(x0)+2h2g′′(ξi)+εi
(c) Using the intuition given by the least-squares estimator, justify the use of the estimator g^N given by
g^N=N1i=1∑NhZi(Yi−g(x0))
(d) Show that
E[∣g^N−g′(x0)∣2]⩽4h2M2+Nh2σ2.
Show that for some choice hN of parameter h, this implies