作者yummy7922 (crucify)
看板R_Language
标题Re: [问题] 分组资料，找符合条件人数
时间Wed Mar 11 22:31:55 2015
※ 引述《yummy7922 (crucify)》之铭言：
: ※ 引述《yummy7922 (crucify)》之铭言：
: : [问题类型]:
: : 程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
: : [软体熟悉度]:
: : 入门(写过其他程式，只是对语法不熟悉)
: : [问题叙述]:
: : 我的资料是一个重复测量的资料，资料中有当期是否使用药物的资料(若有使用，设为1，
: : 没有为0)、还有记录该笔资料为该位病人的第几笔观察值，我希望能算出，药物转换的机
: : 率，例如有多少人第一期使用a药物，但在第二期时转换成使用b药物。
: : 资料中共有6种药物，一位病人最多有48笔观察值，
: : (6x6)种转换可能 x 47(个时间隔) = 1692 个机率值。
: : 想请教各位高手们，该怎麽做比较有效率。
不好意思，我没有说明清楚，
我想算的机率其实是很直观的条件机率，例如：
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表示为该位病人的第几个月的观察值。
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 61.231.220.242
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1426084319.A.DB5.html
1F：推 squallscer: 第一期跟第二期该如何判断呢？ 03/11 23:10
2F：推 celestialgod: 建议以後用修改，增加内容 03/11 23:17
3F：→ celestialgod: 修改请按E 不然会让版上文章变乱 03/11 23:18
4F：→ yummy7922: 我的期别单位是月，也就是duration唷 03/12 00:04
5F：→ yummy7922: 好的，谢谢，因为我不太会用，不好意思 03/12 00:08
6F：→ aaron77217: 你的资料是一次可以吃多种药吗? 03/12 00:27
7F：→ yummy7922: 没错唷，但我先暂时忽略这件事，算一个概略的值 03/12 10:15
8F：→ celestialgod: 那我的code应该不敷使用 03/12 10:59
9F：→ yummy7922: 是的 我昨晚有发现 还在想怎麽做更改 真是谢谢您 03/12 15:36
10F：→ celestialgod: 原本的写法很难改成你需要的... 03/12 17:02