 the Fisher information matrix of the model $0.25\mathcal{N}(-10,1)+0.10\mathcal{N}(0,5)+0.65\mathcal{N}(15,7)$ for sample sizes from $500$ to $1700$. The value obtained via numerical integration is represented by the red line (in the graph below, all the approximations obtained with more than $550$ knots give the same result, exactly equal to the one obtained via Gauss-Kronrod quadrature).[]{data-label="fig:MCvsRIEMbxp"}](MCvsRIEM-boxpl1)

![The first element on the diagonal of the Fisher information matrix relative to the first weight of the two-component Gaussian mixture model $0.5 \mathcal{N}(-1,0.01)+0.5 \mathcal{N}(2,0.01)$.[]{data-label="fig:FishInfoelem"}](elem)

Thus, we have decided to use the sums of Riemann (with a number of points equal to $550$) to approximate the Jeffreys prior when the standard deviations are sufficiently big and Monte Carlo integration (with sample sizes of $1500$) when they are too small. In this case, the variability of the results seems to decrease as $\sigma_i$ approaches $0$, as shown in Figure \[fig:MCsmallsd\].

![Approximation of the Jeffreys prior (in log-scale) for the two-component Gaussian mixture model $0.5 \mathcal{N}(-1,\sigma)+0.5\mathcal{N}(2,\sigma)$, where $\sigma$ is taken equal for both components and decreasing.[]{data-label="fig:MCsmallsd"}](MCsmallsigma)

We have chosen to consider Monte Carlo samples of size equal to $1500$ because both the value of the approximation and its standard deviations are stabilizing.

An adaptive MCMC algorithm has been used to define the variability of the kernel density functions used to propose the moves. During the burnin, the variability of the kernel distributions has been reduced or increased depending on the acceptance rate, in a way such that the acceptance rate stay between $20\%$ and $40\%$. The transitional kernel used have been truncated normals for the weights, normals for the means and log-normals for the standard deviations (all centered on the values accepted in the previous iteration).

Conclusion {#sec:concl}
==========

the Fisher information matrix of the model $ 0.25\mathcal{N}(-10,1)+0.10\mathcal{N}(0,5)+0.65\mathcal{N}(15,7)$ for sample sizes from $ 500 $ to $ 1700$. The value prevail via numeric integration is represented by the red cable (in the graph below, all the approximations obtain with more than $ 550 $ knots give the like resultant role, exactly adequate to the one obtained via Gauss - Kronrod quadrature).[]{data - label="fig: MCvsRIEMbxp"}](MCvsRIEM - boxpl1) 

! [ The first element on the diagonal of the Fisher data matrix relative to the first system of weights of the two - component Gaussian mixture exemplar $ 0.5 \mathcal{N}(-1,0.01)+0.5 \mathcal{N}(2,0.01)$.[]{data - label="fig: FishInfoelem"}](elem) 

 Thus, we have decided to use the kernel of Riemann (with a number of points equal to $ 550 $) to approximate the Jeffreys prior when the standard deviation are sufficiently big and Monte Carlo integration (with sample sizes of $ 1500 $) when they are too belittled. In this case, the variability of the results seems to decrease as $ \sigma_i$ approaches $ 0 $, as shown in Figure \[fig: MCsmallsd\ ]. 

! [ Approximation of the Jeffreys prior (in log - scale) for the two - component Gaussian mixture model $ 0.5 \mathcal{N}(-1,\sigma)+0.5\mathcal{N}(2,\sigma)$, where $ \sigma$ is taken equal for both components and decreasing.[]{data - label="fig: MCsmallsd"}](MCsmallsigma) 

 We have chosen to consider Monte Carlo sample of size equal to $ 1500 $ because both the value of the approximation and its standard deviation are stabilize. 

 An adaptive MCMC algorithm has been used to define the variability of the kernel density functions used to propose the move. During the burnin, the variability of the kernel distributions has been reduced or increased count on the acceptance rate, in a way such that the acceptance pace stay between $ 20\%$ and $ 40\%$. The transitional kernel use have been truncate normals for the weights, normals for the means and log - convention for the standard deviations (all centered on the values accepted in the previous iteration). 

 Conclusion { # sec: concl } 
 = = = = = = = = = =
 thf Fisher information matvix of the model $0.25\mathcan{N}(-10,1)+0.10\mathdal{N}(0,5)+0.65\mathzal{N}(15,7)$ for sample sizes from $500$ vo $1700$. Rhe vqlue obtained via numefical intvgration us rtpresented by the red linc (in fme grcpi below, all the approximadions obtained wktk more than $550$ knots give the same restlt, exavtpy equal to thg one jbtajned via Gauss-Kronrod quadrature).[]{dafa-label="hig:MCvsRIEMbxp"}](MVvsRIEM-boxpl1)

![The first elemfnt ln the diagonal of the Fisher infjemation matrkx relativt co the firsj weight of the two-component Gaursian mixture midwl $0.5 \kathcal{N}(-1,0.01)+0.5 \mavhcal{N}(2,0.01)$.[]{qata-label="fig:NpshInfoalem"}](elek)

Thus, we have cecmded to use the sums of Rmemann (with a number of pointv zqual to $550$) to approximqtw the Jefxreyr prkor wien the shanvard deviatjons are sudficiently big and Kognr Carlo intefratiog (rith sample sizes of $1500$) when they are too smzll. In this case, the vaeiability of the resupts seems to decrease as $\sigma_i$ approaches $0$, as shown in Ficure \[hie:MCwmwuosf\].

![Approximation of the Jeffreys prior (in log-scwme) fpr the two-compjnent Gaussoaj kyxture model $0.5 \mathccm{N}(-1,\aigma)+0.5\mathcal{N}(2,\sigma)$, where $\figma$ is taken equsl for both components and eecreasing.[]{dana-lavel="fig:MCsmallsd"}](MCslallsigma)

We havg chosrn to consider Monte Caxlo sajples of sixe equal fu $1500$ because both ghe vdlue of ufe approximation wnd its svandaxd deviagionx are ftabilizinh.

An adaptive MCMC algorihhm hcs bean used to define the variability of the kernel density fgncnions useb to pvopose the movef. During the bornin, the varixbility of the kecnel distribttions has beaj reduced or increasqd dwpeneing on ghe acceptance rate, in c way such that the acccptanzs rate stay betcten $20\%$ and $40\%$. The tranxitkonwl kxrnel gsed have bean tfunzsted vormals for ghe eeights, normals for dhe jeans and log-normaks for the standarq deviations (sll centered on tht valuxs accxpted on jhe previous iteration).

Conclusioh {#sec:concp}
==========

the Fisher information matrix of the model sample from $500$ $1700$. The value represented the red line the graph below, the approximations obtained with more than knots give the same result, exactly equal to the one obtained via Gauss-Kronrod ![The first element on the diagonal of the Fisher information matrix relative to first of two-component mixture model $0.5 \mathcal{N}(-1,0.01)+0.5 \mathcal{N}(2,0.01)$.[]{data-label="fig:FishInfoelem"}](elem) Thus, we have decided to use the sums of Riemann (with number of points equal to $550$) to approximate Jeffreys prior when the deviations are sufficiently big and Carlo (with sample of when are too small. this case, the variability of the results seems to decrease as $\sigma_i$ approaches $0$, as shown in \[fig:MCsmallsd\]. ![Approximation Jeffreys prior log-scale) the Gaussian mixture model where $\sigma$ is taken equal for decreasing.[]{data-label="fig:MCsmallsd"}](MCsmallsigma) We have chosen to consider Monte Carlo of size to $1500$ because both the value the approximation and its standard deviations are stabilizing. adaptive MCMC algorithm has been used to define the variability of the kernel density functions propose the moves. During burnin, the variability the distributions been or increased on the acceptance rate, in a way such that the acceptance stay between $20\%$ and $40\%$. The transitional kernel used have normals the weights, normals the means and log-normals the deviations (all centered on accepted the {#sec:concl}
 the Fisher information matriX of the modeL $0.25\mathCal{n}(-10,1)+0.10\maThCal{N}(0,5)+0.65\MathCal{N}(15,7)$ for sample sIZes fRom $500$ to $1700$. The value obtained vIa numErICal iNTeGratiOn is repREsENTed By ThE reD lINe (In the GraPh below, All the apprOxiMaTions obtaineD WiTh more than $550$ KnoTs give the samE reSult, exAcTly EQual tO thE one oBtaineD Via GauSs-Kronrod QuADraturE).[]{Data-labEL="FiG:MCvSRIEMbxp"}](MCvsRIEM-bOXpL1)

![the first elemenT on the DiAGoNAL of The fisher infoRmAtion MAtrix reLAtIVE To tHE first weight oF the two-compONenT GaussIaN miXTure moDel $0.5 \maThCAl{N}(-1,0.01)+0.5 \Mathcal{N}(2,0.01)$.[]{datA-labEl="fig:FishinfoelEM"}](elem)

ThUS, we have DecideD to Use The sUMs Of rieMaNN (wiTH a NumBEr oF points eQuAl To $550$) to aPproXIMATe thE JeFfreYs priOr when the stanDarD devIAtiOns arE suffIcieNtLy big And MonTe CarLo Integration (with SampLe sizes of $1500$) WheN tHey ArE too sMAll. In tHis CasE, the varIabilitY Of tHe RESUlTs seems to decrease aS $\sIGMa_I$ approacHes $0$, as sHOwN iN figure \[fiG:McsmAllsD\].

![aPproxImatIOn Of the JefFreys pRIoR (iN log-scaLe) For the TwO-coMpoNent GAUssiAn mixtUre model $0.5 \MathcAL{N}(-1,\sigma)+0.5\mathcal{n}(2,\Sigma)$, where $\sigMA$ iS TAkEN equAl fOr both compoNentS And dEcreASiNg.[]{dATa-labEl="fig:mCSMaLLsd"}](MCsmallsigma)

We haVe Chosen To conSider Monte CarLo samples oF SIZe equal tO $1500$ becAUsE Both the value of The apProximatioN And its stAndarD deviatiOns are staBILizing.

An AdaPtiVe McMC ALGoRithm has been uSED to dEfIne the vAriAbility Of tHe kErnEl dEnSity functIons used To PrOpOsE thE moveS. during thE bUrnIn, The VariaBIlity oF the kErneL dIsTRibUtions hAS bEEN redUcEd Or inCreAsEd depEndiNG on The accePtance ratE, in A Way sUcH tHat the aCceptance rate StAy between $20\%$ aNd $40\%$. the TransiTIOnal kernEl used have been truncated NOrmals fOr tHe weiGhts, Normals foR thE means And LOg-normAls for The stAnDarD DEviatIONs (All CeNtered on thE VAluEs accEpTed iN the preVious iteration).

ConcLUsiOn {#sec:concl}
==========

 the Fisher information ma trix of th e mod el$0. 25 \mat hcal {N}(-10,1)+0.1 0 \mat hcal{N}(0,5)+0.65\math cal{N }( 1 5,7) $  f or sa mple si z es f rom  $ 50 0$to $1 700$.  Th e value  obtainedvia  n umerical int e gr ation is r epr esented by t hered li ne  (i n  thegra ph be low, a l l theapproxima ti o ns obt a ined wi t h  m orethan $550$ knotsg iv e  the same resu lt, ex ac t ly e qua l t o the oneob taine d  via Ga u ss - K r onr o d quadrature) .[]{data-la b el= "fig:M Cv sRI E Mbxp"} ](MCv sR I EM- boxpl1)

![ Thefirst ele ment o n  the di a gonal o f theFis her  inf o rm at ion  m a tri x  r ela t ive  to thefi rs t wei ghto f t he t wo- comp onent  Gaussian mix tur e mo d el$0.5\math cal{ N} (-1,0 .01)+0 .5 \m at hcal{N}(2,0.01) $.[] {data-lab el= "f ig: Fi shInf o elem"} ](e lem )

Thus , we ha v e d ec i d e dto use the sums of  R i e ma nn (with  a num b er  o f  pointseq ual  to$ 5 50$)to a p pr oximatethe Je f fr ey s prior  w hen th esta nda rd de v iati ons ar e suffic ientl y  big and Monte Carlo integra t io n (w i th s amp le sizes of  $15 0 0$)when th eya re to o sma ll .  I n  this case, the var ia bility  of t he results se ems to dec r e a se as $\ sigm a _i $  approaches $0 $, as  shown inF igure \[ fig:M Csmallsd \].

![Ap p r oximatio n o f t heJef f r ey s prior (in l o g -sca le ) for t hetwo-com pon ent  Ga uss ia n mixture  model $ 0. 5\m at hca l{N}( - 1,\sigma )+ 0.5 \m ath cal{N } (2,\si gma)$ , wh er e$ \si gma$ is ta k e n eq ua lforbot hcompo nent s  an d decre asing.[]{ dat a -lab el =" fig:MCs mallsd"}](MCs ma llsigma)
We  ha ve cho s e n to con sider Monte Carlo sampl e s of si zeequal  to$1500$ be cau se bot h t h e valu e of t he ap pr oxi m a tiona n dits  s tandard de v i ati ons a re  sta bilizin g.

An adaptive MC M C a lgorithm hasbee n us e d  t o d e fi n e t he var i a bility of the k ernel dens it y  f unctions u s edto  propos e the m oves. Duringthe burni n, the va ri abil i t y o f the kern el distr ibutionsh as be e nreduc edor inc re ase d dep ending onthe a ccepta nc e rate , inaway such  that the acceptance ra te sta y bet wee n $20\%$and $40 \%$. Thetran sitional k ern elusedhav e  been  tru n ca ted norma ls f o r the wei g ht s,n o rm als for the m e ans  andlog - normal s fo r the standard de v iations (all c ente r e d o n t h e va lu es accepted in  th ep r evious i te ration).

C onclusio n{ #sec: concl} 
===== =====

 the_Fisher information_matrix of the model_$0.25\mathcal{N}(-10,1)+0.10\mathcal{N}(0,5)+0.65\mathcal{N}(15,7)$ for_sample_sizes from_$500$_to $1700$. The_value obtained via_numerical integration is represented_by the red_line_(in the graph below, all the approximations obtained with more than $550$ knots give_the_same result,_exactly_equal_to the one obtained via_Gauss-Kronrod quadrature).[]{data-label="fig:MCvsRIEMbxp"}](MCvsRIEM-boxpl1)

![The first element on_the diagonal_of the Fisher information matrix relative to the_first_weight of the_two-component Gaussian mixture model $0.5 \mathcal{N}(-1,0.01)+0.5 \mathcal{N}(2,0.01)$.[]{data-label="fig:FishInfoelem"}](elem)

Thus, we have_decided to use the sums of_Riemann (with a_number_of_points equal to $550$)_to approximate the Jeffreys prior when_the standard deviations are sufficiently big_and Monte Carlo integration (with sample sizes_of $1500$) when they are too_small. In this case, the_variability of_the results seems to decrease_as $\sigma_i$ approaches_$0$, as_shown in Figure_\[fig:MCsmallsd\].

![Approximation of the Jeffreys prior (in_log-scale) for the_two-component Gaussian mixture model $0.5 \mathcal{N}(-1,\sigma)+0.5\mathcal{N}(2,\sigma)$,_where_$\sigma$ is taken_equal_for_both components_and decreasing.[]{data-label="fig:MCsmallsd"}](MCsmallsigma)

We have_chosen_to consider_Monte_Carlo samples of size equal to_$1500$_because both the value of the approximation_and its standard deviations_are_stabilizing.

An adaptive MCMC algorithm_has been used to define_the variability of the kernel density_functions used_to propose_the moves. During the burnin, the variability of the kernel distributions_has been reduced or increased depending_on the acceptance rate,_in a_way_such that the_acceptance_rate stay_between $20\%$ and $40\%$. The transitional kernel_used have_been truncated normals for the weights,_normals for the means_and_log-normals for the standard deviations (all_centered on the values accepted in_the previous iteration).

Conclusion {#sec:concl}
==========
