在R中绘制二进制协变量响应

Question

我有一个data.frame带有二进制（0或1）指示符列，例如"age_30_35", "age_36_40"等。例如，一个36岁的人将拥有( age_30_35=0, age_36_40=1 ) 。 另外，还有一个参考类别(age_30_35=0, age_36_40=0) 。 这是一个例子：

data<- data.frame(c(1,0,0),c(0,1,0),c(2.5,1.3,4.5))
colnames(data)=c('age_30_35','age_36_40','muscle_strength')

我想对每个类别如何影响连续响应变量“ muscle_strength”进行一些绘制。 例如，每个组的箱形图：参考类别age_30_35 ， age_36_40 。 如果我要制作一个额外的因子列，其中的指标1,2,3代表每个组，我知道该怎么做。 在没有实际创建额外因素的情况下，最轻松的方法是什么？ 特别是，我想要一个可以轻松扩展到许多列的解决方案。

Answer 1

我想您可以通过以下方式创建额外的factor ：

i <- 1:2 #indices of columns with dummies 
f <- factor(
    as.matrix(data[,i]) %*% 1:length(i), 
    labels=c("ref", names(data[,i]))
)
boxplot(data$muscle_strength ~ f)