# The Intertwining Impact of Intraorganizational and Routine Networks on Routine Replication Dynamics: An Agent-Based Model

Authors: Dehua Gao; Flaminio Squazzoni; Xiuquan Deng
Journal: Complexity
 (2018)

---

## Abstract

Organizational routines are means through which organizations can reutilize best practices and so their replication, i.e., duplicating beneficial routines across context, is a key value-creating strategy. However, it is difficult to map network effects on routine replication. Here, we investigated routine replicating dynamics considering two types of network contexts, namely, (1) connections between different (geographically distributed) units in a decentralized organization and (2) the coupling relation between routines, i.e., a bundle of different routines involved in each unit. By considering routine replication as one kind of template-based activities between different units, we examined interrelations between routines with a NK-based fitness landscape model. Our results show that when there is an appropriate level of absorptive capacities (i.e., when organizations are capable of identifying and acquiring externally generated knowledge), there is an optimal combination of these two types of networks, which is beneficial to routine replicating practices and organization adaptation. Furthermore, we also found that intraorganizational variations, including template-duplicating errors and innovative activities, are instrumental to enhance adaptive changes. Our findings suggest measures to control and manage best practice diffusion across organizations.


---

## Body

## 1. Introduction

Recent research has shown that organizational routines, i.e., “repetitive, recognizable pattern of interdependent actions” carried out by multiple individual actors [1], are “ubiquitous” in almost every organization [2]. While certain scholars have focused on the so-called “defensive routines”, i.e., actions or policies that “prevent either individuals or segments of the organization from experiencing embarrassment or threat” [3, 4], and explored their negative effect on organizational learning and change [5–7], there is evidence that organizational routines are a key repository of organizational memory, skills and tacit knowledge which enables both organizational stability and change [1, 8–11] and constitutes organizations’ sustained competitive advantages [12–14]. In this sense, by applying routines representing “best practices” to new environments, organizational core capabilities can be enhanced that in turn further sustain competitive advantages and profitability [15]. This is why the replication of organizational routines is considered as an important value-creating strategy and attracted more increasing attention from both researchers and practitioners [16–18].Firesl and Larty [18] summarized three value-creating mechanisms of replicating routines for which replication is viewed (1) as a strategy for growth; (2) a means to create value by improving organizational efficiency; and (3) an organizational scaffold that may enhance flexibility. These value-creating mechanisms, according to Jenson [19] and D’Adderio [20], are particularly applicable to organizations, such as large or multinational firms, within which multiple similar units perform the same function, or franchise and chain stores, which imply the geographic replication of similar units.According to Szulanski and Jensen [21], replicating a routine involves the creation of another routine that is ‘similar to the original routine in significant respects.’ However, Howard-Grenville [22] pointed out that organizational routines might perform differently depending on specific organizational environments. Given certain context-dependent characteristics, the replication of organizational routines can be understood as a simple knowledge transfer process. Indeed, any match between these newly created routines and the context they were originally coping with imposes what Winter and Szulanski [15] called a “replication dilemma” problem. This means that besides the exact “cloning” of the original routines, these new routines must be to some extent flexible and adaptive to be suitable to the new, in part similar but also different, organizational context [18, 20].Therefore, replication of routines is a means through which organizations can reutilize knowledge that is already in use [21, 23]. Some scholars focused on the importance of replication as a “knowledge transfer” across intraorganizational boundaries [23–26] and considered the intraorganizational transferring process between different units, often across different geographic locations, from a network-based perspective [27–29]. Here, there is the important role of templates during the replicating process of routines [21, 30] as a means to stimulate forward and reverse knowledge flows [18].As suggested by D’Adderio [20] and Friesl and Larty [31], replicating routines is a complex social practice, which involves multiple actors and artifacts. For instance, the capacity of individual actors to understand specific routines is always intertwined with their own behaviors, experiences, and the particular environment within which they live [12]. Therefore, it is likely that individuals have to make context-specific adjustments so as to reproduce “the similar routine” in prospection [11, 32–34]. Unfortunately, as suggested by Friesl and Larty [18, 31], we still know little about the micro processes through which individual actors and actions shape and enact new routines in their new contexts.Furthermore, research on organizational routines has emphasized the connective nature of routines in that organizations and capabilities are better understood as dependent on some “specific networks of routines,” rather than just on bundles of them [11, 35–37]. This is because routines often do not work alone but are interrelated with each other [38]. Secondly, recent findings indicated that the dynamics of several interrelated routines together, which are often what occurs in practice and to some extent determines organizational outcomes, are definitely different from that of individual routines [39]. If so, routine interdependency can be a fundamental source of both organizational stability and flexibility [40]. However, as recently suggested by Feldman et al. [11] and Kremser and Schreyögg [39], more research is needed to understand ‘implications and consequences of routine interrelatedness.’Our paper aims to investigate the replicating dynamics of organizational routines within two types of network contexts. First, we hypothesized connections between (geographically) distributed units of an organization [18, 41, 42]. In a second scenario, we hypothesized a routine-related network – i.e., a network of different routines involved in the whole organization [11]. To do so, we first followed the NK and fitness landscape modeling [43] to depict routine networks and built up the fitness landscape of an organization. Then, we translated this into an agent-based model (from now on, ABM) [44–46] to explore complex behaviors of routine replications.In our model, we focused neither on individual actors nor on artifacts in the enactment of routines in new contexts. Rather, we treated an aggregation of these individuals and artifacts, namely, the unit distributed within an organization, which we called a “unit”, as the basic heterogeneous but interdependent agent. All the units (i.e., agents) randomly interacted with each other, while by using templates specific routines were replicated from unit to unit. Furthermore, we also considered the impact of both internal variability and external dynamics of the organizational environment on routine replications. This included, for example, template-duplicating errors and innovative activities among units and the pace of environmental changes.The rest of the paper is structured as follows. In the second section, we provide a brief literature overview on replication of organizational routines. Further, we also summarized the feasibility and effectiveness of adopting an NK-inspired ABM to examine organizational routines. Then, we present a series of concrete scenarios with combinations of input settings and show our simulation results. Finally, the last section summarizes our main findings, discusses limitations, and elaborates on implications.


## 2. Theoretical Foundations

### 2.1. Replication of Routines within Organizations

Nelson and Winter [8] and Zollo and Winter [47] considered an organization as “a set of interdependent operational and administrative routines.” These different kinds of routines embody an important part of the organization’s productive knowledge [8] and are considered as one of the key sources for sustained competitive advantages [12, 48, 49]. This suggests that sustaining competitive advantages can also depend on the capability of replicating the most efficacious routines from a unit to another, or even from an organization to another, just as efficacious habits are replicated from an individual to another one [17, 50]. However, in some respects, routine replication can be more difficult than the replication of individual habits. First, after discovering an efficacious practice, managers tend to build organizational barriers against replication, so as to exploit the system’s advantage and avoid imitation by others [51]. Second, knowledge embedded in routines is always to some extent not explicit but tacit and enacted in practice [32, 47]. At the same time, either units or organizations have only limited capacities to discern superior practices from their neighbors [50].Furthermore, given the “context-dependent” characteristics of routines [22] and the need for adaptive changes [15, 52], for any organization, routine replication has value when it can support appropriate variations that are instrumental to cope with environmental challenges [15, 18, 20].In this respect, by discussing the connection between intraorganizational replication and interorganizational imitation, Rivkin [51] further pinpointed that replication can be more efficient than imitation in promoting sustained competitive advantage only when the complexity of the task and the environment is moderate. Winter and Szulanski [23] distinguished the term “replication” from similar concepts such as “transfer of knowledge” and “diffusion of innovations” within and/or even between organizations. They suggested that replication involves the creation of “a series of local routines that are quite similar to the original ones in significant respects” [23]. This implies that replication often relates to an organization’s attempts to “reproduce at multiple internal sites the outcome of an existing activity” and so “applies to a broad class of exploitation situations” [8, 23].Since Nelson and Winter [8], many scholars explored the dynamics of replication, such as routine replication barriers [24, 53], the impact of replicator/replicatee relationship, and the role of institutional contexts [54]. Others looked at the movement of employees or independent experts as a carrier of knowledge and experience transfer from one context to another (e.g., Hodgson [50]). If replication is a “knowledge transferring” process across intraorganizational boundaries [24, 25], this process depends on the fact that the ‘core aspects of a routine’, or even “a whole business format”, can be truly communicated from one unit to another, which then transforms related knowledge into practice [18].Research suggested that “templates” are key for replicating organizational routines [21, 52]. Here, the term “template,” according to Nelson and Winter [8], referred to “working examples” that contain some critical aspects of the routine. Given that knowledge embedded in routines is always to some extent not explicit, but tacit and enacted in practice [32, 47], templates define the “details of how the work gets done, in what sequence, and how its various components and subroutines are interconnected” [19]. They can help us codifying and explicitly representing the tacit part of knowledge, thereby providing an effective way to overcome stickiness to the original context of knowledge [21]. Meanwhile, by codifying and explicitly representing knowledge inherent in routines, templates can also be helpful to the replicatee units in discerning and adopting the so-called “best practices” [55]. This requires to look at the relationship between the sending and receiving units [24, 56], as the ability and willingness of replicatees to accept knowledge transferred from the sending unit are also prerequisites for the replicating process of routines [57].


### 2.2. Routine Networks and the NK Fitness Landscape

Organizations can be considered as certain specific aggregations of “bundle of routines” [8, 48, 58, 59], which can form and change through a process of “combinatorics” [12]. However, it is only until recently that scholars began to focus on how routines affect one another in a more complexity-friendly perspective [11, 36, 60]. For instance, Peng et al. [36] emphasized that looking at routine interrelations is key to understand how routines bundled together to generate organizational capabilities – including both incremental and radical innovation. According to Kremser and Schreyögg [39], this interrelatedness can be attributed to the division of labor in that, due to economics of specialization, organizations while dividing work into “multiple, separate routines” are also required to “securing a sufficient integration of efforts” [39]. This suggests the need for looking not only at the division of multiple specialized routines from the complex organizational task, but also at the interrelationship between these routines. This implies that only by specific routine integration organizations can fully accomplish their organizational task.Here, research suggests a holistic perspective as organizational routines are embedded in systems of routines [14, 36, 61]. This means that they are systemically integrated and reinforce each other so giving rise to a collective effect, which is greater than only the single contribution of a simple set of individual components [62]. In this line, Argote and Ingram [25] considered the “network of routines” as an important form of organizational memories and the repository of organizational knowledge. Salvato [63] argued that ‘the interplay of existing routines’, together with individuals’ everyday actions, should be treated as real “unit of analysis” rather than single factors. More recently, Sele and Grand [60] investigated the role of both human and nonhuman actants in the connection between routines and considered the interaction of routines as a source of generativity for the whole organizational system.In this paper, we tried to incorporate some aspects of these networks of routines by using Kauffman’s [43] NK model to formalize interactions among these multiple organizational routines. Being originally developed to analyze biological organizations’ evolutionary problems [64], Kauffman’s NK modeling and fitness landscape theory was then introduced into organizational studies by Levinthal [65]. It eventually became popular in studies on technological innovation [66], strategic planning [67, 68], and organizational learning [69]. NK modeling now can be considered as one of the crucial toolkits for coping with a series of classic organizational problems of search in multidimensional spaces. According to Ganco and Hoetker [68], NK modeling provides a “relatively simple and replicable methodology” addressing organizational problems that cannot be comprehensively understood empirically.


## 2.1. Replication of Routines within Organizations

Nelson and Winter [8] and Zollo and Winter [47] considered an organization as “a set of interdependent operational and administrative routines.” These different kinds of routines embody an important part of the organization’s productive knowledge [8] and are considered as one of the key sources for sustained competitive advantages [12, 48, 49]. This suggests that sustaining competitive advantages can also depend on the capability of replicating the most efficacious routines from a unit to another, or even from an organization to another, just as efficacious habits are replicated from an individual to another one [17, 50]. However, in some respects, routine replication can be more difficult than the replication of individual habits. First, after discovering an efficacious practice, managers tend to build organizational barriers against replication, so as to exploit the system’s advantage and avoid imitation by others [51]. Second, knowledge embedded in routines is always to some extent not explicit but tacit and enacted in practice [32, 47]. At the same time, either units or organizations have only limited capacities to discern superior practices from their neighbors [50].Furthermore, given the “context-dependent” characteristics of routines [22] and the need for adaptive changes [15, 52], for any organization, routine replication has value when it can support appropriate variations that are instrumental to cope with environmental challenges [15, 18, 20].In this respect, by discussing the connection between intraorganizational replication and interorganizational imitation, Rivkin [51] further pinpointed that replication can be more efficient than imitation in promoting sustained competitive advantage only when the complexity of the task and the environment is moderate. Winter and Szulanski [23] distinguished the term “replication” from similar concepts such as “transfer of knowledge” and “diffusion of innovations” within and/or even between organizations. They suggested that replication involves the creation of “a series of local routines that are quite similar to the original ones in significant respects” [23]. This implies that replication often relates to an organization’s attempts to “reproduce at multiple internal sites the outcome of an existing activity” and so “applies to a broad class of exploitation situations” [8, 23].Since Nelson and Winter [8], many scholars explored the dynamics of replication, such as routine replication barriers [24, 53], the impact of replicator/replicatee relationship, and the role of institutional contexts [54]. Others looked at the movement of employees or independent experts as a carrier of knowledge and experience transfer from one context to another (e.g., Hodgson [50]). If replication is a “knowledge transferring” process across intraorganizational boundaries [24, 25], this process depends on the fact that the ‘core aspects of a routine’, or even “a whole business format”, can be truly communicated from one unit to another, which then transforms related knowledge into practice [18].Research suggested that “templates” are key for replicating organizational routines [21, 52]. Here, the term “template,” according to Nelson and Winter [8], referred to “working examples” that contain some critical aspects of the routine. Given that knowledge embedded in routines is always to some extent not explicit, but tacit and enacted in practice [32, 47], templates define the “details of how the work gets done, in what sequence, and how its various components and subroutines are interconnected” [19]. They can help us codifying and explicitly representing the tacit part of knowledge, thereby providing an effective way to overcome stickiness to the original context of knowledge [21]. Meanwhile, by codifying and explicitly representing knowledge inherent in routines, templates can also be helpful to the replicatee units in discerning and adopting the so-called “best practices” [55]. This requires to look at the relationship between the sending and receiving units [24, 56], as the ability and willingness of replicatees to accept knowledge transferred from the sending unit are also prerequisites for the replicating process of routines [57].


## 2.2. Routine Networks and the NK Fitness Landscape

Organizations can be considered as certain specific aggregations of “bundle of routines” [8, 48, 58, 59], which can form and change through a process of “combinatorics” [12]. However, it is only until recently that scholars began to focus on how routines affect one another in a more complexity-friendly perspective [11, 36, 60]. For instance, Peng et al. [36] emphasized that looking at routine interrelations is key to understand how routines bundled together to generate organizational capabilities – including both incremental and radical innovation. According to Kremser and Schreyögg [39], this interrelatedness can be attributed to the division of labor in that, due to economics of specialization, organizations while dividing work into “multiple, separate routines” are also required to “securing a sufficient integration of efforts” [39]. This suggests the need for looking not only at the division of multiple specialized routines from the complex organizational task, but also at the interrelationship between these routines. This implies that only by specific routine integration organizations can fully accomplish their organizational task.Here, research suggests a holistic perspective as organizational routines are embedded in systems of routines [14, 36, 61]. This means that they are systemically integrated and reinforce each other so giving rise to a collective effect, which is greater than only the single contribution of a simple set of individual components [62]. In this line, Argote and Ingram [25] considered the “network of routines” as an important form of organizational memories and the repository of organizational knowledge. Salvato [63] argued that ‘the interplay of existing routines’, together with individuals’ everyday actions, should be treated as real “unit of analysis” rather than single factors. More recently, Sele and Grand [60] investigated the role of both human and nonhuman actants in the connection between routines and considered the interaction of routines as a source of generativity for the whole organizational system.In this paper, we tried to incorporate some aspects of these networks of routines by using Kauffman’s [43] NK model to formalize interactions among these multiple organizational routines. Being originally developed to analyze biological organizations’ evolutionary problems [64], Kauffman’s NK modeling and fitness landscape theory was then introduced into organizational studies by Levinthal [65]. It eventually became popular in studies on technological innovation [66], strategic planning [67, 68], and organizational learning [69]. NK modeling now can be considered as one of the crucial toolkits for coping with a series of classic organizational problems of search in multidimensional spaces. According to Ganco and Hoetker [68], NK modeling provides a “relatively simple and replicable methodology” addressing organizational problems that cannot be comprehensively understood empirically.


## 3. The Model

### 3.1. A fitness Landscape Approach to Routine Networks

Based on Nelson and Winter’s [8] notion of routines as the “genes” of organizations, Ventresca and Kaghan [70] and Yi et al. [71] suggested to apply the NK model to the study of organizational routines. Through the NK-based fitness landscape, we can build up a description of the linkage between the network of routines and the organizational fit with its environment. For example, let us assume that there areN routines involved in a unit, R=r1,r2,…,rN, and for each routine ri∈R  (i = 1, 2, …,N), each owning a number of different states – or in other words, different kinds of approaches for unit members to accomplish the corresponding organizational subtasks Si=s1,s2,…,sA. Thus, according to Kauffman’s [43] theory, this routine’s contribution to the overall unit fitness fiti would depend not only on its own state, but also on the states of some otherK routines RKi=r1i,r2i,…,rKi (where RKi⊂R and ri∉RKi) with which it is interconnected.Kauffman [43] described two main methods for assigning relations between dependent parts – namely nearest neighbor influences and randomly assigned influences. He then proved that the distribution ofK did not dramatically affect the NK model outcomes. In this paper, for simplicity, we adopted the latter strategy and randomly assigned relations between routines [72], as shown in Figure 1.Figure 1
Randomly assigned relations between routines (N = 10, K = 4).Thus, we had the following fitness contribution of each routineri,(1)fiti=fri;r1,r2,…,rK∣r1,r2,…,rK∈RKi,i=1,2,…,NHere, the value of the key parameterK can be from 0 toN – 1. On the one hand,K = 0 represented a situation in which each of theN routines involved in the unit was independent. On the other hand, with the increase of the value ofK (K > 0), the change of one routine would influence the performance of other interrelated routines in that any search for improving the configuration of all of them would be more and more difficult [71]. Finally,K =N – 1 represented another very specific situation in which all routines were interconnected with each other.The fitness value of the unit with its task environment can be calculated as(2)Fitness=1N∑i=1Nfiti=1N∑i=1Nfri;r1,r2,…,rK∣r1,r2,…,rK∈RKiaccording to which the higher the fitness valueFitness, the more adaptive the unit to its living environment.


### 3.2. Agent Design

In our model, we focused on network characteristics of distributed branches/units within an organization. For the sake of simplicity, we assumed that there was a decentralized organization [73] withunit_num = 40 units. All these units G=g1,g2,…,gM were interdependent entities with similar characteristics, facing a recurring but common problem. They were randomly located in a 100 × 100 lattice and interconnecting with each other with a given probabilityrandom_p (0 ≤random_p ≤ 1). The units constituting the whole organization faced the same environment defined by a landscape described in formulas (1) and (2). We assumed that every unit gj∈G (j = 1, 2, …,unit_num) was considered a heterogeneous agent dealing with a series of organizational subtasks OT=ot1,ot2,…,otN and having a correspondingly specific routine ri∈R=r1,r2,…,rN for each task oti∈OT  (i = 1, 2, …,N) [71]. All routines were interrelated with each other, resulting in interrelationships both within and between bundles of routines, as shown in Figure 2.Figure 2
Network of units with different routine networks.Furthermore, we assumed that each routine could generate only two different patterns of action, “0” and “1”, respectively. Although a specific routine may have more than two states, a set of routines could also be interpreted “as subroutines that constitute a routine of more than two states” [71].We assumed that each unit agent had attributes as shown in Table1.Table 1
Definition of agent attributes.
Attributes Definition unitID A unique ID number for agent identification.  unitNeighbors Every unit agent owned a list of neighbors defined by the intraorganizational network structure. Here, we assumed a classic random network following Erdös and Rényi [74], where nodes represented unit agents and linkages between nodes – which represented their neighborhoods – were randomly assigned by the probability value random_p (0 ≤ random_p ≤ 1). Given the bounded rationality characterizing organizational actors [75], we assumed that each unit agent could only interact with and replicate routines from nearest neighbors.  networkOfRoutines Each unit agent had a bundle of routines (i.e., a certain specific vector with “0” or “1” values) to cope with its series of organizational subtasks. However, rather than being only a simple collection, these routines were interrelated with each other [11, 39], so giving rise to a network of routines. We assumed that this network was randomly generated [43, 72], as illustrated in Figure 1.  Fitness The value of this parameter represented the adaptiveness of a unit agent to its living environment. In detail, we considered each routine’s fitness contribution functionf⋅ shown in formula (1) as a random function mapping the vector ri,r1,r2,…,rK (where, r1,r2,…,rK∈RKi and RKi⊂R, ri∉RKi, i=1,2,…,N) with uniformly distributed values fiti (0 ≤ fiti ≤ 1). Then, we assumed that the fitness value of the unit agent was numerically equal to the average of the fit value of all the routines involved (as shown in formula (2)).However, as we mentioned in Section2.1, unit agents often used templates when replicating routines from neighbors. We assumed that a template could not be the exact copy of the network of routines of the replicator unit agent. Indeed, given that at least a part of knowledge embedded in routines is to some extent tacit and enacted in practice [32, 47], these differences mimicked in difficulties for replicatee units to explicitly encode routines into templates.More specifically, we assumed that when replicating activities occurred, the replicatee unit agents obtained a template which was duplicated from the network of routines of the replicator unit agents but with some given probability of errorserror_p (0 ≤error_p ≤ 1). This allowed us to mimic errors that may occur during codifying and explicitly representing process of the original set of organizational routines [18, 25]. However, note that in the model we considered the value oferror_p as constant throughout all simulations. Considering reasons and/or even mechanisms for this kind of errors was out of the scope of our work.Finally, we considered theenvironment as a regime of possibilities and constraints for unit agents while accomplishing their organizational task [22]. This environment determined both the organizational subtask series and the fitness landscape as shown in formulas (1) and (2). In order to simulate environmental dynamics, in our model, we adopted a Boolean variableis_env_change? to represent whether the fitness landscape changed or not. On the one hand, whenis_env_change? =TRUE, the fitness landscape would randomly vary once after every certain number of simulation roundsenv_change_P so as to mimic environmental changes occurred in practice. Further, for anyenv_change_P≥ 1, we have that the lower the value of this variable was, the more rapid environmental changes happened. On the other hand, whenis_env_change? =FALSE, it meant that there was no change of the organizational environment. Thus, any values of the variableenv_change_Pare meaningless in this situation.


### 3.3. Rules of the Replicating Activities

According to Hodgson [50], replicating routines from peers is not easy often also due to organizational inertia (e.g., [71, 76]) and bounded rationality [75]. To mimic this, we assumed that unit agents were prone to replicate “superior practices” from neighbors only when they felt competitive pressures, i.e., they had a significantly lower fitness value than the average of neighbors.Next, atemplate was generated through explicitly codifying the knowledge embedded in the network of routines of the replicator unit agent [19, 21]. In this model, we used a vector of “0” or “1” values as follows: (3)TPL=tpl1,tpl2,…,tplNbeing templates of a copy of the original series of routines, but allowing each of the element values of the vectortplk (0 ≤k ≤N) to vary with a certain probability of errorserror_p (0 ≤error_p ≤ 1). Namely,(4)tplk=1-rk,withaprobabilityerror_prk,withaprobability1-error_pHere, the value of the coefficient rk was either “0” or “1”, which represented the two different states (i.e., patterns of actions) of thekth original routine of the replicator unit agent, and rk∈R=r1,r2,…,rN.Furthermore, given that replicatee unit agents were always constrained by a specific level of absorptive capacity [77, 78], we assumed that they could not adopt the whole template but only a certain part. Here, we used a variableunit_AC to represent the absorptive capacity of unit agents. Its value numerically equaled the number of routines that replicatee unit agents could change at a time according to the template,unit_AC = {1,2,3}. Obviously, the higher the value ofunit_AC, the higher level of the unit agent’s absorptive capacity. In this respect, note that exploring details of agents’ absorptive capacities was beyond the scope of our study. For the sake of simplicity, we assumed that the value ofunit_AC was assigned at the beginning of the simulation and remained fixed throughout our simulations.The replicatee unit agents could randomly select a number of routines (which numerically equaled theunit_AC value), which owned different states with those from the template – if they existed, of course. Then, the replicatee unit agents changed the relating routines’ states (values) according to the template and generated a new list of routinesR’.Finally, when the new collection of routinesR’ was formed, replicatee unit agents checked if they could perform any innovation to match their work of routines (obviously, here, referring to the list of routinesR’) and the specific context in which they were living. This was to mimic the “context-dependent” characteristic of organizational routines [22] as well as the interconnection between routines [15, 18, 20, 39]. For example, replicatee units could revise parts of routines while considering both the embeddedness of the system of routines [11] and certain fitness landscapes – namely, the units’ living environment [22, 42]. By doing these innovative activities, replicatee unit agents transformed the temporary list of routinesR’ into their final network of routinesR with newly updated status.In our model, we assumed that replicatee unit agents could innovate and randomly vary the state of one of their routines with a given probabilityinnovate_p (0 ≤innovate_p ≤ 1) – when this was helpful to improve their fitness values.Figure3 provides a description of the routine replicating mechanisms between different unit agents. Figure 4 shows a stylized description of the procedure in which unit agents involved in the model should execute during the whole simulation process.Figure 3
Mechanisms of replicating routines between unit agents (unit_AC = 2).Figure 4
Flow chart of replicating activities between unit agents.


## 3.1. A fitness Landscape Approach to Routine Networks

Based on Nelson and Winter’s [8] notion of routines as the “genes” of organizations, Ventresca and Kaghan [70] and Yi et al. [71] suggested to apply the NK model to the study of organizational routines. Through the NK-based fitness landscape, we can build up a description of the linkage between the network of routines and the organizational fit with its environment. For example, let us assume that there areN routines involved in a unit, R=r1,r2,…,rN, and for each routine ri∈R  (i = 1, 2, …,N), each owning a number of different states – or in other words, different kinds of approaches for unit members to accomplish the corresponding organizational subtasks Si=s1,s2,…,sA. Thus, according to Kauffman’s [43] theory, this routine’s contribution to the overall unit fitness fiti would depend not only on its own state, but also on the states of some otherK routines RKi=r1i,r2i,…,rKi (where RKi⊂R and ri∉RKi) with which it is interconnected.Kauffman [43] described two main methods for assigning relations between dependent parts – namely nearest neighbor influences and randomly assigned influences. He then proved that the distribution ofK did not dramatically affect the NK model outcomes. In this paper, for simplicity, we adopted the latter strategy and randomly assigned relations between routines [72], as shown in Figure 1.Figure 1
Randomly assigned relations between routines (N = 10, K = 4).Thus, we had the following fitness contribution of each routineri,(1)fiti=fri;r1,r2,…,rK∣r1,r2,…,rK∈RKi,i=1,2,…,NHere, the value of the key parameterK can be from 0 toN – 1. On the one hand,K = 0 represented a situation in which each of theN routines involved in the unit was independent. On the other hand, with the increase of the value ofK (K > 0), the change of one routine would influence the performance of other interrelated routines in that any search for improving the configuration of all of them would be more and more difficult [71]. Finally,K =N – 1 represented another very specific situation in which all routines were interconnected with each other.The fitness value of the unit with its task environment can be calculated as(2)Fitness=1N∑i=1Nfiti=1N∑i=1Nfri;r1,r2,…,rK∣r1,r2,…,rK∈RKiaccording to which the higher the fitness valueFitness, the more adaptive the unit to its living environment.


## 3.2. Agent Design

In our model, we focused on network characteristics of distributed branches/units within an organization. For the sake of simplicity, we assumed that there was a decentralized organization [73] withunit_num = 40 units. All these units G=g1,g2,…,gM were interdependent entities with similar characteristics, facing a recurring but common problem. They were randomly located in a 100 × 100 lattice and interconnecting with each other with a given probabilityrandom_p (0 ≤random_p ≤ 1). The units constituting the whole organization faced the same environment defined by a landscape described in formulas (1) and (2). We assumed that every unit gj∈G (j = 1, 2, …,unit_num) was considered a heterogeneous agent dealing with a series of organizational subtasks OT=ot1,ot2,…,otN and having a correspondingly specific routine ri∈R=r1,r2,…,rN for each task oti∈OT  (i = 1, 2, …,N) [71]. All routines were interrelated with each other, resulting in interrelationships both within and between bundles of routines, as shown in Figure 2.Figure 2
Network of units with different routine networks.Furthermore, we assumed that each routine could generate only two different patterns of action, “0” and “1”, respectively. Although a specific routine may have more than two states, a set of routines could also be interpreted “as subroutines that constitute a routine of more than two states” [71].We assumed that each unit agent had attributes as shown in Table1.Table 1
Definition of agent attributes.
Attributes Definition unitID A unique ID number for agent identification.  unitNeighbors Every unit agent owned a list of neighbors defined by the intraorganizational network structure. Here, we assumed a classic random network following Erdös and Rényi [74], where nodes represented unit agents and linkages between nodes – which represented their neighborhoods – were randomly assigned by the probability value random_p (0 ≤ random_p ≤ 1). Given the bounded rationality characterizing organizational actors [75], we assumed that each unit agent could only interact with and replicate routines from nearest neighbors.  networkOfRoutines Each unit agent had a bundle of routines (i.e., a certain specific vector with “0” or “1” values) to cope with its series of organizational subtasks. However, rather than being only a simple collection, these routines were interrelated with each other [11, 39], so giving rise to a network of routines. We assumed that this network was randomly generated [43, 72], as illustrated in Figure 1.  Fitness The value of this parameter represented the adaptiveness of a unit agent to its living environment. In detail, we considered each routine’s fitness contribution functionf⋅ shown in formula (1) as a random function mapping the vector ri,r1,r2,…,rK (where, r1,r2,…,rK∈RKi and RKi⊂R, ri∉RKi, i=1,2,…,N) with uniformly distributed values fiti (0 ≤ fiti ≤ 1). Then, we assumed that the fitness value of the unit agent was numerically equal to the average of the fit value of all the routines involved (as shown in formula (2)).However, as we mentioned in Section2.1, unit agents often used templates when replicating routines from neighbors. We assumed that a template could not be the exact copy of the network of routines of the replicator unit agent. Indeed, given that at least a part of knowledge embedded in routines is to some extent tacit and enacted in practice [32, 47], these differences mimicked in difficulties for replicatee units to explicitly encode routines into templates.More specifically, we assumed that when replicating activities occurred, the replicatee unit agents obtained a template which was duplicated from the network of routines of the replicator unit agents but with some given probability of errorserror_p (0 ≤error_p ≤ 1). This allowed us to mimic errors that may occur during codifying and explicitly representing process of the original set of organizational routines [18, 25]. However, note that in the model we considered the value oferror_p as constant throughout all simulations. Considering reasons and/or even mechanisms for this kind of errors was out of the scope of our work.Finally, we considered theenvironment as a regime of possibilities and constraints for unit agents while accomplishing their organizational task [22]. This environment determined both the organizational subtask series and the fitness landscape as shown in formulas (1) and (2). In order to simulate environmental dynamics, in our model, we adopted a Boolean variableis_env_change? to represent whether the fitness landscape changed or not. On the one hand, whenis_env_change? =TRUE, the fitness landscape would randomly vary once after every certain number of simulation roundsenv_change_P so as to mimic environmental changes occurred in practice. Further, for anyenv_change_P≥ 1, we have that the lower the value of this variable was, the more rapid environmental changes happened. On the other hand, whenis_env_change? =FALSE, it meant that there was no change of the organizational environment. Thus, any values of the variableenv_change_Pare meaningless in this situation.


## 3.3. Rules of the Replicating Activities

According to Hodgson [50], replicating routines from peers is not easy often also due to organizational inertia (e.g., [71, 76]) and bounded rationality [75]. To mimic this, we assumed that unit agents were prone to replicate “superior practices” from neighbors only when they felt competitive pressures, i.e., they had a significantly lower fitness value than the average of neighbors.Next, atemplate was generated through explicitly codifying the knowledge embedded in the network of routines of the replicator unit agent [19, 21]. In this model, we used a vector of “0” or “1” values as follows: (3)TPL=tpl1,tpl2,…,tplNbeing templates of a copy of the original series of routines, but allowing each of the element values of the vectortplk (0 ≤k ≤N) to vary with a certain probability of errorserror_p (0 ≤error_p ≤ 1). Namely,(4)tplk=1-rk,withaprobabilityerror_prk,withaprobability1-error_pHere, the value of the coefficient rk was either “0” or “1”, which represented the two different states (i.e., patterns of actions) of thekth original routine of the replicator unit agent, and rk∈R=r1,r2,…,rN.Furthermore, given that replicatee unit agents were always constrained by a specific level of absorptive capacity [77, 78], we assumed that they could not adopt the whole template but only a certain part. Here, we used a variableunit_AC to represent the absorptive capacity of unit agents. Its value numerically equaled the number of routines that replicatee unit agents could change at a time according to the template,unit_AC = {1,2,3}. Obviously, the higher the value ofunit_AC, the higher level of the unit agent’s absorptive capacity. In this respect, note that exploring details of agents’ absorptive capacities was beyond the scope of our study. For the sake of simplicity, we assumed that the value ofunit_AC was assigned at the beginning of the simulation and remained fixed throughout our simulations.The replicatee unit agents could randomly select a number of routines (which numerically equaled theunit_AC value), which owned different states with those from the template – if they existed, of course. Then, the replicatee unit agents changed the relating routines’ states (values) according to the template and generated a new list of routinesR’.Finally, when the new collection of routinesR’ was formed, replicatee unit agents checked if they could perform any innovation to match their work of routines (obviously, here, referring to the list of routinesR’) and the specific context in which they were living. This was to mimic the “context-dependent” characteristic of organizational routines [22] as well as the interconnection between routines [15, 18, 20, 39]. For example, replicatee units could revise parts of routines while considering both the embeddedness of the system of routines [11] and certain fitness landscapes – namely, the units’ living environment [22, 42]. By doing these innovative activities, replicatee unit agents transformed the temporary list of routinesR’ into their final network of routinesR with newly updated status.In our model, we assumed that replicatee unit agents could innovate and randomly vary the state of one of their routines with a given probabilityinnovate_p (0 ≤innovate_p ≤ 1) – when this was helpful to improve their fitness values.Figure3 provides a description of the routine replicating mechanisms between different unit agents. Figure 4 shows a stylized description of the procedure in which unit agents involved in the model should execute during the whole simulation process.Figure 3
Mechanisms of replicating routines between unit agents (unit_AC = 2).Figure 4
Flow chart of replicating activities between unit agents.


## 4. Simulation Experiments, Initialization, and Verification of the Model

### 4.1. Simulation Experiments

We constructed our model using NetLogo (https://ccl.northwestern.edu/netlogo/). All the main input settings were shown in Table 2. By applying the unit testing method [79] to each program module and to the whole simulation model, we ensured that our model was comprehensively verified.Table 2
Model input settings.
Input parameters Default settings Notes unit_num 20,40, 60, 80 The number of unit agents involved in the organization  random_p 0.1, 0.3, 0.5, 0.7, 0.9 We adopted a random network model to describe connections between unit agents. Here, the parameter represented the probability of links between two units  N_routine 4, 6, 8, 10 The number of routines involved in the organization  A_routine 2 The number of states for every routine  K_routine 1, 3, …, N_routine – 1 The number of routines that affected one certain routine’s contribution to the unit’s fitness  unit_AC 1, 2, 3 Absorptive capacity of each unit agent  error_p 0, 0.1, 0.3, 0.5, 0.7, 0.9 The probability of errors occurring when replicatee units duplicated routines from a template  innovate_p 0, 0.2, 0.4, 0.6, 0.8 The probability of innovative activities for each unit agent  is_env_change? FALSE, TRUE A Boolean variable representing whether the environmental changes existed or not  env_change_P 200, 1000 The simulation rounds period of environmental changes Note: the underlining values are default settings of the input parameters.The aim of this paper was to investigate the replicating dynamics of routines under the constraints of both the intraorganizational network of unit agents and the routine networks. Furthermore, we considered four other factors, namely, the absorptive capacity of unit agents, the probability of errors occurring during the template-duplicating processes and that of unit agents’ innovative activities, and the dynamic characteristic of the organizational environment.To do so, we first considered one of the simplest cases (“Scenario I”). By assuming that there was no environmental change (in this model, we let the Boolean variableis_env_change? =FALSE, referring to no changes occurred of the fitness landscape) and that there are neither template copying errors nor innovative activities (i.e.,error_p =innovate_p = 0), we purely investigated the influence of both of the two different networks – namely, one is the intraorganizational network of unit agents (determined by the model parameterrandom_p) and the other is the network of routines (represented by the model parameterK_routine) – on routines’ replication dynamics.Then, in “Scenario II”, we added variability. On the one hand, we introduced template-duplicating errors (represented by a series of values of the parametererror_p). On the other hand, we assumed that success probability of innovative activities by different replicatee units (represented by variations of the parameter valuesinnovate_p) could result in variations of the routine system.Next, in “Scenario III”, we introduced uncertainty of the organizational environment, by varying model parameter valuesis_env_change? =TRUE, andenv_change_P, so varying the changing rate of the fitness landscape function. This allowed us to examine routine replication within the dynamic organizational environment.Furthermore, in “Scenario IV”, we combined Scenarios II and III to analyze routine replication in a more complex environment, in which internal sources of variability, i.e., template-duplicating errors and innovative activities, were coupled with external sources, i.e., uncertainty of the organizational environment. It is worth noting that these two sources of variations are typical of complex, dynamic environments in the current business landscape, as suggested by Goldstein, et al. [80].Finally, we developed four aggregate measures to estimate the impact of scenario manipulations. The first one, calledsteadyTime, measured the number of simulation rounds required by the organization system to arrive at a steady state (e.g., every unit agent found the best bundle of routines that matched its own landscape, and thus there was no replication activity left), if such a specific steady state eventually existed, of course. This is crucial to understand the performance of routine replications as it indicates the replicating speed of routines between different units within the organization. The second one is the average fitness values of all the unit agents. This represented the organizational efficiency that every unit could gain by routine replication, so in principle benefiting the whole organization.We then also measured the number of replicating and innovative activities occurred among all the unit agents (here, denoted asreplicate_num andinnovate_num, respectively). When copying routines from their neighbors, unit agents bore some operating cost, i.e., the expenses or resources required to replicating or innovative activities. In general, we hypothesized that the more these activities were performed by unit agents, the higher the operating cost was. In this sense, the number of replicating and innovative activities occurred together could reflect the level of operating cost required by the organization to perform its internal replicating activities.


### 4.2. Initialization and Verification of the Model

For verifying the correctness of the model, we first considered a very basic scenario in which there were neither variations from template-duplicating errors and innovative activities nor changes in the organizational environment. We let the parameter valueunit_num = 20, 40, 60, and 80 and keep other variables in the default settings as shown in Table 2. This led to four corresponding but definitely different parameter configurations. In order to determine the ideal number of simulation run for each parameter configuration, we applied the statistical power analysis suggested by Secchi and Seri [81, 82]. Given the significance level α = 0.01, the power 1 - β = 0.95, and the effect size with a small-levelES = 0.1 according to discussions in [81], we found that 569 runs were sufficient. Then, using the ‘smallest effect size of interest’ (SESOI) approach, we set simulation runs to 600 times per parameter configuration, which ensured an estimated expected power of 0.96.On the other hand, after some pilot runs to study convergence and sensitivity, we decided to let each simulation running for 3000 rounds, each of which representing one routine performance. This allowed us to find a convergence of the routine system to a given steady state, at the same time enabling us to measure one of the most important outputs, i.e., the variablesteadyTime.We then collected output data and calculated the mean values of each variable in each simulation run. Figure5 shows a summary of simulation results. Figure 5(a) indicates that with the increasing of the number of units involved in the organization, there was a continuous variation of the organization’s states, i.e., the values of the three main outputs,steadyTime,replicate_num, andavg_Fitness. It is worth noting here that we did not consider the variableinnovate_num as thatinnovate_p = 0 and no innovative activities occurred. First, when connections between units were kept at a very lower level (random_P = 0.1) and there were only a few units involved (i.e.,unit_num = 20), replicating activities occurred randomly but with a very small mean value (see Figure 5(c)). The organization took a short time to reach a steady state (see Figure 5(b)), and theavg_Fitness values followed a very standard normal distribution with mean value approximately equaling 0.5 (see Figure 5(d)). This corroborated our idea of initializing routine networks of all the units randomly.Figure 5
Simulation outputs of the testing scenario: (a) a sketch of the three output variables; (b) distribution ofsteadyTime; (c) distribution ofreplicate_num; (d) distribution ofavg_Fitness.
(a) (b) (c)  (d)Furthermore, the increasing number of units obviously increased their interactive opportunities and thus led to the emergence of certain specific patterns of replicating activities. More specifically, the distribution of the variablereplicate_num obtained a higher mean value but with a smaller range of variance (see Figure 5(c)). Then, thesteadyTime values increased, and the distribution of the variableavg_Fitness followed a partial peak distribution with a bigger mean value (see Figures 5(b) and 5(d)). Again, all these results confirmed the validity of our logics in building and parameterizing the model.


## 4.1. Simulation Experiments

We constructed our model using NetLogo (https://ccl.northwestern.edu/netlogo/). All the main input settings were shown in Table 2. By applying the unit testing method [79] to each program module and to the whole simulation model, we ensured that our model was comprehensively verified.Table 2
Model input settings.
Input parameters Default settings Notes unit_num 20,40, 60, 80 The number of unit agents involved in the organization  random_p 0.1, 0.3, 0.5, 0.7, 0.9 We adopted a random network model to describe connections between unit agents. Here, the parameter represented the probability of links between two units  N_routine 4, 6, 8, 10 The number of routines involved in the organization  A_routine 2 The number of states for every routine  K_routine 1, 3, …, N_routine – 1 The number of routines that affected one certain routine’s contribution to the unit’s fitness  unit_AC 1, 2, 3 Absorptive capacity of each unit agent  error_p 0, 0.1, 0.3, 0.5, 0.7, 0.9 The probability of errors occurring when replicatee units duplicated routines from a template  innovate_p 0, 0.2, 0.4, 0.6, 0.8 The probability of innovative activities for each unit agent  is_env_change? FALSE, TRUE A Boolean variable representing whether the environmental changes existed or not  env_change_P 200, 1000 The simulation rounds period of environmental changes Note: the underlining values are default settings of the input parameters.The aim of this paper was to investigate the replicating dynamics of routines under the constraints of both the intraorganizational network of unit agents and the routine networks. Furthermore, we considered four other factors, namely, the absorptive capacity of unit agents, the probability of errors occurring during the template-duplicating processes and that of unit agents’ innovative activities, and the dynamic characteristic of the organizational environment.To do so, we first considered one of the simplest cases (“Scenario I”). By assuming that there was no environmental change (in this model, we let the Boolean variableis_env_change? =FALSE, referring to no changes occurred of the fitness landscape) and that there are neither template copying errors nor innovative activities (i.e.,error_p =innovate_p = 0), we purely investigated the influence of both of the two different networks – namely, one is the intraorganizational network of unit agents (determined by the model parameterrandom_p) and the other is the network of routines (represented by the model parameterK_routine) – on routines’ replication dynamics.Then, in “Scenario II”, we added variability. On the one hand, we introduced template-duplicating errors (represented by a series of values of the parametererror_p). On the other hand, we assumed that success probability of innovative activities by different replicatee units (represented by variations of the parameter valuesinnovate_p) could result in variations of the routine system.Next, in “Scenario III”, we introduced uncertainty of the organizational environment, by varying model parameter valuesis_env_change? =TRUE, andenv_change_P, so varying the changing rate of the fitness landscape function. This allowed us to examine routine replication within the dynamic organizational environment.Furthermore, in “Scenario IV”, we combined Scenarios II and III to analyze routine replication in a more complex environment, in which internal sources of variability, i.e., template-duplicating errors and innovative activities, were coupled with external sources, i.e., uncertainty of the organizational environment. It is worth noting that these two sources of variations are typical of complex, dynamic environments in the current business landscape, as suggested by Goldstein, et al. [80].Finally, we developed four aggregate measures to estimate the impact of scenario manipulations. The first one, calledsteadyTime, measured the number of simulation rounds required by the organization system to arrive at a steady state (e.g., every unit agent found the best bundle of routines that matched its own landscape, and thus there was no replication activity left), if such a specific steady state eventually existed, of course. This is crucial to understand the performance of routine replications as it indicates the replicating speed of routines between different units within the organization. The second one is the average fitness values of all the unit agents. This represented the organizational efficiency that every unit could gain by routine replication, so in principle benefiting the whole organization.We then also measured the number of replicating and innovative activities occurred among all the unit agents (here, denoted asreplicate_num andinnovate_num, respectively). When copying routines from their neighbors, unit agents bore some operating cost, i.e., the expenses or resources required to replicating or innovative activities. In general, we hypothesized that the more these activities were performed by unit agents, the higher the operating cost was. In this sense, the number of replicating and innovative activities occurred together could reflect the level of operating cost required by the organization to perform its internal replicating activities.


## 4.2. Initialization and Verification of the Model

For verifying the correctness of the model, we first considered a very basic scenario in which there were neither variations from template-duplicating errors and innovative activities nor changes in the organizational environment. We let the parameter valueunit_num = 20, 40, 60, and 80 and keep other variables in the default settings as shown in Table 2. This led to four corresponding but definitely different parameter configurations. In order to determine the ideal number of simulation run for each parameter configuration, we applied the statistical power analysis suggested by Secchi and Seri [81, 82]. Given the significance level α = 0.01, the power 1 - β = 0.95, and the effect size with a small-levelES = 0.1 according to discussions in [81], we found that 569 runs were sufficient. Then, using the ‘smallest effect size of interest’ (SESOI) approach, we set simulation runs to 600 times per parameter configuration, which ensured an estimated expected power of 0.96.On the other hand, after some pilot runs to study convergence and sensitivity, we decided to let each simulation running for 3000 rounds, each of which representing one routine performance. This allowed us to find a convergence of the routine system to a given steady state, at the same time enabling us to measure one of the most important outputs, i.e., the variablesteadyTime.We then collected output data and calculated the mean values of each variable in each simulation run. Figure5 shows a summary of simulation results. Figure 5(a) indicates that with the increasing of the number of units involved in the organization, there was a continuous variation of the organization’s states, i.e., the values of the three main outputs,steadyTime,replicate_num, andavg_Fitness. It is worth noting here that we did not consider the variableinnovate_num as thatinnovate_p = 0 and no innovative activities occurred. First, when connections between units were kept at a very lower level (random_P = 0.1) and there were only a few units involved (i.e.,unit_num = 20), replicating activities occurred randomly but with a very small mean value (see Figure 5(c)). The organization took a short time to reach a steady state (see Figure 5(b)), and theavg_Fitness values followed a very standard normal distribution with mean value approximately equaling 0.5 (see Figure 5(d)). This corroborated our idea of initializing routine networks of all the units randomly.Figure 5
Simulation outputs of the testing scenario: (a) a sketch of the three output variables; (b) distribution ofsteadyTime; (c) distribution ofreplicate_num; (d) distribution ofavg_Fitness.
(a) (b) (c)  (d)Furthermore, the increasing number of units obviously increased their interactive opportunities and thus led to the emergence of certain specific patterns of replicating activities. More specifically, the distribution of the variablereplicate_num obtained a higher mean value but with a smaller range of variance (see Figure 5(c)). Then, thesteadyTime values increased, and the distribution of the variableavg_Fitness followed a partial peak distribution with a bigger mean value (see Figures 5(b) and 5(d)). Again, all these results confirmed the validity of our logics in building and parameterizing the model.


## 5. Simulation Results

Table2 shows the set of conditions after a series of analyses performed on a testing scenario (as mentioned above) and additional “trial-and-error” experiments. We decided to report our simulation experiments in four subsections corresponding to the four scenarios included in Section 4. To avoid either under- or overpowering our simulation results [81], we let the significance level α = 0.01, the power 1 - β = 0.95, while assuming the effect size with a micro levelES = 0.05 due to a pilot study and previous testing experiments. For each of our four scenarios, we followed a full factorial design and calculated the ideal number of simulation runs (denoted asi_NofRun) needed for each scenario – as shown in Table 3. Furthermore, following the SESOI approach [81, 82], we set certain running times, seer_NofRun, for each parameter configuration, resulting in an estimated expected power (EEP) of 0.96. Moreover, we kept each simulation run at 3000 rounds to ensure that the organization system reached a steady state – if this steady state definitely existed.Table 3
Calculating running times needed for each scenario based on power analysis.
Scenarios Configuration of parameters (CoP) i_NofRun r_NofRun EEP Scenario I(is_env_change? = FALSE) Variables and values: 382 400 0.96 K_routine random_p unit_AC 1 0.1 1 3 0.3 2 5 0.5 3 7 0.7 0.9 CoP:Φ(K_routine) ×  Φ(random_p) ×  Φ(unit_AC) (1)  = 4 × 5 × 3 = 60  Scenario II(is_env_change? = FALSE) Variables and values: 382 400 0.96 K_routine error_p innovate_p 1 0.1 0.2 3 0.3 0.4 5 0.5 0.6 7 0.7 0.9 CoP:Φ(K_routine) ×  Φ(error_p) ×  Φ(innovate_p) = 4 × 5 × 3 = 60  Scenario III(is_env_change? = TRUE) Variables and values: 253 260 0.96 K_routine random_p unit_AC env_change_P 1 0.1 1 200 3 0.3 2 1000 5 0.5 3 7 0.7 0.9 CoP:Φ(K_routine) ×  Φ(random_p) ×  Φ(unit_AC) ×  Φ(env_change_P) = 4 × 5 × 3 × 2 = 120  Scenario IV(is_env_change? = TRUE) Variables and values: 253 260 0.96 K_routine error_p innovate_p env_change_P 1 0.1 0.2 200 3 0.3 0.4 1000 5 0.5 0.6 7 0.7 0.9 CoP:Φ(K_routine) ×  Φ(error_p) ×  Φ(innovate_p) ×  Φ(env_change_P) = 4 × 5 × 3 × 2 = 120 Notes: (1) here, the functionΦ(X) represented the number of different levels involved in the variable X.(1) Routine Replications within a Static Environment.We assumedunit_num = 40 andN_routine = 8 based on the pilot experimental results and assignedrandom_p = 0.1, 0.3, 0.5, 0.7 and 0.9, andK_routine = 1, 3, 5 and 7, to represent different density levels of both the network of units and that of routines, respectively. We letunit_AC = 1, 2 and 3 to represent different levels of the units’ absorptive capacity. For each configuration of these three parameters and with default settings of the other parameters as shown in Table 2 (here, we named this as Scenario I), we collected all the 400 simulation running results and calculated mean values of the three outputs, i.e.,steadyTime,replicate_num, andavg_Fitness but excluded the variableinnovate_num because we removed any possibility of innovative activities due to the parameterinnovate_p = 0 in this scenario.Figure6 shows simulation results of the other three outputs – namely,steadyTime,replicate_num, andavg_Fitenss. We found that the network density of both units and routines had an impact on routine replication processes, so leading to different performance of the organization. However, their effects were also determined by another important variable, namely, the units’ absorptive capacityunit_AC. First, when the parameterunit_AC = 1, units were limited within a very narrow space when searching target neighbors for replication. With this constraint, even increasing the unit network density (random_p from 0.1 to 0.9) could increase the probability for units to link with each other, and all the three outputssteadyTime,replicate_num, andavg_Fitness could increment only minimally. On the other hand, increasing the routine network density (K_routine from 1 to 7) had no significant influence on the outputssteadyTime. However, it improved bothreplicate_num andavg_Fitness with a limited margin – except at a very medium valueK_routine = 3 when the unit network densityrandom_p = 0.7, which resulted in a decrease of theavg_Fitness value (see the left part in Figures 6(a), 6(b), and 6(c)).Figure 6
Simulation outputs within a static environment (Scenario I): (a)steadyTime values; (b)replicate_num values; (c)avg_Fitness values.
(a) (b) (c)Next, when the parameterunit_AC = 2 and 3, we found that the organization system spent a longer time to arrive at a steady state (i.e., a highersteadyTime value, see the medium and right parts in Figure 6(a)). More replicating activities occurred (i.e., a highreplicate_num value, see the medium and right parts in Figure 6(b)), and so a higher level of fitness was reached (i.e., a highavg_Fitness value, see the medium and right parts in Figure 6(c)). In this situation, the lower the value of the variablerandom_pwas, the greater the effects of its increments to the organizational system were.On the other hand, although variation ofK_routine values had a weak effect when the variablerandom_p took low values (e.g.,random_p = 0.1 and 0.3 whenunit_AC = 2; orrandom_p = 0.1 whenunit_AC = 3), achieving an optimal fitness when taking a medium value (here, for example,K_routine = 3 or 5) and the variablerandom_p with higher values (i.e.,random_p = 0.5, 0.7 and 0.9) was extremely beneficial for the organization. Obviously, in this case, reaching a steady state required more time.However, if we increasedK_routine to a significantly higher value (e.g., in our model,K_routine = 7), the system took less time to reach a steady state, but the increasing values ofreplicate_num and the decreasing values ofavg_Fitness indicate that the organization undertook more replication activities. This permitted to overcome resistance against change due to the strong ties between different routines, yet with no obvious benefits for organization adaptiveness.To sum up, if unit agents had only very limited absorptive capacity, variations of both the unit and routine networks’ density do not improve organizational performance. If we consider the cost for unit agents to manage their relationship with partners, the contribution of the replicating process of organizational routines to the organization performance depends on an optimal configuration of the two types of networks – i.e., the network of units in the organization represented by the variablerandom_p and that of routines represented by the variableK_routine.(2) Routine Replications with Template-Duplicating Errors and Units’ Innovative Activities. In this section, we considered two important factors of variation. One is the probability of errors that may occur when we duplicate routines from the template. The other one is the probability of innovation for unit agents while they replicate routines from their partners. We assumed the values of the parametererror_p = 0.1, 0.3, 0.5, 0.7, and 0.9 and another oneinnovate_p = 0.2, 0.4 and 0.6, representing different levels of the two kinds of events’ occurrence probabilities, respectively. According to simulation results of Scenario I (as shown in Figure 6), we selected one of the most representative but not extreme cases by lettingrandom_p = 0.5 andunit_AC = 2 and kept the other default settings as Scenario I. Then, we separately investigated the impact of variation from both template-duplicating errors and innovative activities of unit agents.Figure7 shows that these induced variations created opportunities for replicatee units with significant benefit to the fitness of the whole organization. However, this came to the cost of more replicating and innovative activities required and so more time is required by the organizational system to reach a steady state. First, increasing network density of routines (i.e.,K_routine) meant a tighter coupling relation between these routines. Therefore, the organizations spent more time and were more involved in replicating and innovative activities (see Figures 7(a), 7(b), and 7(c)). However, we found that when the variableerror_p equaled low values (e.g.,error_p = 0.1 and 0.3), a medium value ofK_routine (i.e.,K_routine = 3) still existed which was beneficial to the organization’s adaptiveness (see Figure 7(d)).Figure 7
Simulation results with template-duplicating errors and units’ innovative activities (Scenario II:random_p = 0.5 andunit_AC = 2): (a)steadyTime values; (b)replicate_num values; (c)innovate_num values; (d)avg_Fitness values.
(a) (b) (c)  (d)Second, in case of variation from template-duplicating errors, on the one hand, the higher was the value of the variableerror_p (e.g.,error_p = 0.9), the more randomness occurred during the replicating activities, which could result in higher heterogeneity of routine networks between unit agents. On the other hand, very low values oferror_p (e.g.,error_p = 0.1) caused less fluctuations to the whole routine system. Therefore, these two cases were both beneficial for the organization to economize on the time required to reach a steady state (as shown in Figure 7(a)). Further, increasinginnovate_p values mainly led to more innovative activities among unit agents (see Figure 7(c)). This had a significant influence both onsteadyTime andreplicate_num only whenerror_p equaled low values (see Figures 7(a) and 7(b)). It also influenced theavg_Fitness whenerror_p did not equal medium values – i.e.,error_p = 0.1 or 0.9 (see Figure 7(d)).Furthermore, variations from both template-duplicating errors and unit innovations led to an increase of replicating and innovative activities at the same time reducing average fitness values. However, whenK_routine values increased from 1 to 7, especially whenK_routine = 7, the influence of these variations decreased. In these cases, for differenterror_p values (from 0.1 to 0.9), the three outputssteadyTime,innovate_num, andavg_Fitness tended to converge (see Figures 7(a), 7(c), and 7(d)).(3) Routine Replications in Dynamic Environments.In this section, we investigated routine replicating processes in a dynamic environment by letting the variablesis_env_change? =TRUE andenv_change_P = 200 and 1000, respectively. This was to mimic different changing rates of the organizational environment. For the rest, we kept the other input settings as in Scenario I. In this case, given that environmental changes were pervasive, we expected that the organization system as a whole could not arrive at any steady state. Given that the output valuesteadyTime had no meaning here and that by assuming the variableinnovate_p = 0 the output valuesinnovate_num always equaled 0, we reported here only two simulation outputs includingreplicate_num andavg_Fitness.Figure8 shows that when the fitness landscape changed, we found that the absorptive capacity level of unit agents has a strong effect on routine replicating activities so determining the average fitness value of all the units involved.Figure 8
Simulation results in dynamic environments (Scenario III): (a)replicate_num values; (b)avg_Fitness values.
(a) (b)First, when absorptive capacity was minimal (e.g.,unit_AC = 1), unit agents were very constrained to search for appropriate routine candidates from their neighbors. Independently from the degree of connections between routines (i.e., theK_routine values) or that between unit agents (i.e., therandom_p values), and the changing rate of the organizational environment (i.e., the differentenv_change_P values), the average number of unit agents involved in replicating activities (i.e.,replicate_num) and the average fitness of all the unit agents (i.e.,avg_Fitness) did not have a significant effect. In addition, we also found that environmental changes influenced both these two outputs even more with the connection density between units increased (e.g., increasing therandom_p values from 0.1 to 0.9) (see the left part of Figures 8(a) and 8(b)).Secondly, higher levels of absorptive capacity (e.g.,unit_AC = 3) implied that every unit could easily find its ideal routine targets for replication regardless of the pace of environmental changes. To some extent, this explained why variations ofenv_change_P values had less impact on thereplicate_num values (see the right part of Figure 8(a)). Further, although bothreplicate_num andavg_Fitness were very sensitive to variations of the connection degree between units in case of low values ofrandom_p (e.g.,random_p = 0.1), the variation ofrandom_p had no effect when unit agents were highly connected with each other (i.e.,random_p = 0.5, 0.7 and 0.9).Furthermore, when the variableunit_AC equaled a medium value (e.g.,unit_AC = 2), we found a correlation between the increasing rapidity of environmental changes (e.g., decreasingenv_change_P from 1000 to 200) and less replicating activities (see the middle part of Figure 8(a)). This was especially true when the network density of unit agents and routines was high (e.g.,K_routine = 5 and 7 andrandom_p = 0.7 and 0.9).Finally, Figure8(b) shows that whenunit_AC = 2 and 3, the organization’s fitness (i.e., theavg_Fitness values) decreased when organizational environment varied more rapidly (e.g., decreasingenv_change_P from 1000 to 200) (see the middle and right parts of Figure 8(b)). In addition, we found a positive relationship betweenK_routine andavg_Fitness, except when the unit network density was low (e.g.,random_p = 0.1). In this cases, there was an optimalavg_Fitness whenK_routine had a medium value (i.e.,K_routine = 3), similarly to what occurred in a static environment.(4) Routine Replications with Internally and Externally Coupled Uncertainties.In this section, we combined Scenarios II and III to examine the influence of coupled internal variations and external uncertainties. To do this, we set uprandom_p = 0.5;unit_AC = 2;error_p = 0.1, 0.3, 0.5, 0.7, and 0.9;innovate_p = 0.2, 0.4 and 0.6;is_env_change? =TRUE; andenv_change_P = 200 and 1000, while the other input parameters were the same default settings of scenarios I, II, and III.Figure9 shows that within a dynamic environment, more replicating and innovative activities were required to cope with more coupling of routines (i.e., routines’ connection density parameterK_routine). Moreover, we found that although more template-duplicating errors (i.e., the variableerror_p) had less impact on unit agents’ innovative activities, this caused more randomness to the routines system, and thus required more routine replicating activities. On the other hand, while increasinginnovate_p values could only promote innovative activities of unit agents, this determined less influence of the number of unit agents involving in routine replicating activities (see Figures 9(a) and 9(b)). This was especially true when errors in template-duplication were low (e.g.,error_p = 0.1 and 0.3).Figure 9
Simulation results with internally and externally coupled uncertainties (Scenario IV:random_P = 0.5 andunit_AC = 2): (a)replicate_num values; (b)innovate_num values; (c)avg_Fitness values.
(a) (b)  (c)Secondly, when the organizational environment changed rapidly, density variations in unit agents and routine networks did not have any significant influence on the level of innovative activities, except when less template-duplicating errors (e.g.,error_p = 0.1) were matched with a slowly changing environment (e.g.,env_change_P = 1000). In this situation, any increase of routine coupling (e.g.,K_routine values varied from 1 to 7) led to more innovative activities. However, differences of both replicating and innovative activities caused by variation of eitherK_routine orerror_p values decreased with the increasing rapidity of environmental changes.Figure9(c) shows that environmental changes had a negative effect on organization adaptiveness. Furthermore, in case internal and external variations were not particularly prominent (i.e.,error_p = 0.1 and 0.3, andenv_change_P = 1000), medium values ofK_routine (e.g.,K_routine = 3) were still beneficial to the optimal fitness of the organization. Otherwise, we found a positive relation betweenK_routine andavg_Fitness values whenK_routine = 1, but this relation was not significant whenK_routine value was high (e.g.,K_routine = 5).Furthermore, in a slowly changing environment (e.g.,env_change_P= 1000), with the increase of template-duplicating errors (i.e., increasingerror_p from 0.1 to 0.9), the average fitness values decreased. However, this negative relation decreased when the speed of environmental changes was higher (see the bottom part of Figure 9(c)). On the other hand, no matter whether the organizational environment changed slowly or rapidly, higher values ofinnovate_p only resulted in more innovative activities (i.e., the increase ofinnovate_num values) as a means to cope with the randomness caused by both internal variations and external dynamics, and this did not have significant effects on the outputavg_Fitness (see the medium and right parts of Figure 9(c)).Finally, by comparing the bottom and upper parts of Figures9(a) and 9(c), we found when the organizational environment was moderate (e.g.,env_change_P = 1000), the routine system could reach a higher level of fitness (i.e., theavg_Fitness values) but with less replicating activities in contrast to that within a dramatically changing environment (e.g.,env_change_P = 200). These findings are consistent with Rivkin’s arguments in [51].


## 6. Discussions and Conclusions

In this paper, we considered replication of organizational routines as some specific “template-based” activities between different units within an organization [19, 21, 30] and developed an agent-based model to simulate the routine replicating dynamics from the “bottom-up” approach [83]. We considered two distinguished types of network contexts, the first one referred to connections between (geographically) distributed units [24] and the other one was the coupling relation between different routines [11, 36, 60, 63]. By depicting the former as a random network model [74] and the latter as a NK-based fitness landscape [43], we constructed an agent-based model that allowed us to zoom in the routine replicating process in detail. At the same time, the model permitted to consider the impact of structural factors, such as the two different types of network contexts, the internal variations from both template-duplicating errors and innovative activities of unit agents involved, and the uncertainty of environmental changes (more detailed panel regression results can be found in the Appendix).Our findings showed that in a relatively static environment, a functional replication of routines within a decentralized organization depends on absorptive capacity of unit agents [18, 24]. Only by an appropriate level of absorptive capacity, each unit can have opportunities to replicate the “best practices” from their neighbors. In this case, an optimal combination of the two different types of networks – namely, both the connections between different units and the coupling relations between routines – can be achieved which is beneficial to routine replicating practices that help the organization to adapt to its living environment.Then, we also found that certain degree of internal variations either from template-duplicating errors or innovative activities of units involved could create opportunities for unit agents to develop “performative drifts” of organizational routines [32, 84]. While this is beneficial to the adaptiveness of the whole organization, it also comes at operating expenses. However, it is worth noting that a higher level of variations could be harmful to routine replication dynamics [12, 85].Next, we found that although environmental changes require more routine replicating activities, these processes can be detrimental to the whole organization’s adaptiveness [59, 86]. When the environment is dynamic and uncertain, both routine replicating processes and the performance of the organization are sensitive to smaller connecting density values either of the network of unit agents or that of routines.This implies that increasing connections between different units are not always beneficial to the organization. Managing unit relations always requires balancing costs and benefits [42, 87]. On the other hand, our findings show that the coupling relation of routines is key to cope with the changing pace of the organizational environment [11].Finally, we found that a certain level of variation from either template-duplicating errors or innovative activities can help organizations to cope with environmental changes better [18]. Here, developing platforms to exploit a certain degree of redundancy and uncertainty from the bottom-up can increase the capacity of an organization to leverage on emergent processes and so cope with turbulent environments better (e.g., [88]). Obviously, this depends on the extent of internal variations, as these could increase operating costs without benefiting organizational performance.This said our study has important limitations. First, our agent-based model of routine replication dynamics is based on a series of pioneer work (e.g., [15, 17–21]) that inherited Nelson and Winter’s (1982) seminal book [8]. This implies that we focused on routines that are beneficial for organizations’ survival and development; i.e., by reconfiguring these routines, creating new ones, and using them consistently across organizational sites, organizations can adapt to a changing environment and obtain a high-level of performance [47–49]. However, “defensive routines” as mentioned in [5–7] are beyond the scope of our discussion. Second, we only limited our concentration on the network characteristic of distributed units within a decentralized organization – i.e., the organization within which activities (particularly those regarding decision-making) are distributed or delegated away from a central, authoritative leadership – rather than considering factors such as leadership and culture [73, 80]. These factors are, of course, very crucial ones for organization studies, and we would discuss their effects on routine dynamics in detail in another forthcoming paper. Third, the template-based replication of routines between different units is often linked with some kinds of artifacts – i.e., “physical manifestations” of the performance of organizational routines [89]. These artifacts play an important role in routine dynamics [20, 90, 91]. In this sense, template-duplicating activities can incorporate very complicated interactions between human actors and their material contexts, which are channelled by relating artifacts. Investigating the role of different kinds of artifacts in routine replicating dynamics could complement and enrich our findings. Last but not least, our agent-based model is only atypification and so should be verified with context-specific details [92, 93].



---

*Source: 8496235-2018-11-11.xml*