'YHhddlZddlmZddlmZddZddZedk(r gdZ gdZ d Z d Z e rd D]Z ddlZejZed e \ZZZededeeede de edejez edeeeedeeedeeej)dyyy)N)arma_generate_sample)Armac |gd}|gd}g}g}t|D] }t|||dz|| d}||jz}t|} | j d\} } } } }|j | | j | }tjtj|||z }| ;|j |tjtj| z|j tjtj| ztj|dd|ddftj|tj|fS)zrun Monte Carlo for ARMA(2,2) DGP parameters currently hard coded also sample size `nsample` was not a self contained function, used instances from outer scope now corrected N?ggrg333333?g?)r )rangermeanrfitappend geterrorsnpsqrtdotdiagnan zeros_liker_array)niternsamplearmasigresults results_bse_y2arest2rhohat2acov_x2ainfodictmesgiererr2asige2as \/var/www/html/planif/env/lib/python3.12/site-packages/statsmodels/sandbox/mcevaluate/arma.pymcarma22r+sF z  z GK 5\A !"R c :G89 E bggib17E1B.'8T3x   *e,W45     v0@(AA B   rvv h(?? @A 55ABAB "((7"3RXXk5J JJc0|"tj|jd}tj|j d}t d|j ||}t dt tj||z dzj dt dt ||z j dt dt tj||z dt dt tj||z |z d zdt d t tjtj||z dt d t ||kDj dy) Nr z!fractions of iterations with nansRMSEr rz mean biasz median biaszmedian bias percentdzmedian absolute errorzpositive error fraction) rzerosshapeisnananyprintr rmedianabs)resrtnanrowss r* mc_summaryr:,s& z XXciil #hhsm"G -w||~> wh-C &M "''CFQ;$$Q' () + 3r6--  - "))SVQ     "))SVRKOA &' !" "))BFF3r6N1 %& #$ 38//! r,__main__rrT)g??rr/)rrz Results for Monte Carlotruez nsample =zsigma = zelapsed time for Monte Carloz MC of rho versus truez MC of bse versus zeroz MC of bse versus std) r NNr=)N)numpyrstatsmodels.tsa.arima_processrstatsmodels.tsa.arma_mlerr+r:__name__rrrrun_mcrtimet0r8res_rhores_bser4stdr,r*rKs>) "KH( z.` B BGF ! 0C B#+#3#? B - . &M "I +w C 8 0)$))+b. A + , w # + , w  * + w A /% 0r,