# Twins

30 views

### Ray Kidd

Feb 13, 2022, 5:07:49 AMFeb 13
Hello Cristopher and Group,

I am assisting with an interesting problem.  Perhaps there is already a solution.  I would like comments on my suggestion.

As a result of DNA analysis, two skeletons are identified as twins, 0-2 months old at death.  Radiocarbon determinations are:

R_Date("Skeleton 1", 3347, 22);
R_Date("Skeleton 2", 3258, 23);

Using R_Combine or Combine shows a failed W&W1978 test.  This is hardly surprising as the mean dates are about 100 years apart.

So we might assume that one or the other is wrong ( I favour a mistaken second digit)

The question is, how to identify the one that is wrong and which one is correct.

My solution is to create a single Phase model using R_Simulate for six dates.  I place three before an unbounded Phase  pair of the twins and three following (this is not strictly  necessary, I like things to balance).  For the R_Simulate I make each one the same, i.e. the mean value of the twins R_Dates roughly.  I have used R_Simulate -1600, 25.

The model runs quite well.  Running many times gives a new model at each run, but, here the thing, it consistently shows Skeleton1 with a higher individual agreement index.  From which I take it that Skeleton 1 is the 'good' date and Skeleton 2 should be discarded in any future model.

Plot()
{
Sequence()
{
Boundary("stsim");
Phase("Phsim")
{
R_Simulate("Sim1", -1600, 25);
R_Simulate("Sim2", -1600, 25);
R_Simulate("Sim3", -1600, 25);
Phase("1")
{
R_Date("Skeleton 1P", 3347, 22);
R_Date("Skeleton 2P", 3258, 23);
};
R_Simulate("Sim6", -1600, 25);
R_Simulate("Sim4", -1600, 25);
R_Simulate("Sim5", -1600, 25);
};
Boundary("endsim");
};
};

Best wishes

Ray

### MILLARD, ANDREW R.

Feb 14, 2022, 4:44:39 AMFeb 14

Hello Ray,

Your model compares the two dates to six dates with a true calendar date of -1600. Whichever of the two has a higher likelihood at -1600 will give the better agreement. It tells you nothing about the quality of the dates, just how they compare to the assumptions that went into the model. There’s no way to discriminate between the two dates on the information you have given.

Best wishes

Andrew

--

Dr. Andrew Millard

Associate Professor of Archaeology,

Durham University, UK

Email: A.R.M...@durham.ac.uk

Personal page: https://www.dur.ac.uk/directory/profile/?id=160

Scottish Soldiers Project: https://www.dur.ac.uk/scottishsoldiers

Dunbar 1650 MOOC: https://www.futurelearn.com/courses/battle-of-dunbar-1650

From: 'Ray Kidd' via OxCal <ox...@googlegroups.com>
Sent: 13 February 2022 10:08
Subject: Twins

[EXTERNAL EMAIL]

--
You received this message because you are subscribed to the Google Groups "OxCal" group.
To unsubscribe from this group and stop receiving emails from it, send an email to oxcal+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/oxcal/768128813.667964.1644746863135%40mail.yahoo.com.

### Erik Marsh

Feb 14, 2022, 8:47:59 AMFeb 14
to OxCal
I have to agree with Andrew. This is far from ideal, but perhaps you could still combine them with Ward and Wilson's (1978) sum of squares and use that, despite the failed chi-square test. Doing this would use the logic that that the dates are equally wrong and something in the middle is most likely. That's about all I can see you could try. But it's not a very solid approach – it seems quite plausible, as you suggest, that one date is good and the other isn't.

Sorting this out would require other contextual dates, I think. Or something from the lab on the chemistry of the two samples. Or find the typo in the dates.
Erik

### Ray Kidd

Feb 14, 2022, 10:24:41 AMFeb 14
Hello Erik, Andrew,

I'm aware that it all hinges on the Simulated likelihoods and they change at every run for aa given Date.   If the six simulated likelihoods and the two Radiocarbon likelihoods bias towards one 'lot' of phase calculations then it will favour one Radiocarbon date over the other for that run.  It will continue that bias in the the majority of runs, but not all.

I  think it was just fortuitous that my rough estimate of -1600 happened on a close choice.

R_Combine and Combine, if only one date is correct give a wrong answer (But precisely)

regards

Ray

### Pavol Hnila

Feb 14, 2022, 11:53:33 AMFeb 14