我正在尝试用 Python 重新实现以下代码：http : //www.wordfish.org/uploads/1/2/9/8/12985397/wordfish_1.3.r
但是，我正在努力移植优化（最大化）功能。我制作了一个代码片段并制作了两个应该产生相同结果的简单示例：
回复：
omega = c(0)
y = c(18, 24, 23, 27, 24, 5, 7, 5, 4, 7, 64, 59, 66, 51, 50, 11, 17, 9, 13, 16)
b = c(-0.5047485, -0.6346391, -0.375262, -0.4550382, -0.4839501, 0.3181721, 0.4557756, 0.5618271, 0.3847818, 0.485241, -0.4500023, -0.4224367, -0.461969, -0.5605062, -0.4635604, 0.4995454, 0.5027968, 0.5656253, 0.5506545, 0.5368292 )
psi = c(2.653242, 2.687847, 2.564949, 2.61007, 2.694627, 2.533697, 2.60269, 2.564949, 2.60269, 2.631889, 3.65584, 3.632309, 3.591818, 3.572346, 3.605498, 3.580737, 3.54674, 3.569533, 3.624341, 3.538057)
llik_alpha_1 <- …