cross-reference k parameter in two P_Sequences

9 views
Skip to first unread message

Erik Marsh

unread,
Oct 6, 2025, 6:06:54 PMOct 6
to OxCal
Hi everyone –
I have cores from the same lake with similar dates (but at different depths).I am modeling each as its own P_Sequence but I want to use the same k parameters for both (a possibility mentioned in the manual). This is the code I tried, and it runs, but I the results don't seem right. I can't see the raw data in the second core. I see depth model results, but only for core A. Where can I check to see what k-parameters are being applied to core B? Is that in Model Spec somewhere? I'm guessing I need to code this different – any suggestions?

//core A, depths 1420–1400 cm
P_Sequence("k-parameter",1,2,U(-2,2))
  { ... };
//core B, depths 900–850 cm
P_Sequence("=k-parameter")
  { ... };

Thanks! – Erik


MILLARD, ANDREW R.

unread,
Oct 7, 2025, 4:12:27 AMOct 7
to ox...@googlegroups.com

Hi Erik,

 

I think you can either explicitly model k0, or the modifier of k0.

 

This models k0 as the same for both sequences using the prior you had:

 

Plot()

{

  k0 = 2*10^U(-2,2);

  P_Sequence("Core A",1,"=k0")

  {

   Boundary();

   R_Simulate("DateA1",1100,25)

   {

    z=900;

   };

   R_Simulate("DateA2",1150,25)

   {

    z=850;

   };

   Boundary();

  };

  P_Sequence("Core B",1,"=k0")

  {

   Boundary();

   R_Simulate("DateB1",1200,25)

   {

    z=800;

   };

   R_Simulate("DateB2",1250,25)

   {

    z=750;

   };

   Boundary();

  };

};

 

While this models the modifier. It gives a more familiar result as it is usually the modifier that is plotted when k is allowed to vary.

 

Plot()

{

  k_expr = U(-2,2);

  P_Sequence("Core A",1,2,"=k_expr ")

  {

   Boundary();

   R_Simulate("DateA1",1100,25)

   {

    z=900;

   };

   R_Simulate("DateA2",1150,25)

   {

    z=850;

   };

   Boundary();

  };

  P_Sequence("Core B",1,2,"=k_expr ")

  {

   Boundary();

   R_Simulate("DateB1",1200,25)

   {

    z=800;

   };

   R_Simulate("DateB2",1250,25)

   {

    z=750;

   };

   Boundary();

  };

};

 

 

 

Best wishes 

Andrew 

-- 

Prof. Andrew Millard 

Department of Archaeology,

Durham University, UK 

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

Personal page: https://www.durham.ac.uk/staff/a-r-millard/ 

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

 

Image

 

From: ox...@googlegroups.com <ox...@googlegroups.com> On Behalf Of Erik Marsh
Sent: 06 October 2025 23:07
To: OxCal <ox...@googlegroups.com>
Subject: cross-reference k parameter in two P_Sequences

 

[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 visit https://groups.google.com/d/msgid/oxcal/454922d0-dc4d-4171-b622-c4320231ff46n%40googlegroups.com.

Erik Marsh

unread,
Oct 7, 2025, 9:20:31 AMOct 7
to OxCal
Thanks Andrews, helpful as always! – Erik
Reply all
Reply to author
Forward
0 new messages