Re: [問題] 求取轉移機率矩陣
※ 引述《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
推
03/12 15:33, , 1F
03/12 15:33, 1F
※ 編輯: celestialgod (223.136.172.140), 03/13/2015 15:26:26
討論串 (同標題文章)