Here is an example of something similar that you can use.
dt=Open("$SAMPLE_DATA/Cities.jmp");
Control Chart Builder(
Show Two Shewhart Charts( 0 ),
Variables( Subgroup( :Region ), Y( :"pop- m"n ) ),
Chart( Points( Statistic( "Average" ) ), Limits( Sigma( "Levey Jennings" ) ) ),
Column Switcher( :"pop- m"n, {:"pop- m"n, :POP, :Max deg. F Jan, :OZONE, :CO} ),
SendToReport(
Dispatch(
{},
"Region",
ScaleBox,
{Add Ref Line( 1, "Solid", "Blue", "", 1 ),
Add Ref Line( 2, "Solid", "Blue", "", 1 ),
Add Ref Line( 3, "Solid", "Blue", "", 1 )}
)
)
)