我正在使用omu
软件包进行代谢物分析。我有一个计数数据文件和元数据文件。软件包的omu_summary
函数对计数文件执行t检验,如下所示:
DF_stats <- omu_summary(count_data = DF, metadata = meta.data, numerator = "Summer",
denominator = "Winter", response_variable = "Metabolite", Factor = "Treatment",
log_transform = TRUE, p_adjust = "BH")
执行函数时引发以下错误:
if(stderr <10 * .Machine $ double.eps * max(abs(mx), abs(my)))stop(“数据本质上是恒定的”):缺少值 需要TRUE / FALSE的地方
str(DF)返回:
Classes 'cpd' and 'data.frame': 1340 obs. of 11 variables:
$ Metabolite: chr "Pilocereine" "Quercetin 3-(2G-xylosylrutinoside)" "Cyanidin 3-O-(6''-glucosyl-2''-xylosylgalactoside)" "Troxerutin" ...
$ KEGG : chr "C09609" "C10175" "C08612" "D07180" ...
$ W1 : num 8282 1628 0 238904 42965 ...
$ W2 : num 10520 1528 0 342672 30306 ...
$ S1 : num 12829 14628 0 188096 46352 ...
$ S2 : num 31750 14628 14628 180538 0 ...
$ Class : Factor w/ 13 levels "Alcohols","Antibiotics",..: NA 7 7 NA NA 12 12 7 12 NA ...
$ Subclass_1: Factor w/ 48 levels "Alkaloids","Amines",..: NA 36 36 NA NA 47 3 36 16 NA ...
$ Subclass_2: Factor w/ 201 levels "2-Deoxystreptamines",..: NA 86 86 NA NA 199 55 86 86 NA ...
$ Subclass_3: Factor w/ 327 levels "1-alkyl,2-acylglycerols",..: NA 156 23 NA NA 248 92 23 157 NA ...
$ Subclass_4: Factor w/ 78 levels "Abietane diterpenoids",..: NA 54 54 NA NA 54 54 54 54 NA ...
测试是否在W1,W2,S1和S2上执行,这似乎是数字。我该如何解决这个问题?