我正在尝试使用 RStatix's 执行一系列 T 检验t_test()
,其中因变量在每个测试中都是相同的,并且分组变量会发生变化。我在一个循环中进行这些测试,所以我想用列号而不是列名来选择分组变量。我尝试使用 执行此操作colnames(dataframe)[[columnnumber]]
,但出现以下错误:“无法提取不存在的列”。如何使用列号而不是列名选择分组变量?
下面是一个带有虚构数据框的最小可复制示例;当指示分组变量的名称(性别)时,测试正常工作,但在指示列号时不能正常工作。
library(tidyverse)
library(rstatix)
dat<-data.frame(gender=rep(c("Male", "Female"), 1000),
age=rep(c("Young","Young", "Old", "Old"),500),
tot= round(runif(2000, min=0, max=1),0))
dat %>% t_test(tot ~ gender,detailed=T) ##Works
dat %>% t_test(tot ~ colnames(dat)[[1]],detailed=T) ##Doesn't work