In addition, WOMBAT produces a number of small ‘utility’ files. These serve to monitor progress during estimation, to carry information over to subsequent runs or to facilitate specialised post-estimation calculations by the user.
This file lists the covariance components defined by the model of analysis, together with their running numbers and starting values given. It is written out during the ‘set-up’ phase (see 5.2.3). It can be used to identify the running numbers needed when defining additional functions of covariance components to be evaluated (see 4.11)
This file gives a count of the numbers of repeated records per trait and, if option TSELECT is used, a count of the number of pairs taken at the same time.
Whenever WOMBAT encounters a set of parameters which improves the likelihood, the currently ‘best’ point is written out to the file BestPoint.
The first line of BestPoint gives the following information :
This is followed by the covariance matrices estimated.
N.B.: BestPoint is used in any continuation or post-estimation steps – do not delete is until the analysis is complete !
WOMBAT appends a line of summary information to the file Iterates on completion of an iterate of the AI, PX-EM or EM algorithm. This can be used to monitor the progress of an estimation run – useful for long runs in background mode. Each line gives the following information :
This small file gathers accumulates the number of non-zero elements in the Cholesky factor (or inverse) of the mixed model matrix together with the resulting operation count for the factorisation. This can be used to compare the efficacy of different ordering strategies for a particular analysis. The file contains one line per ordering tried, with the following information :
These files are written out when the AI algorithm is used. After each iterate, they give the average information matrix (not its inverse !) corresponding to the ‘best’ estimates obtained by an AI step, as written out to BestPoint. These can be used to approximate sampling variances and errors of genetic parameters.
N.B.: If the AI iterates are followed by further estimates steps using a different algorithm, the average information matrices given may not pertain to the ‘best’ estimates any longer.
AvInfoParms contains the average information matrix for the parameters estimated. Generally, the parameters are the elements of the leading columns of the Cholesky factors of the covariance matrices estimated. This file is written out for both full and reduced rank estimation.
For full rank estimation, the average information is first calculated with respect to the covariance components and then transformed to the Cholesky scale. Hence, the average information for the covariances is available directly, and is written to the file AvInfoParms.
Both files give the elements of the upper triangle of the symmetric information matrix row-wise. The first line gives the log likelihood value for the estimates to which the matrix pertains – this can be used to ensure corresponding files of estimates and average information are used. Each of the following lines in the file represents one element of the matrix, containing 3 variables :
N.B.: Written out are the information matrices for all parameters. If some parameters (or covariances) are not estimated (such as zero residual covariances for traits measured on different animals), the corresponding rows and columns may be zero.
For random regression analyses, file(s) with the basis functions evaluated for the values of the control variable(s) in the data are written out. These can be used, for example, in calculating covariances of predicted random effects at specific points.
The name of a file is equal to the name of the covariable (or ‘control’ variable), as given in the parameter file (model of analysis part), followed by the option describing the form of basis function (POL, LEG, BSP; see 4.9.2) the maximum number of coefficients, and the extension .baf. The file then contains one row for each value of the covariable, giving the covariable, followed by the coefficients of the basis function.
NB. These files pertain to the random regressions fitted! Your model may contain a fixed regression on the same covariable with the same number of specified regression coefficients, , but with intercept omitted. If so, the coefficients in this file are not appropriate to evaluate the fixed regression curve.
For analyses involving an additional parameter, the values used for the parameter and the corresponding maximum log likelihoods are collected in this file. This is meant facilitate estimation of the parameter through a quadratic approximation of the resulting profile likelihood curve via the run option --quapp. Note that this file is appended to at each run.
If analyses considering a subset of traits are carried out, WOMBAT writes out files EstimSubset.dat (see 7.2.6), to be used as input files in a run with option --itsum. In addition, for each run performed, this file name this appended to SubSetsList. This file contains one line per ‘partial’ run with two entries: the file name (EstimSubset.dat) and a weight given to the corresponding results when combining estimates. The default for the weight is unity.