Plot function for a list of sienaFit objects
Draws a funnel plot for a list of sienaFit objects that all have estimated the same parameter.
funnelPlot(anslist, k, threshold=NULL, origin=TRUE,
plotAboveThreshold=TRUE, verbose=TRUE, ...)anslist |
A list of object of class |
k |
The number of the parameter to be plotted. |
threshold |
threshold for standard errors:
all estimations where the standard error for parameter
|
origin |
Boolean: whether to include the origin in the plot, if all estimates have the same sign. |
plotAboveThreshold |
Boolean: whether to include the estimates
for which the standard error is larger than |
verbose |
Boolean: whether to report in the console all estimates
omitted, because either their standard error is larger than
|
... |
For extra arguments (passed to |
The function funnelPlot plots estimates against standard
errors for a given effect k, with red
reference lines added at the two-sided
significance threshold 0.05. Effects for which a score test
was requested are not plotted (and reported to the console
if verbose).
If not all effects with number k are the same
in all sienaFit objects, a warning is given.
The effect name for the first object is used as the plot title.
Another funnel plot is available as print.sienaMeta.
The two-column matrix of
values of the plotted points is invisibly returned.
Tom Snijders
# A meta-analysis for three groups does not make much sense. # But using three groups shows the idea. Group1 <- sienaDependent(array(c(N3401, HN3401), dim=c(45, 45, 2))) Group3 <- sienaDependent(array(c(N3403, HN3403), dim=c(37, 37, 2))) Group4 <- sienaDependent(array(c(N3404, HN3404), dim=c(33, 33, 2))) dataset.1 <- sienaDataCreate(Friends = Group1) dataset.3 <- sienaDataCreate(Friends = Group3) dataset.4 <- sienaDataCreate(Friends = Group4) OneAlgorithm <- sienaAlgorithmCreate(projname = NULL, nsub=1, n3=50, seed=123) effects.1 <- getEffects(dataset.1) effects.3 <- getEffects(dataset.3) effects.4 <- getEffects(dataset.4) ans.1 <- siena07(OneAlgorithm, data=dataset.1, effects=effects.1, batch=TRUE) ans.3 <- siena07(OneAlgorithm, data=dataset.3, effects=effects.3, batch=TRUE) ans.4 <- siena07(OneAlgorithm, data=dataset.4, effects=effects.4, batch=TRUE) funnelPlot(list(ans.1, ans.3, ans.4), k=2) funnelPlot(list(ans.1, ans.3, ans.4), k=2, origin=FALSE)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.