治疗效果的定量分析

Question

我需要分析一项生物学研究。研究设计相当简单。它包括3组。 1个对照组 和 2 个测试组。两个试验组用相同的药物治疗，但剂量不同。每组都有 大约 14 个观察值，我查看了大约 600 个变量。我已经使用单向方差分析 确定组之间的显着命中。但是，我现在想知道是否有定量 治疗剂量之间的影响。换句话说，当您使用 将测试组 1（较低剂量）与对照组进行比较，或者比较时治疗效果更大 试验组 2（较高剂量）与对照组。我还需要在分析中包含 1 或 2 个协变量。

不幸的是，我不知道 R 中是否有可以用来回答这个问题的统计测试/技术。 非常感谢任何建议。

Syrvn

Answer 1

这听起来像是 multcomp-package 的工作，那里有很多资源。也可以看看this book。

下面是一个将两组的反应与对照组进行比较的示例：

require(multcomp)
mice <- data.frame(group=as.factor(rep(c("C","1","2"),rep(6,3))), 
                   score=c(58, 32, 59, 64, 55, 49, 73, 70, 68, 71, 60, 62, 53, 74, 72, 62, 58, 61))
# reoder factor, so that Control is the 1st level
levels(mice$group) <- c("C", "1", "2")

plot(score ~ group, data = mice)

# Anova
mod <- aov(score ~ group, data = mice)
# Multiple Comparisons with Dunnett contrasts (=Compare to control)
summary(glht(mod, linfct=mcp(group = "Dunnett")))

     Simultaneous Tests for General Linear Hypotheses

Multiple Comparisons of Means: Dunnett Contrasts


Fit: aov(formula = score ~ group, data = mice)

Linear Hypotheses:
           Estimate Std. Error t value Pr(>|t|)  
1 - C == 0   -4.000      4.965  -0.806   0.6427  
2 - C == 0  -14.500      4.965  -2.920   0.0195 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
(Adjusted p values reported -- single-step method)

HTH，