 $$\frac{\partial^2 L}{\partial \mu \, \partial \sigma_{v}} = \frac{\partial^2 L}{\partial \sigma_{v} \, \partial \mu} = -2\sum_i \frac{\sigma_{v} (f'_{{\rm m},j} - \mu)}{(\sigma_{v}^2 + \sigma_{{\rm o}, j}^2)^2},$$ and $$\frac{\partial^2 L}{\partial^2 \sigma_{v}} = \sum_i \frac{(\sigma_{{\rm o}, j}^2 - 3\sigma_{v}^2)\times(f'_{{\rm m},j} - \mu)^2}{(\sigma_{v}^2 + \sigma_{{\rm o}, j}^2)^3} - \frac{(\sigma_{{\rm o}, j}^2 - \sigma_{v}^2)}{(\sigma_{v}^2 + \sigma_{{\rm o}, j}^2)^2}.$$

Then, we computed the inverse of the minus Hessian, $h_{xy} = (-H_{xy})^{-1}$. Finally, and because maximum likelihood theory states that $\sigma_{\theta_{x}}^2 \leq h_{xx}$, we estimated the variances of our inferred parameters as $\sigma_{\mu}^2 = h_{11}$ and $\sigma_{\sigma_{v}}^2 = h_{22}$.

We tested the performance and the limitations of our MLE through synthetic catalogues of merger fractions. We created several sets of 1000 synthetic catalogues, each of them composed by a number $n$ of merger fractions randomly drawn from a log-normal distribution with $\mu_{\rm in} = \log 0.05$ and $\sigma_{v, {\rm in}} = 0.2$, and affected by observational errors $\sigma_{\rm o}$. We explored the $n = 50, 250$ and $1000$ cases for the number of merger fractions, and varied the observational errors from $\sigma_{\rm o} = 0.1$ to 0.5 in 0.1 steps. That is, we explored observational errors in the measurement of the merger fraction from $\Delta \sigma \equiv \sigma_{\rm o}/\sigma_{v} = 0.5$ to 2.5 times the cosmic variance that we want to measure. We checked that the results below are similar for any value of $\sigma_{v, {\rm in}}$. We find that

1.  The median value of the recovered
$ $ \frac{\partial^2 L}{\partial \mu \, \partial \sigma_{v } } = \frac{\partial^2 L}{\partial \sigma_{v } \, \partial \mu } = -2\sum_i \frac{\sigma_{v } (f'_{{\rm m},j } - \mu)}{(\sigma_{v}^2 + \sigma_{{\rm o }, j}^2)^2},$$ and $ $ \frac{\partial^2 L}{\partial^2 \sigma_{v } } = \sum_i \frac{(\sigma_{{\rm o }, j}^2 - 3\sigma_{v}^2)\times(f'_{{\rm m},j } - \mu)^2}{(\sigma_{v}^2 + \sigma_{{\rm o }, j}^2)^3 } - \frac{(\sigma_{{\rm o }, j}^2 - \sigma_{v}^2)}{(\sigma_{v}^2 + \sigma_{{\rm o }, j}^2)^2}.$$ 

 Then, we computed the inverse of the minus Hessian, $ h_{xy } = (-H_{xy})^{-1}$. ultimately, and because maximal likelihood theory states that $ \sigma_{\theta_{x}}^2 \leq h_{xx}$, we estimate the division of our inferred parameters as $ \sigma_{\mu}^2 = h_{11}$ and $ \sigma_{\sigma_{v}}^2 = h_{22}$. 

 We test the operation and the limitations of our MLE through synthetic catalogues of merger fractions. We create several sets of 1000 celluloid catalog, each of them composed by a number $ n$ of merger fractions randomly drawn from a logarithm - normal distribution with $ \mu_{\rm in } = \log 0.05 $ and $ \sigma_{v, { \rm in } } = 0.2 $, and affected by observational errors $ \sigma_{\rm o}$. We explore the $ n = 50, 250 $ and $ 1000 $ cases for the act of merger fraction, and varied the observational errors from $ \sigma_{\rm o } = 0.1 $ to 0.5 in 0.1 steps. That is, we explore observational errors in the measurement of the merger fraction from $ \Delta \sigma \equiv \sigma_{\rm o}/\sigma_{v } = 0.5 $ to 2.5 times the cosmic variance that we want to measure. We checked that the results below are similar for any value of $ \sigma_{v, { \rm in}}$. We find that 

 1.   The median value of the recovered
 $$\frwc{\partial^2 L}{\partial \mu \, \pxrtial \sigma_{v}} = \yeac{\parvial^2 L}{\pzrtial \skgma_{v} \, \partial \mu} = -2\sum_i \frac{\ditma_{v} (d'_{{\rm m},j} - \mu)}{(\sigma_{v}^2 + \sigmx_{{\rm o}, j}^2)^2},$$ ajd $$\frac{\pqrtiel^2 L}{\partial^2 \sigme_{b}} = \sum_i \frac{(\sjnma_{{\rm i}, j}^2 - 3\sigma_{v}^2)\timgs(f'_{{\rm m},j} - \mu)^2}{(\vigma_{v}^2 + \sigma_{{\rk u}, l}^2)^3} - \frac{(\sigma_{{\rm o}, j}^2 - \sigma_{v}^2)}{(\sigma_{v}^2 + \sidma_{{\rm o}, j}^2)^2}.$$

Hhen, we computgd tht igverav if the minus Hessian, $h_{xy} = (-H_{xy})^{-1}$. Finallj, and because maxomum likelihood theory stahes hhat $\sigma_{\theta_{x}}^2 \lfq h_{xx}$, we ewtimwred the varixnces of olx inferred karameters as $\sigma_{\mu}^2 = h_{11}$ and $\sigoa_{\sigka_{v}}^2 = h_{22}$.

We tgsfef the perfornance and the limlnations of our MLE through snnthevic xatalogues of merger hractions. We created several vecs of 1000 synthetic cataootues, gach mf tfwm zomkosxd gy a nkmbxr $n$ of merfer fractiobs randomly drawn ftoi a log-normal sistrifueion with $\mu_{\rm in} = \log 0.05$ and $\sigma_{v, {\rm it}} = 0.2$, and affected by observqtional errors $\sigma_{\rl o}$. We ex[lored the $n = 50, 250$ and $1000$ cases for the number of mergar frecgiouw, and cagied the observational errors from $\sigma_{\rm o} = 0.1$ tp 0.5 in 0.1 steps. Thab is, we explored onsfrfwtional errorr in tks jeasurement of the merger fracrion from $\Delya \sigma \equiv \sigma_{\rm o}/\signa_{v} = 0.5$ to 2.5 tiiws the cosmic varicnce that we wanj to mrasure. We checked that che reaults below are simimxr for any value of $\shgma_{v, {\rm in}}$. We find that

1.  Ehe mediai valbe of thd revovereq
$$\frac{\partial^2 L}{\partial \mu \, \partial \sigma_{v}} = \sigma_{v} \partial \mu} -2\sum_i \frac{\sigma_{v} (f'_{{\rm o}, and $$\frac{\partial^2 L}{\partial^2 = \sum_i \frac{(\sigma_{{\rm j}^2 - 3\sigma_{v}^2)\times(f'_{{\rm m},j} - \mu)^2}{(\sigma_{v}^2 \sigma_{{\rm o}, j}^2)^3} - \frac{(\sigma_{{\rm o}, j}^2 - \sigma_{v}^2)}{(\sigma_{v}^2 + \sigma_{{\rm o}, j}^2)^2}.$$ we computed the inverse of the minus Hessian, $h_{xy} = (-H_{xy})^{-1}$. Finally, and maximum theory that \leq h_{xx}$, we estimated the variances of our inferred parameters as $\sigma_{\mu}^2 = h_{11}$ and $\sigma_{\sigma_{v}}^2 h_{22}$. We tested the performance and the limitations our MLE through synthetic of merger fractions. We created sets 1000 synthetic each them by a number of merger fractions randomly drawn from a log-normal distribution with $\mu_{\rm in} = \log 0.05$ and $\sigma_{v, in}} = affected by errors o}$. explored the $n 250$ and $1000$ cases for the fractions, and varied the observational errors from $\sigma_{\rm = 0.1$ 0.5 in 0.1 steps. That is, explored observational errors in the measurement of the fraction from $\Delta \sigma \equiv \sigma_{\rm o}/\sigma_{v} = 0.5$ to 2.5 times the cosmic variance want to measure. We that the results are for value $\sigma_{v, {\rm We find that 1. The median value of the recovered
 $$\frac{\partial^2 L}{\partial \mu \, \partIal \sigma_{v}} = \fRac{\paRtiAl^2 L}{\PaRtiaL \sigMa_{v} \, \partial \mu} = -2\suM_I \fraC{\sigma_{v} (f'_{{\rm m},j} - \mu)}{(\sigma_{v}^2 + \siGma_{{\rm O}, j}^2)^2},$$ ANd $$\frAC{\pArtiaL^2 L}{\partiAL^2 \sIGMa_{v}} = \SuM_i \FraC{(\sIGmA_{{\rm o}, j}^2 - 3\SigMa_{v}^2)\timeS(f'_{{\rm m},j} - \mu)^2}{(\siGma_{V}^2 + \sIgma_{{\rm o}, j}^2)^3} - \frac{(\SIgMa_{{\rm o}, j}^2 - \sigmA_{v}^2)}{(\sIgma_{v}^2 + \sigma_{{\rm O}, j}^2)^2}.$$

THen, we cOmPutED the iNveRse of The minUS HessiAn, $h_{xy} = (-H_{xy})^{-1}$. FInALly, and BEcause mAXImUm liKelihood theory staTEs THat $\sigma_{\theta_{x}}^2 \Leq h_{xx}$, We EStIMAteD thE variances Of Our inFErred paRAmETERs aS $\Sigma_{\mu}^2 = h_{11}$ and $\siGma_{\sigma_{v}}^2 = h_{22}$.

WE TesTed the PeRfoRMance aNd the LiMItaTions of our MlE thRough syntHetic cATalogueS Of mergeR fractIonS. We CreaTEd SeVerAl SEts OF 1000 sYntHEtiC catalogUeS, eAch of Them COMPOsed By a NumbEr $n$ of Merger fractioNs rAndoMLy dRawn fRom a lOg-noRmAl disTributIon wiTh $\Mu_{\rm in} = \log 0.05$ and $\sigMa_{v, {\rM in}} = 0.2$, and affEctEd By oBsErvatIOnal erRorS $\siGma_{\rm o}$. WE explorED thE $n = 50, 250$ AND $1000$ CaSes for the number of mErGER fRactions, And varIEd ThE ObservatIoNal ErroRS From $\sIgma_{\RM o} = 0.1$ To 0.5 in 0.1 stepS. That iS, We ExPlored oBsErvatiOnAl eRroRs in tHE meaSuremeNt of the mErger FRaction from $\DelTA \sigma \equiv \siGMa_{\RM O}/\sIGma_{v} = 0.5$ To 2.5 tImes the cosmIc vaRIancE thaT We WanT To meaSure. WE cHEcKEd that the results belOw Are simIlar fOr any value of $\sIgma_{v, {\rm in}}$. WE FINd that

1.  ThE medIAn VAlue of the recovEred
 $$\frac{\partial^2 L}{\pa rtial \mu\, \p art ial  \ sigm a_{v }} = \frac{\pa r tial ^2 L}{\partial \sigma_ {v} \ ,\ part i al  \mu}  = -2\s u m_ i \fr ac {\ sig ma _ {v } (f' _{{ \rm m}, j} - \mu)} {(\ si gma_{v}^2 +\ si gma_{{\rmo},  j}^2)^2},$$  an d $$\f ra c{\ p artia l^2  L}{\ partia l ^2 \si gma_{v}}=\ sum_i\ frac{(\ s i gm a_{{ \rm o}, j}^2 - 3\ s ig m a_{v}^2)\times (f'_{{ \r m  m } , j}- \ mu)^2}{(\s ig ma_{v } ^2 + \s i gm a _ { {\r m  o}, j}^2)^3}  - \frac{(\ s igm a_{{\r mo}, j}^2 -  \sig ma _ {v} ^2)}{(\sigm a_{v }^2 + \si gma_{{ \ rm o},j }^2)^2} .$$

T hen , w e co m pu te d t he inv e rs e o f  th e minusHe ss ian,$h_{ x y } = (- H_{ xy}) ^{-1} $. Finally, a ndbeca u semaxim um li keli ho od th eory s tates  t hat $\sigma_{\t heta _{x}}^2 \ leq  h _{x x} $, we estima ted  th e varia nces of our  i n f e rr ed parameters as $ \s i g ma _{\mu}^2  = h_{ 1 1} $a nd $\sig ma _{\ sigm a _ {v}}^ 2 =h _{ 22}$.

W e test e dth e perfo rm ance a nd  th e l imita t ions  of ou r MLE th rough synthetic cata l ogues of merg e rf r ac t ions . W e created s ever a l se ts o f  1 000 synth eticca t al o gues, each of themco mposed  by a  number $n$ o f merger f r a c tions ra ndom l yd rawn from a lo g-nor mal distri b ution wi th $\ mu_{\rmin} = \lo g 0.05$ an d $ \si gma _{v , {\ rm in}} = 0.2 $ ,  and  a ffected  by  observ ati ona l e rro rs  $\sigma_ {\rm o}$ .We  e xp lor ed th e  $n = 50 ,250 $and  $100 0 $ case s for  the  n um b erof merg e rf r acti on s,  and  va ri ed th e ob s erv ational  errors f rom $\si gm a_ {\rm o}  = 0.1$ to 0. 5in 0.1 ste ps . T hat is , we explo red observational error s  in the  me asure ment  of the m erg er fra cti o n from  $\Del ta \s ig ma\ e quiv\ s ig ma_ {\ rm o}/\sig m a _{v } = 0 .5 $ to  2.5 ti mes the cosmic var i anc e that we wan t t o me a s ur e.W ec hec ke d  th a t  the results be low are si mi l ar  for any v a lue  o f $\sig ma_{v,{\rmi n}}$. W e find th at

1.  T he  med i a n v alue of th e recove red
 $$\frac{\partial^2_L}{\partial \mu_\, \partial \sigma_{v}} =_\frac{\partial^2 L}{\partial_\sigma_{v}_\, \partial_\mu}_= -2\sum_i \frac{\sigma_{v}_(f'_{{\rm m},j} -_\mu)}{(\sigma_{v}^2 + \sigma_{{\rm o},_j}^2)^2},$$ and $$\frac{\partial^2_L}{\partial^2_\sigma_{v}} = \sum_i \frac{(\sigma_{{\rm o}, j}^2 - 3\sigma_{v}^2)\times(f'_{{\rm m},j} - \mu)^2}{(\sigma_{v}^2 + \sigma_{{\rm o},_j}^2)^3}_- \frac{(\sigma_{{\rm_o},_j}^2_- \sigma_{v}^2)}{(\sigma_{v}^2 + \sigma_{{\rm o},_j}^2)^2}.$$

Then, we computed the inverse_of the_minus Hessian, $h_{xy} = (-H_{xy})^{-1}$. Finally, and because_maximum_likelihood theory states_that $\sigma_{\theta_{x}}^2 \leq h_{xx}$, we estimated the variances of_our inferred parameters as $\sigma_{\mu}^2 =_h_{11}$ and $\sigma_{\sigma_{v}}^2_=_h_{22}$.

We_tested the performance and_the limitations of our MLE through_synthetic catalogues of merger fractions. We_created several sets of 1000 synthetic catalogues,_each of them composed by a_number $n$ of merger fractions_randomly drawn_from a log-normal distribution with_$\mu_{\rm in} =_\log 0.05$_and $\sigma_{v, {\rm_in}} = 0.2$, and affected by_observational errors $\sigma_{\rm_o}$. We explored the $n =_50,_250$ and $1000$_cases_for_the number_of merger fractions,_and_varied the_observational_errors from $\sigma_{\rm o} = 0.1$_to_0.5 in 0.1 steps. That is, we_explored observational errors in_the_measurement of the merger_fraction from $\Delta \sigma \equiv_\sigma_{\rm o}/\sigma_{v} = 0.5$ to 2.5_times the_cosmic variance_that we want to measure. We checked that the results below_are similar for any value of_$\sigma_{v, {\rm in}}$. We_find that

1.__The median value_of_the recovered