Nothing
set.seed(12344)
# suppressWarnings(RNGversion("3.4.0"))
net <- rbernoulli(4, .3)
model <- ergmito(net ~ edges)
gof0 <- gof_ergmito(model)
gof1 <- gof_ergmito(model, sim_ci = TRUE)
# Alternative calculation
prb <- exact_loglik(
x = model$formulae$stats_statmat[[1]],
params = coef(model),
stats_weights = model$formulae$stats_weights[[1]],
stats_statmat = model$formulae$stats_statmat[[1]]
)
prb <- exp(prb)
# Adds up to one
expect_equal(sum(prb*model$formulae$stats_weights[[1]]), 1)
min_max <- range(model$formulae$stats_statmat[[1]])
# Quantiles
ord <- order(model$formulae$stats_statmat[[1]])
cumprb <- cumsum(prb[ord]*model$formulae$stats_weights[[1]][ord])
# 90% CI
idx <- c(max(which(cumprb < .05)), min(which(cumprb > .95)))
expect_equivalent(
gof0$ci[[1]][, c("lower-q", "upper-q")],
model$formulae$stats_statmat[[1]][ord][idx]
)
# expect_equivalent(
# gof1$ci[[1]][, c("lower-q", "upper-q")],
# model$formulae$stats_statmat[[1]][ord][idx],
# tol = 1
# )
expect_equivalent(gof0$ci[[1]][, c("lower-p", "upper-p")], cumprb[idx])
set.seed(1244)
# suppressWarnings(RNGversion("3.4.0"))
net <- rbernoulli(4, .3)
model <- ergmito(net ~ edges + ttriad)
gof0 <- gof_ergmito(model)
# Alternative calculation
prb <- exact_loglik(
x = model$formulae$stats_statmat[[1]],
params = coef(model),
stats_weights = model$formulae$stats_weights[[1]],
stats_statmat = model$formulae$stats_statmat[[1]]
)
prb <- exp(prb)
# Adds up to one
expect_equal(sum(prb*model$formulae$stats_weights[[1]]), 1)
min_max <- range(model$formulae$stats_statmat[[1]])
# Quantiles
for (k in 1L:2L) {
ord <- order(model$formulae$stats_statmat[[1]][, k])
S <- sort(unique(model$formulae$stats_statmat[[1]][ord, k]))
# Aggregating probs
cumprb <- numeric(length(S))
for (i in seq_along(S)) {
ids <- which(model$formulae$stats_statmat[[1]][,k] == S[i])
cumprb[i] <- sum(model$formulae$stats_weights[[1]][ids]*prb[ids])
}
cumprb <- cumsum(cumprb)
# 90% CI
idx <- c(max(1, which(cumprb < .05)), min(which(cumprb > .95)))
expect_equivalent(
gof0$ci[[1]][, c("lower-q", "upper-q")][k,],
S[idx]
)
expect_equivalent(
gof0$ci[[1]][k, c("lower-p", "upper-p")],
cumprb[idx])
}
set.seed(1823)
net <- rbernoulli(c(3,4,5,3,4))
ans <- ergmito(net ~ ttriad + edges + mutual)
gof <- gof_ergmito(ans)
expect_output(print(gof), "Exact")
expect_silent(plot(gof))
set.seed(184223)
net <- rbernoulli(c(3,4,5,3, 5, 3))
ans <- ergmito(net ~ ttriad + edges)
# Adding a term
expect_silent(gof <- gof_ergmito(ans, GOF = ~ mutual + balance))
expect_output(print(gof), "Exact")
expect_silent(plot(gof))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.