2017-08-02 136 views
1

我想使用包msm来做一个基本的马尔可夫模型,并且事情工作正常,直到我突然开始接收下面的错误代码。我不知道为什么它突然开始抛出这个,因为它在早期工作得很好,我不认为我已经改变了任何东西。错误代码似乎是指向线性代数库,但我不知道如何处理它到底做...r msm BLAS/LAPACK例程'DGEBAL'给出了错误代码-3

Error in balance(baP$z, "S") : 
    BLAS/LAPACK routine 'DGEBAL' gave error code -3 

的代码如下:

statesDistMatrix2 <- matrix(c(.1,0,0,.1), nrow = 2, ncol = 2) 

    msm1 <- msm(error ~ stop_datetime, subject = TRIP_ID, data = train_245_mk, 
     qmatrix = statesDistMatrix2, control=list(fnscale=5000,maxit=500)) 

回答

0

从这个文件http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.230.5929&rep=rep1&type=pdf约LAPACK

似乎(部分错误指示灯及警告) “错误或警告由例程检测:INFO 如果INFO = -i,第i个参数有一个非法值”。

还是从这个文档看,你的案例中似乎是问题的第三个参数必须是真实的。有可能在过程的某个时刻出现一些复杂的价值。确保此参数仅具有实际值。有时,某些操作可能会导致如1 + 0.000000001j这样的结果,并且这是有问题的,请尝试通过仅传递参数的实部来调用该函数。

希望这有助于您。