医学研究的R语言应用(4):
方差分析与多重比较
1. 研究实例
在评价某临床新药耐受性及安全性的2a期临床试验中,对符合纳入标准的30名健康自愿者随机分为3组,每组10名。各组注射剂量分别为0.5U、1U、2U,观察48小时部分凝血活酶时间(s),试问不同剂量的部分凝血活酶时间有无不同?数据库见time48.sav(SPSS格式)
2. 数据导入和基本分析
#加载必要的包
library(foreign) # 第1篇介绍
library(nortest) # 第1篇介绍
library(CBCgrps) # 第2篇介绍
library(tidyverse) # 第1篇介绍
library(car) #用于方差齐性检验
#导入数据
e4=read.spss("time48.sav") #读取数据集
e4
# 浏览数据变量结构
str(e4)
#利用函数整体进行快速统计描述
fun1
median
mean
IQR
sd
y5
return(y5)
}
fun1(e4$time,e4$group)
# 此处也可以利用dplyr函数进行