updateTheta {RSiena} | R Documentation |
A function to update the initial values of theta, and a function to update an effects object.
Description
updateTheta
copies the final values of any matching selected effects
from a sienaFit
object to a Siena effects object.
updateSpecification
includes in a Siena effects object
a set of effects that are included in another effects object.
Usage
updateTheta(effects, prevAns, varName=NULL)
updateSpecification(effects.to, effects.from,
effects.extra=NULL, name.to=NULL, name.from=NULL)
Arguments
effects |
Object of class |
prevAns |
|
varName |
Character string or vector of character strings;
is this is not |
effects.to |
Object of class |
effects.from |
Object of class |
effects.extra |
Object of class |
name.to |
Character string, name of dependent variable in |
name.from |
Character string, name of dependent variable in
|
Details
The initial values of any selected effects
in the input effects object which match an effect estimated in
prevAns
will be updated by updateTheta
.
If the previous run was conditional, the estimated rate parameters for
the dependent variable on which the run was conditioned are added to
the final value of theta. If varName
is not NULL
,
this update is restricted to effects for the dependent variable/s
specified by varName
.
By updateSpecification
, the effects included in effects.from
are also included in effects.to
; if name.to
and/or
name.from
is specified, this is restricted to effects
for those dependent variables.
If effects.from
contains interaction effects, the corresponding
main effects will be looked for in effects.from
; if they are not
found there, they will be looked for in effects.extra
.
It is not guaranteed that this will be successful.
For effects.extra
, it is best to use an effects object constructed for
the same data set as effects.from
, and by the same version of RSiena.
Value
Updated effects object.
Note
Using updateTheta
explicitly before calling siena07
rather than using it via the argument
prevAns
of siena07
will not permit the use of
the previous derivative matrix.
In most cases, using siena07
with prevAns
will be more
efficient.
Author(s)
Ruth Ripley, Tom A.B. Snijders
References
See https://www.stats.ox.ac.uk/~snijders/siena/
See Also
Examples
## For updateTheta:
mynet1 <- sienaDependent(array(c(tmp3, tmp4), dim=c(32, 32, 2)))
mydata <- sienaDataCreate(mynet1)
myeff1 <- getEffects(mydata)
myeff1 <- includeEffects(myeff1, transTrip)
myeff1 <- includeInteraction(myeff1, recip, inPop)
myalgorithm <- sienaAlgorithmCreate(nsub=1, n3=100, projname=NULL)
ans <- siena07(myalgorithm, data=mydata, effects=myeff1, batch=TRUE)
ans$theta
(myeff <- updateTheta(myeff1, ans))
##
## For updateSpecification:
myeff2 <- getEffects(mydata)
myeff2 <- includeEffects(myeff2, outAct)
updateSpecification(myeff2, myeff1)