测试按组的差异

Question

我有一个由x个变量组成的数据集，var1，var2，...，varx 我感兴趣的是变量1，var1，一个有三个级别的因子;“a”，“b”，“c”，以及它如何受到其他变量的影响。

我用

做了一些描述性统计

describeBy(dataset, group=var1)

现在我想查看每个变量，看看不同grops中的varj的意思是“a”，“b”，“c”。

编辑1： 很抱歉没有像我想象的那样清晰...... 其实我的问题是这个，因为我知道我想要使用双面t检验

t.test(varj, alternative="two.sided",conf.level = 0.95,
   subset=var1)

我确实得到了一个结果，但只是针对varj的平均值，而不是针对不同组的差异。 我做错了什么？

编辑2： 我觉得我有点厌倦了。

但这是我想要做的事情

anova(lm(var1 ~ varj))

似乎很好。

Answer 1

这听起来像单因素方差分析？？

试试这个：

model1 = aov(var1 ~ varj)

summary(model1)

如果你想知道哪些配对/群体彼此显着不同，你可以使用例如Bonferroni

pairwise.t.test(var1, varj, p.adj = "bonf")