故事源于在某项目中,我需要统计0.75和1两个剂量组AE of intrest的发生例数,最后是得到下面这张表:

其实,一般情况下,这里有一个很简便的做法就是用table函数:
> diabetes <- c("Type1", "Type2", "Type1", "Type1","Type1")
> status <- c("Poor", "Improved", "Excellent", "Poor","Poor")
> patientdata <- data.frame(patientID, age, diabetes, status)
> patientdata

> table(patientdata$diabetes, patientdata$status)
Excellent Improved Poor
Type1 1 0 3
Type2 0 1 0
但是呢,我当时的情况和上面这个例子中的数据一样,ID不是唯一的,也就是说一个受试者同一个Dose和不良反应会有多行数据,这里直接用table()会把ID相同的受试者的两行当成两个人算。
无奈之下,我只能筛选出剂量为特定值,不良反应为特定值的数据集,然后求该数据集的ID(factor类)列有几中元素,在R中的表达是这样的:用unique()和lengh()
