python: invalid value encountered in divide以及invalid value encountered in double_scalars报错
运行命令python eqtl_prepare_expression.py data.tpm.gct data.reads_count.gct --tpm_threshold 0.1 --count_threshold 2 --sample_frac_threshold 0.2 --normalization_method tmm --output data.txt
时出现了报错“invalid value encountered in divide”以及“invalid value encountered in double_scalars”的报错,问chatgpt,给了我以下的解决方案:
按照上述方案,去除无效值(如NaN或Inf):
cleaned_data <- data %>%
mutate(across(col4:coln,
~ as.numeric(.) %>% # 将列转换为 double 类型
na_if(Inf) %>% # 将 Inf 替换为 NA
na_if(NaN))) # 将 NaN 替换为 NA
再次运行命令python eqtl_prepare_expression.py data.tpm.gct data.reads_count.gct --tpm_threshold 0.1 --count_threshold 2 --sample_frac_threshold 0.2 --normalization_method tmm --output data.txt
,没啥卵用~~,还是出现相同的报错,说明不是无效值导致的。
再次检查数据,发现是某列均是0导致运行出现的报错。那就去除全是0的列:
tr1=tr %>% select(where(is.numeric)) %>% select_if(~ sum(.) != 0)
再次运行命令,总算正常了~
本文来自博客园,作者:橙子牛奶糖(陈文燕),转载请注明原文链接:https://www.cnblogs.com/chenwenyan/p/18474015
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 10年+ .NET Coder 心语 ── 封装的思维:从隐藏、稳定开始理解其本质意义
· 地球OL攻略 —— 某应届生求职总结
· 提示词工程——AI应用必不可少的技术
· Open-Sora 2.0 重磅开源!
· 字符编码:从基础到乱码解决
2022-10-18 LDSC 运行 Partitioned Heritability 报错:ValueError: operands could not be broadcast together with shapes (1,210) (1,200)
2020-10-18 本周最新文献速递20201018
2015-10-18 1-3