※ 引述《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表示为该位病人的第几个月的观察值。