Implements the maximum likelihood algorithm of Ludwig and Titterington (1994)
titterington(x)a [9 x n] matrix with the following columns:
X, sX, Y, sY, Z, sZ, rhoXY, rhoXZ, rhoYZ.
a four-element list of vectors containing:
4-element vector c(a,b,A,B) where a is
the intercept of the X-Y regression, b
is the slope of the X-Y regression, A
is the intercept of the X-Z regression, and
B is the slope of the X-Z regression.
[4 x 4]-element covariance matrix of par
the mean square of the residuals (a.k.a `reduced Chi-square') statistic
Ludwig, K.R. and Titterington, D.M., 1994. Calculation of \(^{230}\)Th/U isochrons, ages, and errors. Geochimica et Cosmochimica Acta, 58(22), pp.5031-5042.