Re: [問題] 求取轉移機率矩陣

看板R_Language作者celestialgod (攸藍)時間9年前 (2015/03/12 14:17)推噓1(1推 0噓 0→)

留言1則, 1人參與討論串3/3 (看更多)

※ 引述《yummy7922 (crucify)》之銘言： : ※ 引述《yummy7922 (crucify)》之銘言： : 不好意思，我沒有說明清楚， : 我想算的機率其實是很直觀的條件機率，例如： : P(第二期使用b藥物|第一期使用a藥物) : = (第一期使用a藥且第二期使用b藥的人數)/(第一期使用a藥的人數) : 但是我不知道該怎麼在分組之後，還能夠給定條件， : 計算出第一期使用a藥第二期使用b藥的人數。 : 資料大概是長這樣： : id M_01 M_02 M_03 M_04 M_05 M_06 duration : ace01 1 1 ..... 1 : ace01 0 1 2 : ace01 0 0 3 : swe02 0 1 1 : swe02 1 1 2 : kue03 0 1 1 : kue03 1 0 2 : kue03 0 1 3 : kue03 1 0 4 : kue03 1 1 5 : 其中，M_01表示有沒有服用第一種藥，共有六種藥， : duration表示為該位病人的第幾個月的觀察值。 Updata a faster way to do. code: http://pastebin.com/Eg7vn058 > transitMatrix_eachTime[[1]] # 第二期的轉移矩陣 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.1579643 0.1652346 0.1625909 0.1711831 0.1784534 0.1645737 [2,] 0.1698612 0.1692003 0.1639128 0.1625909 0.1771315 0.1573034 [3,] 0.1635638 0.1775266 0.1569149 0.1675532 0.1682181 0.1662234 [4,] 0.1661085 0.1654387 0.1634293 0.1654387 0.1694575 0.1701273 [5,] 0.1722746 0.1648721 0.1561238 0.1641992 0.1709287 0.1716016 [6,] 0.1675862 0.1620690 0.1648276 0.1655172 0.1765517 0.1634483 (1,1)位置就是上一期服dose 1，下一期也是dose 1的條件機率 (1,2)位置就是上一期服dose 1，下一期服dose 2的條件機率，剩下依此類推 PS: 因此，轉移矩陣的每一個列，其和為1 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 223.136.172.140 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1426141048.A.600.html