KraljicMatrix Package
Purchasing strategies are a significant contributing factor to the overall performance of a firm. Not only has it been argued that transitioning purchasing from a tactical to a strategic role provides sustainable competitive advantage, compelling evidence supports the relationship between strategic purchasing and supplier integration, performance, power, and total cost of ownership among other benefits. Consequently, it has become widely accepted that establishing a strategic purchasing approach is imperative for today’s firm.
Towards this goal, the most established portfolio model in literature and practice is Kraljic’s model^{1} known as the Kraljic Portfolio Matrix (KPM). However, one of the primary weaknesses of the KPM is the qualitative nature of the model, which results in a subjective method for weighting and positioning suppliers or commodities in the quadrants of the KPM. Thus, Jason Freels and I developed the kraljicMatrix
package to be used as a quantitative framework to objectively position purchases in the KPM matrix along with identifying which products and services to begin strategic sourcing efforts with. This post provides some explanation behind the functionality of kraljicMatrix
, which is largely based on research illustrated in this Kraljic Matrix paper produced by two other authors and myself.
The Kraljic Portfolio Matrix
Kraljic (1983) introduced the first comprehensive portfolio approach for purchasing and supply management. Although other models have been developed, Kraljic’s KPM approach has become the dominant methodology built upon in literature and the foundation for purchasing strategies for many organizations. Kraljic’s approach leverages a portfolio matrix that classifies products on the basis of two dimensions: the external dimension (‘Supply Risk’) concerns the factors regarding suppliers and supply market, while the internal dimension (‘Profit Impact’) relates to the importance and profit impact of a given product. Each dimension needs to be assessed against a number of variables where an overall classification score (‘low’ and ‘high’) is established. The result is a matrix and a classification in four categories: noncritical, leverage, bottleneck, and strategic items as illustrated in Figure 1. This categorization allows commodities to be classified in a way that minimizes the supply risk while maximizing purchasing power and profits of the firm. Furthermore, positioning commodities in this matrix form not only enables managers to visually see tradeoffs between items, but the classification also allows firms to develop distinctive management approaches for each category.^{2}
A Quantitative Approach for Positioning
A significant question in applying the KPM is how to categorize purchased products and services within the matrix. Organizations vary in how they measure supply risk and profit impact. Furthermore, some organizations may choose to measure their purchases against different dimensions; however, the fundamental purpose of the matrix remains. For our purpose, we assume that the organization has developed a means to condense their measurement of each dimension to a single value. For example, many organizations use an index (i.e. IBISWorld Buyer Power Score) to measure one of the dimensions. Or some organizations develop an indexed value function that generates a single value score for many attributes (i.e. profit impact can be a function of volume purchased, expected growth in demand, percent of total purchase cost, impact on product quality and business growth, etc.). However, once you have a single value that represents each dimension, subjectivity still largely drives how they are positioned in the KPM. The kraljicMatrix
package was designed to assist with this concern and the examples that follow walk you through how to implement kraljicMatrix
functions.
Primary functions in kraljicMatrix
There are eight primary functions provided by kraljicMatrix
as follows:
Function  Purpose 

SAVF_score 
Computes a utility score based on an exponential single attribute value function 
SAVF_preferred_rho 
Computes the preferred rho that minimizes the squared error between subject matter inputs and exponentially fitted utility scores 
SAVF_plot_rho_error 
Plots the squared error terms for the rho search space to illustrate the preferred rho that minimizes the squared error between subject matter desired values and exponentially fitted scores 
SAVF_plot 
Plots the single attribute utility curve along with the subject matter desired values for comparison 
MAVF_score 
Computes the multiattribute value score based on x and y attribute utility scores and their respective weights 
MAVF_sensitivity 
Computes summary statistics for multiattribute value scores for x and y given a range of swing weights for each attribute 
kraljic_quadrant 
Identifies the Kraljic purchasing matrix quadrant for each product or service based on the attribute utility scores of \code{x} and \code{y} 
kraljic_matrix 
Plots each product or service in the Kraljic purchasing matrix based on the attribute value score of x and y 
Furthermore, an example dataset psc
is provided in the kraljicMatrix
package. This data contains 200 product and service contracts (PSC). Each PSC has an x attribute (i.e. supply risk) score from 1 (worst) to 5 (best) and y attribute (i.e. profit impact) score from 1 (worst) to 10 (best).
psc
## # A tibble: 200 × 3
## PSC x_attribute y_attribute
## <chr> <int> <int>
## 1 D233 3.01 4.84
## 2 F352 4.34 5.64
## 3 T713 3.37 4.30
## 4 K833 2.67 5.53
## 5 Q121 3.48 4.33
## 6 C791 3.32 7.32
## 7 Y207 3.48 5.42
## 8 W439 2.47 3.35
## 9 N290 1.66 4.02
## 10 C251 1.00 7.47
## # ... with 190 more rows
Implementation of kraljicMatrix
The x and y attributes are simply evaluation measures. They enable each product and service to obtain a score for each dimension being measured. For example, the x attribute score (15 in .01 increments) could be the IBISWorld Buyer Power Score measuring supply market complexity. However, to plot these attributes on the KPM matrix we need to normalize the value scores such that the values are between 01. To do this we can use an exponential single attribute value function (SAVF). For example, let represent the normalized value of the x attribute such that and are the lowest and highest preferred value of attribute x respectively. Thus, and . Consequently, let be the SAVF of exponential form whereby each is an input and is the exponential constant for :
However, prior to applying the SAVF to our x and y attributes we must first identify the appropriate $\rho$ value. The benefit of applying an exponential SAVF is that it can take on many forms of increasing rates, along with aligning to a linear value function. Consequently, if certain x attribute values are valued more than other values an exponential SAVF will capture this utility curve. To identify the appropriate exponential rate, subject matter expert (SME) inputs are typically evaluated and an exponential rate that most closely matches the preffered values provided by the SMEs is chosen. Thus, let’s assume for our given x attribute the SME inputs suggest that x attribute values of 3, 4, & 5 provide a utility score of .75, .90 & 1.0 respectively (this represents a decreasing rate of return utility curve). Knowing that our x attribute is bounded between 1 and 5 we can search for a rho value between 01 that provides the best fit utility function using the SAVF_preferred_rho
function.
SAVF_preferred_rho(desired_x = c(3, 4, 5),
desired_v = c(.8, .9, 1),
x_low = 1,
x_high = 5,
rho_low = 0,
rho_high = 1)
## [1] 0.6531
Thus, we can see that provides the best fit exponential SAVF. We can illustrate this two ways. First, we can use SAVF_plot
to plot the single attribute utility curve compared to the subject matter desired values.
SAVF_plot(desired_x = c(3, 4, 5),
desired_v = c(.8, .9, 1),
x_low = 1,
x_high = 5,
rho = 0.6531)
We can also visualize the errors of the search space with SAVF_plot_rho_error
, which plots the squared error terms for all values within the search space to illustrate the preferred rho that minimizes the squared error between subject matter desired values and exponentially fitted scores.
SAVF_plot_rho_error(desired_x = c(3, 4, 5),
desired_v = c(.75, .9, 1),
x_low = 1,
x_high = 5,
rho_low = 0,
rho_high = 1)
Once we’ve identified the preferred value, we can now apply the exponential SAVF with SAVF_score
to normalize our attributes based on our utility curve.
# using dplyr to add a new variable while preserving existing data
library(dplyr)
# here we are assuming we found the appropriate rho value for the y attribute using
# the same process as mentioned above
psc < psc %>%
mutate(x_SAVF_score = SAVF_score(x_attribute, 1, 5, .653),
y_SAVF_score = SAVF_score(y_attribute, 1, 10, .70))
psc
## # A tibble: 200 × 5
## PSC x_attribute y_attribute x_SAVF_score y_SAVF_score
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 D233 3.01 4.84 0.7887459 0.9336977
## 2 F352 4.34 5.64 0.9573299 0.9629164
## 3 T713 3.37 4.30 0.8495938 0.9023958
## 4 K833 2.67 5.53 0.7165401 0.9598009
## 5 Q121 3.48 4.33 0.8655080 0.9044624
## 6 C791 3.32 7.32 0.8419735 0.9898314
## 7 Y207 3.48 5.42 0.8655080 0.9564360
## 8 W439 2.47 3.35 0.6659448 0.8084720
## 9 N290 1.66 4.02 0.3778582 0.8808636
## 10 C251 1.00 7.47 0.0000000 0.9910284
## # ... with 190 more rows
Now that we have the normalized x and y attribute utility scores we can proceed with plotting each PSC within the Kraljic matrix with kraljic_matrix
.
kraljic_matrix(psc, x_SAVF_score, y_SAVF_score)
This illustrates that most of our PSCs fall in the “Leverage” (upper left) quadrant while a few fall in the “Strategic” (upper right) and “Noncritical” (lower left) quadrants and no PSCs fall in the “Bottleneck” quadrant. Keep in mind that each category benefits from a different strategic sourcing approach. So decisionmakers benefit from understanding specifically which products and services align to each so that they can coordinate the appropriate sourcing strategy for that particular product or service. We can easily do this with the kraljic_quadrant
function.
psc %>%
mutate(quadrant = kraljic_quadrant(x_SAVF_score, y_SAVF_score))
## # A tibble: 200 × 6
## PSC x_attribute y_attribute x_SAVF_score y_SAVF_score quadrant
## <chr> <dbl> <dbl> <dbl> <dbl> <chr>
## 1 D233 3.01 4.84 0.7887459 0.9336977 Leverage
## 2 F352 4.34 5.64 0.9573299 0.9629164 Leverage
## 3 T713 3.37 4.30 0.8495938 0.9023958 Leverage
## 4 K833 2.67 5.53 0.7165401 0.9598009 Leverage
## 5 Q121 3.48 4.33 0.8655080 0.9044624 Leverage
## 6 C791 3.32 7.32 0.8419735 0.9898314 Leverage
## 7 Y207 3.48 5.42 0.8655080 0.9564360 Leverage
## 8 W439 2.47 3.35 0.6659448 0.8084720 Leverage
## 9 N290 1.66 4.02 0.3778582 0.8808636 Strategic
## 10 C251 1.00 7.47 0.0000000 0.9910284 Strategic
## # ... with 190 more rows
Lastly, it is important to keep in mind that decisionmakers may weight the importance of each attribute differently. Consequently, due to certain market environments, decisionmakers may weight the x attribute (i.e. supply risk) of greater importance than the y attribute (i.e. profit impact). Thus, we can prioritize PSCs based on this preference by applying a multiattribute value function (MAVF) with swing weights. Swing weight values for x and y attributes ( and respectively) are typically elicited from SMEs. This allows for calculation of the interaction swing weight . Thus, we can calculate the MAVF as outlined by Keeney and Raiffa (1993)^{3}:
Thus, we can apply the MAVF_score
function to compute the multiattribute value score based on x
and y
attribute utility scores and their respective swing weights. So if through discussions with decisionmakers we identify swing weight values of 0.65 and 0.35 for the x and y attributes respectively, we can obtain the computed MAVF score for each PSC:
psc %>%
mutate(MAVF = MAVF_score(x_SAVF_score, y_SAVF_score, 0.65, 0.35))
## # A tibble: 200 × 6
## PSC x_attribute y_attribute x_SAVF_score y_SAVF_score MAVF
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 D233 3.01 4.84 0.7887459 0.9336977 0.8394790
## 2 F352 4.34 5.64 0.9573299 0.9629164 0.9592852
## 3 T713 3.37 4.30 0.8495938 0.9023958 0.8680745
## 4 K833 2.67 5.53 0.7165401 0.9598009 0.8016814
## 5 Q121 3.48 4.33 0.8655080 0.9044624 0.8791420
## 6 C791 3.32 7.32 0.8419735 0.9898314 0.8937237
## 7 Y207 3.48 5.42 0.8655080 0.9564360 0.8973328
## 8 W439 2.47 3.35 0.6659448 0.8084720 0.7158293
## 9 N290 1.66 4.02 0.3778582 0.8808636 0.5539101
## 10 C251 1.00 7.47 0.0000000 0.9910284 0.3468599
## # ... with 190 more rows
This allows us to quickly dissect our PSCs. For example, if decisionmakers are most concerned with the “Leverage” quadrant but want to assess the top 10 PSCs based on the decisionmakers preferences of the attributes we can efficiently make this assessment. This identifies the top 10 PSCs that are most likely to benefit from a strategic sourcing approach specifically designed for “Leverage” PSCs.
psc %>%
mutate(MAVF = MAVF_score(x_SAVF_score, y_SAVF_score, 0.65, 0.35),
quadrant = kraljic_quadrant(x_SAVF_score, y_SAVF_score)) %>%
filter(quadrant == "Leverage") %>%
top_n(10, wt = MAVF)
## # A tibble: 10 × 7
## PSC x_attribute y_attribute x_SAVF_score y_SAVF_score MAVF
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Y357 5.00 6.55 1.0000000 0.9812541 0.9934389
## 2 E103 4.46 7.71 0.9665150 0.9927003 0.9756799
## 3 C432 4.65 6.05 0.9796633 0.9726273 0.9772007
## 4 P402 5.00 5.82 1.0000000 0.9675243 0.9886335
## 5 Q255 4.95 6.41 0.9973714 0.9791345 0.9909885
## 6 H426 5.00 5.58 1.0000000 0.9612468 0.9864364
## 7 E908 4.75 7.11 0.9859554 0.9879299 0.9866465
## 8 X634 5.00 5.19 1.0000000 0.9485047 0.9819766
## 9 O288 5.00 5.00 1.0000000 0.9409177 0.9793212
## 10 V870 5.00 5.88 1.0000000 0.9689357 0.9891275
## # ... with 1 more variables: quadrant <chr>
And finally, since our swing weight inputs are subjective in nature we may wish to perform a senstivity analysis on these swing weights to see their impact on MAVF scores. The MAVF_sensitivity
function executes a sensitivity analysis by performing a Monte Carlo simulation with 1000 trials for each product or service (row). Each trial randomly selects a weight from a uniform distribution between lower and upper bound swing weight parameters and calculates the multattribute utility score. From these trials, summary statistics for each product or service (row) are calculated and reported for the final output.
MAVF_sensitivity(psc,
x = x_SAVF_score,
y = y_SAVF_score,
x_wt_min = .55,
x_wt_max = .75,
y_wt_min = .25,
y_wt_max = .45) %>%
select(PSC, starts_with("MAVF"))
## # A tibble: 200 × 8
## PSC MAVF_Min MAVF_1st_Q MAVF_Median MAVF_Mean MAVF_3rd_Q
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 D233 0.8152 0.8293 0.8393 0.8395 0.8490
## 2 F352 0.9519 0.9570 0.9594 0.9593 0.9614
## 3 T713 0.8469 0.8610 0.8683 0.8680 0.8746
## 4 K833 0.7722 0.7881 0.8015 0.8018 0.8151
## 5 Q121 0.8593 0.8728 0.8793 0.8791 0.8851
## 6 C791 0.8775 0.8861 0.8938 0.8938 0.9016
## 7 Y207 0.8812 0.8906 0.8972 0.8974 0.9035
## 8 W439 0.6772 0.7021 0.7163 0.7158 0.7290
## 9 N290 0.4958 0.5265 0.5539 0.5542 0.5812
## 10 C251 0.2478 0.2980 0.3461 0.3475 0.3974
## # ... with 190 more rows, and 2 more variables: MAVF_Max <dbl>,
## # MAVF_Range <dbl>
Additional Resources
This package originated from the research illustrated in this Kraljic Matrix paper. Please refer to it for more details regarding Kraljic purchasing matrix literature and the single and multiattribute value functions discussed in this vignette.

Kraljic, P. (1983). Purchasing must become supply management. Harvard Business Review, 61(5), 109117. ↩

See Gelderman & Semeijn (2006) and Gelderman & Van Weele (2003) for more discussion regarding strategic management approaches. ↩

Keeney, R., & Raiffa, H. (1993). Decisions with multiple objectives: preferences and value tradeoffs. Cambridge University Press. ↩