1

我有一个 500 行 38 列的基因表达数据框。我想对第 435 行的特定基因的两组（第 1-27 列和第 28-39 列）执行 t.test 。这是我的 R 代码

x <- df[435,1:27]
y <- df[435,28:38]
t.test(x,y)

x 的输出是

  G1    G2    G3    G4    G5    G6    G7    G8    G9    G10

435 -0.19524 -0.13841 -0.05966 -0.0742 -0.4371 0.05094 0.45927 0.04794 -0.495 0.39688

   G11    G12    C13    G14    G15   G16    G17    G18    G19    G20 

 -0.3524 0.04483 -0.25422 0.2458 -0.49332 -0.18619  -0.56592 -0.13109 -0.103 -0.57249
   G21    G22    G23    G24    G25    G26    G27

435 -0.26623 0.02744 0.07891 -0.04138 -0.35251 0.17627 0.90422

y 的输出是

    G28     G29      G30    G31     G32     G33     G34     G35

435 -0.20156 -0.56749 0.1803 0.47034 0.43657 0.31767 0.29962 0.62897

G36 G37 G38

0.32455 0.22439 -0.07216

在运行 t.test 后给我一条消息

if (stderr < 10 * .Machine$double.eps * max(abs(mx), abs(my))) stop("数据基本上是恒定的") 中的错误：需要 TRUE/FALSE 的缺失值另外：警告消息：1：在 mean.default(x) 中：参数不是数字或逻辑：返回 NA 2：在 mean.default(y) 中：参数不是数字或逻辑：返回 NA

但是，如果我以这种方式输入数据，则 t.test 可以执行

x <- c(-0.19524,-0.13841,-0.05966,-0.0742,-0.4371,0.05094,0.45927,0.04794,-0.495,0.39688,-0.3524,0.04483,-0.25422,0.2458,-0.49332,-0.18619,-0.56592,-0.13109,-0.103,-0.57249,-0.26623,0.02744,0.07891,-0.04138,-0.35251,0.17627,0.90422)

y <- c(-0.20156,-0.56749,0.1803,0.47034,0.43657,0.31767,0.29962,0.62897,0.32455,0.22439,-0.07216)

如何修改我的 R 代码以便 t.test 可以在 RStudio 中执行？谢谢

4

0 回答 0