leifeld
commented
2 years ago Thanks. I can confirm there is something off. Here is a minimal example:
library("network")
library("btergm")
set.seed(5)
networks <- list()
for (i in 1:10) { # create 10 random networks with 10 actors
mat <- matrix(rbinom(100, 1, .25), nrow = 10, ncol = 10)
diag(mat) <- 0 # loops are excluded
nw <- network(mat) # create network object
networks[[i]] <- nw # add network to the list
}
covariates <- list()
for (i in 1:10) { # create 10 matrices as covariate
mat <- matrix(rnorm(100), nrow = 10, ncol = 10)
covariates[[i]] <- mat # add matrix to the list
}
fit <- btergm(networks ~ edges + istar(2) + edgecov(covariates), R = 100)
checkdegeneracy(fit)The part that causes the problem is the following command in Line 119 in the checkdegeneracy.R file:
target.stats[[index]] <- summary(statnet.common::filter_rhs.formula(form),
response = NULL)I will look into this next time I work on the package.
daevra
Hello,
I'm currently working on some research using the bterg model, and am encountering an issue when I try to run a degeneracy check on my btergm object. I've tried finding a solution elsewhere, but to no avail. I'm fairly new to programming so I'm not entirely sure if the problem is on my end or with the program itself. So far, my code runs smoothly for estimating the btergm and for checking gof. However, once I attempt to run the degeneracy check on the btergm object using the code in the documentation:
checkdegeneracy( object, nsim = 1000, MCMC.interval = 1000, MCMC.burnin = 10000, verbose = FALSE )
I receive the following error code:
Is there something that I'm missing? I've double checked to make sure that the object is a btergm object and it is. I've also tried with different data and still experienced the same issue.
Thank you so much for the help and your time!