The PRIIP (packaged products) regulation prescribes Monte Carlo bootstrapping simulation for calculation of VaR for products of category III (non-linearly leveraged products). The idea is based on Geometric Brownian Motion. The bank shall for each product realize Monte Carlo simulation of future returns and disclose the total-life 97.5% VaR.
There is one example published (https://ec.europa.eu/info/system/files/risk-section-kid-11072016_en.pdf). This example follows this logic:
1) History of 9 daily log-returns is available, with the sample mean of $\mu$=-0,00227252 and standard deviation $\sigma$=0,011551773.
2) Bootstrapping of the returns is performed via random sampling out of these 9 historical log-returns. N=10 000 paths, each T=10 days long are boostrapped. Total return for each path is calculated as sum of sampled log-returns.
3) The path responding to 97.5% quantile is selected, with total log-return of $r_{bootstrapped}^{97,5\%}$=0,029342.
4) The risk-risk free rate of $r_{free}$=1.2% per annum is assumed.
Now the intermediary result is clear, the VaR log-return $r_{bootstrapped}^{97,5\%}$=0,029342 is however not risk-neutral. The regulation prescribes to derive a risk neutral return, based on the risk-free rate $r_{free}$ and historical mean $\mu$ and sample deviation $\sigma$.
How could the risk-neutral return be calculated ?
The example states deemed result of $r_{neutral}^{97,5\%}$=0,055447241. I am however not able to derive to this results. The regulation states formula
$Return = E[Return_{risk-neutral}] - E[Return_{measured}] - 0.5\sigma^2T$
which is very unclear. The calculation I have tried (assuming 261 business days, which is however information not given directly in the example) is based on the logic, that we have to replace the historical mean with risk-free rate. I tried to do it in following way, not sure if correct (?):
$r_{neutral}^{97,5\%}=r_{bootstrapped}^{97,5\%}+(r_{free}*T/261-0,5*\sigma^2*T-\mu*T)=0,029342+(0,012*10/261-0,5*0,011551773^2*10-(-0,00227252*10)) )= 0,05185976$
Is this approach correct or how should one derive risk-neutral return from the bootstrapping Monte Carlo simulation ?