Question

我看了this SO question和that one，但仍然无法解决我的问题。我有以下data.table，其中只包含我的data.table的总列和行的几个。

library(data.table)
structure(list(Patient = c("MB108", "MB108", "MB108", "MB108", 
"MB108", "MB108", "MB108", "MB108", "MB108", "MB108"), Visit = c(1, 
1, 1, 1, 9, 9, 9, 9, 12, 12), Stimulation = c("NC", "SEB", "PPD", 
"E6C10", "NC", "SEB", "PPD", "E6C10", "NC", "SEB"), `CD38   ` = c(83.3, 
63.4, 83.2, 91.5, 90.9, 70.9, 71, 88.4, 41.7, 47.9)), .Names = c("Patient", 
"Visit", "Stimulation", "CD38   "), class = c("data.table", "data.frame"
), row.names = c(NA, -10L), .internal.selfref = <pointer: 0x102806578>)

我想在第4列进行t.test，当访问次数为1时，访问次数为9次。我检查了NAs以及两列的长度。

感谢您的帮助！

      #na.omit(boolean_dt3)
      #print(length(unlist(boolean_dt3[Visit== 1,4, with = FALSE])))
      #print(length(unlist(boolean_dt3[Visit== 9,4, with = FALSE])))

wilcox.test( unlist(boolean_dt3[Visit== 1,4, with = FALSE])~ unlist(boolean_dt3[Visit== 9,4, with = FALSE]) , paired = T, correct=FALSE)

Answer 1

我刚认为,而不是~适用于我的问题。

Answer 2

以下是column 4

对Value分组进行wilcoxon测试的方法

 library(dplyr)
 wilcox.test( filter(df, Visit==1)$CD38, filter(df, Visit==9)$CD38, paired=TRUE)

Answer 3

尝试这个表述：

wilcox.test(numeric_var ~ two_level_group_var)

R：wilcoxon测试错误：分组因子必须正好有2个级别

3 个答案: