R(rstudio)语言协方差分析

使用不同剂量的药物对小鼠进行处理怀孕时间对产下小鼠体重的影响协方差分析

1.载入multcomp包

> library(multcomp)

载入需要的程辑包：mvtnorm

载入需要的程辑包：survival

载入需要的程辑包：TH.data

载入需要的程辑包：MASS

载入程辑包：‘TH.data’

The following object is masked from ‘package:MASS’:

    geyser

 

 2.使用multcomp包中的litter数据集

> litter

  3.统计分组情况

> table(litter$dose)

4.分组统计平均数

> attach(litter)

> aggregate(weight, by=list(dose), FUN=mean)

 结果分析：不同处理组体重不同

5.定义fit变量进行方差分析

> fit <- aov(weight ~gesttime+dose,data=litter)  #gesttime为协变量

> summary(fit)

            Df Sum Sq Mean Sq F value  Pr(>F)   

gesttime     1  134.3  134.30   8.049 0.00597 **

dose         3  137.1   45.71   2.739 0.04988 *

Residuals   69 1151.3   16.69                   

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

 

 结果分析：小鼠怀孕时间与小鼠体重相关    药物剂量对小鼠体重有影响

6.plot绘制方差结果

> plot(fit)

Hit <Return> to see next plot:

Hit <Return> to see next plot:

Hit <Return> to see next plot:

Hit <Return> to see next plot:

>