Выборка из модели autologisti c с использованием rautologisti c в пакете ngspatial R - PullRequest
Новая
       52

Выборка из модели autologisti c с использованием rautologisti c в пакете ngspatial R

0 голосов
/ 13 июля 2020

Я пытаюсь сгенерировать случайный dr aws из центрированной модели autologisti c, используя функцию rautologistic в пакете ngspatial R. Код возвращает результат для определенных значений параметров (например, rho=0.2), но не работает (без сообщения об ошибке) для другого значения (rho=0.3).

Любые советы о том, почему нельзя использовать случайную выборку будет очень признателен, или как вызвать сообщение об ошибке, которое нарушает работу функции. Спасибо! 

library("igraph")
library("ngspatial")

n <- 100 # sample size
set.seed(1234)

# simulate graph
g <- sample_pa(n, power = 1, m = n/8, directed = FALSE)
a_mtx <- as.matrix(as_adjacency_matrix(g)); rm(g)
i_perm <- sample(n,n,replace=FALSE) ## randomly permute indices
a_mtx <- a_mtx[i_perm,i_perm]

# covariates
p <- 2
l <- matrix(rnorm(n=n*p),nrow=n,ncol=p)

gamma <- rep(1/p,p)
rho <- 0.2
stm <- tryCatch(
  out <- rautologistic(X=cbind(1,l), A=a_mtx, theta=c(0,gamma,rho)),
  error=function(cond) return(NA))

rho <- 0.3
stm <- tryCatch(
  out <- rautologistic(X=cbind(1,l), A=a_mtx, theta=c(0,gamma,rho)),
  error=function(cond) return(NA))
...