Neural Codes Constructs Based on Combinatorial Design

Abstract

Neuroscience (also known as neurobiology) is a science that studies the structure, function, development, pharmacology and pathology of the nervous system. In recent years, C. Cotardo has introduced coding theory into neuroscience, proposing the concept of combinatorial neural codes. And it was further studied in depth using algebraic methods by C. Curto. In this paper, we construct a class of combinatorial neural codes with special properties based on classical combinatorial structures such as orthogonal Latin rectangle, disjoint Steiner systems, groupable designs and transversal designs. These neural codes have significant weight distribution properties and large minimum distances, and are thus valuable for potential applications in information representation and neuroscience. This study provides new ideas for the construction method and property analysis of combinatorial neural codes, and enriches the study of algebraic coding theory.

Share and Cite:

Huang, J. (2025) Neural Codes Constructs Based on Combinatorial Design. Applied Mathematics, 16, 42-60. doi: 10.4236/am.2025.161003.

1. Introduction

In 1948, Shannon pioneered the theory of coding in interference channels. Based on Shannon theory, scholars have designed many error correcting codes with efficient decoding algorithms. Information theory has had a big impact on theoretical neuroscience. Mathematical coding theory received a lot of attention more than a decade ago [1] [2], and has been fueled by research on positional cell neurons. O’keefe and Dostrovsky’s discovery of positional cells was a major breakthrough. O’keefe was awarded both the 2014 Nobel Prize in Medicine and the Nobel Prize in Physiology. Position cells encode spatial information about an organism’s surroundings by firing precisely when the organism is located in the appropriate position field [3] [4]. A positional field is a collection of several convex open sets that may overlap each other, and the firing activity of a group of neurons produces a set of common firing patterns that can be represented as binary vectors or codewords. Each codeword represents a set of neurons stimulated at the same time, and the set of codewords on n neurons C is called a combinatorial neural code [5].

Based on the connection between neuroscience and algebraic coding, one can use coding to study neuroscience [5]-[7]. The earliest work originated in 2013, when C. Curto et al. defined a discrete receptive field code and modeled it as a binary code C { 0,1 } n . This led to the study of neural coding using algebraic and combinatorial methods [8] [9]. Let the n neurons in the brain be respectively { 1,2,3,,n } . Each neuron i has its receptive field S i , X= i1 n S i , Per stimulus xX . We got a code c( x )=( x 1 , x 2 , x n ) F 2 n , among others 1in . When x i S i , that is x stimulates neuron i , x i =1 , when x i S i . That is, neuron i is not stimulated by x , x i =0 of which zero codeword c( x )=( 0,0,,0 ) is the absence of a stimulus acting on the neuron. Let Σ be the set consisting of K stimulus. A subset C of the K codewords in F 2 n is called a combinatorial neural code (CN code for short) [10]. As shown in Figure 1, x is a stimulus and the codeword for x is c( x )=( 00110 ) .

Figure 1. c(x) = (00110).

Stimulus x codewords for c( x )=( x 1 , x 2 ,, x n ) F 2 n , but in practice, the received vectors may be c ( x ) =( x 1 , x 2 ,, x n ) F 2 n . It deviates from c( x )=( x 1 , x 2 ,, x n ) F 2 n . That is to say, there is x i S i , but it’s not being fired. Or x i S i , but shows that it is exciting. In reference [10], the following ideal mathematical model is proposed.

1) x i =1 and x i =0 probability of occurrence is q< 1 2 , this means that x i S i but the probability that neuron i is not stimulated is q .

2) x i =0 and x i =1 probability of occurrence is p and 0pq , this means that x i S i but the probability that neuron i is not stimulated is p .

From c( x ) to c ( x ) , consider the binary asymmetric memoryless channel transmission of the following Figure 2. If q=p , then the channel is a binary symmetric, but in general in neurology q<p , the channel is an asymmetric binary channel [11] [12]. Out of a desire to revert to the received codeword c ( x ) to the codeword c( x ) , Therefore, since the 1950’s, research on error-correcting codes in binary asymmetric channels has been initiated [12]-[14]. Further proposed new “Matched metrics” for asymmetric δ r ( x,y ) [15] to replace hamming distance, which is not really a distance.

Figure 2. Binary asymmetric memoryless channel.

The aim of this paper is to construct a special class of combinatorial neural codes with special properties using classical combinatorial structures such as Orthogonal Latin rectangle, disjoint Steiner systems, Group Divisible designs, and Transversal designs which have significant weight distribution properties and large minimum distances, and which are of greater value for applications in information representation and neuroscience. Relatively speaking the weight distribution of this paper is very special and the minimum distance is also reached to the maximum, the specific code words of the neural code have indicated, which is more conducive to the research of neuroscience and coding theory.

The paper is structured as follows, and in Section 2, the relevant basics and definitional citations are explained in detail; in Section 3, it is shown that combinatorial neural codes with weight 4 are constructed by orthogonal Latin rectangle, disjoint Steiner systems, Groupable Divisible and a class of combinatorial neural codes are constructed based on the transversal design; and finally, the conclusion of the paper tells about the main conclusions of the paper.

2. Preparatory Knowledge

This section introduces the concepts and properties related to antichain codes [16], it is a combinatorial neural code. The underlying lemmas used to construct the antichain code are described below.

Let A 2 [ n ] be a family of sets, if any distinct x,yA , with xy , A is said to be an antichain. For example, layer

( [ n ] k )={ x[ n ]:| x |=k } , (1)

For all 0kn form an antichain. From the classical theory of Sperner [17] it is known that the antichain in 2 [ n ] is at most ( n n 2 ) .

A family of vectors B { 0,1 } n is called a code with distance r , if any two vectors in B have a Hamming distance of at least r . Identifying { 0,1 } n and 2 [ n ] , we call a set family A a distance- r codes if the size of the symmetric difference xy of any two distinct x,yA is always at least r . In the literature [16], for any fixed rN , if A 2 [ n ] is an antichain and a distance- ( 2r+1 ) code, then | A |=Ο( 2 n n r 1 2 ) .

Definition 1 [18]: For y=( y 1 , y 2 ,, y n ) and x=( x 1 ,x,, x n ) F 2 n . Let

d 10 ( y,x )=#{ i|1in,( y i , x i )=( 1,0 ) }

d 01 ( y,x )=#{ i|1in,( y i , x i )=( 0,1 ) }

The asymmetric discrepancy of y and x is defined as

δ( y,x )= δ r ( y,x )=r d 10 ( y,x )+ d 01 ( y,x ) (2)

For a combinatorial neural code C F 2 n , the constituent parameters are ( n,K,δ( C ) ) , where size of C is K=| C |2 , and the minimum discrepancy is

δ( C )= δ r ( C )=min{ δ r ( c, c ):c, c C,c c }

For p=q , we get r=1 , and δ 1 ( y,x )= d 10 ( y,x )+ d 01 ( y,x )= d H ( y,x ) is the Hamming distance. Hence δ 1 ( C )= d H ( C ) , as the minimum Hamming distance of C. It is easy to find δ( y,x )0 and δ( y,x )=0 if and only if y=x . If p<q ( r>1 ) usually δ r ( y,x ) and δ r ( x,y ) are not intended to be the same, but are satisfied by the triangle inequality [4].

According to the definition of combinatorial neural codes, it is easy to know that the antichain code is a class of neural codes, so the antichain code with distance of constitutes a combinatorial neural code with parameters

( n,| A |=Ο( 2 n n r 1 2 ), δ r =4r+1 ) (3)

Definition 2 [19]: A m×n rectangle G , whose elements consist of a set of s elements such that each of them occurs only once in each row and column is said to be a Latin rectangle; L 1 and L 2 are two Latin rectangles consisting of a basis set of S . G is said to be an orthogonal Latin rectangle if there exists a unique ( i,j )m×n for each pair ( x 1 , x 2 )S×S such that L 1 ( i,j )= x 1 , L 2 ( i,j )= x 2 , then G is said to be an orthogonal Latin rectangle.

Let the family of subsets V ={ F 1 , F 2 ,, F s },s2 of [ V ]={ 0,1,2,,v } , V is said to be k -regular if i , | F i |=k , otherwise it is not. If  i,j[ s ] , | F i F j |t1 is satisfied, it is called t -intersecting. R. Gabrys et al. [19] called this subset family V ={ F 1 , F 2 ,, F s } , s2 ,and the sum of the number of positive and negative labels in each subset is exactly equal to { 1,0,+1 } by setting a label L on the elements of [ V ] , some of which are labelled as P + , and some of which are labelled as P , family of sets is the family of extremal balanced sets denoted as ( t,k,v ) , A( t,k,v ) denotes the number of subsets satisfying ( t,k,v ) .

Lemma 1 [20]: If t<k<v , where p + =| P + |, p =| P |, p + > k+1 2 ,v>2 , then

A( t,k,v )<S( t,k,v )

where S( t,k,v ) denotes the Steiner system: the set X has v points, the set of subsets of X of size k is called a block, and t points of X are in exactly one block. p + denotes the number of elements denoted as P + , p denotes the number of elements denoted as P

Lemma 2 [20]: Extremal balanced systems based on Latin rectangle constructions. If v8 , then

A( t,k,v )= v 2 v 2 2 . (4)

3. Main Results

3.1. Based on Orthogonal Latin Rectangle

Let v be an even, make a token for each element of V as follows P + ={ 0,1,, v 2 1 } , P ={ v 2 , v 2 +1,,v1 } , p + =| P + | , p =| P | then v=2 p + =2 p , r i,j denotes the element of Latin rectangle, i is the row indicator i P + and j is the column indicator j P . If the element of a Latin rectangle is equal to its row or column indicator, such a point is called a fixed point.

The Latin rectangle constructed in this paper cannot contain fixed points, if they contain fixed points that the elements in the corresponding set have the same contradict the definition of the set, two Latin rectangles without fixed points if they satisfy the definition of orthogonal Latin rectangle, the orthogonal Latin rectangle formed also do not contain fixed points. Latin rectangle without fixed points can form a Steiner system with parameters ( 2,3,v ) and also form an antichain. Their collection forms an antichain code, which results in a combinatorial neural code with special parameters.

The following shows the construction of an antichain based on orthogonal Latin rectangle with parameters t=2,k=4 .

If the first row of the Latin rectangle is all P + or P , then after the cycle will inevitably appear in an element of row will be equal to its row index or column index, so the first row needs to contain both P + and P , in order to facilitate the construction of the Latin rectangle is now divided into two sub-rectangle, represents the Latin rectangle, p + means that all the element of the positive sub-Latin rectangle are all from the P + , p means that all the negative sub-Latin rectangle are all from the P , and its rows and columns to do a substitution. It is easy to find that if p + and p do not have fixed points it means that does not have fixed points either, if r i,j =k p + then i P + , j P , mod p + for any element, and if r i,j =k p then i P , j P + , mod2 p + for any element.

The row index and column index of Latin moments and the corresponding elements in the above can form a ternary array. For example, r i,j =k can form a ternary array of { i,j,k } . It is not hard to see that knowing two of these three elements tells us whether r i,j is an element in a positive sub-Latin rectangle or a negative sub-Latin rectangle. For the convenience of representation, the row index i and the corresponding element k form a pair { i,k } , and the set formed by such a pair of positive sub-Latin rectangles is called the set of positive sub-Latin rectangles, and the set formed by such a pair of negative sub-Latin rectangles is called the set of negative sub-Latin rectangles.

In order for the set to which the Latin rectangle are pairs to satisfy the definition of set and t-intersectivity, the following three conditions are satisfied on top of the above:

Condition1 (Update rule): If r i,j =k then r i1,j+1 =k1,i P + ,j P .

Condition2 (Symmetry condition): If r i,j =h p + then r h,j =i p + ; if r i,j =h p then r h,j =i p (Note the row and column substitution).

Condition3 (Regularity condition): Different pairs { i,h } and { i , h } are belonging to positive sub-Latin rectangle negative sub-Latin rectangle then there is i i h h .

Note: Consider condition 1, e.g. r 2,0 =5 r 1,1 =4 and it can be found that the elements of the Latin rectangle are completely determined by the first column, so that a good choice of the first column constructs the Latin square; Condition 2, e.g. r 0,8 =7 and r 7,8 =0 form the same triad ( 0,7,8 ) . The symmetry condition groups the elements of the Latin rectangle into pairs that result in identical elements in which exactly half of them are identical.

It can be found that the Latin rectangles as long as the construction of the first column in accordance with the update rule can be constructed throughout the Latin rectangle, but the Latin rectangle also to meet the inside of all the elements are not duplicated, so the construction of the first column to meet the first column cannot be duplicated meta, but also note that the peer cannot be duplicated meta, and the condition of the 2 is to meet the requirements of this. For example: if the pairs { 6,2 } and { 4,0 } are part of the set of positive sub-Latin rectangle, r 6,0 =2 r 5,3 =1 r 4,4 =0 and r 4,0 =0 , the presence of two 0 in the same rows contradicts the definition of a Latin rectangle. Thus, the regularity condition ensures that there are no same elements in each row.

Construct 1: p + = p is a multiple of 4 if v=4m .

The positive sub-Latin rectangle set of the first Latin rectangle 1 has elements { i,i1 } , where i[ p + 4 1 ] ; the negative sub-Latin rectangle set has elements { i,i1 } , where i[ p + 4 , p + 4 +1,, p + 2 1 ] ;

The positive sub-Latin rectangle of the second Latin rectangle 2 has elements { i,i1 } , where i[ p + 4 , p + 4 +1,, p + 2 1 ] ; the negative sub-Latin rectangle has element { i,i1 } , where i[ p + 4 1 ] .

Orthogonal Latin rectangle ( x 1 , x 2 )v×v are composed and ( i,j ) p + × p , 1 ( i,j )= x 1 , 2 ( i,j )= x 2 .

Lemma 3: The Latin rectangle in Construction 1 has no fixed points, nor does the orthogonal Latin rectangle contain fixed points.

Proof: it is only necessary to prove that one Latin rectangle has no fixed point, and the second Latin rectangle is proved using the same steps.

1) There is no such thing as an element appearing twice in the first column, and the subsequent columns are sequentially minus one from the first, so there is also no such thing as an element existing twice, i.e., there are no identical elements in the same column.

2) Let different points { i,j },{ i , j } such that i i j j , by construction 1 j=i1 , j = i 1 thus i i =i1( i 1 )=( i i ) , i.e. i i =0 or i i = p + 2 , obviously i i =0 needs to be rounded off since i is not the same as i , so i i = p + 2 . From Construction 1 it follows that if i, i P + , max( i i )= p + 2 2,min( i i )=1 .

If i, i P ,

maxi i = p + 2,min( i i )= p + 2 2

Neither satisfies i i = p + 2 , so only one of i  and i   belongs to P and one to P + , which contradicts that both have to belong to positive sub-Latin rectangle or both have to belong to negative sub-Latin rectangle, i.e. there is no r i,j = r i , j , and so the peers don’t have the same element.

3) Let r i,j =i r i+1,j1 =i+1 r i+j,0 =i+j , contradictory to Construct 1. For r i,j =j the case is the same as for r i,j =i , i.e. it is proved that there is no fixed point.

The two Latin rectangles have no fixed points, i.e., neither do their mutually orthogonal Latin rectangle. #

Example 1: Let v=16, p + = p =8 , constructed Latin rectangles 1 and 2 (Figure 3 and Figure 4).

Figure 3. 2 .

Figure 4. 1 .

Constructed orthogonal Latin rectangles (Figure 5).

Figure 5. Orthogonal Latin rectangle.

Figure 6. Rounding off the orthogongal Latin rectangle.

Rounding off the Latin rectangle after repeating 3 elements (Figure 6).

The orthogonal Latin rectangle forms six different quaternions:

{0, 7, 8, 15}, {1, 5, 9, 14}, {0, 3, 10, 13}, {1, 6, 8, 13}, {1, 4, 9, 12}, {2, 5, 8, 11}.

These seven different quaternion feet do not contain an inclusion relation, i.e., they form a set that also forms an antichain:

{0, 7, 8, 15}, {1, 5, 9, 14}, {0, 3, 10, 13}, {1, 6, 8, 13}, {1, 4, 9, 12}, {2, 5, 8, 11},

Based on the relationship between antichain codes and combinatorial neural codes, this results in a combinatorial neural code with parameters

( n=16,K=6, δ r =3r+3 )

and a weight of 4, a character in a binary code

(10000001100000001) (0100010001000010) (1001000000100100) (0100001010000100) (0100100001001000) (0010010010010000).

Consider below the case where p + 2 is odd and v=4m .

Construct 2: v=4m , p + 2 is odd.

The elements of the set of positive sub-Latin rectangle of the first Latin rectangle 1 are { i,i1 } at i[ P + 2 4 1 ] , at i{ P + 2 4 , P + 2 4 +1,, P + 2 2 } for { i,i2 } ; the element of the set of negative sub-Latin rectangle is { P + 2 1, P + 2 4 1 } at i= P + 2 1 .

The elements of the set of positive sub-Latin rectangle of the second Latin matrix 2 are { P + 2 1, P + 2 4 1 } at i= P + 2 1 ; the elements of the set of negative sub-Latin rectangle are { i,i1 } at i[ P + 2 4 1 ] and at i{ P + 2 4 , P + 2 4 +1,, P + 2 2 } for { i,i2 } .

The orthogonal Latin rectangle ( x 1 , x 2 )v×v are composed and ( i,j ) p + × p , 1 ( i,j )= x 1 , 2 ( i,j )= x 2 .

Lemma 4: The Latin rectangle in Construction 2 has no fixed points, nor does the orthogonal Latin rectangle contain fixed points.

Proof: It is only necessary to prove that one Latin rectangle has no fixed point, and the second Latin rectangle is proved using the same steps.

1) There is no such thing as an element appearing twice in the first column, and the subsequent columns are in order minus one from the first, so there can be no such thing as an element appearing twice, i.e., there are no identical elements in the same column.

2) Let  { i,j } and { i , j } be two different elements belonging to the set of positive Latin rectangle, and let i i= j j .

a) When j=i1 , j = i 1 ,

i i =j j =i1( i 1 )=( i i ) ,

We get i i = p + 2 , by construction 2 if

i, i [ P + 2 4 1 ]max( i i )= p + 2 3,min( i i )=1

so i i { 1,2,, p + 2 3 } , a contradiction.

b) When j=i1 , j = i 2 ,

i i =j j =i1( i 2 )=( i i )+12( i i )=1 ,

Obviously not available, since P + is even, and in the case of mod P + , P + +1 being odd cannot equal 2( i i ) , a contradiction.

c) When j=i2 , j = i 1 ,

i i =j j =i2( i 1 )=( i i )1 ,

the same obviously cannot be obtained.

d) When j=i2 , j = i 2 ,

i i =j j =i1( i 1 )=( i i ) ,

We get i i = p + 2 , by construction 2 it follows that if

i, i { P + 2 4 , P + 2 4 +1,, P + 2 2 } , i i { p + 2 +3,, p + 2 } ,

a contradiction, that satisfies regularity for positive sub-Latin rectangle, discussed below for negative sub-Latin rectangle:

i i = P + 2 1( P + 2 4 1 )= 3 p + 2 4 j j = P + 2 4 1( p + 2 1 )= 23 p + 4

so i i j j , satisfying regularity.

3) Let r i,j =i r i+1,j1 =i+1 r i+j,0 =i+j , contradiction with Construct 2. For r i,j =j the case is the same as for r i,j =i , i.e. it is proved that there is no fixed point.

There is no fixed point for two Latin rectangle, i.e. there is no fixed point for their mutually orthogonal Latin rectangle. #

Example 2: Let v=20, p + 2 =5 , constructed Latin rectangles 1 and 2 (Figure 7 and Figure 8).

Figure 7. 1 .

Figure 8. 2 .

Constructed orthogonal Latin rectangles (Figure 9)

Figure 9. Orthogonal Latin rectangle.

Rounding off the Latin rectangle after repeating 3 elements. (Figure 10)

Figure 10. Rounding off the orthogonal Latin rectangle.

Orthogonal Latin rectangle forms an array of 12 distinct quaternions:

{0, 9, 10, 19}, {0, 7, 11, 18}, {0, 4, 12, 16}, {0, 3, 14, 17}, {1, 8, 10, 17}, {1, 5, 11, 15}, {2, 6, 10, 14}, {3, 5, 10, 12}, {2, 4, 11, 13}, {6, 7, 12, 13}, {0, 2, 16, 18}, {1, 9, 14, 16}.

These 12 different quaternions arrays do not have an inclusion relationship therefore the set formed constitutes an antichain:

{{0, 9, 10, 19}, {0, 7, 11, 18}, {0, 4, 12, 16}, {0, 3, 14, 17}, {1, 8, 10, 17}, {1, 5, 11, 15}, {2, 6, 10, 14}, {3, 5, 10, 12}, {2, 4, 11, 13}, {6, 7, 12, 13}, {0, 2, 16, 18}, {1, 9, 14, 16}}.

This results in a parameter of ( n=20,k=12, δ r =3r+3 ) and a combined neural code with a weight of 4. Its corresponding codeword is:

(1000000001100000001) (10000001000100000010) (100010000000100001000) (10010000000000100100) (01000000101000000100) (01000100000100010000) (00100010001000100000) (00010100001001000000) (00101000000101000000) (00000011000011000000) (10100000000000001010) (01000000010000101000).

Consider the case v=4m+1 below.

Construction 3: One column on top of construction 12 constitutes [ 2m ]×[ 2m+1 ] Latin rectangle, the first Latin rectangle adds a column r i,2m+1 = p + 2 +i p + , and the second Latin rectangle adds a column r i,2m+1 = p + 2 i p .

Lemma 5: Construction 3 is an orthogonal Latin rectangle without fixed points.

Proof: the first Latin rectangle obviously occurs at most once in the same column of elements, r 0,0 , r 0,1 ,, r 0,2m+1 the elements belonging to p + in this row are all odd, and the one element added is even. For row i , if i is even, i1 is odd, and p + 2 +i is even; if i is odd, i1 is even, and p + 2 +i is odd; therefore the same element cannot occur in the same row. The same proof for the second Latin rectangle, i.e., that there are no fixed points for either Latin rectangle, and no fixed points for the orthogonal Latin rectangle.

Consider below v=4m+2, p + = p =2m+1 .

Construct 4: v=4m+2, p + = p =2m+1 , construct orthogonal Latin rectangle of [ 2m+1 ]×[ 2m+1 ] .

The element of the set of positive sub-Latin rectangle of the first Latin moment 1 is { i,i1 } at i[ P + 1 2 1 ] the elements of the set of negative sub-Latin rectangle are { P + 1 2 , P + 1 2 1 } at i= P + 1 2

The elements of the set of positive sub-Latin rectangle of the second Latin rectangle 2 are { P + 1 2 , P + 1 2 1 } for i= P + 1 2 ; the elements of the set of negative sub-Latin rectangle are { i,i1 } for i[ P + 1 2 1 ] .

Lemma 6: Construction 4 is a Latin rectangle without fixed points, and orthogonal Latin rectangle does not contain fixed points.

Proof: it is sufficient to prove that one Latin rectangle has no fixed point, and the second Latin rectangle is proved using the same steps.

1) There is no such thing as an element appearing twice in the first column, and the subsequent columns are sequentially minus one from the first column, so there can be no such thing as an element appearing twice, i.e., there are no identical elements in the same column.

2) Let { i,j } and { i , j } be two different elements belonging to the set of positive Latin rectangle, and let i i= j j , By the construction 4 know

i i= j j= i 1( i1 )=( i i ) ,

that is, i i=0 or 2( i i )= p + obviously cannot be satisfied, a contradiction; and belongs to the negative sub-Latin rectangle of the set of a column of a row of only one element therefore the peer will not be duplicated, that is, to satisfy the regularity.

3) Let r i,j =i r i+1,j1 =i+1 r i+j,0 =i+j , Contradiction with Construct 4. For r i,j =j the case is the same as for r i,j =i , it is proved that there is no fixed point.

The two Latin rectangles have no fixed points, i.e., neither do their mutually orthogonal Latin rectangle.

Example 3: v=10, p + = p =5 , rounding off the Latin rectangle after repeating 3 elements. (Figure 11)

Figure 11. Rounding off the orthogonal Latin rectangle.

The orthogonal Latin rectangle forms 3 distinct quaternions:

{0, 4, 9, 5}, {0, 2, 6, 8}, {1, 3, 7, 5}.

These 3 different quaternions do not have an inclusion relationship therefore the set formed constitutes an antichain:

{{0, 4, 9, 5}, {0, 2, 6, 8}, {1, 3, 7, 5}}

This results in a parameter of ( n=10,k=3, δ r =3r+3 ) and the combined neural code with weight 4 has the corresponding codeword:

(1001100001) (1010001010) (0101010100)

Example 4: v=14, p + = p =7 , the orthogonal Latin rectangle forms 3 distinct quaternion arrays (Figure 12):

Figure 12. Rounding off the orthogonal Latin rectangle.

The orthogonal Latin rectangle forms 6 distinct quaternions:

{0, 7, 6, 13}, {0, 4, 8, 12}, {0, 2, 11, 9}, {1, 5, 7, 11}, {1, 3, 8, 10}, {2, 4, 7, 9}.

These six different quaternions do not have an inclusion relationship therefore the set formed constitutes an antichain:

{{0, 7, 6, 13}, {0, 4, 8, 12}, {0, 2, 11, 9}, {1, 5, 7, 11}, {1, 3, 8, 10}, {2, 4, 7, 9}}.

This results in a parameter of ( n=14,k=6, δ r =3r+3 ) and the combined neural code with weight 4 has the corresponding codeword:

(10100000010100) (01000101000100) (01010000101000) (00101001010000) (10000011000001) (10001000100010).

Consider the case v=4m+3 below.

Construct 5: One column on top of construct 4 constitutes the Latin rectangle of [ 2m+1 ]×[ 2m+2 ] ,

the first Latin rectangle adds a column r i,2m+2 = p + 1 2 +i p ,

the second Latin rectangle adds a column r i,2m+2 = p + 1 2 +i p + .

Lemma 7: Construction 5 is an orthogonal Latin rectangle without fixed points.

Proof: two proofs of similarity of Latin rectangle. The following only proves the first one. Obviously the same column element occurs at most once, r 0,0 , r 0,1 ,, r 0,2m+1 is only one element belonging to p in this row and it is odd, the added one is an element belonging to p and it is even. For row i , if i is even, i1 is odd and p + 1 2 +i is even; if i is odd, i1 is even and p + 1 2 +i is odd, so there can be no identical element in the row.

3.2. Construction Based on Disjoint Steiner Systems

F is a family of sets, If the set in F can be split into F 1 , F 2 ,, F n satisfying that ( V, F i ) is a ( t ,k,v ) system 1in , then ( V,F ) is t -partitionable. ( t,k,v ) In layman’s terms, this means that each set of size k in the family of sets formed by v elements has at most t1 elements that intersect every two sets.

Let ( V 1 , F 1 ) be t 1 -partitionable, with partition classes F 1 1 , F 2 1 ,, F n 1 , then ( V 2 , F 2 ) is t 2 -partitionable, with partition classes F 1 2 , F 2 2 ,, F n 2 , assuming that V 1 and V 2 do not intersect, then a new system can be formed, blocks are

{ AB:A F i 1 ,D F i 2 ,1in }

Lemma 8 [16] [19]: The above V 1 and V 2 form a new system with parameters

( t=max( k 2 + t 1 , k 1 + t 2 ), k 1 + k 2 , v 1 + v 2 ) . (5)

Proof: the new system contains v= v 1 + v 2 points, each block has k= k 1 + k 2 points, and the block

{ AB:A F i 1 ,B F i 2 ,1in } .

Let A 1 , A 2 F i 1 , B 1 , B 2 F i 2 then

t=max( ( A 1 B 1 )( A 2 B 2 ) ) .

And because the value of max( ( A 1 B 1 )( A 2 B 2 ) ) is only in two cases: the first one where B 1 is the same as B 2 ,

A 1 A 2 = t 1 1 ;

the second where A 1 is the same as A 2 ,

B 1 B 2 = t 2 1 , with | B i |= k 2 ,| A i |= k 1 ,

Since,

max( ( A 1 B 1 )( A 2 B 2 ) )=max( k 2 + t 1 , k 1 + t 2 )1 ,

So,

t=max( ( A 1 B 1 )( A 2 B 2 ) ) ,

That is, the parameters of the system constituting the new

( t=max( k 2 + t 1 , k 1 + t 2 ), k 1 + k 2 , v 1 + v 2 ) .

Construct 6: Based on 1-partitionable and 0-partitionable constructs ( 2,4,v ) .

The 1-factorization is ( 2,4,2m ) system 1-partitionable with 2m1 classes, the 1-factorization can be expressed as ( 1,2,2m ) and 2m1 points can form ( 1,1,2m1 ) with one point as a block of districts, which is 0-partitionable. By Lemma 6 a new system parameter can be formed as ( 2,3,4m1 ) . Obviously the number of block of ( 2,3,4m1 ) is m( 2m1 )=2 m 2 m , and 2 m 2 m points different from the previous ones are added to each block of ( 2,3,4m1 ) respectively to form a new system parameter as ( 2,4,2 m 2 +3m1 ) .

( 2,4,2 m 2 +3m1 ) forms 2 m 2 m quaternion, each of which has at most one intersecting element and no containment relation, so the set family formed by these quaternion is an antichain, i.e., it can be constituted as a 4-weighted combinatorial neural code with the parameters ( n=2 m 2 +3m1,k=2 m 2 m, δ r =3r+3 ) .

Theorem 1: Construct 6 forms a combinatorial neural code with parameter

( n=2 m 2 +3m1,k=2 m 2 m, δ r =3r+3 ) (6)

3.3. The Based Group Divisible Design Construct (2, 4, v), v Is a Multiple of 12

Group Divisible Design (GD design) [18]: with X as a given set of positive integers, K and M as a given set of positive integers, let D=( V,G,A ) be a finite associative structure, where V is a v-element set, G constitutes a division of V , the elements of A are called the block, and the elements of G are called the group, if the following conditions are satisfied:

  • For any BA , there is | B |K ;

  • For any GG , there is | G |M ;

  • For any BA and any GG , then | BG |1 ;

  • Any pair of elements of V belonging to different groups is contained in exactly λ block at the same time.

Then D is said to be a Group Divisible Design or GD design, denoted GD( K,X,M;V ) .

Theorem 2: GD( 4,1,3;v ) , v is a multiple of 12,

| B( 4,1;v+1 ) |= v 3 + v( v3 ) 12 (7)

Example 7: Let v={ 1,2,,12 } , the groups in GD( 4,1,3;v ) are {1, 2, 3}, {4, 5, 6}, {7, 8, 9}, {10, 11, 12}; blocks: {1, 4, 7, 10}, {1, 5, 8, 11}, {1, 6, 9, 12}, {2, 4, 8, 12}, {2, 5, 9, 10}, {2, 6, 7, 11}, {3, 4, 9, 11}, {3, 5, 7, 12}, {3, 6, 8, 10}. Adding a point 13 to each of these groups results in 4 new blocks: {1, 2, 3, 13}, {4, 5, 6, 13}, {7, 8, 9, 13}, {10, 11, 12, 13}. Together with the above, the 9 blocks form a B( 4,1;13 ) , | B( 4,1;13 ) |=13 . The above quaternion array forms a combinatorial neural code with parameters ( n=13,K=13, δ r =3r+3 ) with the following structure:

(1001001001000) (1000100100100) (1000010010010) (0101000100010) (0100100011000) (0100011000100) (0010101000010) (0011000010100) (1110000000001) (0010010101000) (0001110000001) (0000001110001) (0000000001111)

3.4. Based on the Transversal Design Construction (t, k, v), v = km, m, k Is Even

Transversal design [18]: i) The set V has kv elements called points; ii) Partition of V into sets { G i :i[ k ] } , each G i contains v points called groups; iii) A set B of k -subsets is called a block if it satisfies the following two conditions: a. Each group and each block intersect by exactly 1 point; b. t -subsets of V appear in only one block or two and more points are in one group but not all points.

Theorem 3: TD( t,k,v ) , v=km , m,k is even, then

| ( k1,k,v )= m k1 + 2 k 2 ×( m 2 )× ( m 2 ) k2 2 | (8)

Proof: v=km divide the v points into k groups, each group has m points, will be  v points with two-dimensional label [ k ]×[ m ] , based on the transversal design there are already m k1 blocks, because it is determined that one is selected for each group, and the last is determined by determining the k1 in front of it. Below add new block which is formed by

{ a 1 , a 2 ,, a k1 , a k } ,

where { a 1 , a 2 } Φ i g 1 , { a 3 , a 4 } Φ i g 2 , , { a k1 , a k } Φ i g k 2 , and

g 1 [ 0,1 ], g 2 [ 2,3 ],, g k 2 [ k2,k1 ], Φ ( g ) ={ Φ 1 ( g ) , Φ 2 ( g ) ,, Φ m1 ( g ) }

is a 1-factor decomposition of m elements in each group. There are a total of k groups for g[ k1 ] because g consists of { g 1 , g 2 ,, g k 2 } , and each g i has two choices, i.e., there are 2 k 2 choices for g . If { g 1 , g 2 ,, g k 2 } is fixed, { a 1 , a 2 } are related so there are ( m 2 ) choices, and { a 3 , a 4 } have m 2 choices, and similarly { a k1 , a k } also have m 2 choices, then { a 1 , a 2 ,, a k1 , a k } has a total of 2 k 2 ×( m 2 )× ( m 2 ) k2 2 choices. So

| ( k1,k,v ) |= m k1 + 2 k 2 ×( m 2 )× ( m 2 ) k2 2

TD( t,k,v ) constructs a family of sets of set size k and the sets in the family intersect at most t1 identically, which do not have an inclusion relation, thus constituting a combinatorial neural code. TD( k1,k,v ) constitutes a combinatorial neural code with parameters

( n=v,K= m k1 + 2 k 2 ×( m 2 )× ( m 2 ) k2 2 , δ r =2r+2 )

Since there are four choices for it follows that the total number of added blocks is 48.

Example 8: Suppose that is a transversal design whose with points of the form Z 4 ×{ 0,1,2,3 } , and blocks

{ { ( i 0 ,0 ),( i 1 ,1 ),( i 2 ,2 ),( i 3 ,3 ) }: i 0 , i 1 , i 2 , i 3 Z 4 , i 3 = i 0 + i 1 + i 2 }

Since there are four choices for each of i 0 , i 1 , i 2 it follows that we have 4 3 =64 blocks.

The factorization of interest reads as

Φ g ={ Φ 1 ( g ) , Φ 2 ( g ) , Φ 3 ( g ) } ={ { { ( 0,g ),( 1,g ) },{ ( 2,g ),( 3,g ) } }, { { ( 0,g ),( 2,g ) },{ ( 1,g ),( 3,g ) } }, { { ( 0,g ),( 3,g ) },{ ( 1,g ),( 2,g ) } } }

Where g{ 0,1,2,3 } , for a g 1 { 0,1 } and a g 2 { 2,3 } , we obtain 12 blocks:

{ { ( 0, g 1 ),( 1, g 1 ),( 0, g 2 ),( 1, g 2 ) } ,{ ( 0, g 1 ),( 1, g 1 ),( 2, g 2 ),( 3, g 2 ) }, { ( 2, g 1 ),( 3, g 1 ),( 0, g 2 ),( 1, g 2 ) },{ ( 2, g 1 ),( 3, g 1 ),( 2, g 2 ),( 3, g 2 ) }, { ( 0, g 1 ),( 2, g 1 ),( 0, g 2 ),( 2, g 2 ) },{ ( 0, g 1 ),( 2, g 1 ),( 1, g 2 ),( 3, g 2 ) }, { ( 1, g 1 ),( 3, g 1 ),( 0, g 2 ),( 2, g 2 ) },{ ( 1, g 1 ),( 3, g 1 ),( 1, g 2 ),( 3, g 2 ) }, { ( 0, g 1 ),( 3, g 1 ),( 0, g 2 ),( 3, g 2 ) },{ ( 0, g 1 ),( 3, g 1 ),( 1, g 2 ),( 2, g 2 ) }, { ( 1, g 1 ),( 2, g 1 ),( 0, g 2 ),( 3, g 2 ) }, { ( 1, g 1 ),( 2, g 1 ),( 1, g 2 ),( 2, g 2 ) } }

Since there are four choices for g 1 , g 2 it follows that the total number of added blocks is 48. There are 112 blocks in total, i.e., they constitute a combinatorial neural code with parameters ( n=16,K=112, δ r =2r+2 ) .

4. Conclusion

This paper focuses on constructing a class of combinatorial neural codes using orthogonal Latin rectangle, transversal designs, and group divisible designs. Section 1 introduces the research background and properties of combinatorial neural codes. Section 2 introduces antichain codes as a class of combinatorial neural codes, and the associated lemmas. Section 3 constructs combinatorial neural codes with parameters ( n,k, δ r =3r+3 ) using the method of orthogonal Latin rectangle and it has all 4-weight distributions; then a 4-weight code with parameters ( n=2 m 2 +3m1,K=2 m 2 m, δ r =3r+3 ) is constructed by using the disjointed Steiner system; and then a 4-weight code with parameters ( n,k= n 3 + n( n3 ) 12 , δ r =3r+3 ) 4-weighted code; and finally the combined neural code with parameters ( n,k, δ r ) is constructed based on the transversal design.

Acknowledgements

The authors gratefully acknowledge the discussion with Tang Chunming, Southwest Jiaotong University.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1] Atick, J.J. (1992) Could Information Theory Provide an Ecological Theory of Sensory Processing? Network: Computation in Neural Systems, 3, 4-44.
https://doi.org/10.3109/0954898X.2011.638888
[2] Borst, A. and Theunissen, F.E. (1999) Information Theory and Neural Coding. Nature Neuroscience, 2, 947-957.
https://doi.org/10.1038/14731
[3] O’Keefe, J. and Dostrovsky, J. (1971) The Hippocampus as a Spatial Map. Preliminary Evidence from Unit Activity in the Freely-Moving Rat. Brain Research, 34, 171-175.
https://doi.org/10.1016/0006-8993(71)90358-1
[4] Cotardo, G. and Ravagnani, A. (2022) Parameters of Codes for the Binary Asymmetric Channel. IEEE Transactions on Information Theory, 68, 2941-2950.
https://doi.org/10.1109/tit.2022.3147593
[5] Chen, A., Frick, F. and Shiu, A. (2019) Neural Codes, Decidability, and a New Local Obstruction to Convexity. SIAM Journal on Applied Algebra and Geometry, 3, 44-66.
https://doi.org/10.1137/18m1186563
[6] Gross, E., Obatake, N. and Youngs, N. (2018) Neural Ideals and Stimulus Space Visualization. Advances in Applied Mathematics, 95, 65-95.
https://doi.org/10.1016/j.aam.2017.10.002
[7] Curto, C., Itskov, V., Veliz-Cuba, A. and Youngs, N. (2013) The Neural Ring: An Algebraic Tool for Analyzing the Intrinsic Structure of Neural Codes. Bulletin of Mathematical Biology, 75, 1571-1611.
https://doi.org/10.1007/s11538-013-9860-3
[8] Curto, C., Itskov, V., Morrison, K., Roth, Z. and Walker, J.L. (2013) Combinatorial Neural Codes from a Mathematical Coding Theory Perspective. Neural Computation, 25, 1891-1925.
https://doi.org/10.1162/neco_a_00459
[9] Borst, A. and Theunissen, F.E. (1999) Information Theory and Neural Coding. Nature Neuroscience, 2, 947-957.
https://doi.org/10.1038/14731
[10] Curto, C., Itskov, V., Morrison, K., Roth, Z. and Walker, J.L. (2013) Combinatorial Neural Codes from a Mathematical Coding Theory Perspective. Neural Computation, 25, 1891-1925.
https://doi.org/10.1162/neco_a_00459
[11] Fu, F.-W., Ling, S. and Xing, C. (2003) New Lower Bounds and Constructions for Binary Codes Correcting Asymmetric Errors. IEEE Transactions on Information Theory, 49, 3294-3299.
https://doi.org/10.1109/tit.2003.820028
[12] Firer, M. and Walker, J.L. (2016) Matched Metrics and Channels. IEEE Transactions on Information Theory, 62, 1150-1156.
https://doi.org/10.1109/tit.2015.2512596
[13] Kløve, T. (1981) Errors Correcting Codes for Asymmetric Channel. Department of Informatics, University of Bergen.
[14] Qureshi, C.M. (2019) Matched Metrics to the Binary Asymmetric Channels. IEEE Transactions on Information Theory, 65, 1106-1112.
https://doi.org/10.1109/tit.2018.2885782
[15] Cotardo, G. and Ravagnani, A. (2020) Parameters of Combinatorial Neural Codes. arXiv: 2011.101v1.
[16] Gunby, B., He, X., Narayanan, B., et al. (2022) Antichain Codes. arXiv: 2212.08406.
[17] Sperner, E. (1928) Ein Satz über Untermengen einer endlichen Menge. Mathematische Zeitschrift, 27, 544-548.
https://doi.org/10.1007/bf01171114
[18] Zhang, X., Jing, X. and Feng, K. (2023) Optimal Combinatorial Neural Codes with Matched Metric δr: Characterization and Constructions. IEEE Transactions on Information Theory, 69, 5440-5448.
https://doi.org/10.1109/TIT.2023.3266010
[19] Shen, H. (2008) Combinatorial Design Theory. 2nd Edition, Shanghai Jiao Tong University Press.
[20] Gabrys, R., Dau, H., Colbourn, C.J. and Milenkovic, O. (2019) Set-Codes with Small Intersections and Small Discrepancies. 2019 IEEE International Symposium on Information Theory (ISIT), Paris, 7-12 July 2019, 2359-2363.
https://doi.org/10.1109/isit.2019.8849651.

Copyright © 2025 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.