Updates and changes to WOMBAT

This page tries to keep track of what has been changed in WOMBAT and when — including `bug' fixes.

John K. noted that fitting correlated random effects did not work for effects with the GIN option; fixed.

Changed rules for sire model analyses: require “subject” identifier in data & model part of parameter file.

Romain D. reported a problem for an analysis with many dependencies among fixed effects – increased limit and improved fault trapping for exceeding the maximum number.

Added checks for sorting of data file; new SPECIAL option CLONES. Analyses with very low (<10% but >0%) proportions of repeated records for a trait are no longer accepted.

Changes made last year caused WOMBAT to no longer recognise the DIAG option for random effect covariance matrices; fixed. Option ––choozhz was only recognised for animal model analyses; now enabled for sire models as well.

Christie I. reported that the COVZER option did not work properly for multivariate analyses with numbers of covariates differing between traits. Klint G. found that the Windows versions had problems if data and pedigree file were the same. Both fixed.

Added facility to read in “pre-defined” inbreeding coefficients from the pedigree file.

Fixed bug which affected convergence of reduced rank analyses under a sire model with known maternal grand sires.

Fixed bug which caused WOMBAT to ignore some of the effects chosen to be set to zero through the ZEROUT directive, and improved `automatic' identification of dependencies among fixed effects somewhat.

WOMBAT now stops if it finds repeated records for a trait for less than 20% of individuals, under the assumption that this is due to careless data preparation. For the rare cases were this is a planned scenario, this can be overridden by a new 'special' option (see section 4.8 of the manual). In addition, WOMBAT now checks for duplicate identities in the pedigree file and stops if these give different parental identities.

When writing out solutions for random effects for multivariate analyses using the “PC” option, WOMBAT now produces additional files with corresponding solutions transformed `back' to the original scale.

WOMBAT now checks for values of zero for covariables fitted. If the data contains such values, it expects a 'special' option (see section 4.8 of the manual) telling it whether to ignore such records or treat them as valid values for the covariable.

Added facilities to a) fix selected covariance matrices at their `input' values, b) replace the “1” in design matrices (standard uni- and multivariate analyses only) by a value (weight) obtained from the data file.

Warren S. pointed out that WOMBAT would not allow correlated random effects in conjunction with a general relationship matrix. Changed - WOMBAT now allows both NRM and GIN options for such effects.

Added option −−mmeout — this is like −−solvit but writes out the complete mixed model equations to file, including solutions and mapping of effects fitted to equation numbers. This is the option for folks who want to calculate linear contrasts, F-ratios etc. for fixed effects.

Martin L. reported that WOMBAT would complain about an animal occurring as sire and dam for a sire model when a bull occurred as sire and MGS. Fixed.

New run option: −−solvit provides a BLUP analysis solving the mixed model equations iteratively.

Pia B. reported that WOMBAT had trouble coping with blank lines in the pedigree file. Fixed.

Added option for `stand-alone' matrix inversion.

Have changed (made faster) computation of average information matrix in reduced rank analyses - Cholesky factor of numerator relationship matrix no longer required or written out (Files cholnrm*.bin).

David H. reported that the Windows versions failed for multivariate, reduced rank analyses while the Linux versions did not. This turned out to be due to a minor bug in the g95 compiler. Both compiler and Windows versions have been fixed.

New features have been added :

1. WOMBAT can now carry out analyses involving subsets of traits in multivariate analyses, without the need to extracts subsets of the data file or renumber traits.
2. An option is provided to write out the parameter files for bivariate analyses of all pairs of traits, using the parameter file for the complete analysis as input
3. When encountering such analyses of subsets, abbreviated files of partial results are written out automatically. These can be used as input files when pooling results.
4. WOMBAT can now combine results from multiple analyses of subsets of traits into positive definite, pooled covariance matrices using a (weighted) iterative summation of expanded part matrices.

• Andrew S. reported that the Windows versons of WOMBAT crash when attempting to write out a “Curve_” file for a nested covariable - this seems to be a problem associated with a '*' in the file name.

• Occasionally WOMBAT hangs up or crashes during the set-up phase. This occurs mainly for very large data sets or models with a dense coefficient matrix in the mixed model equations, and can be attributed to 'failure' of the routine used to re-order the mixed model equations.
  WOMBAT employs three different public domain routines to carry out the ordering. Unfortunately, none of them is `bomb-proof'. In particular, the routine from the MeTiS package (which is used by default for large analyses) is prone to problems and does not provide the possibility of a graceful exit when this occurs. If this happens, please try one of the other ordering strategies or other MeTiS options.