Ошибка, говоря, что у меня есть значение NA / NaN / Inf, когда кажется, что нет ни одного. Ошибка в hclustfun (distc): NA / NaN / Inf при вызове сторонней функции (аргумент 11) - PullRequest
Новая
       36

Ошибка, говоря, что у меня есть значение NA / NaN / Inf, когда кажется, что нет ни одного. Ошибка в hclustfun (distc): NA / NaN / Inf при вызове сторонней функции (аргумент 11)

0 голосов
/ 31 октября 2018

У меня есть файл, который выглядит так: 

Taxa,D0,D1,...
1,0,10,...
20,0,0,...
...

Файл доступен здесь: https://www.dropbox.com/s/lvtxztz3kypnxre/Taxa_Conf.csv?dl=0

Когда я запускаю скрипт по 

$CONCOCT/scripts/ConfPlot.R -c Taxa_Conf.csv -o Taxa_Conf.pdf

Я получаю следующую ошибку: 

Error in hclustfun(distc) : NA/NaN/Inf in foreign function call (arg 11)
Calls: heatmap.2 -> hclustfun
Execution halted

где ConfPlot.R выглядит как 

#!/usr/bin/Rscript

#load libraries
library(gplots)
library(getopt)

spec = matrix(c('verbose','v',0,"logical",'help','h',0,"logical",'confile','c',1,"character",'ofile','o',1,"character"),byrow=TRUE,ncol=4)

opt=getopt(spec)

if( !is.null(opt$help)) {
cat(getopt(spec, usage=TRUE)); 
q(status=1);
}

confFile <- opt$confile

Conf <- read.csv(confFile,header=TRUE,row.names=1)
Conf.t <- t(Conf)
ConfP <- Conf.t/rowSums(Conf.t)

crp <- colorRampPalette(c("blue","red","orange","yellow"))(100)

pdf(opt$ofile)

heatmap.2 (as.matrix(t(ConfP)),col=crp,trace="none",dendrogram="none",Rowv=FALSE,lwid = c(1.25,4.5),lhei = c(1.25,4.5),cexRow=0.5)

dev.off()

После прочтения Ошибка при создании тепловых карт - NA / NaN / Inf при вызове сторонней функции (arg 11) , я пытался 

> a=read.csv("/home/mathed/Simulation/custom_desman/1/Strains/Simulation/Metabat/Taxa_Conf.csv",header=TRUE,row.names=1)
> any(is.na(a))
[1] FALSE

Кажется, нет значения NA. ConfPlot.R работает для других моих файлов. В чем проблема?? Кажется, что нет значений na, nan или inf. Заранее спасибо

1 Ответ

0 голосов
/ 31 октября 2018

Проверьте эту строку: 

ConfP <- Conf.t/rowSums(Conf.t)

Должно быть, по существу, деление на ноль: 

> 0/0
[1] NaN

> 1/0
[1] Inf
...