求解：错误于numericDeriv(form[[3L]], names(ind), env) : 在计算模型的时候产...

fenghua

我运行如下的命令，出现错误， 求解
> M.gen.diss <- read.table(file = 'D:/DResD-T/GenDistance.csv', sep = ',', header = T, row.names = 1, fill = T)
> M.gen.diss <- as.dist(M.gen.diss)
> sample.points <- read.table(file = 'D:/DResD-T/GeoDistance.csv', sep = ',', header = T, row.names = 1, fill = T)
> names(sample.points) <- c('x', 'y')
> M.geo.dist <- dist(sample.points[,1:2])
> IBD.data <- data.frame(cbind(as.vector(M.geo.dist), as.vector(M.gen.diss)))
> names(IBD.data) <- c('geo.dist', 'gen.diss')
> IBD.model.prelim <- NLSstAsymptotic(sortedXyData(x = IBD.data[1:1000,1], y = IBD.data[1:1000,2]))
错误于numericDeriv(form[[3L]], names(ind), env) :
  在计算模型的时候产生了缺省值或无限值
此外: 警告信息：
In log(-coef(lm(log(abs(y - rt)) ~ x, data = xy))[2L]) : 产生了NaNs

fenghua

请高手帮忙啊！

biostatistic

警告说的很明显啊
此外: 警告信息：
In log(-coef(lm(log(abs(y - rt)) ~ x, data = xy))[2L]) : 产生了NaNs

你数据中产生 NAN值了
也就是说log(abs(y - rt))中 abs(y-rt)的值 有一些是0