If your datasets are small, then it is not difficult to "roll your own" using optim(). If you look at http://www.qimr.edu.au/davidD/sib-pair.R you will find such a routine at lines 1511-1562. This calls kinship.rel() that produces NRM etc.
-- | David Duffy (MBBS PhD) ,-_|\ | email: [EMAIL PROTECTED] ph: INT+61+7+3362-0217 fax: -0101 / * | Epidemiology Unit, Queensland Institute of Medical Research \_,-._/ | 300 Herston Rd, Brisbane, Queensland 4029, Australia GPG 4D0B994A v ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
