R中的分层聚类

Question

使用hclust功能时出现以下错误？我怎么解决这个问题？ 我使用的是Windows 7和2.12.3版本的R。

n_seq <- 250                                                   
mat <- matrix(NA, ncol=n_seq, nrow=n_seq)    
for (idx in 1:n_seq) 
{mat[idx,idx] <- 0.0}         
for(idx in 1:(n_seq-1) )
{intemp <- read.xls("C:// clustal.xls", sheet = idx ); 
mat[(1+idx):n_seq,idx] <- intemp[1:(n_seq-idx), 11]}

fit <- hclust(as.dist(mat), method="single")

Error in hclust(as.dist(mat), method = "single") : 
NA/NaN/Inf in foreign function call (arg 11)

请帮我解决这个问题。

Answer 1

出现此错误消息是因为距离矩阵as.dist(mat)中包含错误值（NA，NaN或Inf）。如果查看hclust的代码，外部函数调用的第11个参数（即通过编译代码）是距离矩阵的值，这就是R抱怨的内容。默认情况下，外部函数调用不接受NA，NaN或Inf。有关详细信息，请参阅?.Fortran。