通常错误是初始设置时矩阵是非正定的，比如在OLS估计时矩阵未满秩就会报这样的错误。如模型两个或多个变量是共线性。可以尝试分段运行，看哪里报错再确定怎么修改

我感觉是你自己的数据在变成矩阵运算时，矩阵不是正定矩阵了。这个应该是自己数据本身的问题

不知道 Matrix not positive提示的地方是在哪？是在运算过程中还是计算参数标准差的地方？  无论出现在哪，跟数据没有关系，采用极大似然类的参数估计方法经常会遇到这个问题，原因在于 给出的初始值不太理想。可以试着换换不同的初始值。

建议你把数据、代码放上来，以便群友们用 Gauss 跑一下、才能知道你哪里出问题了；你现在这样，别人怎么帮你呢？