Directory Contents  3 Oct 95|  send questions or comments via e-mail to
------------------          |  bkirk@sbc.edu
 
text of a Ph. D. thesis, "Constructing Confidence Regions for a
	Composition of Tectonic Plate Rotations under Heteroscedasticity", 
	by B. Kirkwood
	filenames: text.ps (the whole thesis)
        or, in pieces, chap1.ps, chap2.ps, chap3.ps, chap4.ps, chap5.ps, app.ps



postscript files for figures in the thesis
		filenames:  fig*.plot


FORTRAN code for constructing the confidence regions described in the thesis:

add2.f    composition of two rotations, unequal kappas, using fits from 
	hellinger1 and all the data for both rotations

addplus.f   constructs conservative 95% confidence regions for a composition of
        two rotations using only summary statistics, not full data.  Valid
	if any of the rotations is a triple junction estimate employing a 
	single kappa but not valid with a triple-junction, unequal-kappas
	estimate.  Output from addplus can be used by addplus again.

add3.f   composition of three rotations, unequal kappas, full data

trikap.f  finds a triple junction fit assuming 3 different kappas on the
	branches

add2tj.f  confidence region construction for combining a triple junction,
	single kappa estimate with another rotation having a different kappa,
	using all data for both estimates
	(minor revision 10 Oct 1999)

addtj4k.f  combines a triple junction, 3 kappa estimate with another rotation,
	using all data for both estimates
	(minor revision 10 Oct 1999)

averotkap.f   averages independent estimates of the same rotation, allowing for
        unequal kappas, full data

averotplus.f   conservative 95% confidence region for the average of independent
	estimates of the same rotation, allowing for unequal kappas; uses only
	summary statistics provided by hellinger1

hellinger1.f   estimates a single tectonic plate rotation

hellinger3.f   estimates rotations at triple junctions

invrot.f   inverts a rotation using output from hellinger1

addrot.f  combines two rotations assuming equal kappas

conreg.f  produces a file of points on the boundary of the region of axes in a
	rotation confidence region.  Uses a rotation fit file from any of the
	programs above.  

subroutine libraries:

	numlib.f      
	(minor revision 10 Oct 1999)

	numlib2.f  (LAPACK routines)

	dtrans.f

	teclib.f 
	(minor revision 10 Oct 1999)

	bingham.f  (needed only by conreg.f)


pval.f    tests whether a rotation is the Identity

royer2.f  combines independent estimates of the same rotation (assumes large 
    samples or equal kappa's).  This program uses a first order approximation
    unlike averotkap.f and averotplus.f which use corrections to the first
    order approximation.