Plots a "missmodel" object produced by data augmentation

DESCRIPTION:

Produces time series plots of parameter iterates which can help diagnose convergence of a data augmentation algorithm.

USAGE:

plot.missmodel(x, select = F, variables = (all), ask = T, ...) 

REQUIRED ARGUMENTS:

x
class "missmodel" object produced by a data augmentation algorithm.

OPTIONAL ARGUMENTS:

select
logical flag: if TRUE, you may select the variables (parameters) to plot. If FALSE, variables indicated by the variables argument are plotted.
variables
vector of indices of variables to plot. By default all are plotted.
ask
logical flag: if TRUE, you must hit the return key between pages. This allows you to inspect each graph when the graphics device does not stop automatically between pages.

SIDE EFFECTS:

For all or selected parameters: produces x-y plots of parameter vs. iteration number.

DETAILS:

This function is a method for the generic function plot for class "missmodel". It can be invoked by calling plot for an object of the appropriate class, or directly by calling plot.missmodel .

When select=T, rather than produce all plots sequentially, plot.missmodel displays a menu listing all the plots that can be produced. Corresponding to each plot is a number. You may supply one or more numbers. See examples below. To obtain a pause between pages of plots, set ask=T.

SEE ALSO:

, , .

EXAMPLES:

crime.DA <- daLoglin(crime, margins = count~Visit.1:Visit.2,  
                 control = list(save = 101:500)) 
plot(crime.DA) 
# plot the last 2 parameters 
plot(crime.DA, variables = 3:4) 
# plot the first 3 parameters, interactively
plot(crime.DA, select = T) 
Select Variable Indices (or 0 to exit): 
 1 Visit.1=1;Visit.2=1 
 2 Visit.1=2;Visit.2=1 
 3 Visit.1=1;Visit.2=2 
 4 Visit.1=2;Visit.2=2 
? 1:3 
# plot the first 2 parameters, then the fourth 
plot(crime.DA, select = T) 
Select Variable Indices (or 0 to exit): 
 1 Visit.1=1;Visit.2=1 
 2 Visit.1=2;Visit.2=1 
 3 Visit.1=1;Visit.2=2 
 4 Visit.1=2;Visit.2=2 
? c(1,2) 
Select Variable Indices (or 0 to exit): 
 1 Visit.1=1;Visit.2=1 
 2 Visit.1=2;Visit.2=1 
 3 Visit.1=1;Visit.2=2 
 4 Visit.1=2;Visit.2=2 
? 4