Fit a linear model to the residuals in a gofN object.
This non-method runs a properly weighted linear model on the raw residuals of a gofN simulation for a multi-network ERGM fit.
lm.gofN(formula, data, ...)
formula: an lm-style formula. See Details for interpretation.data: a gofN object....: additional arguments to lm(), excluding weights.A list of lm objects, one for each element of the vector on the LHS.
The formula's RHS is evaluated in an environment comprising the network statistics used in the gofN() call (which refer to the raw residuals for the corresponding statistic) and the network attributes.
The LHS is handled in a nonstandard manner, designed to make it easier to reference the usually lengthy network statistics: first, it is evaluated in the formula's environment. If the evaluation is successful and the result is numeric, these numbers are used as indices of the statistics in the gofN object to use on the RHS. If it is a character vector, it is treated as names of these statistics.
data(samplk) # Add time indices: samplk1 %n% "t" <- 1 samplk2 %n% "t" <- 2 samplk3 %n% "t" <- 3 monks <- Networks(samplk1, samplk2, samplk3) fit <- ergm(monks~N(~edges+nodematch("group"))) fit.gof <- gofN(fit) # GOF = original model # Is there a time effect we should incorporate? fit.gof.lm <- lm.gofN((1:2)~t, data=fit.gof) lapply(fit.gof.lm, summary)
gofN() and related methods.