simstats0c {RSiena} | R Documentation |
The function to be called as 'FRAN'. Calls compiled C++.
simstats0c(z, x, INIT = FALSE, TERM = FALSE, initC=FALSE, data=NULL, effects=NULL, fromFiniteDiff=FALSE, profileData=FALSE, prevAns=NULL, returnDeps=FALSE)
z |
control object, passed in automatically in Siena07 |
x |
model object, passed in automatically in Siena07 |
INIT |
if TRUE, do initial processing. May be required to set up z |
TERM |
if TRUE, do end processing. |
initC |
if TRUE, call is to setup the data and model in C++. For use with multiple processes only |
data |
A siena object |
effects |
list of data frames as returned by getEffects |
fromFiniteDiff |
Boolean used during calculation of derivatives by finite differences. Not for user use. |
profileData |
boolean to force dumping of the data for profiling
with sienaProfile.exe |
prevAns |
An object of class "sienaFit", from which scaling information (derivative matrix and standard deviation of the deviations) will be extracted along with the latest version of the parameters which will be used as the initial values, unless the model requests the use of standard initial values. Phase 1 will then be omitted. |
returnDeps |
Boolean. Whether to return the simulated networks in Phase 3. Currently only returns the first network. |
This should be used for the element FRAN of the model object. The arguments with no defaults must be passed in on the call to siena07.
Depends on the call. If INIT
or initC
or
TERM
are true, returns z
, the control object. Otherwise,
returns a list containing:
fra |
Simulated statistics |
sc |
Scores with which to calculate the derivative (not phase 2 or if using finite differences or maximum likelihood) |
dff |
For maximum likelihood, 2nd deriv?? |
ntim |
For conditional processing, time taken |
feasible |
currently set to TRUE |
OK |
could be set to FALSE if serious error has occurred |
nets |
A list of simulated networks, as a list of edgelists for
each period. Only if returnDeps is TRUE. |
Ruth Ripley
See http://www.stats.ox.ac.uk/~snijders/siena/
mynet1 <- sienaNet(array(c(tmp3,tmp4),dim=c(32,32,2))) mydata <- sienaDataCreate(mynet1) myeff<- getEffects(mydata) myeff[myeff$effectName=='transitive triplets'& myeff$type=='eval','include']<- TRUE mymodel<- model.create(findiff=TRUE, fn = simstats0c) ans<- siena07(mymodel, data=mydata, effects=myeff, batch=TRUE)