matrix

Calculate a matrix of the specified type from input coordinates.

matrix [out <filename>] [start <#>] [stop|end <#>] [offset <#>] [name <name>] [ byatom | byres [mass] | bymask [mass] ] [ ired [order <#>] ] [ {distcovar | idea} <mask1> ] [ {dist | correl | covar | mwcovar} <mask1> [<mask2>] ] [ dihcovar dihedrals <dataset arg> ]

[out <filename>] If specified, write matrix to <filename>.
[start <#>] [stop|end <#>] [offset <#>] Start, stop, and offset frames to use (as a subset of all frames read in).
[name <name>] Name of the matrix dataset (for referral in subsequent analysis).
byatom Write results by atom (default). This is the sole option for covar, mwcovar, and ired.
byres Write results by calculating an average for each residue (mass weighted if mass is specified).
bymask Write average over <mask1>, and if <mask2> is specified <mask1> x <mask2> and <mask2> as well (mass weighted if mass is specified).

Calculate matrix of the specified type from input coordinate frames:
dist <mask1> [<mask2>] Distance matrix (default).
correl <mask1> [<mask2>] Correlation matrix (aka dynamic cross correlation).
covar <mask1> [<mask2>] Coordinate covariance matrix.
mwcovar <mask1> [<mask2>] Mass-weighted coordinate covariance matrix.
distcovar <mask1> Distance covariance matrix.
idea <mask1> Isotropically Distributed Ensemble Analysis matrix.
ired [order <#>] Isotropic Reorientational Eigenmode Dynamics matrix with Legendre polynomials of specified order (default 1). IRED vectors must have been specified previously with ’vector ired’.
dihcovar dihedrals <dataset arg> Dihedral covariance matrix. Dihedral data sets must have been previously defined with e.g. dihedral or multidihedral commands or read in externally with readdata and marked as dihedrals.

Matrix dimensions will be of the order of N x M for dist, correl, idea, and ired, 2N x 2N for dihcovar, 3N x 3M for covar and mwcovar, and N(N-1) x N(N-1) / 4 for distcovar (with N being the number of data sets in the case of ired and dihcovar and the number of atoms in <mask1> otherwise, and M being the number of atoms in <mask2> if specified or <mask1> otherwise). No mask is required for ired; the matrix will be made up of previously defined

IRED vectors (see the vector command). Similarly no mask is required for dihcovar; dihedral data sets must have been previously defined. Only one mask can be used with distcovar and idea matrices (i.e. they can be symmetric only), otherwise one or two masks can be used (for symmetric and full matrices respectively). If two masks are specified the number of atoms covered by mask1 must be greater than or equal to the number of atoms covered by mask2, and on output <mask1> corresponds to columns while <mask2> corresponds to rows.

As a simple example, a distance matrix of all CA atoms is generated and output to distmat.dat.

matrix dist @CA out distmat.dat