'YHh# x d dl Zd dlZd dlmZ d dlmZmZ d dl Z ej g dg dg dg dg dg d g d
g d g d g d
g dg dg dg dg dgg d Z e dd D ]6 Z
ej e j dde
f e j dde
f< 8 d Zd Zd Zd Zd Zd Zd Zd Zd Zy) N)_MultivariateOLS)assert_array_almost_equal
assert_raises)MorphineNg{Gz?g?皙?{Gz?)r r {Gz? Q?r
r
)r r Q ?ffffff?g Q? Q?)r r g(\?g=
ףp=?gffffff?r )r Yr
ףp=
?p=
ף?gQ ?)placebor r r r r )r r 皙?r r r )r r Q?gףp=
?gףp=
?g)\(?)r r r g?g\(\?g333333?)r r r g(\?g Q ?g?) Trimethaphanr r g
ףp=
@g{Gz @gGz?)r r r r r r )r r r r r r )r r r r Q?r )r r r r r r )DrugDepleted
Histamine0
Histamine1
Histamine3
Histamine5)columns c t j dt }|j | }|j }g dg dg dg dg}t |d d j |d g d
g d g d g d
g}t |d d j |d g dg dg dg dg}t |d d j |d g dg dg dg dg}t |d d j |d y)a Testing within-subject effect interact with 2 between-subject effect
Compares with R car library Anova(, type=3) output
Note: The test statistis Phillai, Wilks, Hotelling-Lawley
and Roy are the same as R output but the approximate F and degree
of freedoms can be different. This is due to the fact that this
implementation is based on SAS formula [1]
.. [*] https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_introreg_sect012.htm
CHistamine0 + Histamine1 + Histamine3 + Histamine5 ~ Drug * Depletedmethodg~c? r E+K@1$~?g#?r% r r& r' gˏP B@r% r r& r' Interceptstatr decimal)gTD~? (@g!i
@gstaŕ?)gI"?r. ,@g勤zf@g^Ȯ?)gJ @r. :^ܡ@gXV@g L? @r/ g @g0GĶ?)ghw?r7 r0 g_X0@gro?)gݭ
@r7 r1 gYua?@g6E6F?)g>F@r3 r2 g{|}J@g r헑r?z
Drug:DepletedNr from_formuladatafitmv_testr values)r# modras m/var/www/html/planif/env/lib/python3.12/site-packages/statsmodels/multivariate/tests/test_multivariate_ols.pycompare_r_output_dogs_datarB s
'
'M C v A A ? ? ? ? AA a nV 4 ; ;QJ L L L L NA ai / 6 61E 5 5 5 5 7A a
mF 3 : :AqI = = = = ?A a 0 8 ? ?AN c 4 t d t d y )Nsvdr" pinv)rB rC rA test_glm_dogs_examplerH G s e, f-rC c t j dt } | j }|j ddgd gg }g dg dg dg dg}t |d d j |d g d
}g d }|j d ||gg }g d
g dg dg}t |d d |d g dg dg dg}t |d d j |d y )Nr! r*
hypothesesr$ r( r) r+ r r, )r* zDrug[T.Trimethaphan]zDrug[T.placebo])r r r r@ r r r r r r rM r r r r r r rM r r r
contrast_L
)r rM r r )r r rM r )r r r rM transform_M)r r9 r: r; r< r r= T)r> r? r1r@ LMs rA test_specify_L_M_by_stringrW L s
'
'M C A
k]D A B CB ? ? ? ? AA bof 5 < Histamine0 + Histamine1 + Histamine3 + Histamine5 ~ Drug * dup)r: copyr r9 r
ValueErrorr; )data1r> s rA "test_independent_variable_singularr] d sf IIKE=E%L
'
'H
C *cgg&
'
'H
C *sww'rC c t j dt } | j d }|j }t j dt d \ }}t j g dg }t || j d }|j d|d gg }t |d d j |d d j d
t |j |j j d }|j d|d gg }t |d d j |d d j d
t j g d g d
g }|j d|d gg }t || j d }|j d|d gg }t |d d j |d d j d
t |j |j j d }|j d|d gg }t |d d j |d d j d
y )Nr! rE r" dataframe) return_typerL r* rJ r+ r r, rN rO r ) r r9 r: r; r< patsy dmatricesnparrayr r= )r> r? r0endogexogrU rT s rA test_from_formula_vs_no_formularh q s4
'
'M C u A
B//M +'KE4 $%&A%)))7A
Q 5 6 7B bof 5 < < of 5 < r? rU rV s rA test_L_M_matrices_1D_arrayrj s
'
'M C u A
#$A*aiifa5F4GH
$%&A
#$A*aiifa^4DErC c t j dt } | j d }|j }g dg dg dg dg}t |d d j |d
y ) Nz@Histamine0 + Histamine1 + Histamine3 + Histamine5 ~ 0 + DepletedrE r" )gŏ1w-!_?r7 g 4@g0*K@ )gB>٬?r7 g 6@gŏ1w-a:@rl )g{xX@r7 g[Ӽ](@gN@G]@rl )g $(~QW@r3 g &@g^)p@rl r r+ r% r, r8 )r> r? re r@ s rA test_exog_1D_arrayrm se
'
'J C u A
B 3 3 5 5 7A bnV 4 ; ;QJrC c L t t t j dt y )NzHistamine0 ~ 0 + Depleted)r r[ r r9 r: rG rC rA test_endog_1D_arrayro s * . ; ;#T+rC c z t j dt } | j d }t j
g dg dg }d }t j
g dg dg }|j
d|||fg }g d
g d g d g d
g}t |d d j |d |j ddd y )Nr! rE r" )r g333333?g?g?g ?r
)r g @g @gffffff
@g @g@)rM r r% ) r r. test1rJ )g q?r7 r/ gjM @gMbP?)g?r7 r0 gzv@gMb?)gjqD3@r7 gI&@gJ ++@g-C 6Z?)gFX2@r3 r2 gǺ
@@g-C 6?r+ r% r, T)show_contrast_Lshow_transform_Mshow_constant_C)
r r9 r: r; rc rd r< r r= summary)r> r? rU rV Cre r@ s rA test_affine_hypothesisrz s
'
'M C u A
..0 1A A
, !A
Aq 1 2 3B 2 2 3 3 5A bk& 1 8 8!QGJJtd # %rC )numpyrc pandaspd)statsmodels.multivariate.multivariate_olsr
numpy.testingr r ra DataFramer: rangeilogilocrB rH rW r] rh rj rm ro rz rG rC rA r s F B r||::;::999::A>>>> @ 9 :&
q! .A bffTYYq!t_-DIIadO.&OR.
G0
(!DF F
K+%rC