tests/nolatent.R

#  File tests/nolatent.R in package latentnet, part of the
#  Statnet suite of packages for network analysis, https://statnet.org .
#
#  This software is distributed under the GPL-3 license.  It is free,
#  open source, and has the attribution requirements (GPL Section 7) at
#  https://statnet.org/attribution .
#
#  Copyright 2003-2022 Statnet Commons
################################################################################
library(latentnet)
library(statnet.common)

data(sampson)

monks.nmr<-ergmm(samplike~nodematch("group")+rreceiver)
mcmc.diagnostics(monks.nmr)

print(summary(monks.nmr))
# Should produce a meaningful error message.
print(try(plot(monks.nmr)))

opttest({
  
monks.dnmr<-ergmm(samplike~nodematch("group",diff=TRUE)+rreceiver)
mcmc.diagnostics(monks.dnmr)
print(summary(monks.dnmr))

monks.dnmr2<-ergmm(samplike~nodematch("group",diff=TRUE)+rreceiver,prior=monks.dnmr$prior)
if(!all.equal(monks.dnmr2$prior,monks.dnmr$prior)) stop("Prior specification problem!")

monks.nm<-ergmm(samplike~nodematch("group"))
mcmc.diagnostics(monks.nm)
print(summary(monks.nm))

monks.dnm<-ergmm(samplike~nodematch("group",diff=TRUE))
mcmc.diagnostics(monks.dnm)
print(summary(monks.dnm))

# tests importing of ergm terms with local variable as inputs
set.seed(1)
nw <- samplike
n <- network.size(samplike)
covar <- matrix(rbinom(n^2, 1, 0.2), nrow=n)
covar.nw <- network(covar)
test3 <- ergmm(nw ~ euclidean(d = 2) + edgecov(covar)) 
test4 <- ergmm(nw ~ euclidean(d = 2) + edgecov(covar.nw)) 


}, "Some non-latent-space")

Try the latentnet package in your browser

Any scripts or data that you put into this service are public.

latentnet documentation built on May 11, 2022, 5:16 p.m.