Python / Numpy / Scipy，使用分组协变量求解非线性最小二乘？

Question

基本上，我试图让这个功能发生：

我在哪里解决问题。 gamma，alpha和x都来自数据。

最初，我只使用了汇总统计量均值（xi / gamma_i），这意味着可以预先计算出总和中的所有内容，然后我只会向非线性优化器提供一个简单的np数组...但现在没有办法预先计算摘要统计数据，因为当f响应alpha_i时f变化时，beta会如何影响f。因此，我不确定如何展示该阵列。是否有可能将这些协变量作为列表（numpy Objects）嵌入仍然呈现一个numpy数组，然后在剩余函数中解压缩列表？我是以错误的方式解决这个问题吗？