Maritime Science Life

Calculations with context

Crisp DEMATEL influence model

Normalize perceived direct influence by the maximum row/column sum; inspect convergence checks and total influence. The result does not establish empirical causality.

Scientific methods · transparent local models

Crisp DEMATEL: direct ratings, feedback and total influence

Follow a declared influence model, and verify that its feedback sum exists.

This tool studies a small, complete matrix of perceived direct influences. It uses one explicitly cited max-row/max-column normalization and refuses to manufacture a total-relation matrix when convergence or numerical conditioning fails. Every input, matrix operation and graphical threshold is visible.

Interactive workbench

Use 2–6 factors with a complete square matrix. Ratings are 0 or 10⁻¹²–4; the diagonal must be exactly zero. The fixed 0–4 perceived-influence scale allows fractional averaged ratings. No blank cell is interpreted as zero. τ lies in [0,10⁹] and filters only the drawing.

Complete input matrix
FactorF1F2F3F4

Example inputs are ready. Calculate to inspect your current inputs.

All calculations and CSV creation run locally in your browser. No server, file upload or external runtime service is used.

Method, assumptions and interpretation

DEMATEL calculation and verification sequenceDirect ratings A are divided by the declared s to obtain D. Total T is computed only after checking matrix-power contraction, inverse conditioning and residuals. Failed checks stop the calculation.Direct judgmentsAD = A / sDeclared normalizerVerify first‖Dᵏ‖∞ < 1κ∞(I−D)ResidualsT = D(I−D)⁻¹R ± CNo total influence without convergenceFailed check → stop; no hidden damping
Arrows describe a computational workflow; they are not causal evidence about a maritime system.

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

Symbols and units
SymbolDefinition and unit
n; A=[aᵢⱼ]Factor count; complete direct-influence rating matrix on the declared 0–4 scale
s; DMaximum row/column sum; normalized direct matrix (dimensionless)
I; T=[tᵢⱼ]Identity matrix; total-relation matrix (dimensionless model units)
ρ(D); kSpectral 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.

Complete input matrix
↓ / →PlanningSparesWork readinessExecution
Planning0310
Spares1021
Work readiness0102
Execution1000
Raw influence sums and normalization
FactorRaw row sumRaw column sum
Planning42
Spares44
Work readiness33
Execution13

s = 4

Normalized direct-influence matrix D
↓ / →PlanningSparesWork readinessExecution
Planning00.750.250
Spares0.2500.50.25
Work readiness00.2500.5
Execution0.25000
System matrix I − D
↓ / →PlanningSparesWork readinessExecution
Planning1-0.75-0.250
Spares-0.251-0.5-0.25
Work readiness0-0.251-0.5
Execution-0.25001
Inverse (I − D)⁻¹
↓ / →PlanningSparesWork readinessExecution
Planning1.6115107911.4964028781.1510791370.9496402878
Spares0.6906474821.7841726621.0647482010.9784172662
Work readiness0.37410071940.63309352521.4100719420.8633093525
Execution0.40287769780.37410071940.28776978421.237410072
Total-relation matrix T, including feedback diagonal
↓ / →PlanningSparesWork readinessExecution
Planning0.61151079141.4964028781.1510791370.9496402878
Spares0.6906474820.78417266191.0647482010.9784172662
Work readiness0.37410071940.63309352520.41007194240.8633093525
Execution0.40287769780.37410071940.28776978420.2374100719
Factor summaries from the full T matrix
FactorOutgoing RIncoming CProminence R + CRelation R − CδᵢInterpretation
F1: Planning4.2086330942.0791366916.2877697842.1294964031.86166095e-12Net transmitter in supplied model
F2: Spares3.5179856123.2877697846.8057553960.23021582732.015024323e-12Net transmitter in supplied model
F3: Work readiness2.280575542.9136690655.194244604-0.63309352521.537893828e-12Net receiver in supplied model
F4: Execution1.3021582733.0287769784.330935252-1.7266187051.282288206e-12Net receiver in supplied model
Convergence and numerical verification
CheckValueAcceptance rule
Power k2k ≤ 65536
‖Dᵏ‖∞ upper bound0.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-16Raw residual
Scaled total-equation residual3.536768535e-17≤ 10⁻¹²
Smallest absolute pivot0.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.

Threshold-filtered total-influence networkUp to six factors use a fixed circular layout. Only off-diagonal total influences strictly greater than the display threshold are drawn. Arrows go from row factor to column factor; the following table provides full values.F1 → F2: 1.496402878F1 → F3: 1.151079137F1 → F4: 0.9496402878F2 → F1: 0.690647482F2 → F3: 1.064748201F2 → F4: 0.9784172662F3 → F2: 0.6330935252F3 → F4: 0.8633093525F1F2F3F4τ = 0.5; displayed edges: 8Direction: row factor to column factor
Circular placement is display-only. The diagonal is not drawn but remains in every calculation. F1 = Planning; F2 = Spares; F3 = Work readiness; F4 = Execution.
Prominence–relation planeHorizontal position is R plus C; vertical position is R minus C. Outside the numerical guard, above zero denotes a net transmitter and below zero a net receiver in the model. Factor codes match the network figure; the summary table gives all point coordinates.0-2.6618701.9056115-1.3309353.811223005.71683451.33093527.62244602.6618705F1F2F3F4Prominence R + CRelation R − C
Points summarize perceived influence ratings; they do not prove empirical causation. Grey points have an unresolved relation sign within the numerical guard.
Displayed directed edges only
FromToTᵢⱼ
F1: PlanningF2: Spares1.496402878
F1: PlanningF3: Work readiness1.151079137
F1: PlanningF4: Execution0.9496402878
F2: SparesF1: Planning0.690647482
F2: SparesF3: Work readiness1.064748201
F2: SparesF4: Execution0.9784172662
F3: Work readinessF2: Spares0.6330935252
F3: Work readinessF4: Execution0.8633093525
Display-only threshold audit
τEdges drawnCalculation matrix
012Full T; unchanged
0.58Full T; unchanged
1.4964028780Full 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.

Related context

The method explanation and worked example are on this page. The articles below provide additional context.

All calculators