“…We label these groups 1, …, κ and let group g have n (g) experimental units within it. We model the responses of experimental units by the block network model (BNM), which is an extension of the LNM (Parker et al, 2016) and is described by the equation where i = 1, 2, …, ; j = 1, 2, …, n (i) , y ij is the continuous response from unit j in the ith block receiving the treatment s = r(ij) ∈ {1, …, m}, μ represents the response for a baseline treatment or (unit) average, b i is the effect of block i, r(ij) is the (direct) treatment effect, A {ij,gh} / A = A {ij,gh} with j, h ∈ V is the adjacency matrix indicating the edge between units j and h belonging to blocks i and g ∈ {1, 2, …, κ} respectively, r(gh) is the network effect (neighbour or indirect treatment effect) and ij are the errors, which we assume to be independent and identically distributed with mean 0 and constant variance 2 . To overcome the model overparametrisation requires imposing some constraints, otherwise the normal equations have an infinite number of solutions and our parameters cannot be uniquely estimated.…”