The linear regression model assumes that the errors are independently and identically normally distributed. It also assumes that the range of the function is negative infinity to positive infinity. You could clamp the model prediction.
temp = model( X );
temp = Maximum( 0, temp );
Y = Minimum( temp, 1 );
Model(X) is the saved regression model.
Alternatively, you could model Logit( Y ); This treats Y as a response with a range of [0,1]. Your model will be saved with a series of column formulas. The one for probability is your response.