View source: R/sienaDataCreate.r
sienaDataCreate | R Documentation |
Creates a Siena data object from input dependent variables (networks and possibly behavioural variables), covariates, and composition change objects.
sienaDataCreate(..., nodeSets=NULL, getDocumentation=FALSE)
... |
objects of class |
nodeSets |
list of Siena node sets. Default is the single node
set named |
getDocumentation |
Flag to allow documentation of internal functions, not for use by users. |
The function checks that the objects fit, that there is at least one dependent
variable, and adds various attributes to each variable describing the
data. If there is more than one nodeSet they must all be specified.
Function print01Report
will give a basic description of the data
object and is a check useful, e.g., for diagnosing problems.
An object of class siena
which is designed to be used in a siena
model fit by siena07
. The components of the object are:
nodeSets |
List of node sets involved |
observations |
Integer indicating number of waves of data |
depvars |
List of networks and behavior variables |
cCovars |
List of constant covariates |
vCovars |
List of changing covariates |
dycCovars |
List of constant dyadic covariates |
dyvCovars |
List of changing dyadic covariates |
compositionChange |
List of composition change objects corresponding to the node sets |
Ruth Ripley
See https://www.stats.ox.ac.uk/~snijders/siena/
sienaDependent
, coCovar
, varCovar
,
coDyadCovar
, varDyadCovar
,
sienaNodeSet
, sienaCompositionChange
,
sienaGroupCreate
,
sienaDataConstraint
, sienaNodeSet
,
print01Report
mynet <- sienaDependent(array(c(s501, s502, s503), dim=c(50, 50, 3)))
mybeh <- sienaDependent(s50a, type="behavior")
mydata <- sienaDataCreate(mynet, mybeh)
# This gives the same result as
mydata <- sienaDataCreate(list(mynet, mybeh))
## And for a two-mode network
mynet1 <- sienaDependent(array(c(s501, s502), dim=c(50, 50, 2)), nodeSet="senders")
senders <- sienaNodeSet(50, nodeSetName="senders")
receivers <- sienaNodeSet(30, nodeSetName="receivers")
mynet2 <- sienaDependent(array(c(s501[,1:30], s502[,1:30]), dim=c(50, 30, 2)),
nodeSet=c("senders", "receivers"))
(mydata <- sienaDataCreate(mynet1, mynet2, nodeSets=list(senders, receivers)))
## Not run:
print01Report(mydata, modelname = "mydescription")
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.