Блог

Я пытаюсь выполнить начальную загрузку модели пространственного подсчета, используя spatcounts-package в R. Однако повторная выборка с заменой может привести к «островам», потому что иногда соседи не будут выбраны. Я полагаю, что это вызывает сообщение об ошибке: Error: NA/NaN/Inf in foreign function call (arg 1) . Может ли кто-нибудь подтвердить это, если это действительно так? Если да, то есть ли способ обойти это, чтобы я все еще мог выполнять загрузочную загрузку?

Пожалуйста, найдите ниже код для воссоздания проблемы, используя пример набора данных пакета spatcounts.

Спасибо вам за ваше время, и я надеюсь, что кто-нибудь сможет мне немного помочь.

 ####START CODE#####
set.seed(987654321)

library(spatcounts)

AllData <- cbind(sim.Yin, sim.region, sim.fm.X)
colnames(AllData)[1:2] <- c("Yin", "Region")

idx <- sample(1:nrow(AllData), 100, replace=TRUE)
newdata.df <- AllData[idx,]
newdata.df <- newdata.df[order(newdata.df$Region),]
X <- as.data.frame(newdata.df[,3:4])
region <- as.data.frame(newdata.df$Region)
colnames(region) <- "V1"
Yin <- as.data.frame(newdata.df$Yin)

temp.idx <- sort(unique(newdata.df$Region))
g.temp <- sim.gmat[temp.idx,temp.idx]

TotalN <- rowSums(g.temp) ##CHECK IF THERE ARE ZEROS PRESENT IN THIS VECTOR.
##OTHERWISE RUN AGAIN UNTIL AT LEAST ONE ZERO

n.temp <- sim.nmat[temp.idx,]
n.temp$V2 <- temp.idx[match(n.temp$V2, temp.idx)]
n.temp$V3 <- temp.idx[match(n.temp$V3, temp.idx)]
n.temp$V4 <- temp.idx[match(n.temp$V4, temp.idx)]
n.temp$V5 <- temp.idx[match(n.temp$V5, temp.idx)]
n.temp$V6 <- TotalN

n.temp[is.na(n.temp)] <- 0

Yin.NB <- est.sc(Yin, ~ X[,1]   X[,2] -1, 
region, model="NB", g.temp, n.temp, totalit=10) ##ERROR