
    !gx$                         d dl Z d dlZe j                  j	                  e j                  j                  e            Z G d d      Z G d d      Z	 G d d      Z
y)    Nc                       e Zd ZdZddZy)ARLagResultsz
    Results are from R vars::VARselect for sunspot data.

    Comands run were

    var_select <- VARselect(SUNACTIVITY, lag.max=16, type=c("const"))
    c                 p    |dk(  r1g d}t        j                  |      j                  ddd      | _        y y )Nconst)@g3;?@gjXYI@gd>X@gp8@g#c(@ga@g479@gsԹq@gqk@gԀ@g޼@g&&]}q@gÎ@g#@g@g,ɔq@ghc'㛕@gNR@g=i@g6q@g	Cy@gT㎧@gjF@g}ؔcq@gf59V@gUo@g瓍@gTp@g&*@g.YX@gon3@gVo@gG@g*[xC"@g?e^o@gkg	n@g!ZͲ@gGҽ\.@g'ˇ-@gQ	Hn@gY@@gops:@g@g^|n@g,.ϟ@gokF@gxY=Ӫ@g];\n@g%p)
@g?(R@g
@g03Ln@gZ@gO/W[@g\&E@g*#a0o@g٥ @gFzlb@gen@gRUo@g%mQ@g1D_g@gu@g.c+o@   F)order)npasarrayreshapeic)selftyper   s      e/var/www/dash_apps/app1/venv/lib/python3.12/site-packages/statsmodels/tsa/tests/results/results_ar.py__init__zARLagResults.__init__   s8    7?%B. jjn,,Q#,>DG1     N)r   __name__
__module____qualname____doc__r    r   r   r   r      s    ?r   r   c                       e Zd ZdZddZy)ARResultsOLSzy
    Results of fitting an AR(9) model to the sunspot data.

    Results were taken from Stata using the var command.
   c                    d| _         |r[g d| _        g d| _        g d| _        d| _        d| _        d| _        t        j                  j                  t        d      }t        j                  |      }|d d	d
f   }|d	d df   }~|| _        || _        |dd  | _        |d d | _        t        j"                  |dd  |d d f      | _        | j$                  | _        t        j"                  |d   |f      | _        t        j"                  ||d d f      | _        || _        t        j"                  |d   |d d f      | _        t        j"                  |d   |d d f      | _        y |sbg d| _        g d| _        g d| _        d| _        d| _        d| _        d| _        t        j                  j                  t        d      }t        j                  |      }|d d	d
f   }|d	d df   }~|| _        || _        |dd  | _        |d d | _        t        j"                  |dd  |d d f      | _        | j$                  | _        t        j"                  |d   |f      | _        t        j"                  ||d d f      | _        || _        t        j"                  |d   |d d f      | _        t        j"                  |d   |d d f      | _        y y )Ng     r@)
gn
@gHm?gޖNB`ٿg\B)QſgG,?gX&J9gzhkt?ggj֩?g(/Cgs29?)
gKN@gＰ?gb+X
c?g%?g7'_"?g1
?gmw	?g%)Bѐ?g78%D?gȍ?)
gvN@g$M.?g,Ŷ?gMu?g$jj?go?g*Pd?g6O[?g?gj"?g3A.@g\}rm@g)<NzAROLSConstantPredict.csv,  r      [         \   r            )	g="?gĦXAٿg$:n=Ŀg_ZF?g3#gϐ6?g{b?g-wlUgM:?)	gXx6}?g ?gyhY?g|eoE?gO]cS?gj!B?gT<?gOԭ?g^I)?)	gno?gQ7?gsϴ?g$ܠ?g?gI`s?g`H?gG?gŀ6	?gm .@gl_l@g*{]g*AuM n@zAROLSNoConstantPredict.csv)avobsparams	bse_stata	bse_gretlrmsefpellfospathjoincur_dirr   loadtxtFVOLSnneg1start0FVOLSnneg1start9FVOLSnneg1start100FVOLSn200start0hstackFVOLSn200start200FVOLSn200startneg109FVOLSn100start325FVOLSn301start9FVOLSdefaultFVOLSn4start312FVOLSn15start312sigma)r   constantfilenamepredictresultsfvpvs         r   r   zARResultsOLS.__init__3   s   
!DKDN&DN )DI(DH)DH ww||G-GHHZZ1Na(Ba(B %'D!$&D!&(jD##%ds8D %'YY56
BvJ/G%HD"(,(>(>D%%'YY2|%<D"#%99b"Ra&\#:D  "D#%99bfb!f-=#>D $&IIr"vr#2w.?$@D!"DKDNDN *DI*DJ(DH
 )DHww||G-IJHZZ1Na(Ba(B %'D!$&D!&(jD##%ds8D %'YY56
BvJ/G%HD"(,(>(>D%%'YY2|%<D"#%99b"Ra&\#:D  "D#%99bfb!f-=#>D $&IIr"vr#2w.?$@D!i r   NTr   r   r   r   r   r   -   s    
wAr   r   c                       e Zd ZdZddZy)ARResultsMLEzt
    Results of fitting an AR(9) model to the sunspot data using exact MLE.

    Results were taken from gretl.
    c                 R   d| _         |rt        j                  j                  t        d      }t        j                  j                  t        d      }t        j                  |d      }|d d df   }t        j                  |dd      }|d d | _        |d	d | _	        |d
d | _
        |d d | _        |dd  | _        |dd  | _        |d	d | _        |d d | _        |dd | _        |dd | _        |d d df   | _        |d d df   | _        |d d df   | _        |d d df   | _        y y )Nr   zARMLEConstantPredict.csvz#results_ar_forecast_mle_dynamic.csv,)	delimiterr   )rI   skip_headeri5  	   d         i4  i6  i.  r   i9  r#   r$   r      )r&   r-   r.   r/   r0   r   r1   
genfromtxtFVMLEdefaultFVMLEstart9end308FVMLEstart100end308FVMLEstart0end200FVMLEstart200end334FVMLEstart308end334FVMLEstart9end309FVMLEstart0end301FVMLEstart4end312FVMLEstart2end7fcdynfcdyn2fcdyn3fcdyn4)r   r?   r@   	filename2rA   rC   	dynamicpvs          r   r   zARResultsMLE.__init__   s:   
 ww||G-GHHW%JLIZZC@N1%Bi3ANI !#4CD%'#YD"')#c{D$%'XD"')#$xD$')#$xD$%'#YD"%'XD"%'#YD"#%a7D "1a4DJ#AqD/DK#AqD/DK#AqD/DK r   NrD   r   r   r   r   rF   rF      s    
(r   rF   )r-   numpyr   r.   dirnameabspath__file__r0   r   r   rF   r   r   r   <module>re      sN    	 
''//"''//(3
4"? "?J}A }A@. .r   