You could use 'Fit Special' in 'Fit Y By X', then the approach of @txnelson:
NamesDefaultToHere(1);
// Example data
dt = Open("$SAMPLE_DATA/Big Class.jmp");
// Plot data
biv = dt << Bivariate(Y( :height ), X( :weight ), Invisible);
// Fit Special
myIntercept = 0;
biv << Fit Special( Intercept( myIntercept ));
// Save formula to table
biv << (Curve[1] << savePredicteds);
// Close Bivariate
Report(biv) << closeWindow;
// Get the formula
form = Column(dt, "Predicted Height") << getFormula;
// Express the formula as a function of x so 'YFunction()' will work
SubstituteInto(form, Expr(:weight), Expr(x));
// Use Graph Builder
gb = dt << Graph Builder(
Show Control Panel( 0 ),
Variables( X( :weight ), Y( :height ) ),
Elements( Points( X, Y, Legend( 3 ) ) )
);
// Add the graphics script
gbRep = Report(gb);
addScript = Expr(gbRep[FrameBox(1)] << addGraphicsScript(PenColor("Red"); YFunction(formTBD, x)));
SubstituteInto(addScript, Expr(formTBD), NameExpr(form));
addScript;