ergm icon indicating copy to clipboard operation
ergm copied to clipboard

After dropping extreme values, MPLE and MLE likelihoods report different numbers of observations.

Open krivit opened this issue 7 months ago • 0 comments

suppressPackageStartupMessages(library(ergm))
data(sampson)

samplike.m <- as.matrix(samplike, matrix.type="adjacency")
samplike.m[4:10,4:10] <- 0

mple <- ergm(samplike~edges+edgecov(samplike.m))
#> Observed statistic(s) edgecov.samplike.m are at their greatest attainable values. Their coefficients will be fixed at +Inf.
mcmc <- ergm(samplike~edges+edgecov(samplike.m), control=control.ergm(force.main=TRUE, MCMLE.maxit=10))
#> Observed statistic(s) edgecov.samplike.m are at their greatest attainable values. Their coefficients will be fixed at +Inf.
waldo::compare(logLik(mple), logLik(mcmc), tolerance = 0.01)
#> `attr(old, 'nobs')`: 237
#> `attr(new, 'nobs')`: 306
#> 
#> `attr(old, 'br')` is absent
#> `attr(new, 'br')` is a list

Created on 2024-07-23 with reprex v2.1.1

NB: Here, the br attribute being different is as intended. The problem is with nobs.

krivit avatar Jul 23 '24 12:07 krivit