治疗效果的定量分析

Question

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

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

Syrvn

Answer 1

这听起来像是一个multcomp-package的工作，那里有很多资源。另请查看this本书。

以下是比较两组对照组的反应的例子：

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，