// Competing Causes Interactive Survival Plot
x=1; y=1;
shape1 = 2.4; 	scale1 = log10(46); 
shape2 = 0.6; 	scale2 = log10(60); 

survival = expr((1-Weibull Distribution(x, shape1, 10^scale1, 0))*(1-Weibull Distribution(x, shape2, 10^scale2, 0)));

gBox = Graph Box(FrameSize(500,300),Double Buffer,
	XScale(0,200),yScale(0,1),
	pen color("green");pen size(3);
	YFunction(survival,x),
	pen color("blue");pen size(1);
	YFunction(1-Weibull Distribution(x, shape1, 10^scale1, 0),x),
	pen color("red");pen size(1);
	YFunction(1-Weibull Distribution(x, shape2, 10^scale2, 0),x));

plot = gBox[FrameBox(1)];
NewWindow("Competing Causes Survival Function",
	v list box (gBox,
		h list box (textbox("  Scale (Blue)  "),
			slider box (0,4, scale1,plot<<reshow()),
			textbox("  Shape (Blue)  "),
			slider box (0.1,3, shape1,plot<<reshow())),
		h list box (textbox("  Scale (Red)  "),
			slider box (0,4, scale2,plot<<reshow()),
			textbox("  Shape (Red)  "),
			slider box (0.1,3, shape2,plot<<reshow()))
		)
	);