The complex model is
The model takes into account the successive increase of GM parameter along with the increase of source size, non-linear source influence, non-elastic damping, non-linear geometric scattering and influence of site effect at the station. The model can be customised by selecting any parts.
The ordinary one-stage regression method is used.
This application is prepared in 3 steps and the model can be improved in two of them.
Step by Step
- Preparation of the data - GM Parameters Catalog is needed.
To prepare GM Parameters Catalog 2 files are needed: the 1st input file is a Seismic Catalog and the 2nd input file is Ground Motion Catalog from the same Episode (AH EPISODES). Thus, choose a Catalog and GM Catalog from a selected episode. Then run Ground Motion Parameters Catalog builder, and as a result GM Parameters Catalog will be created.
The part of the GM Parameters Catalog can be used with Catalog filter. It is worth to remove weak ground motion eg. PGA <=0.03 m/s2.
- Selection of the model parameters (STEP 1 of analysis)
After the User adds the Application into his/her personal workspace, the following window appear on the screen (Figure 1), and the user is now requested to fill the fields shown in Figure 1:
GM Parameters: The user may choose from all ground motion parameters which were available in GM Catalog.
Source characteristic parameter: The user may choose from all source parameters which were available in Catalog. The model may contain M, M2, both or none of them. Recommended is to select both and check their statistical significance in results. When coefficients or one of them is insignificant, they can be removed from the model.
Path characteristic: The model may contain R, logR, both or none of them. Additionally, epicentral distance may be calculated with fixed depth ('h' parameter in the model). 'h' parameter can be estimated (then option Estimate depth with model should be chosen) or can be typed (the value is given in kilometres). Recommended is to select estimation of fix depth with model and choose both coefficients and check their statistical significance in results. When coefficients or one of them is insignificant, they can be removed from the model.
Site characteristic: Site characteristic depends on station position.
Figure 1. Input window of Ground Motion Prediction Equations: GMPE calculation
After defining the parameters mentioned above, the user shall click on the button (Figure 1), and the calculations are performed. The output will be created soon.
The output results include:gmpe_model_report
Figure 2. Ground Motion Prediction Equations model report.
R2 - (R-squared value) - coefficient of determination is the proportion of the variance in the dependent variable that is predictable from the independent variable(s). It shows in the percentage of which the model explains the variability in the response.
Root Mean Square Error (RMSE) - is a frequently used measure of the differences between values (sample or population values) predicted by a model or an estimator and the values observed.
p - probability o F-statistic for testing statistical significance of the model.
SEE chart (when 'h' is chosen as one of coefficient)
Figure 3. SEE plot.
gmpe_dataset.mat - mat file with chosen parameters based on GM Parameters Catalog for which GMPE is calculated,
gmpe_model_residuals.mat - mat file with residuals of the model,
gmpe_model_step1.mat - mat file with results of 1. step.
- Improving the model:
According to the result - the mainly statistical significance of coefficients and its physical consistency - the model can be improved. To improve the model, click on the .
In this case, the 'R' coefficient of distance will be removed because it is non-negative.
Figure 4. Input window of GMPE recalculation.
Figure 5. Result of recalculated model.
- As an additional result, after improving the model, the comparison of previous vs current model parameter appears.
The new (recalculated) model is better when 'SEE' and 'p' are smaller, and 'R2' and 'F' are bigger than in the previous model.
This model does not need improving. All coefficient are proper and statistically significant. The calculation can be done as many times until the results are satisfactory.
- Residual Analysis:
After receiving the satisfying model, the residual analysis can be performed (STEP2 of analysis). To start choose .
Figure 6. Residuals plots.
Figure 7. Input window of remove residuals part.
Figure 8. Results of the model after the removing outlieres.
Depending on the result - especially when scattering of residuals is evident - the analysis can be repeated. To repeat choose . Recalculation can be done as many times until the results are satisfactory.
After receiving the satisfying model, the final model can be performed (STEP3 of analysis). To calculate the last model with additional plots, please select .
The final output results include:
gmpe_model_raport, gmpe_model_dataset, figures (Normal probability plot of residuals, box plot - residuals according to the station, Mw vs residuals, r vs residuals).
There is also a possibility to configure selected figures based on used and estimated parameters.