Method, assumptions and interpretation
1. Define the factors and the meaning of a rating
A factor should have a stable definition shared by everyone supplying a judgment. In the fictional engine-room maintenance example, planning, spares, work readiness and execution are four distinct aspects of a process. Entry aᵢⱼ means the assessed direct influence from row factor i to column factor j. It is not the reverse influence, a correlation coefficient, a measured probability or a time-series causal effect. State the period and operating context before collecting ratings. Duplicate or overlapping factors can create artificial feedback. This workbench accepts a completed single matrix; it does not collect or validate expert opinions.
2. Keep the scale and the diagonal policy explicit
The declared scale is 0 to 4: zero means no perceived direct influence, followed by low, medium, high and very high judgments. Fractional values are allowed for an already-aggregated assessment; every cell is still required. Blank is not synonymous with no influence. The direct diagonal must be exactly zero, a modeling policy that excludes direct self-influence at entry. Indirect paths can nevertheless return to their starting point, so the total-relation diagonal is generally positive. Those feedback terms are retained in every row and column sum even though self-loops are omitted from the drawing.
3. Normalize with one named rule
Compute every row sum and every column sum of A. The divisor s is the maximum of all of them; D=A/s when s is positive. This is the max-row/max-column version specified in the cited 2011 paper, not the row-only alternative. Multiplying every rating by the same positive constant leaves D unchanged, provided the new ratings remain within the accepted scale. Consequently the normalized output does not preserve the absolute overall strength of uniformly rescaled judgments. Scaling one entry or one row can change s and therefore every normalized entry. The raw matrix and s are always shown to make this dependency inspectable.
4. Ask whether feedback actually decays
D describes direct propagation, D² describes two-step propagation and later powers describe longer walks, including repeated visits. The series D+D²+… exists when the spectral radius is below one. Merely bounding row and column sums by one does not guarantee that condition: the two-factor matrix with ratings 4 in both directions normalizes to an undiminished cycle. Its inverse formula is invalid. This implementation seeks an outward-rounded upper bound on the infinity norm of a power Dᵏ that is strictly below one. Such contraction certifies decay of subsequent blocks of powers for the stored matrix.
5. Verify the inverse before reporting T
The check examines powers 1, 2, 4 and so on up to 65536, with at most seventeen norm checks. It is a bounded numerical acceptance policy, so a valid but very slowly contracting model can be refused. After the convergence gate, pivoted elimination computes (I−D)⁻¹. Pivots below 10⁻¹², infinity-norm condition numbers above 10⁸, inverse residuals above 10⁻⁹ or scaled total-equation residuals above 10⁻¹² stop the calculation. No epsilon is added to the diagonal and no damping factor is inserted into the ratings. If the model is rejected, any revision must be explicit and justified by the person responsible for it.
6. Interpret totals without claiming causal discovery
The total matrix is T=D(I−D)⁻¹. Row sum Rᵢ is the model’s total outgoing influence from factor i; column sum Cᵢ is the total incoming influence to that factor. Prominence R+C summarizes involvement in the supplied network. Relation R−C describes a net transmitter or net receiver within that network. These names summarize an algebraic consequence of judgments; they do not show that changing a factor will cause a measured outcome. Near-zero relation signs can be numerical artifacts, especially when the inverse is ill-conditioned. A conditioning- and scale-aware guard marks these cases as numerically unresolved; raw values remain visible. This guard is a conservative display policy, not a confidence interval or rigorous error bound. A large prominence is not a maintenance priority, an importance weight, a risk estimate or a budget allocation. External evidence and a separate decision framework are needed for those uses.
7. Keep the graph separate from the calculation
The directed network displays an off-diagonal edge only when its total influence is strictly greater than τ. A high threshold gives a simpler drawing and may hide substantial connections; a low threshold can create a dense picture. The full T matrix, row and column totals, prominence and relation always use every entry regardless of τ. Node positions are a fixed circular layout, not a proximity measure. The second figure locates factors by prominence and relation, and its horizontal zero line separates positive and negative net relations. Both figures have complete numeric tables so the drawing is not the only way to inspect the result.
8. Handle absence of information honestly
If every entry is zero, s is zero and A/s is undefined. The explicit no-information branch reports zero D, zero T, identity inverse and zero summaries instead of dividing by zero or inventing a weak network. This branch does not imply that the real system has no relationships; it only records what the submitted matrix contains. The worked example is entirely synthetic and demonstrates a valid contracting system. Loading it restores all factor labels, the objective, sixteen ratings and τ=0.5. Export records the scale, normalizer, all intermediate matrices, checks and displayed edges, so a later reader can distinguish the model from its presentation.
s = max(maxᵢ Σⱼ aᵢⱼ, maxⱼ Σᵢ aᵢⱼ); D = A/s when s > 0
T = D + D² + D³ + … = D(I − D)⁻¹, only when ρ(D) < 1
Contraction certificate: ‖Dᵏ‖∞ ≤ upward-rounded bound < 1
κ∞(I−D) = ‖I−D‖∞ × ‖(I−D)⁻¹‖∞
Rᵢ = Σⱼ tᵢⱼ; Cᵢ = Σⱼ tⱼᵢ; prominence = Rᵢ+Cᵢ; relation = Rᵢ−Cᵢ
Numerical relation guard: δᵢ = max(10⁻¹², 128 ε κ∞(I−D) max(1,Rᵢ+Cᵢ)); |Rᵢ−Cᵢ| ≤ δᵢ is unresolved
Display edge i→j iff i ≠ j and tᵢⱼ > τ; all tᵢⱼ remain in calculations
| Symbol | Definition and unit |
|---|---|
| n; A=[aᵢⱼ] | Factor count; complete direct-influence rating matrix on the declared 0–4 scale |
| s; D | Maximum row/column sum; normalized direct matrix (dimensionless) |
| I; T=[tᵢⱼ] | Identity matrix; total-relation matrix (dimensionless model units) |
| ρ(D); k | Spectral radius; power used for the contraction certificate |
| ‖M‖∞; κ∞ | Maximum absolute row sum; infinity-norm condition number |
| Rᵢ; Cᵢ | Total outgoing and incoming influence, including diagonal feedback |
| ε; δᵢ | Machine epsilon 2⁻⁵²; conservative numerical relation-sign guard |
| τ | Display-only threshold in the same model units as T |
A complete synthetic maritime example
All names, ratings and performance figures below were invented for this lesson. They are not observations from a vessel, port or company.
Decision or study question: Explore fictional engine-room maintenance influences.
Factor summaries from the full T matrix
These are consequences of supplied perceived influence ratings. A positive relation is not proof of empirical causality, and prominence is not a safety priority or an importance weight.
| ↓ / → | Planning | Spares | Work readiness | Execution |
|---|---|---|---|---|
| Planning | 0 | 3 | 1 | 0 |
| Spares | 1 | 0 | 2 | 1 |
| Work readiness | 0 | 1 | 0 | 2 |
| Execution | 1 | 0 | 0 | 0 |
| Factor | Raw row sum | Raw column sum |
|---|---|---|
| Planning | 4 | 2 |
| Spares | 4 | 4 |
| Work readiness | 3 | 3 |
| Execution | 1 | 3 |
s = 4
| ↓ / → | Planning | Spares | Work readiness | Execution |
|---|---|---|---|---|
| Planning | 0 | 0.75 | 0.25 | 0 |
| Spares | 0.25 | 0 | 0.5 | 0.25 |
| Work readiness | 0 | 0.25 | 0 | 0.5 |
| Execution | 0.25 | 0 | 0 | 0 |
| ↓ / → | Planning | Spares | Work readiness | Execution |
|---|---|---|---|---|
| Planning | 1 | -0.75 | -0.25 | 0 |
| Spares | -0.25 | 1 | -0.5 | -0.25 |
| Work readiness | 0 | -0.25 | 1 | -0.5 |
| Execution | -0.25 | 0 | 0 | 1 |
| ↓ / → | Planning | Spares | Work readiness | Execution |
|---|---|---|---|---|
| Planning | 1.611510791 | 1.496402878 | 1.151079137 | 0.9496402878 |
| Spares | 0.690647482 | 1.784172662 | 1.064748201 | 0.9784172662 |
| Work readiness | 0.3741007194 | 0.6330935252 | 1.410071942 | 0.8633093525 |
| Execution | 0.4028776978 | 0.3741007194 | 0.2877697842 | 1.237410072 |
| ↓ / → | Planning | Spares | Work readiness | Execution |
|---|---|---|---|---|
| Planning | 0.6115107914 | 1.496402878 | 1.151079137 | 0.9496402878 |
| Spares | 0.690647482 | 0.7841726619 | 1.064748201 | 0.9784172662 |
| Work readiness | 0.3741007194 | 0.6330935252 | 0.4100719424 | 0.8633093525 |
| Execution | 0.4028776978 | 0.3741007194 | 0.2877697842 | 0.2374100719 |
| Factor | Outgoing R | Incoming C | Prominence R + C | Relation R − C | δᵢ | Interpretation |
|---|---|---|---|---|---|---|
| F1: Planning | 4.208633094 | 2.079136691 | 6.287769784 | 2.129496403 | 1.86166095e-12 | Net transmitter in supplied model |
| F2: Spares | 3.517985612 | 3.287769784 | 6.805755396 | 0.2302158273 | 2.015024323e-12 | Net transmitter in supplied model |
| F3: Work readiness | 2.28057554 | 2.913669065 | 5.194244604 | -0.6330935252 | 1.537893828e-12 | Net receiver in supplied model |
| F4: Execution | 1.302158273 | 3.028776978 | 4.330935252 | -1.726618705 | 1.282288206e-12 | Net receiver in supplied model |
| Check | Value | Acceptance rule |
|---|---|---|
| Power k | 2 | k ≤ 65536 |
| ‖Dᵏ‖∞ upper bound | 0.93750000000000056 | < 1 |
| κ∞(I−D) | 10.41726619 | ≤ 10⁸ |
| ‖(I−D)B−I‖∞ | 4.718447855e-16 | ≤ 10⁻⁹ |
| ‖B(I−D)−I‖∞ | 4.440892099e-16 | ≤ 10⁻⁹ |
| ‖T(I−D)−D‖∞ | 3.330669074e-16 | Raw residual |
| Scaled total-equation residual | 3.536768535e-17 | ≤ 10⁻¹² |
| Smallest absolute pivot | 0.8081395349 | ≥ 10⁻¹² |
B = (I−D)⁻¹. Scaled residual = raw residual / max(1, ‖T‖∞‖I−D‖∞ + ‖D‖∞). Outward rounding is used only for the contraction check; D and A are not altered.
Substitute the numbers
s = max(4, 4, 3, 1, 2, 4, 3, 3) = 4
D₁₂ = 3 / 4 = 0.75
T₁₁ = 0 × 1.611510791 + 0.75 × 0.690647482 + 0.25 × 0.3741007194 + 0 × 0.4028776978 = 0.6115107914
R₁ = 0.6115107914 + 1.496402878 + 1.151079137 + 0.9496402878 = 4.208633094
C₁ = 0.6115107914 + 0.690647482 + 0.3741007194 + 0.4028776978 = 2.079136691
R₁ + C₁ = 4.208633094 + 2.079136691 = 6.287769784; R₁ − C₁ = 2.129496403
Subscripts follow the input factor order; numbers are rounded for display.
| From | To | Tᵢⱼ |
|---|---|---|
| F1: Planning | F2: Spares | 1.496402878 |
| F1: Planning | F3: Work readiness | 1.151079137 |
| F1: Planning | F4: Execution | 0.9496402878 |
| F2: Spares | F1: Planning | 0.690647482 |
| F2: Spares | F3: Work readiness | 1.064748201 |
| F2: Spares | F4: Execution | 0.9784172662 |
| F3: Work readiness | F2: Spares | 0.6330935252 |
| F3: Work readiness | F4: Execution | 0.8633093525 |
| τ | Edges drawn | Calculation matrix |
|---|---|---|
| 0 | 12 | Full T; unchanged |
| 0.5 | 8 | Full T; unchanged |
| 1.496402878 | 0 | Full T; unchanged |
Tables display up to 10 significant digits as formatting, not guaranteed accuracy. Computation and CSV use unrounded JavaScript numbers. Numerical tolerance is not an operational acceptance criterion.
Exact result for this example: T = (1/139) × [[85,208,160,132], [96,109,148,136], [52,88,57,120], [56,52,40,33]].
What the result does and does not establish
- Direct-influence ratings are subjective inputs. The words transmitter and receiver refer only to this supplied model, not established empirical causes.
- Uniform rescaling cancels during normalization. The method cannot compare absolute overall influence between uniformly scaled matrices.
- The max-row/max-column normalizer does not guarantee convergence. Singular, nearly singular and numerically unverified inputs are rejected.
- The 0–4 scale, zero diagonal, bounded dimensions and numerical gates are declared implementation policies. No fuzzy, revised or damped DEMATEL variant is silently substituted.
- Correlation, experiments, domain evidence and temporal mechanisms require separate examination. An influence picture alone cannot approve a vessel operation or select a safe maintenance strategy.
Does raising the graph threshold reduce total influence?
No. It removes visible arrows only. For the exact example, τ=0.5 shows eight off-diagonal edges. The sixteen entries of T and all four factor summaries remain unchanged, including diagonal feedback.
Primary sources and implementation scope
Chen & Tzeng (2011), Assessment Model for Improving Educational Curriculum Materials Based on the DANP Technique with Grey Relational Analysis, CSEDU, 299–308. p.302, (1)–(2): max-row/max-column normalization and total relation; row/column summaries on the same page. Only its crisp DEMATEL component is implemented here, not DANP or grey relational analysis.
Lee et al. (2013), Revised DEMATEL: Resolving the Infeasibility of DEMATEL, Applied Mathematical Modelling 37, 6746–6757. Cited for the warning that convergence cannot be assumed; this tool does not implement the paper’s revised normalization.
Original teaching text and synthetic examples. Numerical limits, tie handling and the no-information branch are explicit implementation policies. dematel v1.0.0.