M/Ph# f 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 ]-Z
ej e j dde
f e j dde
f< .d Zd Zd Zd Zd Zd Zd Zd Zd ZdS ) 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 }| | }| }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 dS )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 t/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/statsmodels/multivariate/tests/test_multivariate_ols.pycompare_r_output_dogs_datarB s
'M C v A A ? ? ? ? ? ? ? ? ? ? ? ? AA a nV 4 ;QJJJJ L L L L L L L L L L L L NA ai / 61EEEE 5 5 5 5 5 5 5 5 5 5 5 5 7A a
mF 3 :AqIIII = = = = = = = = = = = = ?A a 0 8 ?ANNNNNN c F t d t d d S )Nsvdr" pinv)rB rC rA test_glm_dogs_examplerH G s* e,,,, f------rC c t j dt } | }| ddgd gg }g dg dg dg dg}t |d d j |d g d
}g d }| d ||gg }g d
g dg dg}t |d d |d g dg dg dg}t |d d j |d d S )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 sO
'M C A
k]D A B C CB ? ? ? ? ? ? ? ? ? ? ? ? AA bof 5 Histamine0 + Histamine1 + Histamine3 + Histamine5 ~ Drug * dup)r: copyr r9 r
ValueErrorr; )data1r> s rA "test_independent_variable_singularr] d sz IIKKE=E%L
'H
C *cg&&&
'H
C *sw'''''rC c t j dt } | d }| }t j dt d \ }}t j g dg }t || d }| d|d gg }t |d d j
|d d j
d
t |j
|j
d }| d|d gg }t |d d j
|d d j
d
t j g d g d
g }| d|d gg }t || d }| d|d gg }t |d d j
|d d j
d
t |j
|j
d }| d|d gg }t |d d j
|d d j
d
d S )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 s
'M C u A
B/M +' ' 'KE4 $$$%&&A%%)))77A
Q 5 6 7 7B bof 5 < of 5 r? rU rV s rA test_L_M_matrices_1D_arrayrj s
'M C u A
###$$A*aifa5F4GHHHH
$$$%&&A
###$$A*aifa^4DEEEEEErC c t j dt } | d }| }g dg dg dg dg}t |d d j |d
d S ) 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 s
'J C u A
B 3 3 3 3 3 3 5 5 5 5 5 5 7A bnV 4 ;QJJJJJJrC c R t t t j dt d S )NzHistamine0 ~ 0 + Depleted)r r[ r r9 r: rG rC rA test_endog_1D_arrayro s+ * . ;#T+ + + + +rC c t j dt } | d }t j g dg dg }d }t j g dg dg }| d|||fg }g d
g d g d g d
g}t |d d j |d | ddd d S )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 1 1A A
,,, ,, ! !A
Aq 1 2 3 3B 2 2 2 2 2 2 3 3 3 3 3 3 5A bk& 1 8!QGGGGJJtd # % % % % %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 F F F F F B B B B B B B B r|::::::;;;::::::999999999::::::AAA>>>>>>>>>>>> @ 9 9 9 : : :&
q! . .A bfTYqqq!t_--DIaaadOO&O &O &OR. . .
G G G0
(
(
(!D !D !DF F F F
K
K
K+ + +% % % % %rC