Plot a matrix of transfer functions
# S3 method for tfam
plot(x, xvar = NULL, yvar = NULL, mar = c(1.1, 1.1, 0.1, 0.1), ...)
an object of class 'tfam' (transfer function analysis matrix)
created using tfam_lambda
or tfam_inertia
.
(optional) the variables to plot on the x and y axes. May
be "p"
, "lambda"
or "inertia"
. Defaults to
xvar="p"
and yvar="lambda"
for objects created using
tfam_lambda
, and xvar="p"
and yvar="inertia"
for
objects created using tfam_inertia
.
the margin limits on the plots: see par
plot.tfam
plots matrices of transfer functions (class
tfam
) created using tfam_lambda
or
tfam_inertia
. The plot is laid out to correspond with
the nonzero entries of the matrix used to generate the transfer functions,
for easy visual comparison of how perturbation affects different matrix
elements.
Constructor functions: tfam_lambda
, tfam_inertia
# Create a 3x3 matrix
( A <- matrix(c(0,1,2,0.5,0.1,0,0,0.6,0.6), byrow=TRUE, ncol=3) )
# Calculate the matrix of transfer functions using default arguments
( tfmat<-tfam_lambda(A) )
# Plot the matrix of transfer functions
plot(tfmat)
# Create an initial stage structure
( initial <- c(1,3,2) )
# Calculate the matrix of transfer functions for inertia and
# specified initial stage structure using default arguments
( tfmat2<-tfam_inertia(A,vector=initial) )
# Plot the result (defaults to inertia ~ p)
plot(tfmat2)
# Plot inertia ~ lambda
plot(tfmat2, xvar="lambda", yvar="inertia")
Run the code above in your browser using DataLab