@@ -1111,24 +1111,30 @@ def test_polyagamma_moment(h, z, size, expected):
1111
1111
(np .array ([0.25 , 0.25 , 0.25 , 0.25 ]), 1 , None , np .array ([1 , 0 , 0 , 0 ])),
1112
1112
(np .array ([0.3 , 0.6 , 0.05 , 0.05 ]), 2 , None , np .array ([1 , 1 , 0 , 0 ])),
1113
1113
(np .array ([0.3 , 0.6 , 0.05 , 0.05 ]), 10 , None , np .array ([4 , 6 , 0 , 0 ])),
1114
- (np .array ([[0.3 , 0.6 , 0.05 , 0.05 ], [0.25 , 0.25 , 0.25 , 0.25 ]]), 10 , None , np .array ([[4 , 6 , 0 , 0 ], [4 , 2 , 2 , 2 ]])),
1115
1114
(
1116
- np .array ([[0.25 , 0.25 , 0.25 , 0.25 ], [0.26 , 0.26 , 0.26 , 0.22 ]]),
1117
- np . array ([ 1 , 10 ]),
1118
- None ,
1119
- np .array ([[1 , 0 , 0 , 0 ], [2 , 3 , 3 , 2 ]])
1115
+ np .array ([[0.3 , 0.6 , 0.05 , 0.05 ], [0.25 , 0.25 , 0.25 , 0.25 ]]),
1116
+ 10 ,
1117
+ None ,
1118
+ np .array ([[4 , 6 , 0 , 0 ], [4 , 2 , 2 , 2 ]]),
1120
1119
),
1121
1120
(
1122
- np .array ([0. 26 , 0.26 , 0.26 , 0.22 ]),
1123
- np .array ([1 , 10 ]),
1124
- None ,
1125
- np .array ([[1 , 0 , 0 , 0 ], [2 , 3 , 3 , 2 ]])
1121
+ np .array ([[ 0.25 , 0.25 , 0.25 , 0.25 ], [ 0. 26 , 0.26 , 0.26 , 0.22 ]]),
1122
+ np .array ([1 , 10 ]),
1123
+ None ,
1124
+ np .array ([[1 , 0 , 0 , 0 ], [2 , 3 , 3 , 2 ]]),
1126
1125
),
1127
1126
(
1128
- np .array ([[0.25 , 0.25 , 0.25 , 0.25 ], [0.26 , 0.26 , 0.26 , 0.22 ]]),
1129
- np .array ([1 , 10 ]),
1130
- 2 ,
1131
- np .full ((2 , 2 , 4 ), [[1 , 0 , 0 , 0 ], [2 , 3 , 3 , 2 ]])),
1127
+ np .array ([0.26 , 0.26 , 0.26 , 0.22 ]),
1128
+ np .array ([1 , 10 ]),
1129
+ None ,
1130
+ np .array ([[1 , 0 , 0 , 0 ], [2 , 3 , 3 , 2 ]]),
1131
+ ),
1132
+ (
1133
+ np .array ([[0.25 , 0.25 , 0.25 , 0.25 ], [0.26 , 0.26 , 0.26 , 0.22 ]]),
1134
+ np .array ([1 , 10 ]),
1135
+ 2 ,
1136
+ np .full ((2 , 2 , 4 ), [[1 , 0 , 0 , 0 ], [2 , 3 , 3 , 2 ]]),
1137
+ ),
1132
1138
],
1133
1139
)
1134
1140
def test_multinomial_moment (p , n , size , expected ):
0 commit comments