R语言的rbind函数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

> table(dose)
  dose
  0   5  50 500
20  19  18  17

> contrast<-rbind("no drug vs. drug"=c(3,-1,-1,-1))  #解释是设定第一组和其他三组的均值进行比较，不懂为什么是c(3，-1，-1，-1)
求大神赐教！！！！！！！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：rbind bind IND bin R语言 contrast

不太理解你的问题。 一般地，c()是定义list或说向量的意思，而rbind()常用于将向量/list按行组成矩阵matrix

hyu9910 发表于 2016-1-6 20:09
不太理解你的问题。 一般地，c()是定义list或说向量的意思，而rbind()常用于将向量/list按行组成矩阵matrix

这是单因素协方差分析的一个例子，来自于multcomp包中的litter数据集，怀孕小鼠
被分为四个小组，每个小组接受不同剂量（0、5、50或500）的药物处理。产下幼崽的体重均值
为因变量，怀孕时间为协变量。利用multcomp包来检验用户自定义的均值假设。假定对未用药条件（第一组）与其他三种用药条件（后三组）影响是否不同感兴趣，代码是：
> library(multcomp)
>contrast<-rbind("no drug vs.drug"=c(3,-1,-1,-1))
>summary(glht(fit,linfct=mcp(dose=contrast)))

         Simultaneous Tests for General Linear Hypotheses

Multiple Comparisons of Means: User-defined Contrasts


Fit: aov(formula = weight ~ gesttime + dose)

Linear Hypotheses:
                      Estimate Std. Error t value Pr(>|t|)  
no drug vs. drug == 0    8.284      3.209   2.581    0.012 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Adjusted p values reported -- single-step method)
我不懂的是为什么c(3,-1,-1,-1)表示第一组和后三组的比较，如果我想表示前两组和后两组的比较应该是什么呢？

fresher@_@ 发表于 2016-1-7 08:51
这是单因素协方差分析的一个例子，来自于multcomp包中的litter数据集，怀孕小鼠
被分为四个小组，每个小 ...

没有做你的项目，建议你结合R命令的定义去领会哈

hyu9910 发表于 2016-1-7 10:34
没有做你的项目，建议你结合R命令的定义去领会哈

嗯呢，还是谢谢你！

fresher@_@ 发表于 2016-1-7 15:06
嗯呢，还是谢谢你！

请问你的问题解决了吗？

这里的多重比较函数mcp的参数由rbind给出，3表示该剂量组跟其他几组对比，这里是跟其他三组对比，-1表示该组为对照组，0表示该组不参与比较；dose=contrast表示以剂量分组进行多重比较。

dazzlingpuck 发表于 2017-9-23 17:49
这里的多重比较函数mcp的参数由rbind给出，3表示该剂量组跟其他几组对比，这里是跟其他三组对比，-1表示该组 ...

厉害了！

dazzlingpuck 发表于 2017-9-23 17:49
这里的多重比较函数mcp的参数由rbind给出，3表示该剂量组跟其他几组对比，这里是跟其他三组对比，-1表示该组 ...

非常感谢，get~

dazzlingpuck 发表于 2017-9-23 17:49
这里的多重比较函数mcp的参数由rbind给出，3表示该剂量组跟其他几组对比，这里是跟其他三组对比，-1表示该组 ...

厉害了！！