profileLikelihoods {RSienaTest} | R Documentation |
Calculates the profile likelihood for a Siena model for one or two parameters
profileLikelihoods(resp, x, data, effects, i, j = NULL, gridl = c(0.8, 1.2), seqlen = 5, maxit = 2, method = "BFGS", trace = 0, nIter = 100, ...)
resp |
Value returned from a call to |
x |
A |
data |
A |
effects |
A |
i |
Required. Index of parameter to maximise over. |
j |
Optional. Index of second parameter to maximise over. |
gridl |
Multiplier used to create grid for evaluation |
seqlen |
Number of points for evaluation. |
maxit |
Number of iterations of the optimizer |
method |
Method for |
trace |
trace argument for |
nIter |
Number of samples to use |
... |
Other arguments to be passed to the initialize call. |
A sienaFit
object, with added component zz
containing the points used for the plot.
Produces a plot.
Not for general use
Ruth Ripley
See http://www.stats.ox.ac.uk/~snijders/siena/
## Not run: mynet1 <- sienaDependent(array(c(tmp3, tmp4), dim=c(32,32,2))) mydata <- sienaDataCreate(mynet1) myeff <- getEffects(mydata) profileLikelihoods(list(theta=c(5,-2, 2)), sienaModelCreate(maxlike=TRUE), mydata, myeff, i=1) ## End(Not run)