

 APACE20140106
 Various performance enhancements (by precomputing a number of matrix inversions, a substantial speedup was possible for the estimation of the parameters). Note that imagewise inference requires computation of the likelihood ratio statistic, which remains the slow (relatively); setting NoImg=1 will only provide the summary measures but will offer the fastest possible performance.
 APACE20140104
 Updated to "LRSD_ACE_FAM.m" have better handling of exceptional badly behaved elements/voxels.
 APACE20140102
 More normalisation options added, both for usual heritability and aggregate heritability ("Steve's method"). In particular, the default suggested "WK_ParallelOutline.m" how runs aggregate heritability two ways (with and without variance normalisation) and APACEsummary.m has been updated to show these results. Also minor changes to "LRSD_ACE_FAM.m" to suppress occasional warnings of rank deficient warnings.
 APACE20131227
 APACEsummary.m now gracefully works with all possible configurations of ACEfit_Par.NoImg and ACEfit_Par.Vs.ClustInf (i.e. when imagewise computation and/or clustersize inference are/are not used).
 APACE20131224
 Steve's estimate was noisily reported by APACEsummary (due to the point estimate being estimated as the center of the CI instead of directly saved). This is now fixed.
 APACEsummary.m now reads the Voxelwise FWE Pvalue from Pvals_MaxLRT_h2.mat instead of recomputing it.
 APACE20131220
 FDR maps no longer thresholded. The FDR image, ACE_A_LRT_vox_FDRP.nii, was unintentionally thresholded before writing out the image, giving apparently empty (i.e. allzero) images if no voxels were significant at FDRalpha=0.05. Now the FDR image shows, as intended, all FDRcorrected Pvalues after a log10 transformation.
 Documentation update: log10 Pvalue images. The documentation (README_APACE_outputs.docx) failed to mention that the images reporting Pvalues were actually log10 Pvalues, with the exception of the two asymptotic Pvalue image (see next). The documentation has been updated.
 Asymptotic Pvalues images now log10 P: Inconsistent with all other Pvalue images, the ACE_A_LRT_vox_Pasym.nii & ACE_C_LRT_vox_Pasym.nii images were written as Pvalues, instead of log10 Pvalues. Now these images are log10 P.
 Parallel output name change. Previously, parallel results were written out with spaces in the names, e.g. "ACEfit_Parallel_ 3.mat". This was not intended, and has now been fixed, e.g. "ACEfit_Parallel_0003.mat". It should be backwards compatible (i.e. read the results with spaces); let me know if not.
 Voxelwise FWE Pvalue was reported in the PDF figure but not saved; now it is saved in Pvals_MaxLRT_h2.mat.
