SlideShare a Scribd company logo
A technique to construct linear trend free fractional factorial design using some linear codes
IJSM
A technique to construct linear trend free fractional
factorial design using some linear codes
Poonam Singh1
, Puja Thapliyal2*
, Veena Budhraja3
1
Department of Statistics,Faculty of Mathematical Sciences, University of Delhi,Delhi.110007
2*
Department of Statistics, Faculty of Mathematical Sciences, University of Delhi,Delhi.110007
3
Department of Statistics, Sri Venkateswara College, University of Delhi
Co-author: Email: veena_budhraja@yahoo.co.in, pbs_93@yahoo.co.in
*Corresponding author email: pujathapliyal98@gmail.com
Sometimes, an experimental situation may arise where experimental units are following
smooth pattern of trend over time and space. In that situation, to eliminate the effect of
uncontrollable variables that are correlated with the time, we use systematic run order
instead of randomising the run order. An ordering of treatments, thus obtained is known as
Trend free design. This article presents a method for constructing trend free fractional
factorial design using parity check matrix of a linear code. The method provides a
systematic approach to construct fractional and blocked fractional factorial design with
trend free main effects and some two- factor interactions.
Keywords: Time count , Linear Codes, Reed Muller Codes, BCH codes,Generalised Foldover Scheme
INTRODUCTION
Fractional Factorial experiments, especially those of exploratory type are often conducted sequentially in a specific run
order. In some situations, the experimentation scheme may suffer from some undesirable effects of factors. For
instance, in many industrial and agricultural experiments, treatments are applied to experimental units sequentially in
time or space, where there may be unknown or uncontrollable variables influencing the experimental process that are
highly correlated with the order in which the observations are obtained. Hence, the treatments which are not being
influenced by the effect of the trend under study is required. Such type of designs are known as trend free designs.
The need of trend free designs was mentioned by Joiner and Campbell (1976) who described an experiment in which
the measurement drift with time due to the build –up of carbon in a spectrophotometer and also by Freeny and Lai
(1997), in an experiment taken from the electronic industry in which a photolithographic polisher showed the tendency
drift lower through time .
The study of trend free designs was begun by Cox (1951) and has been addressed by different authors like Hill
(1960), Daniel and Wilcoxon(1966), Draper and Stoneman (1968), for the construction of trend free design. The
important work of Daniel and Wilcoxon (1966) was on constructing trend resistant contrasts for 2
n
design and using
the contrasts to generate a new run order. Cheng and Jacroux (1988) generalised Daniel and Wilcoxon work for 2
n
designs. Coster and Cheng (1988) introduced a generalized foldover method for constructing trend free run order
from a sequence of generators in a simple manner. Further, Bailey, Cheng and Kipnis (1992) extended the method
to general symmetric and asymmetrical factorial designs. Adekeys and Kuner(2006) compared the run order of 2
k-p
design in the presence of time trend. The two level trend robust fractional factorial designs were analysed by Mee
and Romanova(2010). In Recent years, Hilow(2012) developed an assignment procedures of constructing minimum
cost trend free two level fractional factorial designs.Some trend free designs in blocked fractional factorial designs
were being studied and investigated by Wu and Soushan(2013).
In this article, we present a method to construct the generator matrix for the required designs using the parity check
matrix of a linear code. Section 2 gives the preliminaries required. A brief introduction of coding theory is given in
section 3. Section 4 presents the construction technique to generate the designs. A method using some linear codes
is described and illustrated through examples in section 5. The section 6 presents the construction for trend free
blocked fractional factorial designs using the linear codes, mentioned in previous section.
International Journal of Statistics and Mathematics
Vol. 3(1), pp. 073-081, February, 2016. © www.premierpublishers.org, ISSN: 2375-0499x
Research Article
A technique to construct linear trend free fractional factorial design using some linear codes
Singh et al. 073
PRELIMINARIES
Fractional Factorial Design: Let D = (qr
n-p
) denote a q
-p
fraction of the q
n
factorial design in q
r
incomplete blocks. Let N
= q
n-p
be the number of treatment combinations, K=q
n-p-r
be the size of each block and B = q
r
be the number of
blocks.The linear model for the q
n
factorial design is given as
(1)
Where y is the N × 1 column vector of observations , X is the N× n design matrix of known constants; N is the n×1
column vector of regression coefficients and e is the N×1 column vector of random errors with zero means and
variance σ
2.
. When block effects are considered, the model is changed to
(2)
with B= is N×B block matrix with block size K and γ is the corresponding (K ×1) block .
Definition1: Let denotes the ordered vector of observations, for
x = 0,1,2,... be the N × 1 vector of trend coefficients and let ui be the contrast for main effect Ai ; i=1,2,...n, in the
run order. Then the quantity is known as the time count for the main effect Ai ..
A necessary and sufficient condition for a main effect contrast u to be trend free is that
(1.1)
In general, an N × 1 vector u is called trend free if equation(1.1) holds.
When the assumption of trend effect of different factors of the design is also considered, then this assumption changes
the model (2) to
(3)
Where , t is the (K ×1) linear vector and is the trend effect coefficient.The model (3) is a linear model
for factorial designs with trend. The first K rows of correspond to K treatment combinations in the
principal block, the next K rows to the treatment combinations in the second block and so on.
Definition 2: A run order is optimal for the estimation of the factor effects of interest in the presence of nuisance v-
degree polynomial trend iff
(1.2)
Where is an Nn matrix of factor effect coefficients and T is the Nv matrix of polynomial trend coefficients. If
equation(1.2) is satisfied then the run order is said to be v-trend free. If x is any column of nd t is any
column of T then the usual inner product is the time count between x and t. Criterion given in equation (1.2)
states that all the time count are zero for optimal run order.
Coster and Cheng(1988) introduced Generalised Foldover Scheme to generate systematic run order for full and
fractional factorial plans. We give here the definition of Generalised Foldover Scheme (Coster and Cheng(1988)).
Definition 3: Generalised Foldover Scheme (GFS): For a factorial design D= (qr
n-p
) with n factors A1, A2,…..An let G
be the generator matrix of design given by
A technique to construct linear trend free fractional factorial design using some linear codes
Int. J. Stat. Math. 074
Let be a 1×n matrix of zeros. Then the run order of design D is produced by the GFS with respect to the generator
sequence is given by where
for i = 1,2,…,n-p. In the run order , the principal block consists of the first q
n-p-r
runs; the second block consist of
the next q
n-p-r
runs, and so on. For a factorial design D1=(qr
n-p
), let be the n factors and G be the
generator matrix of design D1 given by Coster and Cheng (1988) derived following conditions for -trend free effects in
GFS. These conditions involve the generator matrix.
1. The main effect of a given factor is trend free if the corresponding letter appears at least ( +1) times in the
generator sequence.
2. A 2-factor interaction is -trend free if and only if there are at least ( +1) generator each of which exactly one of
the two factors appears at non-zero level.
For linear trend free design, the above conditions can be described in the following properties:
Property 1: For any factor , if there are at least two non zero elements , then all the main effect
components of factor are linear trend free.
Property 2:Forany factor ,if there are at least two pairs ,
such that one element is zero and the other element is non zero, then all components of the interaction are
linear trend free.
LINEAR CODES
A linear [n,k,d]q
code C over GF(q), where q is prime or prime power, n is the length, k is the dimension and d is the
minimum distance, is a k-dimensional subspace of the n-dimensional vector space V(n,q) over GF(q). The dual code
C
of an [n,k,d]q
code C is C
= { v V(n,q)/ v.w=0 for all w C}. This is an [n, n-k, d
]q
code and an (n-k)  n generator
matrix H of C
is called a parity check matrix of C. If the generator matrix is given in the standard form, a
corresponding parity check matrix is given as
H = [
T
A  knI  ]
Any d
-1 columns in generator matrix G of C are linearly independent and any d-1 columns in parity check matrix
H are linearly independent.
FACTORIAL DESIGNS WITH SOME LINEAR TREND FREE EFFECTS
GFS provides a technique to construct linear trend free factional factorial designs using generator matrix, but there is
no general method to construct the generator matrix. We use the parity check matrix of a linear [n,k,d]q
code to obtain
the set of generators to construct the desired designs.
Method of Construction
1. Consider the parity check matrix Hn-k×n of a linear [n,k,d]q code.
2. Partition the matrix H in two submatrix F1 and F2, where F2 consists of the vectors with weight (w) one.
3. Let m denotes the number of columns in F2. If m <d-1, then select any d-1 columns from H to form a matrix M of
order n-k × d-1 otherwise select M such that there is at least one column from F2. The order of M is n-k × d-1.
4. Write the transpose of M and the transpose be denoted by MT of order (d-1×n-k). Delete the column with all zero
entries and repeated columns to form the matrix Gα of order ,say (d-1)×l
5. Retain Gα , if the number of columns (l ) in Gα>d-1 Otherwise go to step 3.
6. Apply GFS on Gα to obtain the fractional factorialdesign q
l-r
; where 1≤ r ≤ (n-k)-(d-1) & (d-1)+1 ≤l ≤ n-k .
The resultant fractional factorial designis such that some of the main effects and two factor interactions are linear
trend free.
A technique to construct linear trend free fractional factorial design using some linear codes
Singh et al. 075
Remark
1. The total number of possible choices of M matrices is
2. In case d-1 = l in step 5, we get the full factorial design q
l
with some main effects and two factor
interactions Linear trend free.
3. The design of trend freeness of the factors depends upon the weight (w) of the corresponding
column in Gα.
4. All the main effects will be trend free if the weight(w) of all columns in Gα is at least two.
Above method of construction can be summarized in following theorem. The proof of the theorem can be
obtained from the sequence of steps given in the construction method. The proof(*) is available with the authors.
*
Theorem1: Existence of a linear [n, k, d]q code implies the existence of q
l-r
{1≤ r ≤ (n-k)-(d-1) & (d-1)+1≤l≤ n-k
},fractional factorial design and q
l
( l=d-1) full factorial design in which all main effects and some of the two-factor
interactions are linear trend free, where (d-1) is number of linearly independent columns in matrix H and l is the
number of factors/columns in Gα.
LINEAR TREND FREE FRACTIONAL FACTORIAL DESIGNS USING DIFFERENT CODES
We consider different types of linear codes and generate linear trend free fractional factorial designs using the method
of construction described in previous section. For details see Hedayat et.al (1999).
a) Reed Muller Codes
The r
th
order binary Reed-Muller code R(r,a) of length n = 2
a
, for 0  r  a, is the set of all vectors f, where f(i1,…..ia) is a
Boolean function which is a polynomial of degree at most r. For any a and any r, 0 ≤ r ≤ a, there is a binary r
th
order RM
code R(r,a) with the following properties:
Length n = 2
a
, dimension k = +…+ and minimum distance 2
a-r
. The parity check matrix of R(r,a) code is
the generator matrixof its dual code. The dual of RM(r,a) is RM(a-r-1,a) code.
Example 1: Consider a parity check matrix of RM(2,4)
H=
Any three columns in parity check matrix H matrix are linearly independent. Here d-1=3, and m =5 (in F2) and we get
the designs and corresponding generator matrices depending on the columns selected. According to the selection of
columns ,we have the following cases :
i)11 columns with none of the column having weight one and there are possible ways of selection where
none of the column is with weight (w) one and the corresponding factorial designs will have all main effects linear trend
free. All generated matrices and the table contains generator matrices, their corresponding designs, trend free effects
and defining relations are available with th authors.
A technique to construct linear trend free fractional factorial design using some linear codes
Int. J. Stat. Math. 076
ii) For m=5 and d-1=3, n-m=11.The possible choices of columns selected for generator matrix are
+ + . Suppose we choose the three columns 4
th
, 5
th
and 6
th
to form the matrix
MRM = On transposing we get matrix as GRM
in which all three rows d-1 =3 are independent. here d-1 =3 and l=5, Using the method of GFS on generator
matrix we get a resolution III, fractional factorial design with defining relation I =A1A3A5 = A1A2A3A4= A2A4A5.
Using the properties 1and 2 we observe that main effect component of A1, A3, A5 and the two factor interaction A1A2,
A1A3 and A2A4 and A3A4 are linear trend free.
Cyclic Codes
A linear code over GF(q) is said to be cyclic if whenever (c0
,c1
,…..,cn-2
,cn-1
) is a codeword so also is (c1
,c2
,…..,cn-1
,c0
).
Cyclic arrays can be described by a single generating vector z = ( z0
z1
……….zn-1
) such that the generator matrix consists
of this vector and its first ( k-1 ) cyclic shifts. The generating vector z is represented by a polynomial z(x) = z0 + z1X + …+zn-
1X
n-1
which is called a generator polynomial for the code. If a code is cyclic, so is its dual, and the generator
polynomial of its dual can be obtained by the following result given in Macwilliam and Sloane (1977).
Theorem 2: If C is a cyclic code of length n over GF(q), with generator polynomial z(x), then the dual code C
is also
cyclic and has generator polynomial
where z
*
(x) = X
deg.z
z(x
-1
) is reciprocal polynomial to z(x).
Example 2: Let = be the generator polynomial for
[15,7,5]2 code. Then the generator polynomial for the dual of this code is given as
h (x) = x
7
( 1+ x
-4
+x
-6
+x
-7
) = 1+x+x
3
+x
7
Hence the parity check matrix of [15,7,5]2 code is obtained by writing the coefficients and giving the cyclic shift to the
coefficients,as given below
Hcy= =
Any four columns in the H matrix are linearly independent, d-1 = 4 and m=5, n-m= 15-5 = 10 column with none of
them of weight(w) one. Then , there are 210 possible selection of such sets of columns from the matrix H. Following
the method of construction, we generated all generator matrices from the set of selected columns and listed
corresponding fractional factorial designs with their trend free effects and defining relations {available with the
authors}.The result of selecting columns 4
th
, 5
th
, 6
th
and 7
th
gives the matrix Mcy(8×4
Mcy(8×4=
A technique to construct linear trend free fractional factorial design using some linear codes
Singh et al. 077
Deleting column of zeroes and transposing
MT=
we get generator matrix as Gcy of order( (d-1)× l)
Gcy =
where l =7 and d-1 =4 . Using GFS, a 2
7-3
fractional factorial design is obtained with some main effects A3, A4, A5,
A6 and two factor interactions A1A2, A2A3, A2A6, A2A7, A3A4, A3A5, A3A6, A3A7, A5A4, A6A4, A5A6, A5A7 are linear trend
free.
BCH Codes
The BCH codes over GF(q) of length n = q
m
-1 and designed distance δ is the largest possible cyclic code having
zeroes α
b
,
,
α
b+1
,…,α
b+δ-2
where α є GF( q
m
) is the primitive n
th
root of unity, b is a non negative integer and m is the
multiplicative order of q mod n. The parity check matrix of a BCH code with b=1 is given by
C=


























1n2222
1n5255
1n3233
1n2
)(.........)()(1
......
......
......
)(.........)(1
)(.........)(1
.........1
where each entry is replaced by the corresponding binary m-tuple.
Example 3.: The parity check matrix of [15,5,7]2 code is
H=
and we observe here that d-1=6 columns in above matrix are linearly independent. Since the matrix cannot be
partitioned into the sub matrices ,we consider all 15 columns to get the design and they are
Suppose we select columns 8
th
,9
th
,10
th
,11
th
,12
th
and 13
th
from above matrix H.The matrix obtained
by selection is
M =
A technique to construct linear trend free fractional factorial design using some linear codes
Int. J. Stat. Math. 078
On transposing and deleting columns of zeroes we get matrix GBCH as
GBCH =
Applying GFS on generator matrix with l = 10 and d-1 = 6 we get the design 2
10-4
fractional factorial design with all
main effects linear trend free. Some of the designs that can be generated by selecting different columns of the
generator matrix GBCH of the code [15,5,7]2 along with their defining relations are available with the author.
Ternary Golay Code
The Golay code were discovered by M.J.E.Golay in late 1940’s. The (unextended) Golay code are examples of perfect
codes. A q-nary code that attains the hamming ( or sphere packing) bound i.e. the one which has
codewords, is said to be perfect code. Consider the ternary Golay Code [11,6,5]3
over a ternary alphabet, the relative distance of the codes is as large as it possibly can be for a ternary code, and it
satisfies Hamming bound and is therefore a perfect ternary Golay [11,6,5]3 code. For perfect code its dual distance is
same as its covering radius. In terms of design the strength is same as the estimation index of an orthogonal array
obtained using this code. Golay codes are unique in the sense that binary or ternary codes with same parameters, can
be shown to be equivalent to them.
Example 4: Consider the parity check matrix of [11,6,5]3 and partition it into two submatrices as given in the method in
section 4, we get the matrix H as
H=
Any four columns in matrix H5×11 are linearly independent d-1 = 4 and m = 5 . Thus,
15 sets of column in which none of the columns has weight (w) one. Whereas the total possible selection of columns
to formmatrix M is 325.
Suppose we choose columns 1
st
,2
nd
,3
rd
and 4
th
from matrix H, we get the matrix
M = .
Next , we generate the matrices of order (4×5) that ensures the linear trend freeness of main effects and following
the steps of construction method , we get 3
5-1
fractional factorial design.
and then 3
5-1
fractional factorial design I= A1A2
2
A3A4
2
is generated using GFS on matrix Ggolay. The generated
design has all its main effects linear trend free with respect to the properties [1-3].
A technique to construct linear trend free fractional factorial design using some linear codes
Singh et al. 079
BLOCKING IN FRACTIONAL FACTORIAL DESIGN WITH SOME LINEAR TREND FREE EFFECTS.
When the block size is smaller than the number of treatment combinations in any factorial experiment, the technique of
blocking is used to carry out the analysis. The factorial/fractional factorial experiment thus obtained is known as
Blocked fractional factorial experiment. When we go for blocking of fractional factorial design, the block structure
affects the linear trend-freeness of the effects. We state here the result in continuation with the properties [1-3] . The
first h=n-p-r generates the principle block and let zv , h+1 ≤ v ≤ n-p be the generators of other blocks. Then the
following holds:
Property 4: For any given v, h+1 ≤ v ≤ n-p and factor say A1 , suppose z1v ≠ 0. Then
(a) All (q-1) main effect components are linear trend-free.
(b) If ziv = 0 , 2 ≤ i ≤ n-p, then all components of A1 × Ai interaction are linear trend-free.
(c)
The method described is used to construct Blocked fractional factorial designs with some linear trend free effects.
Reed Muller Code: We use here the generator matrix constructed in Example 1 by selecting columns 8
th
,9
th
and 10
th
We get 2
4-1
fractional factorial design, here q=2, n=3, p=1, r1=1, h=n-p-r1=1 independent generators. Using GFS, h=2
generators forms the principle block and remaining generators form the contents of the other block in which main
effects A and C are linear trend free. The confounded effect of the design is ABC.
Table 1. 24-1
fractional factorial design with
resolution III, I=ABC
Cyclic Code: We consider here generator matrix constructed in example 2 and by selecting 1
st
, 2
nd
,5
th
and 9
th
columns
We get Gcy=
Using GFS 2
5-1
fractional factorial design is obtained n=5, p1=1, s=2, r1=1, h = n-p-r1 = 5-1-1 = 3 independent
generators, generates the principal block and other generates the other block contents. Table 1 gives the design
generated.
Table 2 : 25-1
blocked fractional factorial design with resolution IV, I = ABCE
BCH Code: Consider the generator matrix in Example 3 by selecting columns 1
st
, 4
th
, 5
th
, 6
th
, 7
th
and 8
th
as
GBCH =
2
9-3
fractional factorial design is obtained, here n=9, p=3 , r1=1 and h = n-p-r1 = 9-3-1= 5 independent generators
forms the principal block and remaining generators the other block. Thus, Blocked 2
9-3
fractional factorial.design with
defining relation I = DEGJ= ABDEH = BCDF is obtained. Table 3 displays the first 16 runs of the design
constituting the principle block and the next 16 the other block and so on.
Block 1: (1) , acd , abd, bc
Block 2: c , ac, ab, bcd
Block 1: (1) , ab, bc, ac, acd, bcd, abd, d
Block 2: ce, abce, be, ae, ade, bde, abcde, cde
A technique to construct linear trend free fractional factorial design using some linear codes
Int. J. Stat. Math. 080
Table 3. 29-3
blocked fractional factorial design with resolution IV, I= DEFJ = ABDEH=BCDF
Block 1: (1), aej, dfhj, adefh, abefgh, bfghj, abdegj, bdg, bcej, abc, cdefh, abcdfhj, acfghj, cdhj, acdeh, cf,
acefj, abcdefj, bcdfg, abcegh, bcghj, bdeh, bdhj, bef, abf, adfg, defgj, aghj, egh.
Block 2: abdgh, bdeghj, abfgj, befg, def, adfj, ehj, ah, acdeghj, bcfg, acefg, cfgj, bcdfj, abcdef, bch, abcehj, abcgj,
abcdfgh, bcdefghj, cefhj, acfh, bcd, acdj, aeg, gj, adefg. dfgh, bfh, abefhj, bdj, abde.
Ternary Golay Code: Consider the generator matrix from Example 4 by selecting columns 1
st
,2
nd
, 3
rd
and 4
th
of
matrix H as
Ggolay =
3
5-1
fractional factorial design is obtained. Here n= 5, p1=1 and q=3 , r1=1 ,h=n-p-r1=3 independent generators are
obtained. Thus, principal block is generated by first three generators and remaining two generates the content of the
other block. First 3
3
treatment combinations in 3
5-1
fractional factorial design (given in Table 4 ) forms the principle
block and next 3
3
combinations forms the contents of second block and the remaining the third block contents. The
factors are denoted by A,B,C,D and E.
Table 4: 35-1
blocked fractional factorial design with resolution IV, I= ABCDE
Block 1: (1), abcde, a
2
b
2
c
2
d
2
e
2
, abc
2
d
2
, a
2
b
2
e, cde
2
, a
2
b
2
cd, c
2
d
2
e, abe
2
, ab
2
ce
2
, a
2
c
2
d , bd
2
e, a
2
d
2
e
2
, bc,
ab
2
c
2
de, bc
2
de
2
, ab
2
d, a
2
ce, a
2
bc
2
e, b
2
de
2
, acd
2
, b
2
cd
2
e, ac
2
d
2
, b
2
ce, ade, a
2
bcd
2
e
2
, b
2
d
2
.
Block 2: a
2
bde
2
, b
2
cd
2
, ac
2
e, b
2
c
2
e
2
, ad, a
2
bcd
2
e, acd
2
e
2
, a
2
bc
2
, b
2
de, cde, abc
2
d
2
e
2
, a
2
b
2
, abe, a
2
b
2
cde
2
, c
2
d
2
,
a
2
b
2
c
2
d
2
e, e
2
, abcd, ab
2
c
2
d, a
2
d
2
e, bce
2
, a
2
c, ab
2
cd, ab
2
d
2
e
2
, bd
2
, ab
2
de, a
2
c
2
de
2
.
Block 3: ab
2
d
2
e, a
2
ce
2
, bc
2
d, a
2
c
2
de, bd
2
e
2
, ab
2
c, bce, ab
2
c
2
de
2
, a
2
d
2
, a
2
bcd
2
, b
2
d
2
e, ade, b
2
d, acd
2
e, a
2
bc
2
e
2
, ac
2
,
a
2
bde, b
2
cd
2
e
2
, c
2
d
2
e
2
, ab, a
2
b
2
cde, abcde
2
, a
2
b
2
c
2
d
2
, e, a
2
b
2
e
2
, cd, abc
2
d
2
e.
CONCLUSION
A technique to construct fractional factorial designs with factors at q level, q>2, with some linear trend free effects
using parity check/generator matrix of linear code is developed. Since the generator matrix is not unique in nature,
therefore fractional factorial designs with linear trend free effects with same/ different resolution can be easily
constructed.
REFERENCES
Adekeys K, Kunert J (2006). On the Comparison of Run Orders of Unreplicated 2
k-p
designs in the presence of a
Time Trend. Metrika. 63: 257-269.
Bailey RA, Cheng CS, Kipnis P (1992). Construction of trend resistant factorial designs. Statist. Sinica. 2:393-
411.
Betsumiya K, Harada M (2001). Classification of Formally Self-Dual Even Codes of length to 16. Designs, Codes
and Cryptography. 23: 325-332.
Cheng CS, Jacroux M. (1988).On the construction of trend free run-order of two –level factorial designs. Journal of
American statistical association. 83:1152-1158.
Cox DR (1951). Some systematic experimental designs. Biometrica. 38: 312- 323.
Coster DC, Cheng CS (1998). Minimum cost trend free run orders of fractional factorial design. The Annals of
Statistics. 16:1188-1205.
Daniel C, Wilcoxon F (1966). Factorial 2
p-q
plans robust against linear and quadratic trends.Technometrics. 8: 259-
278.
Draper NR, Stoneman DM (1968). Factor changes and linear trends in eight run two level factorial designs.
Technometrics. 10: 301-311.
A technique to construct linear trend free fractional factorial design using some linear codes
Singh et al. 081
Freeny AE, Lai WYC (1997). Planarization by chemical Mechanical polishing: A rateand uniformity study. Statistical
case studies for industrial Process Improvement. V. Czitrom and P.D. Spagon eds Philadelphia, Alexandria: ASA-
SIAM, pp. 251-264.
Hedayat, AS, Sloane NJA, Stuffken J (1999). Orthogonal Arrays: Theory and applications. New York. US: Springer
and Verlag . 61-67, 87-108, 247-288.
Hill HH (1960).Experimental Designs to adjust for time trends. Technometrics. 2: 67-82.
Hilow H (2012). Minimum cost linear trend free fractional factorial designs. Journal of statistical theory and
practice. 6: 580-589.
Joiner BL, Campbell C (1976). Designing experiments when run order is important. Technometrics. 18: 249-259.
Mac Williams FJ, Sloane NJA (1977). The Theory of Error-Correcting Codes. North- Holland. Amesterdam: Oxford.
Mee R, Romanora A (2010). Constructing and Analyzing Two-level Trend Robust designs. Quality Enigeneering.
22: 306-316.
Wang PC, Wu S (2013). Trend Free Designs in Blocked Fractional Factorial Experiments. Communications in
Statistics- Theory and Methods. 42(10):1870-1878.
Accepted December 31, 2015.
Citation: Singh P, Thapliyal P, Budhraja V (2016). A Technique to Construct Linear Trend Free Fractional Factorial
Design Using Some Linear Codes. International Journal of Statistics and Mathematics 3(1): 073-081.
Copyright: © 2016 Singh et al. This is an open-access article distributed under the terms of the Creative Commons
Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the
original author and source are cited.

More Related Content

PPT
DESIGN AND ANALYSIS OF ALGORITHMS
PPTX
Daa unit 3
PPTX
Divide and Conquer - Part II - Quickselect and Closest Pair of Points
RTF
Design and Analysis of algorithms
PPTX
Data Structures- Hashing
PDF
Ch24 efficient algorithms
DESIGN AND ANALYSIS OF ALGORITHMS
Daa unit 3
Divide and Conquer - Part II - Quickselect and Closest Pair of Points
Design and Analysis of algorithms
Data Structures- Hashing
Ch24 efficient algorithms

What's hot (20)

PDF
Algorithm chapter 2
PPT
5.5 back tracking 02
PDF
Design and analysis of algorithm
DOCX
CS 6402 – DESIGN AND ANALYSIS OF ALGORITHMS questio
PDF
Backtracking & branch and bound
PPTX
Spline interpolation numerical methods presentation
PDF
01. design & analysis of agorithm intro & complexity analysis
PPT
Divide and conquer
PPT
Design and Analysis of Algorithms
PPTX
Lecture2b algorithm
PPT
Spline Interpolation
PPTX
Asymptotic Notation and Data Structures
PDF
Jmestn42351212
PDF
Analysis and design of algorithms part 4
PPTX
Mathematical Analysis of Recursive Algorithm.
PPTX
PPT
CS8461 - Design and Analysis of Algorithms
PDF
Digit Factorial Chains .(Euler Problem -74) (Matlab Programming Solution)
PPT
CS8451 - Design and Analysis of Algorithms
PPT
Analysis of Algorithm
Algorithm chapter 2
5.5 back tracking 02
Design and analysis of algorithm
CS 6402 – DESIGN AND ANALYSIS OF ALGORITHMS questio
Backtracking & branch and bound
Spline interpolation numerical methods presentation
01. design & analysis of agorithm intro & complexity analysis
Divide and conquer
Design and Analysis of Algorithms
Lecture2b algorithm
Spline Interpolation
Asymptotic Notation and Data Structures
Jmestn42351212
Analysis and design of algorithms part 4
Mathematical Analysis of Recursive Algorithm.
CS8461 - Design and Analysis of Algorithms
Digit Factorial Chains .(Euler Problem -74) (Matlab Programming Solution)
CS8451 - Design and Analysis of Algorithms
Analysis of Algorithm
Ad

Viewers also liked (15)

PPTX
Toma de decisiones
PPTX
Prezentacja1
PDF
Nuage Wordcreate 5
PDF
La amistad 2
DOCX
Trabajo practico n 2
DOCX
Definiciones
DOCX
Strategy Inventory for Language Learning (SILL) survey
PPTX
Vibration Controller | Electrodynamic Vibration System
PPTX
Presentation1
PDF
Lab 1: Modifying a module's attributes and viewing its history
PDF
Portforio works 02_haseko
PDF
Closet Works Beautifully Designed Spaces
PDF
Real Data For Real Change - SQL Power Brochure
PPS
Gedachtegang van een_man
Toma de decisiones
Prezentacja1
Nuage Wordcreate 5
La amistad 2
Trabajo practico n 2
Definiciones
Strategy Inventory for Language Learning (SILL) survey
Vibration Controller | Electrodynamic Vibration System
Presentation1
Lab 1: Modifying a module's attributes and viewing its history
Portforio works 02_haseko
Closet Works Beautifully Designed Spaces
Real Data For Real Change - SQL Power Brochure
Gedachtegang van een_man
Ad

Similar to A technique to construct linear trend free fractional factorial design using some linear codes (20)

PPTX
Design of Engineering Experiments Part 5
PDF
Comparison of the optimal design
PPTX
Optimization technique in pharmaceutical product development.pptx
PPT
factorial design
PPT
Fractional factorial design tutorial
PPT
Design of experiments
PPTX
factorial design.pptx
PDF
RM_05_DOE.pdf
PPTX
516factorialdesign-161119134528-converted.pptx
PPTX
Grds international conference on pure and applied science (9)
PPTX
Factorial design ,full factorial design, fractional factorial design
PPTX
Factorial Design (M. Pharmacy- 1st Semester)
PDF
23 industrial engineering
PPTX
Response Surface.pptx
PDF
Javier Garcia - Verdugo Sanchez - Six Sigma Training - W3 Fractional Factoria...
PDF
Information Matrices and Optimality Values for various Block Designs
PPTX
Drawing Cumulative Frequency Histogram Polygon
PPTX
Chetan dhal-Optimization techniques in pharmaceutics, formulation and processing
PPTX
Factorial design \Optimization Techniques
Design of Engineering Experiments Part 5
Comparison of the optimal design
Optimization technique in pharmaceutical product development.pptx
factorial design
Fractional factorial design tutorial
Design of experiments
factorial design.pptx
RM_05_DOE.pdf
516factorialdesign-161119134528-converted.pptx
Grds international conference on pure and applied science (9)
Factorial design ,full factorial design, fractional factorial design
Factorial Design (M. Pharmacy- 1st Semester)
23 industrial engineering
Response Surface.pptx
Javier Garcia - Verdugo Sanchez - Six Sigma Training - W3 Fractional Factoria...
Information Matrices and Optimality Values for various Block Designs
Drawing Cumulative Frequency Histogram Polygon
Chetan dhal-Optimization techniques in pharmaceutics, formulation and processing
Factorial design \Optimization Techniques

More from Premier Publishers (20)

PDF
Evaluation of Agro-morphological Performances of Hybrid Varieties of Chili Pe...
PDF
An Empirical Approach for the Variation in Capital Market Price Changes
PDF
Influence of Nitrogen and Spacing on Growth and Yield of Chia (Salvia hispani...
PDF
Enhancing Social Capital During the Pandemic: A Case of the Rural Women in Bu...
PDF
Impact of Provision of Litigation Supports through Forensic Investigations on...
PDF
Improving the Efficiency of Ratio Estimators by Calibration Weightings
PDF
Urban Liveability in the Context of Sustainable Development: A Perspective fr...
PDF
Transcript Level of Genes Involved in “Rebaudioside A” Biosynthesis Pathway u...
PDF
Multivariate Analysis of Tea (Camellia sinensis (L.) O. Kuntze) Clones on Mor...
PDF
Causes, Consequences and Remedies of Juvenile Delinquency in the Context of S...
PDF
The Knowledge of and Attitude to and Beliefs about Causes and Treatments of M...
PDF
Effect of Phosphorus and Zinc on the Growth, Nodulation and Yield of Soybean ...
PDF
Influence of Harvest Stage on Yield and Yield Components of Orange Fleshed Sw...
PDF
Performance evaluation of upland rice (Oryza sativa L.) and variability study...
PDF
Response of Hot Pepper (Capsicum Annuum L.) to Deficit Irrigation in Bennatse...
PDF
Harnessing the Power of Agricultural Waste: A Study of Sabo Market, Ikorodu, ...
PDF
Influence of Conferences and Job Rotation on Job Productivity of Library Staf...
PDF
Scanning Electron Microscopic Structure and Composition of Urinary Calculi of...
PDF
Gentrification and its Effects on Minority Communities – A Comparative Case S...
PDF
Oil and Fatty Acid Composition Analysis of Ethiopian Mustard (Brasicacarinata...
Evaluation of Agro-morphological Performances of Hybrid Varieties of Chili Pe...
An Empirical Approach for the Variation in Capital Market Price Changes
Influence of Nitrogen and Spacing on Growth and Yield of Chia (Salvia hispani...
Enhancing Social Capital During the Pandemic: A Case of the Rural Women in Bu...
Impact of Provision of Litigation Supports through Forensic Investigations on...
Improving the Efficiency of Ratio Estimators by Calibration Weightings
Urban Liveability in the Context of Sustainable Development: A Perspective fr...
Transcript Level of Genes Involved in “Rebaudioside A” Biosynthesis Pathway u...
Multivariate Analysis of Tea (Camellia sinensis (L.) O. Kuntze) Clones on Mor...
Causes, Consequences and Remedies of Juvenile Delinquency in the Context of S...
The Knowledge of and Attitude to and Beliefs about Causes and Treatments of M...
Effect of Phosphorus and Zinc on the Growth, Nodulation and Yield of Soybean ...
Influence of Harvest Stage on Yield and Yield Components of Orange Fleshed Sw...
Performance evaluation of upland rice (Oryza sativa L.) and variability study...
Response of Hot Pepper (Capsicum Annuum L.) to Deficit Irrigation in Bennatse...
Harnessing the Power of Agricultural Waste: A Study of Sabo Market, Ikorodu, ...
Influence of Conferences and Job Rotation on Job Productivity of Library Staf...
Scanning Electron Microscopic Structure and Composition of Urinary Calculi of...
Gentrification and its Effects on Minority Communities – A Comparative Case S...
Oil and Fatty Acid Composition Analysis of Ethiopian Mustard (Brasicacarinata...

Recently uploaded (20)

PPTX
1st Inaugural Professorial Lecture held on 19th February 2020 (Governance and...
PDF
Trump Administration's workforce development strategy
PDF
GENETICS IN BIOLOGY IN SECONDARY LEVEL FORM 3
PPTX
Final Presentation General Medicine 03-08-2024.pptx
PDF
1_English_Language_Set_2.pdf probationary
PPTX
Chinmaya Tiranga Azadi Quiz (Class 7-8 )
PPTX
Orientation - ARALprogram of Deped to the Parents.pptx
PDF
LNK 2025 (2).pdf MWEHEHEHEHEHEHEHEHEHEHE
PDF
A GUIDE TO GENETICS FOR UNDERGRADUATE MEDICAL STUDENTS
PDF
RTP_AR_KS1_Tutor's Guide_English [FOR REPRODUCTION].pdf
PPTX
CHAPTER IV. MAN AND BIOSPHERE AND ITS TOTALITY.pptx
PDF
احياء السادس العلمي - الفصل الثالث (التكاثر) منهج متميزين/كلية بغداد/موهوبين
PPTX
Onco Emergencies - Spinal cord compression Superior vena cava syndrome Febr...
PDF
A systematic review of self-coping strategies used by university students to ...
PPTX
Digestion and Absorption of Carbohydrates, Proteina and Fats
PDF
OBE - B.A.(HON'S) IN INTERIOR ARCHITECTURE -Ar.MOHIUDDIN.pdf
PDF
IGGE1 Understanding the Self1234567891011
PDF
Chinmaya Tiranga quiz Grand Finale.pdf
PPTX
Unit 4 Skeletal System.ppt.pptxopresentatiom
PPTX
202450812 BayCHI UCSC-SV 20250812 v17.pptx
1st Inaugural Professorial Lecture held on 19th February 2020 (Governance and...
Trump Administration's workforce development strategy
GENETICS IN BIOLOGY IN SECONDARY LEVEL FORM 3
Final Presentation General Medicine 03-08-2024.pptx
1_English_Language_Set_2.pdf probationary
Chinmaya Tiranga Azadi Quiz (Class 7-8 )
Orientation - ARALprogram of Deped to the Parents.pptx
LNK 2025 (2).pdf MWEHEHEHEHEHEHEHEHEHEHE
A GUIDE TO GENETICS FOR UNDERGRADUATE MEDICAL STUDENTS
RTP_AR_KS1_Tutor's Guide_English [FOR REPRODUCTION].pdf
CHAPTER IV. MAN AND BIOSPHERE AND ITS TOTALITY.pptx
احياء السادس العلمي - الفصل الثالث (التكاثر) منهج متميزين/كلية بغداد/موهوبين
Onco Emergencies - Spinal cord compression Superior vena cava syndrome Febr...
A systematic review of self-coping strategies used by university students to ...
Digestion and Absorption of Carbohydrates, Proteina and Fats
OBE - B.A.(HON'S) IN INTERIOR ARCHITECTURE -Ar.MOHIUDDIN.pdf
IGGE1 Understanding the Self1234567891011
Chinmaya Tiranga quiz Grand Finale.pdf
Unit 4 Skeletal System.ppt.pptxopresentatiom
202450812 BayCHI UCSC-SV 20250812 v17.pptx

A technique to construct linear trend free fractional factorial design using some linear codes

  • 1. A technique to construct linear trend free fractional factorial design using some linear codes IJSM A technique to construct linear trend free fractional factorial design using some linear codes Poonam Singh1 , Puja Thapliyal2* , Veena Budhraja3 1 Department of Statistics,Faculty of Mathematical Sciences, University of Delhi,Delhi.110007 2* Department of Statistics, Faculty of Mathematical Sciences, University of Delhi,Delhi.110007 3 Department of Statistics, Sri Venkateswara College, University of Delhi Co-author: Email: veena_budhraja@yahoo.co.in, pbs_93@yahoo.co.in *Corresponding author email: pujathapliyal98@gmail.com Sometimes, an experimental situation may arise where experimental units are following smooth pattern of trend over time and space. In that situation, to eliminate the effect of uncontrollable variables that are correlated with the time, we use systematic run order instead of randomising the run order. An ordering of treatments, thus obtained is known as Trend free design. This article presents a method for constructing trend free fractional factorial design using parity check matrix of a linear code. The method provides a systematic approach to construct fractional and blocked fractional factorial design with trend free main effects and some two- factor interactions. Keywords: Time count , Linear Codes, Reed Muller Codes, BCH codes,Generalised Foldover Scheme INTRODUCTION Fractional Factorial experiments, especially those of exploratory type are often conducted sequentially in a specific run order. In some situations, the experimentation scheme may suffer from some undesirable effects of factors. For instance, in many industrial and agricultural experiments, treatments are applied to experimental units sequentially in time or space, where there may be unknown or uncontrollable variables influencing the experimental process that are highly correlated with the order in which the observations are obtained. Hence, the treatments which are not being influenced by the effect of the trend under study is required. Such type of designs are known as trend free designs. The need of trend free designs was mentioned by Joiner and Campbell (1976) who described an experiment in which the measurement drift with time due to the build –up of carbon in a spectrophotometer and also by Freeny and Lai (1997), in an experiment taken from the electronic industry in which a photolithographic polisher showed the tendency drift lower through time . The study of trend free designs was begun by Cox (1951) and has been addressed by different authors like Hill (1960), Daniel and Wilcoxon(1966), Draper and Stoneman (1968), for the construction of trend free design. The important work of Daniel and Wilcoxon (1966) was on constructing trend resistant contrasts for 2 n design and using the contrasts to generate a new run order. Cheng and Jacroux (1988) generalised Daniel and Wilcoxon work for 2 n designs. Coster and Cheng (1988) introduced a generalized foldover method for constructing trend free run order from a sequence of generators in a simple manner. Further, Bailey, Cheng and Kipnis (1992) extended the method to general symmetric and asymmetrical factorial designs. Adekeys and Kuner(2006) compared the run order of 2 k-p design in the presence of time trend. The two level trend robust fractional factorial designs were analysed by Mee and Romanova(2010). In Recent years, Hilow(2012) developed an assignment procedures of constructing minimum cost trend free two level fractional factorial designs.Some trend free designs in blocked fractional factorial designs were being studied and investigated by Wu and Soushan(2013). In this article, we present a method to construct the generator matrix for the required designs using the parity check matrix of a linear code. Section 2 gives the preliminaries required. A brief introduction of coding theory is given in section 3. Section 4 presents the construction technique to generate the designs. A method using some linear codes is described and illustrated through examples in section 5. The section 6 presents the construction for trend free blocked fractional factorial designs using the linear codes, mentioned in previous section. International Journal of Statistics and Mathematics Vol. 3(1), pp. 073-081, February, 2016. © www.premierpublishers.org, ISSN: 2375-0499x Research Article
  • 2. A technique to construct linear trend free fractional factorial design using some linear codes Singh et al. 073 PRELIMINARIES Fractional Factorial Design: Let D = (qr n-p ) denote a q -p fraction of the q n factorial design in q r incomplete blocks. Let N = q n-p be the number of treatment combinations, K=q n-p-r be the size of each block and B = q r be the number of blocks.The linear model for the q n factorial design is given as (1) Where y is the N × 1 column vector of observations , X is the N× n design matrix of known constants; N is the n×1 column vector of regression coefficients and e is the N×1 column vector of random errors with zero means and variance σ 2. . When block effects are considered, the model is changed to (2) with B= is N×B block matrix with block size K and γ is the corresponding (K ×1) block . Definition1: Let denotes the ordered vector of observations, for x = 0,1,2,... be the N × 1 vector of trend coefficients and let ui be the contrast for main effect Ai ; i=1,2,...n, in the run order. Then the quantity is known as the time count for the main effect Ai .. A necessary and sufficient condition for a main effect contrast u to be trend free is that (1.1) In general, an N × 1 vector u is called trend free if equation(1.1) holds. When the assumption of trend effect of different factors of the design is also considered, then this assumption changes the model (2) to (3) Where , t is the (K ×1) linear vector and is the trend effect coefficient.The model (3) is a linear model for factorial designs with trend. The first K rows of correspond to K treatment combinations in the principal block, the next K rows to the treatment combinations in the second block and so on. Definition 2: A run order is optimal for the estimation of the factor effects of interest in the presence of nuisance v- degree polynomial trend iff (1.2) Where is an Nn matrix of factor effect coefficients and T is the Nv matrix of polynomial trend coefficients. If equation(1.2) is satisfied then the run order is said to be v-trend free. If x is any column of nd t is any column of T then the usual inner product is the time count between x and t. Criterion given in equation (1.2) states that all the time count are zero for optimal run order. Coster and Cheng(1988) introduced Generalised Foldover Scheme to generate systematic run order for full and fractional factorial plans. We give here the definition of Generalised Foldover Scheme (Coster and Cheng(1988)). Definition 3: Generalised Foldover Scheme (GFS): For a factorial design D= (qr n-p ) with n factors A1, A2,…..An let G be the generator matrix of design given by
  • 3. A technique to construct linear trend free fractional factorial design using some linear codes Int. J. Stat. Math. 074 Let be a 1×n matrix of zeros. Then the run order of design D is produced by the GFS with respect to the generator sequence is given by where for i = 1,2,…,n-p. In the run order , the principal block consists of the first q n-p-r runs; the second block consist of the next q n-p-r runs, and so on. For a factorial design D1=(qr n-p ), let be the n factors and G be the generator matrix of design D1 given by Coster and Cheng (1988) derived following conditions for -trend free effects in GFS. These conditions involve the generator matrix. 1. The main effect of a given factor is trend free if the corresponding letter appears at least ( +1) times in the generator sequence. 2. A 2-factor interaction is -trend free if and only if there are at least ( +1) generator each of which exactly one of the two factors appears at non-zero level. For linear trend free design, the above conditions can be described in the following properties: Property 1: For any factor , if there are at least two non zero elements , then all the main effect components of factor are linear trend free. Property 2:Forany factor ,if there are at least two pairs , such that one element is zero and the other element is non zero, then all components of the interaction are linear trend free. LINEAR CODES A linear [n,k,d]q code C over GF(q), where q is prime or prime power, n is the length, k is the dimension and d is the minimum distance, is a k-dimensional subspace of the n-dimensional vector space V(n,q) over GF(q). The dual code C of an [n,k,d]q code C is C = { v V(n,q)/ v.w=0 for all w C}. This is an [n, n-k, d ]q code and an (n-k)  n generator matrix H of C is called a parity check matrix of C. If the generator matrix is given in the standard form, a corresponding parity check matrix is given as H = [ T A  knI  ] Any d -1 columns in generator matrix G of C are linearly independent and any d-1 columns in parity check matrix H are linearly independent. FACTORIAL DESIGNS WITH SOME LINEAR TREND FREE EFFECTS GFS provides a technique to construct linear trend free factional factorial designs using generator matrix, but there is no general method to construct the generator matrix. We use the parity check matrix of a linear [n,k,d]q code to obtain the set of generators to construct the desired designs. Method of Construction 1. Consider the parity check matrix Hn-k×n of a linear [n,k,d]q code. 2. Partition the matrix H in two submatrix F1 and F2, where F2 consists of the vectors with weight (w) one. 3. Let m denotes the number of columns in F2. If m <d-1, then select any d-1 columns from H to form a matrix M of order n-k × d-1 otherwise select M such that there is at least one column from F2. The order of M is n-k × d-1. 4. Write the transpose of M and the transpose be denoted by MT of order (d-1×n-k). Delete the column with all zero entries and repeated columns to form the matrix Gα of order ,say (d-1)×l 5. Retain Gα , if the number of columns (l ) in Gα>d-1 Otherwise go to step 3. 6. Apply GFS on Gα to obtain the fractional factorialdesign q l-r ; where 1≤ r ≤ (n-k)-(d-1) & (d-1)+1 ≤l ≤ n-k . The resultant fractional factorial designis such that some of the main effects and two factor interactions are linear trend free.
  • 4. A technique to construct linear trend free fractional factorial design using some linear codes Singh et al. 075 Remark 1. The total number of possible choices of M matrices is 2. In case d-1 = l in step 5, we get the full factorial design q l with some main effects and two factor interactions Linear trend free. 3. The design of trend freeness of the factors depends upon the weight (w) of the corresponding column in Gα. 4. All the main effects will be trend free if the weight(w) of all columns in Gα is at least two. Above method of construction can be summarized in following theorem. The proof of the theorem can be obtained from the sequence of steps given in the construction method. The proof(*) is available with the authors. * Theorem1: Existence of a linear [n, k, d]q code implies the existence of q l-r {1≤ r ≤ (n-k)-(d-1) & (d-1)+1≤l≤ n-k },fractional factorial design and q l ( l=d-1) full factorial design in which all main effects and some of the two-factor interactions are linear trend free, where (d-1) is number of linearly independent columns in matrix H and l is the number of factors/columns in Gα. LINEAR TREND FREE FRACTIONAL FACTORIAL DESIGNS USING DIFFERENT CODES We consider different types of linear codes and generate linear trend free fractional factorial designs using the method of construction described in previous section. For details see Hedayat et.al (1999). a) Reed Muller Codes The r th order binary Reed-Muller code R(r,a) of length n = 2 a , for 0  r  a, is the set of all vectors f, where f(i1,…..ia) is a Boolean function which is a polynomial of degree at most r. For any a and any r, 0 ≤ r ≤ a, there is a binary r th order RM code R(r,a) with the following properties: Length n = 2 a , dimension k = +…+ and minimum distance 2 a-r . The parity check matrix of R(r,a) code is the generator matrixof its dual code. The dual of RM(r,a) is RM(a-r-1,a) code. Example 1: Consider a parity check matrix of RM(2,4) H= Any three columns in parity check matrix H matrix are linearly independent. Here d-1=3, and m =5 (in F2) and we get the designs and corresponding generator matrices depending on the columns selected. According to the selection of columns ,we have the following cases : i)11 columns with none of the column having weight one and there are possible ways of selection where none of the column is with weight (w) one and the corresponding factorial designs will have all main effects linear trend free. All generated matrices and the table contains generator matrices, their corresponding designs, trend free effects and defining relations are available with th authors.
  • 5. A technique to construct linear trend free fractional factorial design using some linear codes Int. J. Stat. Math. 076 ii) For m=5 and d-1=3, n-m=11.The possible choices of columns selected for generator matrix are + + . Suppose we choose the three columns 4 th , 5 th and 6 th to form the matrix MRM = On transposing we get matrix as GRM in which all three rows d-1 =3 are independent. here d-1 =3 and l=5, Using the method of GFS on generator matrix we get a resolution III, fractional factorial design with defining relation I =A1A3A5 = A1A2A3A4= A2A4A5. Using the properties 1and 2 we observe that main effect component of A1, A3, A5 and the two factor interaction A1A2, A1A3 and A2A4 and A3A4 are linear trend free. Cyclic Codes A linear code over GF(q) is said to be cyclic if whenever (c0 ,c1 ,…..,cn-2 ,cn-1 ) is a codeword so also is (c1 ,c2 ,…..,cn-1 ,c0 ). Cyclic arrays can be described by a single generating vector z = ( z0 z1 ……….zn-1 ) such that the generator matrix consists of this vector and its first ( k-1 ) cyclic shifts. The generating vector z is represented by a polynomial z(x) = z0 + z1X + …+zn- 1X n-1 which is called a generator polynomial for the code. If a code is cyclic, so is its dual, and the generator polynomial of its dual can be obtained by the following result given in Macwilliam and Sloane (1977). Theorem 2: If C is a cyclic code of length n over GF(q), with generator polynomial z(x), then the dual code C is also cyclic and has generator polynomial where z * (x) = X deg.z z(x -1 ) is reciprocal polynomial to z(x). Example 2: Let = be the generator polynomial for [15,7,5]2 code. Then the generator polynomial for the dual of this code is given as h (x) = x 7 ( 1+ x -4 +x -6 +x -7 ) = 1+x+x 3 +x 7 Hence the parity check matrix of [15,7,5]2 code is obtained by writing the coefficients and giving the cyclic shift to the coefficients,as given below Hcy= = Any four columns in the H matrix are linearly independent, d-1 = 4 and m=5, n-m= 15-5 = 10 column with none of them of weight(w) one. Then , there are 210 possible selection of such sets of columns from the matrix H. Following the method of construction, we generated all generator matrices from the set of selected columns and listed corresponding fractional factorial designs with their trend free effects and defining relations {available with the authors}.The result of selecting columns 4 th , 5 th , 6 th and 7 th gives the matrix Mcy(8×4 Mcy(8×4=
  • 6. A technique to construct linear trend free fractional factorial design using some linear codes Singh et al. 077 Deleting column of zeroes and transposing MT= we get generator matrix as Gcy of order( (d-1)× l) Gcy = where l =7 and d-1 =4 . Using GFS, a 2 7-3 fractional factorial design is obtained with some main effects A3, A4, A5, A6 and two factor interactions A1A2, A2A3, A2A6, A2A7, A3A4, A3A5, A3A6, A3A7, A5A4, A6A4, A5A6, A5A7 are linear trend free. BCH Codes The BCH codes over GF(q) of length n = q m -1 and designed distance δ is the largest possible cyclic code having zeroes α b , , α b+1 ,…,α b+δ-2 where α є GF( q m ) is the primitive n th root of unity, b is a non negative integer and m is the multiplicative order of q mod n. The parity check matrix of a BCH code with b=1 is given by C=                           1n2222 1n5255 1n3233 1n2 )(.........)()(1 ...... ...... ...... )(.........)(1 )(.........)(1 .........1 where each entry is replaced by the corresponding binary m-tuple. Example 3.: The parity check matrix of [15,5,7]2 code is H= and we observe here that d-1=6 columns in above matrix are linearly independent. Since the matrix cannot be partitioned into the sub matrices ,we consider all 15 columns to get the design and they are Suppose we select columns 8 th ,9 th ,10 th ,11 th ,12 th and 13 th from above matrix H.The matrix obtained by selection is M =
  • 7. A technique to construct linear trend free fractional factorial design using some linear codes Int. J. Stat. Math. 078 On transposing and deleting columns of zeroes we get matrix GBCH as GBCH = Applying GFS on generator matrix with l = 10 and d-1 = 6 we get the design 2 10-4 fractional factorial design with all main effects linear trend free. Some of the designs that can be generated by selecting different columns of the generator matrix GBCH of the code [15,5,7]2 along with their defining relations are available with the author. Ternary Golay Code The Golay code were discovered by M.J.E.Golay in late 1940’s. The (unextended) Golay code are examples of perfect codes. A q-nary code that attains the hamming ( or sphere packing) bound i.e. the one which has codewords, is said to be perfect code. Consider the ternary Golay Code [11,6,5]3 over a ternary alphabet, the relative distance of the codes is as large as it possibly can be for a ternary code, and it satisfies Hamming bound and is therefore a perfect ternary Golay [11,6,5]3 code. For perfect code its dual distance is same as its covering radius. In terms of design the strength is same as the estimation index of an orthogonal array obtained using this code. Golay codes are unique in the sense that binary or ternary codes with same parameters, can be shown to be equivalent to them. Example 4: Consider the parity check matrix of [11,6,5]3 and partition it into two submatrices as given in the method in section 4, we get the matrix H as H= Any four columns in matrix H5×11 are linearly independent d-1 = 4 and m = 5 . Thus, 15 sets of column in which none of the columns has weight (w) one. Whereas the total possible selection of columns to formmatrix M is 325. Suppose we choose columns 1 st ,2 nd ,3 rd and 4 th from matrix H, we get the matrix M = . Next , we generate the matrices of order (4×5) that ensures the linear trend freeness of main effects and following the steps of construction method , we get 3 5-1 fractional factorial design. and then 3 5-1 fractional factorial design I= A1A2 2 A3A4 2 is generated using GFS on matrix Ggolay. The generated design has all its main effects linear trend free with respect to the properties [1-3].
  • 8. A technique to construct linear trend free fractional factorial design using some linear codes Singh et al. 079 BLOCKING IN FRACTIONAL FACTORIAL DESIGN WITH SOME LINEAR TREND FREE EFFECTS. When the block size is smaller than the number of treatment combinations in any factorial experiment, the technique of blocking is used to carry out the analysis. The factorial/fractional factorial experiment thus obtained is known as Blocked fractional factorial experiment. When we go for blocking of fractional factorial design, the block structure affects the linear trend-freeness of the effects. We state here the result in continuation with the properties [1-3] . The first h=n-p-r generates the principle block and let zv , h+1 ≤ v ≤ n-p be the generators of other blocks. Then the following holds: Property 4: For any given v, h+1 ≤ v ≤ n-p and factor say A1 , suppose z1v ≠ 0. Then (a) All (q-1) main effect components are linear trend-free. (b) If ziv = 0 , 2 ≤ i ≤ n-p, then all components of A1 × Ai interaction are linear trend-free. (c) The method described is used to construct Blocked fractional factorial designs with some linear trend free effects. Reed Muller Code: We use here the generator matrix constructed in Example 1 by selecting columns 8 th ,9 th and 10 th We get 2 4-1 fractional factorial design, here q=2, n=3, p=1, r1=1, h=n-p-r1=1 independent generators. Using GFS, h=2 generators forms the principle block and remaining generators form the contents of the other block in which main effects A and C are linear trend free. The confounded effect of the design is ABC. Table 1. 24-1 fractional factorial design with resolution III, I=ABC Cyclic Code: We consider here generator matrix constructed in example 2 and by selecting 1 st , 2 nd ,5 th and 9 th columns We get Gcy= Using GFS 2 5-1 fractional factorial design is obtained n=5, p1=1, s=2, r1=1, h = n-p-r1 = 5-1-1 = 3 independent generators, generates the principal block and other generates the other block contents. Table 1 gives the design generated. Table 2 : 25-1 blocked fractional factorial design with resolution IV, I = ABCE BCH Code: Consider the generator matrix in Example 3 by selecting columns 1 st , 4 th , 5 th , 6 th , 7 th and 8 th as GBCH = 2 9-3 fractional factorial design is obtained, here n=9, p=3 , r1=1 and h = n-p-r1 = 9-3-1= 5 independent generators forms the principal block and remaining generators the other block. Thus, Blocked 2 9-3 fractional factorial.design with defining relation I = DEGJ= ABDEH = BCDF is obtained. Table 3 displays the first 16 runs of the design constituting the principle block and the next 16 the other block and so on. Block 1: (1) , acd , abd, bc Block 2: c , ac, ab, bcd Block 1: (1) , ab, bc, ac, acd, bcd, abd, d Block 2: ce, abce, be, ae, ade, bde, abcde, cde
  • 9. A technique to construct linear trend free fractional factorial design using some linear codes Int. J. Stat. Math. 080 Table 3. 29-3 blocked fractional factorial design with resolution IV, I= DEFJ = ABDEH=BCDF Block 1: (1), aej, dfhj, adefh, abefgh, bfghj, abdegj, bdg, bcej, abc, cdefh, abcdfhj, acfghj, cdhj, acdeh, cf, acefj, abcdefj, bcdfg, abcegh, bcghj, bdeh, bdhj, bef, abf, adfg, defgj, aghj, egh. Block 2: abdgh, bdeghj, abfgj, befg, def, adfj, ehj, ah, acdeghj, bcfg, acefg, cfgj, bcdfj, abcdef, bch, abcehj, abcgj, abcdfgh, bcdefghj, cefhj, acfh, bcd, acdj, aeg, gj, adefg. dfgh, bfh, abefhj, bdj, abde. Ternary Golay Code: Consider the generator matrix from Example 4 by selecting columns 1 st ,2 nd , 3 rd and 4 th of matrix H as Ggolay = 3 5-1 fractional factorial design is obtained. Here n= 5, p1=1 and q=3 , r1=1 ,h=n-p-r1=3 independent generators are obtained. Thus, principal block is generated by first three generators and remaining two generates the content of the other block. First 3 3 treatment combinations in 3 5-1 fractional factorial design (given in Table 4 ) forms the principle block and next 3 3 combinations forms the contents of second block and the remaining the third block contents. The factors are denoted by A,B,C,D and E. Table 4: 35-1 blocked fractional factorial design with resolution IV, I= ABCDE Block 1: (1), abcde, a 2 b 2 c 2 d 2 e 2 , abc 2 d 2 , a 2 b 2 e, cde 2 , a 2 b 2 cd, c 2 d 2 e, abe 2 , ab 2 ce 2 , a 2 c 2 d , bd 2 e, a 2 d 2 e 2 , bc, ab 2 c 2 de, bc 2 de 2 , ab 2 d, a 2 ce, a 2 bc 2 e, b 2 de 2 , acd 2 , b 2 cd 2 e, ac 2 d 2 , b 2 ce, ade, a 2 bcd 2 e 2 , b 2 d 2 . Block 2: a 2 bde 2 , b 2 cd 2 , ac 2 e, b 2 c 2 e 2 , ad, a 2 bcd 2 e, acd 2 e 2 , a 2 bc 2 , b 2 de, cde, abc 2 d 2 e 2 , a 2 b 2 , abe, a 2 b 2 cde 2 , c 2 d 2 , a 2 b 2 c 2 d 2 e, e 2 , abcd, ab 2 c 2 d, a 2 d 2 e, bce 2 , a 2 c, ab 2 cd, ab 2 d 2 e 2 , bd 2 , ab 2 de, a 2 c 2 de 2 . Block 3: ab 2 d 2 e, a 2 ce 2 , bc 2 d, a 2 c 2 de, bd 2 e 2 , ab 2 c, bce, ab 2 c 2 de 2 , a 2 d 2 , a 2 bcd 2 , b 2 d 2 e, ade, b 2 d, acd 2 e, a 2 bc 2 e 2 , ac 2 , a 2 bde, b 2 cd 2 e 2 , c 2 d 2 e 2 , ab, a 2 b 2 cde, abcde 2 , a 2 b 2 c 2 d 2 , e, a 2 b 2 e 2 , cd, abc 2 d 2 e. CONCLUSION A technique to construct fractional factorial designs with factors at q level, q>2, with some linear trend free effects using parity check/generator matrix of linear code is developed. Since the generator matrix is not unique in nature, therefore fractional factorial designs with linear trend free effects with same/ different resolution can be easily constructed. REFERENCES Adekeys K, Kunert J (2006). On the Comparison of Run Orders of Unreplicated 2 k-p designs in the presence of a Time Trend. Metrika. 63: 257-269. Bailey RA, Cheng CS, Kipnis P (1992). Construction of trend resistant factorial designs. Statist. Sinica. 2:393- 411. Betsumiya K, Harada M (2001). Classification of Formally Self-Dual Even Codes of length to 16. Designs, Codes and Cryptography. 23: 325-332. Cheng CS, Jacroux M. (1988).On the construction of trend free run-order of two –level factorial designs. Journal of American statistical association. 83:1152-1158. Cox DR (1951). Some systematic experimental designs. Biometrica. 38: 312- 323. Coster DC, Cheng CS (1998). Minimum cost trend free run orders of fractional factorial design. The Annals of Statistics. 16:1188-1205. Daniel C, Wilcoxon F (1966). Factorial 2 p-q plans robust against linear and quadratic trends.Technometrics. 8: 259- 278. Draper NR, Stoneman DM (1968). Factor changes and linear trends in eight run two level factorial designs. Technometrics. 10: 301-311.
  • 10. A technique to construct linear trend free fractional factorial design using some linear codes Singh et al. 081 Freeny AE, Lai WYC (1997). Planarization by chemical Mechanical polishing: A rateand uniformity study. Statistical case studies for industrial Process Improvement. V. Czitrom and P.D. Spagon eds Philadelphia, Alexandria: ASA- SIAM, pp. 251-264. Hedayat, AS, Sloane NJA, Stuffken J (1999). Orthogonal Arrays: Theory and applications. New York. US: Springer and Verlag . 61-67, 87-108, 247-288. Hill HH (1960).Experimental Designs to adjust for time trends. Technometrics. 2: 67-82. Hilow H (2012). Minimum cost linear trend free fractional factorial designs. Journal of statistical theory and practice. 6: 580-589. Joiner BL, Campbell C (1976). Designing experiments when run order is important. Technometrics. 18: 249-259. Mac Williams FJ, Sloane NJA (1977). The Theory of Error-Correcting Codes. North- Holland. Amesterdam: Oxford. Mee R, Romanora A (2010). Constructing and Analyzing Two-level Trend Robust designs. Quality Enigeneering. 22: 306-316. Wang PC, Wu S (2013). Trend Free Designs in Blocked Fractional Factorial Experiments. Communications in Statistics- Theory and Methods. 42(10):1870-1878. Accepted December 31, 2015. Citation: Singh P, Thapliyal P, Budhraja V (2016). A Technique to Construct Linear Trend Free Fractional Factorial Design Using Some Linear Codes. International Journal of Statistics and Mathematics 3(1): 073-081. Copyright: © 2016 Singh et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are cited.