M/PhddlZddlmZddlmZddZddZed kr gd Z gd Z d Z d Z e rdD]Z ddlZejZede \ZZZededeeede de edejez edeeeedeeedeeeddSdSdS)N)arma_generate_sample)Arma ?c ,|gd}|gd}g}g}t|D])}t|||dz|| d}||z}t|} | d\} } } } }|| | | }tjtj |||z }| =||tjtj | z|tj tj | z+tj |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?ggr g333333?g?r)r )rangermeanrfitappend geterrorsnpsqrtdotdiagnan zeros_liker_array)niternsamplearmasigresults results_bse_y2arest2rhohat2acov_x2ainfodictmesgiererr2asige2as c/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/statsmodels/sandbox/mcevaluate/arma.pymcarma22r-s z    z   GK 5\\AA !"R c : :G899 E bggiib17E1B1B.'8T3x     **e,,W455     v0@0@(A(AA B B B B   rv h(?(?? @ @ @ @ 5ABBABB "(7"3"3RXk5J5J JJc|tj|jd}tj|d}t d|||}t dt tj||z dzdt dt ||z dt dt tj||z dt dt tj||z |z d zdt d t tjtj ||z dt d t ||kddS) Nr z!fractions of iterations with nansRMSEr rz mean biasz median biaszmedian bias percentdzmedian absolute errorzpositive error fraction) rzerosshapeisnananyprintrrmedianabs)resrtnanrowss r, mc_summaryr<,s z Xcil # #hsmm""G -w||~~>>> wh-C &MMM "'CFQ;$$Q'' ( ())) + 3r6--   - ")SVQ          ")SVRKOA & &''' !""" ")BF3r6NN1 % %&&& #$$$ 38//!  r.__main__r r T)g?rr r1)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)rrNNr)N)numpyrstatsmodels.tsa.arima_processrstatsmodels.tsa.arma_mlerr-r<__name__rrrrun_mcrtimet0r:res_rhores_bser6stdr.r,rKs>>>>>>)))))) "K"K"K"KH( z.`   B   BGF 0! 0 0C KKKB#+8#3#?#?#? B E- . . . E&MMM E"III E+w C 8 8 8 E0)$)++b. A A A E+ , , , Jw # # # E+ , , , Jw    E* + + + Jw A / / / /e~00 0 0r.