我让机器跑了一下，看来20万次不够，我用3条CHAIN，每条CHAIN不同的初值，75万次（每条CHAIN25万次）算出个基本收敛的结果。由于CHAIN数大于1，可以用BGR RATIO来看是否收敛（BGR接近于1表示收敛），算到75万次，基本收敛，结果如下（其实，除了muy以外的参数早就收敛，10万次足矣，但是就是这个muy收敛非常慢，你看了结果会发现其他参数跟你那个结果差不多，唯有muy差别很大，75万次算出来那个MEAN是负数）：
        mean        sd        MC_error        val2.5pc        median        val97.5pc
k        0.3789        0.05451        6.878E-4        0.2785        0.3766        0.4934       
kt        1.472E-4        2.09E-5        2.561E-7        1.09E-4        1.462E-4        1.908E-4
lambda        0.003693        0.002599        3.025E-6        4.482E-4        0.003109        0.01025
mu        9.924E-5        4.769E-4        4.301E-6        -8.072E-4        8.497E-5        0.001072
muy        -2.423        9.849        0.2522        -21.2        -2.358        17.65
rho        -0.9659        0.007285        5.123E-5        -0.9782        -0.9666        -0.9498
sigv2        0.001012        9.768E-5        6.388E-7        8.387E-4        0.001006        0.001221
sigy2        5.047        2.857        0.03457        1.968        4.325        12.46