Example for Probit Analysis

An engineer of aircraft windshields wants to investigate how well the windshields can withstand projectile impacts at a range of velocities. The engineer subjects a random sample of windshields to projectiles at one of eight velocities and records whether the windshields withstood the impact.

The engineer performs probit analysis to determine the range of velocities at which a certain percentage of the windshields will break when subjected to the projectile impact.

  1. Open the sample data, WindshieldStress.MTW.
  2. Choose Stat > Reliability/Survival > Probit Analysis.
  3. Select Response in event/trial format.
  4. In Number of events, enter Breaks.
  5. In Number of trials, enter N.
  6. In Stress (stimulus), enter Stress.
  7. From Assumed distribution, select Normal.
  8. Click OK.

Interpret the results

To evaluate the distribution fit, the engineer uses a significance level of 0.1. The goodness-of-fit p-values (0.977 and 0.975) are greater than the significance level, and the points on the probability plot fall along an approximate straight line. Therefore, the engineer can assume that the normal distribution model provides a good fit for the data.

To evaluate significant effects, the engineer uses a significance level of 0.05. Because the p-value for Stress (0.000) is less than the significance level (0.05), the engineer concludes that the velocity of the projectile does have a statistically significant effect on whether or not the windshield breaks.

The table of percentiles indicates that the engineer can be 95% confident that 1% of the windshields will fail at a velocity between 300.019 mph and 501.649 mph.

Distribution:   Normal

Response Information

VariableValueCount
BreaksEvent37
  Non-event52
NTotal89
Estimation Method: Maximum Likelihood

Regression Table

VariableCoefStandard
Error
ZP
Constant-6.203761.06565-5.820.000
Stress0.00895960.00156155.740.000
Natural       
Response0     
Log-Likelihood = -38.516

Goodness-of-Fit Tests

MethodChi-SquareDFP
Pearson1.1997260.977
Deviance1.2285860.975

Parameter Estimates



Standard
Error
95.0% Normal CI
ParameterEstimateLowerUpper
Mean692.41618.3649656.421728.410
StDev111.61219.451879.3167157.058

Table of Percentiles



Standard
Error
95.0% Fiducial CI
PercentPercentileLowerUpper
1432.76745.8542300.019501.649
2463.19241.0355345.266525.291
3482.49638.0450373.838540.427
4497.01835.8391395.242551.902
5508.83034.0781412.585561.304
6518.88432.6067427.289569.364
7527.69931.3403440.133576.480
8535.59230.2277451.589582.896
9542.77129.2352461.967588.771
10549.37928.3398471.482594.217
20598.48022.4304540.595636.280
30633.88619.4337587.639669.400
40664.13918.1881624.815700.723
50692.41618.3649656.409733.152
60720.69219.8068685.039768.545
70750.94522.4716713.104808.979
80786.35126.5977743.723858.524
90835.45333.3805783.926929.497
91842.06034.3538789.210939.174
92849.23935.4233794.925949.712
93857.13236.6126801.183961.326
94865.94837.9558808.140974.328
95876.00239.5048816.041989.192
96887.81441.3455825.2801006.70
97902.33543.6350836.5851028.27
98921.63946.7171851.5351057.03
99952.06551.6465874.9541102.50