+ ## 6 可选：持续上调/下调基因 + coef <- tfit$sig.models + up <- rownames(coef)[coef$time >...
🚨 错误信息
+ ## 6 可选：持续上调/下调基因
+ coef <- tfit$sig.models
+ up <- rownames(coef)[coef$time > 0 & coef$time < 0.05 & is.na(coef$time2)]
+ down <- rownames(coef)[coef$time < 0 & coef$time < 0.05 & is.na(coef$time2)]
+ write_csv(tibble(gene = up), str_c(prefix, "_up_genes.csv"))
+ write_csv(tibble(gene = down), str_c(prefix, "_down_genes.csv"))
+ }
====== running Control ======
Rows: 13692 Columns: 13
── Column specification ──────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): symble
dbl (12): 6hC1, 6hC2, 6hC3, 24hC1, 24hC2, 24hC3, 72hC1, 72hC2, 72hC3, 120hC1, 120hC2, 120hC3
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
Rows: 12 Columns: 4
── Column specification ──────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (2): sample, group
dbl (2): time, replicate
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
[1] "fitting gene 100 out of 13692"
[1] "fitting gene 200 out of 13692"
[1] "fitting gene 300 out of 13692"
[1] "fitting gene 400 out of 13692"
[1] "fitting gene 500 out of 13692"
[1] "fitting gene 600 out of 13692"
[1] "fitting gene 700 out of 13692"
[1] "fitting gene 800 out of 13692"
[1] "fitting gene 900 out of 13692"
[1] "fitting gene 1000 out of 13692"
[1] "fitting gene 1100 out of 13692"
[1] "fitting gene 1200 out of 13692"
[1] "fitting gene 1300 out of 13692"
[1] "fitting gene 1400 out of 13692"
[1] "fitting gene 1500 out of 13692"
[1] "fitting gene 1600 out of 13692"
[1] "fitting gene 1700 out of 13692"
[1] "fitting gene 1800 out of 13692"
[1] "fitting gene 1900 out of 13692"
[1] "fitting gene 2000 out of 13692"
[1] "fitting gene 2100 out of 13692"
[1] "fitting gene 2200 out of 13692"
[1] "fitting gene 2300 out of 13692"
[1] "fitting gene 2400 out of 13692"
[1] "fitting gene 2500 out of 13692"
[1] "fitting gene 2600 out of 13692"
[1] "fitting gene 2700 out of 13692"
[1] "fitting gene 2800 out of 13692"
[1] "fitting gene 2900 out of 13692"
[1] "fitting gene 3000 out of 13692"
[1] "fitting gene 3100 out of 13692"
[1] "fitting gene 3200 out of 13692"
[1] "fitting gene 3300 out of 13692"
[1] "fitting gene 3400 out of 13692"
[1] "fitting gene 3500 out of 13692"
[1] "fitting gene 3600 out of 13692"
[1] "fitting gene 3700 out of 13692"
[1] "fitting gene 3800 out of 13692"
[1] "fitting gene 3900 out of 13692"
[1] "fitting gene 4000 out of 13692"
[1] "fitting gene 4100 out of 13692"
[1] "fitting gene 4200 out of 13692"
[1] "fitting gene 4300 out of 13692"
[1] "fitting gene 4400 out of 13692"
[1] "fitting gene 4500 out of 13692"
[1] "fitting gene 4600 out of 13692"
[1] "fitting gene 4700 out of 13692"
[1] "fitting gene 4800 out of 13692"
[1] "fitting gene 4900 out of 13692"
[1] "fitting gene 5000 out of 13692"
[1] "fitting gene 5100 out of 13692"
[1] "fitting gene 5200 out of 13692"
[1] "fitting gene 5300 out of 13692"
[1] "fitting gene 5400 out of 13692"
[1] "fitting gene 5500 out of 13692"
[1] "fitting gene 5600 out of 13692"
[1] "fitting gene 5700 out of 13692"
[1] "fitting gene 5800 out of 13692"
[1] "fitting gene 5900 out of 13692"
[1] "fitting gene 6000 out of 13692"
[1] "fitting gene 6100 out of 13692"
[1] "fitting gene 6200 out of 13692"
[1] "fitting gene 6300 out of 13692"
[1] "fitting gene 6400 out of 13692"
[1] "fitting gene 6500 out of 13692"
[1] "fitting gene 6600 out of 13692"
[1] "fitting gene 6700 out of 13692"
[1] "fitting gene 6800 out of 13692"
[1] "fitting gene 6900 out of 13692"
[1] "fitting gene 7000 out of 13692"
[1] "fitting gene 7100 out of 13692"
[1] "fitting gene 7200 out of 13692"
[1] "fitting gene 7300 out of 13692"
[1] "fitting gene 7400 out of 13692"
[1] "fitting gene 7500 out of 13692"
[1] "fitting gene 7600 out of 13692"
[1] "fitting gene 7700 out of 13692"
[1] "fitting gene 7800 out of 13692"
[1] "fitting gene 7900 out of 13692"
[1] "fitting gene 8000 out of 13692"
[1] "fitting gene 8100 out of 13692"
[1] "fitting gene 8200 out of 13692"
[1] "fitting gene 8300 out of 13692"
[1] "fitting gene 8400 out of 13692"
[1] "fitting gene 8500 out of 13692"
[1] "fitting gene 8600 out of 13692"
[1] "fitting gene 8700 out of 13692"
[1] "fitting gene 8800 out of 13692"
[1] "fitting gene 8900 out of 13692"
[1] "fitting gene 9000 out of 13692"
[1] "fitting gene 9100 out of 13692"
[1] "fitting gene 9200 out of 13692"
[1] "fitting gene 9300 out of 13692"
[1] "fitting gene 9400 out of 13692"
[1] "fitting gene 9500 out of 13692"
[1] "fitting gene 9600 out of 13692"
[1] "fitting gene 9700 out of 13692"
[1] "fitting gene 9800 out of 13692"
[1] "fitting gene 9900 out of 13692"
[1] "fitting gene 10000 out of 13692"
[1] "fitting gene 10100 out of 13692"
[1] "fitting gene 10200 out of 13692"
[1] "fitting gene 10300 out of 13692"
[1] "fitting gene 10400 out of 13692"
[1] "fitting gene 10500 out of 13692"
[1] "fitting gene 10600 out of 13692"
[1] "fitting gene 10700 out of 13692"
[1] "fitting gene 10800 out of 13692"
[1] "fitting gene 10900 out of 13692"
[1] "fitting gene 11000 out of 13692"
[1] "fitting gene 11100 out of 13692"
[1] "fitting gene 11200 out of 13692"
[1] "fitting gene 11300 out of 13692"
[1] "fitting gene 11400 out of 13692"
[1] "fitting gene 11500 out of 13692"
[1] "fitting gene 11600 out of 13692"
[1] "fitting gene 11700 out of 13692"
[1] "fitting gene 11800 out of 13692"
[1] "fitting gene 11900 out of 13692"
[1] "fitting gene 12000 out of 13692"
[1] "fitting gene 12100 out of 13692"
[1] "fitting gene 12200 out of 13692"
[1] "fitting gene 12300 out of 13692"
[1] "fitting gene 12400 out of 13692"
[1] "fitting gene 12500 out of 13692"
[1] "fitting gene 12600 out of 13692"
[1] "fitting gene 12700 out of 13692"
[1] "fitting gene 12800 out of 13692"
[1] "fitting gene 12900 out of 13692"
[1] "fitting gene 13000 out of 13692"
[1] "fitting gene 13100 out of 13692"
[1] "fitting gene 13200 out of 13692"
[1] "fitting gene 13300 out of 13692"
[1] "fitting gene 13400 out of 13692"
[1] "fitting gene 13500 out of 13692"
[1] "fitting gene 13600 out of 13692"
错误于match.arg(method):
'arg' should be one of “holm”, “hochberg”, “hommel”, “bonferroni”, “BH”, “BY”, “fdr”, “none”
🤖 AI解决方案
错误于match.arg(method):
'arg' should be one of “holm”, “hochberg”, “hommel”, “bonferroni”, “BH”, “BY”, “fdr”, “none”- 错误翻译：
Error in match.arg(method):
'arg' 应该是 "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none" 中的一个
在使用
p.adjust 函数（或其他需要 method 参数的函数）时，传入了一个非法的多重检验校正方法。检查调用
p.adjust 或类似函数的地方，确保 method 参数的值是上述列表中的一个有效选项。例如，如果你想使用 Benjamini-Hochberg 方法，应该写 method = "BH"。1. 在代码中硬编码可用的
method 选项列表，并在用户输入或使用这些选项前进行校验。2. 为
method 参数提供一个默认值，该默认值必须是有效的选项。