Question

我有"age_30_35", "age_36_40"二进制（0或1）指标列，如( age_30_35=0, age_36_40=1 )等。例如，36岁的人会(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')。这是一个例子：

age_30_35

我想制作一些关于每个类别如何影响连续响应变量＆＃34; muscle_strength的图表。＆＃34;例如，每个组的框图：参考类别age_36_40，{ "respuesta": true }。我知道怎么做，如果我要制作一个额外的因子列，指示符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)

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

1 个答案: