-1:h}^m.$$ and $$\label{eq:pZhat}
\widehat{\frac{Z_{H}}{Z_{0}}} = \prod_{h = 1}^H \widehat{\frac{Z_{h}}{Z_{h-1}}} = \prod_{h = 1}^H \sum_{m = 1}^M {W_h^m}{\widetilde{w}}_{h-1:h}^m$$ is an unbiased estimator of $Z_H/Z_0$. Now, $Z_H =  p({\boldsymbol{Y}})$ and $Z_0 = 1 $.

Another estimate is given by the path sampling identity. Indeed, under non-restrictive regularity assumptions, the following equality holds: $$\label{eq:Int}
\log p(Y) - \log Z_0 = \int_{0}^{1} \mathbb{E}_{p_{\rho}}\left[\frac{d \log \gamma_{\rho}(\cdot)}{d \rho} \right]d\rho$$ where $\gamma_\rho(\theta) = {{\widetilde{p}}_{{\boldsymbol{Y}}}}(\theta)^{1-\rho}(\ell({\boldsymbol{Y}}| \theta) \pi(\theta))^{\rho}$, $p_\rho()$ is the associated probability density distribution and, in our geometric path sampling: $$\frac{d \log \gamma_{\rho}(\cdot)}{d \rho} = \log \frac{\ell({\boldsymbol{Y}}| {\theta}) \pi({\theta})}{{{\widetilde{p}}_{{\boldsymbol{Y}}}}({\theta})} = \log \alpha(\theta)$$ An elementary trapezoidal scheme and Monte Carlo approximations of the expectations involved in lead to the following approximation of the marginal likelihood: $$\label{eq:pZhat2}
\widehat{\widehat{\log p({\boldsymbol{Y}})}} = \sum_{h = 1}^H \frac{\rho_h - \rho_{h-1}}{2} (U^M_h + U^M_{h-1})$$ where $U^M_h = \widehat{ \mathbb{E}}_{p_{\rho_h}}\left[ \log \alpha(\theta)\right] = \sum_{m = 1}^M W_{h}^m \log \alpha^m_{h}
 $.

Note that as suggested in [@Zhou2016], we noticed on simulation studies that the two estimators behave similarly in our examples. A precise
-1: h}^m.$$ and $ $ \label{eq: pZhat } 
 \widehat{\frac{Z_{H}}{Z_{0 } } } = \prod_{h = 1}^H \widehat{\frac{Z_{h}}{Z_{h-1 } } } = \prod_{h = 1}^H \sum_{m = 1}^M { W_h^m}{\widetilde{w}}_{h-1: h}^m$$ is an unbiased estimator of $ Z_H / Z_0$. Now, $ Z_H =   p({\boldsymbol{Y}})$ and $ Z_0 = 1 $. 

 Another estimate is give by the way sampling identity. Indeed, under non - restrictive regularity assumptions, the follow equality hold: $ $ \label{eq: Int } 
 \log p(Y) - \log Z_0 = \int_{0}^{1 } \mathbb{E}_{p_{\rho}}\left[\frac{d \log \gamma_{\rho}(\cdot)}{d \rho } \right]d\rho$$ where $ \gamma_\rho(\theta) = { { \widetilde{p}}_{{\boldsymbol{Y}}}}(\theta)^{1-\rho}(\ell({\boldsymbol{Y}}| \theta) \pi(\theta))^{\rho}$, $ p_\rho()$ is the associated probability density distribution and, in our geometric way sampling: $ $ \frac{d \log \gamma_{\rho}(\cdot)}{d \rho } = \log \frac{\ell({\boldsymbol{Y}}| { \theta }) \pi({\theta})}{{{\widetilde{p}}_{{\boldsymbol{Y}}}}({\theta }) } = \log \alpha(\theta)$$ An elementary trapezoidal scheme and Monte Carlo approximations of the expectation involved in lead to the keep up approximation of the borderline likelihood: $ $ \label{eq: pZhat2 } 
 \widehat{\widehat{\log p({\boldsymbol{Y } }) } } = \sum_{h = 1}^H \frac{\rho_h - \rho_{h-1}}{2 } (U^M_h + U^M_{h-1})$$ where $ U^M_h = \widehat { \mathbb{E}}_{p_{\rho_h}}\left [ \log \alpha(\theta)\right ] = \sum_{m = 1}^M W_{h}^m \log \alpha^m_{h } 
  $. 

 notice that as suggested in [ @Zhou2016 ], we noticed on model studies that the two estimators act similarly in our examples. A precise
-1:h}^m.$$ wnd $$\label{eq:pZhat}
\widehat{\fvac{Z_{H}}{Z_{0}}} = \prod_{h = 1}^H \widehav{\frac{Z_{h}}{A_{h-1}}} = \prod_{f = 1}^H \sum_{m = 1}^M {W_h^m}{\widetilde{w}}_{h-1:h}^l$$ us an unbiased estimator of $Z_H/Z_0$. Now, $D_H =  p({\bolesymuol{Y}})$ and $Z_0 = 1 $.

Anovger estliate ls girei by the path ssmpling idantity. Indeed, gnaex non-restrictive regularity assumptijns, the flllowing equalyty njlds: $$\label{eq:Int}
\log p(Y) - \log Z_0 = \int_{0}^{1} \mathgb{E}_{p_{\rho}}\neft[\frac{d \log \gamma_{\rho}(\cdot)}{d \rho} \right]d\rjo$$ wjere $\gamma_\rho(\theta) = {{\widetilde{k}}_{{\golqwymbol{Y}}}}(\theta)^{1-\fho}(\ell({\boldsymbol{Y}}| \thetz) \pi(\theta))^{\rho}$, $p_\rho()$ is the associxted 'robability dwnsljy distributmon anq, in our geometric pdth samlling: $$\frac{d \lon \gamka_{\ryo}(\cdot)}{d \rho} = \log \frac{\xll({\boldsymbol{Y}}| {\theta}) \pi({\theta})}{{{\whdztilde{p}}_{{\boldsymbol{Y}}}}({\thera})} = \log \alpva(\thdra)$$ Xn tlekehtary hra'ezoidal scgeme and Mobte Carlo approximauiogw of the expedtatiogs involved in lead to the following apprmxijation of the marginal oikelihood: $$\label{eq:pZhwt2}
\widehat{\ridehat{\log p({\boldsymbol{Y}})}} = \sum_{h = 1}^H \frac{\rho_h - \rho_{h-1}}{2} (U^K_h + U^J_{f-1})$$ wkcve $U^O_y = \widehat{ \mathbb{E}}_{p_{\rho_h}}\left[ \log \alpha(\theta)\right] = \sim_{k = 1}^M W_{h}^m \log \ajpha^m_{h}
 $.

Note tjay as suggested in [@Zhou2016], se noticed on simupation ftudiws that tre teo estimators behave similaely in our eqampoes. A precise
-1:h}^m.$$ and $$\label{eq:pZhat} \widehat{\frac{Z_{H}}{Z_{0}}} = \prod_{h = = = 1}^H = 1}^M {W_h^m}{\widetilde{w}}_{h-1:h}^m$$ $Z_H/Z_0$. $Z_H = p({\boldsymbol{Y}})$ $Z_0 = 1 Another estimate is given by the sampling identity. Indeed, under non-restrictive regularity assumptions, the following equality holds: $$\label{eq:Int} \log - \log Z_0 = \int_{0}^{1} \mathbb{E}_{p_{\rho}}\left[\frac{d \log \gamma_{\rho}(\cdot)}{d \rho} \right]d\rho$$ where $\gamma_\rho(\theta) = \theta) $p_\rho()$ the probability density distribution and, in our geometric path sampling: $$\frac{d \log \gamma_{\rho}(\cdot)}{d \rho} = \log \frac{\ell({\boldsymbol{Y}}| \pi({\theta})}{{{\widetilde{p}}_{{\boldsymbol{Y}}}}({\theta})} = \log \alpha(\theta)$$ An elementary trapezoidal scheme Monte Carlo approximations of expectations involved in lead to following of the likelihood: \widehat{\widehat{\log = \sum_{h = \frac{\rho_h - \rho_{h-1}}{2} (U^M_h + U^M_{h-1})$$ where $U^M_h = \widehat{ \mathbb{E}}_{p_{\rho_h}}\left[ \log \alpha(\theta)\right] = \sum_{m = 1}^M \log \alpha^m_{h} that as in we on simulation studies two estimators behave similarly in our
-1:h}^m.$$ and $$\label{eq:pZhat}
\widehat{\fRac{Z_{H}}{Z_{0}}} = \prod_{H = 1}^H \widEhaT{\frAc{z_{h}}{Z_{h-1}}} = \Prod_{H = 1}^H \sum_{m = 1}^M {W_h^m}{\wideTIlde{W}}_{h-1:h}^m$$ is an unbiased estimaTor of $z_H/z_0$. now, $Z_h =  P({\bOldsyMbol{Y}})$ anD $z_0 = 1 $.

ANOTheR eStImaTe IS gIven bY thE path saMpling idenTitY. INdeed, under noN-ReStrictive rEguLarity assumpTioNs, the fOlLowINg equAliTy holDs: $$\labeL{Eq:Int}
\lOg p(Y) - \log Z_0 = \iNt_{0}^{1} \MAthbb{E}_{P_{\Rho}}\left[\FRAc{D \log \Gamma_{\rho}(\cdot)}{d \rho} \rIGhT]D\rho$$ where $\gamma_\Rho(\theTa) = {{\WIdETIldE{p}}_{{\bOldsymbol{Y}}}}(\ThEta)^{1-\rhO}(\Ell({\boldSYmBOL{y}}| \thETa) \pi(\theta))^{\rho}$, $p_\Rho()$ is the assOCiaTed proBaBilITy densIty diStRIbuTion and, in ouR geoMetric patH samplINg: $$\frac{d \LOg \gamma_{\Rho}(\cdoT)}{d \rHo} = \lOg \frAC{\eLl({\BolDsYMboL{y}}| {\tHetA}) \Pi({\tHeta})}{{{\wideTiLdE{p}}_{{\bolDsymBOL{y}}}}({\ThetA})} = \loG \alpHa(\theTa)$$ An elementarY trApezOIdaL scheMe and montE CArlo aPproxiMatioNs Of the expectatioNs inVolved in lEad To The FoLlowiNG approXimAtiOn of the MarginaL LikElIHOOd: $$\Label{eq:pZhat2}
\widehaT{\wIDEhAt{\log p({\boLdsymbOL{Y}})}} = \SuM_{H = 1}^H \frac{\rhO_h - \Rho_{H-1}}{2} (U^M_h + u^m_{H-1})$$ wherE $U^M_h = \WIdEhat{ \mathBb{E}}_{p_{\rhO_H}}\lEfT[ \log \alpHa(\Theta)\rIgHt] = \sUm_{m = 1}^m W_{h}^m \lOG \alpHa^m_{h}
 $.

NoTe that as SuggeSTed in [@Zhou2016], we notICed on simulatiON sTUDiES thaT thE two estimatOrs bEHave SimiLArLy iN Our exAmpleS. A PReCIse
-1:h}^m.$$ and $$\label{eq :pZhat}
\w ideha t{\ fra c{ Z_{H }}{Z _{0}}} = \prod _ {h =  1}^H \widehat{\frac{Z _{h}} {Z _ {h-1 } }}  = \p rod_{h=  1 } ^ H \ su m_ {m=1 }^ M {W_ h^m }{\wide tilde{w}}_ {h- 1: h}^m$$ is an un biased est ima tor of $Z_H/ Z_0 $. Now ,$Z_ H  =  p ({\ bolds ymbol{ Y }})$ a nd $Z_0 =  1 $.

An o ther es t i ma te i s given by the pa t hs ampling identi ty. In de e d, u nde r n on-restric ti ve re g ularity as s u m pti o ns, the follo wing equali t y h olds:$$ \la b el{eq: Int}\l o g p (Y) - \logZ_0= \int_{0 }^{1}\ mathbb{ E }_{p_{\ rho}}\ lef t[\ frac { d\l og\g a mma _ {\ rho } (\c dot)}{d\r ho } \ri ght] d \ r h o$$whe re $ \gamm a_\rho(\theta ) =  {{\ w ide tilde {p}}_ {{\b ol dsymb ol{Y}} }}(\t he ta)^{1-\rho}(\e ll({ \boldsymb ol{ Y} }|\t heta) \pi(\t het a)) ^{\rho} $, $p_\ r ho( )$ i s  t he associated prob ab i l it y densit y dist r ib ut i on and,in  ou r ge o m etric  pat h  s ampling:  $$\fr a c{ d\log \g am ma_{\r ho }(\ cdo t)}{d \rho } = \l og \frac {\ell ( {\boldsymbol{Y } }| {\theta})\ pi ( { \t h eta} )}{ {{\widetild e{p} } _{{\ bold s ym bol { Y}}}} ({\th et a }) }  = \log \alpha(\the ta )$$ An  elem entary trapez oidal sche m e and Mont e Ca r lo approximations  of t he expecta t ions inv olved  in lead  to the f o l lowing a ppr oxi mat ion o fthe marginall i keli ho od: $$\ lab el{eq:p Zha t2} 
\w ide ha t{\wideha t{\log p ({ \b ol ds ymb ol{Y} } )}} = \s um _{h  =  1} ^H \f r ac{\rh o_h -  \rh o_ {h - 1}} {2} (U^ M _h +  U^M _{ h- 1})$ $ w he re $U ^M_h = \ widehat { \mathbb {E} } _{p_ {\ rh o_h}}\l eft[ \log \al ph a(\theta)\ ri ght ] = \s u m _{m = 1} ^M W_{h}^m \log \alpha^ m _{h}
 $ .
Notethat  as sugge ste d in [ @Zh o u2016] , we n otice dons i mulat i o nstu di es that th e two  esti ma tors  behave  similarly in oure xam ples. A preci se
-1:h}^m.$$ and_$$\label{eq:pZhat}
\widehat{\frac{Z_{H}}{Z_{0}}} =_\prod_{h = 1}^H \widehat{\frac{Z_{h}}{Z_{h-1}}}_= \prod_{h_=_1}^H \sum_{m_=_1}^M {W_h^m}{\widetilde{w}}_{h-1:h}^m$$ is_an unbiased estimator_of $Z_H/Z_0$. Now, $Z_H_=  p({\boldsymbol{Y}})$_and_$Z_0 = 1 $.

Another estimate is given by the path sampling identity. Indeed, under_non-restrictive_regularity assumptions,_the_following_equality holds: $$\label{eq:Int}
\log p(Y) -_\log Z_0 = \int_{0}^{1} \mathbb{E}_{p_{\rho}}\left[\frac{d_\log \gamma_{\rho}(\cdot)}{d_\rho} \right]d\rho$$ where $\gamma_\rho(\theta) = {{\widetilde{p}}_{{\boldsymbol{Y}}}}(\theta)^{1-\rho}(\ell({\boldsymbol{Y}}| \theta) \pi(\theta))^{\rho}$,_$p_\rho()$_is the associated_probability density distribution and, in our geometric path sampling:_$$\frac{d \log \gamma_{\rho}(\cdot)}{d \rho} = \log_\frac{\ell({\boldsymbol{Y}}| {\theta}) \pi({\theta})}{{{\widetilde{p}}_{{\boldsymbol{Y}}}}({\theta})}_=_\log_\alpha(\theta)$$ An elementary trapezoidal_scheme and Monte Carlo approximations of_the expectations involved in lead to_the following approximation of the marginal likelihood:_$$\label{eq:pZhat2}
\widehat{\widehat{\log p({\boldsymbol{Y}})}} = \sum_{h = 1}^H_\frac{\rho_h - \rho_{h-1}}{2} (U^M_h +_U^M_{h-1})$$ where_$U^M_h = \widehat{ \mathbb{E}}_{p_{\rho_h}}\left[ \log_\alpha(\theta)\right] = \sum_{m_= 1}^M_W_{h}^m \log \alpha^m_{h}
_$.

Note that as suggested in [@Zhou2016],_we noticed on_simulation studies that the two estimators_behave_similarly in our_examples._A_precise