Expander walk sampling: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Citation bot
m [344]+: issue.
 
en>ChrisGualtieri
m Uses: General Fixes using AWB
Line 1: Line 1:
Myrtle Benny is how I'm called and I really feel comfy when people use the complete title. Hiring is my profession. California is our beginning location. To do aerobics is a factor that I'm completely addicted to.<br><br>my web page :: [http://xrambo.com/user/NEme xrambo.com]
In [[statistics]], a '''central composite design''' is an experimental design, useful in [[response surface methodology]], for building a second order (quadratic) model for the [[response variable]] without needing to use a complete three-level [[factorial experiment]].  
 
After the designed experiment is performed, [[linear regression]] is used, sometimes iteratively, to obtain results. Coded variables are often used when constructing this design.
 
==Implementation==
 
The design consists of three distinct sets of experimental runs:
 
# A [[factorial design|factorial]] (perhaps [[Fractional factorial design|fractional]]) design in the factors studied, each having two levels;
# A set of ''center points'', experimental runs whose values of each factor are the medians of the values used in the factorial portion. This point is often replicated in order to improve the precision of the experiment;
# A set of ''axial points'', experimental runs identical to the centre points except for one factor, which will take on values both below and above the median of the two factorial levels, and typically both outside their range. All factors are varied in this way.
 
==Design matrix ==
The design matrix for a central composite design experiment involving ''k'' factors is derived from a matrix, '''d''', containing the following three different parts corresponding to the three types of experimental runs:
#The matrix '''F''' obtained from the factorial experiment.  The factor levels are scaled so that its entries are coded as +1 and &minus;1. 
#The matrix '''C''' from the center points, denoted in coded variables as (0,0,0,...,0), where there are ''k'' zeros.
#A matrix '''E''' from the axial points, with 2''k'' rows. Each factor is sequentially placed at ±α and all other factors are at zero. The value of α is determined by the designer; while arbitrary, some values may give the design desirable properties. This part would look like:
 
:<math> \bold E = \begin{bmatrix}
  \alpha  & 0 & 0 &  \cdots  &  \cdots  &  \cdots  & 0  \\
  { - \alpha } & 0 & 0 &  \cdots  &  \cdots  &  \cdots  & 0  \\
  0 & \alpha  & 0 &  \cdots  &  \cdots  &  \cdots  & 0  \\
  0 & { - \alpha } & 0 &  \cdots  &  \cdots  &  \cdots  & 0  \\
    \vdots  & {} & {} & {} & {} & {} & \vdots  \\
  0 & 0 & 0 & 0 &  \cdots  &  \cdots  & \alpha  \\
  0 & 0 & 0 & 0 &  \cdots  &  \cdots  & { - \alpha }  \\
\end{bmatrix}. </math>  
 
Then '''d''' is the vertical concatenation:
:<math> \bold d = \begin{bmatrix} \bold F \\ \bold C \\ \bold E
                  \end{bmatrix}. </math>
 
The design matrix '''X''' used in linear regression is the horizontal concatenation of a column of 1s (intercept), '''d''', and all elementwise products of a pair of columns of '''d''':
 
:<math>\bold X = \begin{bmatrix} \bold 1 & \bold d & \bold d(1)\times\bold d(2) & \bold d(1)\times\bold d(3) & \cdots & \bold d(k-1)\times\bold d(k) & \bold d(1)^2 &\bold d(2)^2 &\cdots & \bold d(k)^2 \end{bmatrix}, </math>
 
where '''d'''(''i'') represents the ''i''<sup>th</sup> column in '''d'''.  
=== Choosing &alpha; ===
 
There are many different methods to select a useful value of α. Let ''F'' be the number of points due to the factorial design and ''T'' = 2''k'' + ''n'', the number of additional points, where ''n'' is the number of central points in the design. Common values are as follows (Myers, 1971):
#'''Orthogonal design:''': <math>\alpha  = (Q\times F/4)^{1/4}\,\!</math>, where <math> Q = (\sqrt{F + T}  -\sqrt{F})^2 </math>;
#'''Rotatable design''': α = ''F''<sup>1/4</sup> (the design implemented by [[MATLAB]]’s ''ccdesign'' function).
 
==References==
Myers, Raymond H. ''Response Surface Methodology''. Boston: Allyn and Bacon, Inc., 1971
{{Experimental design}}
{{Statistics}}
[[Category:Design of experiments]]
[[Category:Educational research]]
[[Category:Social sciences methodology]]
[[Category:Evaluation methods]]
[[Category:Mathematical optimization]]

Revision as of 00:45, 24 October 2013

In statistics, a central composite design is an experimental design, useful in response surface methodology, for building a second order (quadratic) model for the response variable without needing to use a complete three-level factorial experiment.

After the designed experiment is performed, linear regression is used, sometimes iteratively, to obtain results. Coded variables are often used when constructing this design.

Implementation

The design consists of three distinct sets of experimental runs:

  1. A factorial (perhaps fractional) design in the factors studied, each having two levels;
  2. A set of center points, experimental runs whose values of each factor are the medians of the values used in the factorial portion. This point is often replicated in order to improve the precision of the experiment;
  3. A set of axial points, experimental runs identical to the centre points except for one factor, which will take on values both below and above the median of the two factorial levels, and typically both outside their range. All factors are varied in this way.

Design matrix

The design matrix for a central composite design experiment involving k factors is derived from a matrix, d, containing the following three different parts corresponding to the three types of experimental runs:

  1. The matrix F obtained from the factorial experiment. The factor levels are scaled so that its entries are coded as +1 and −1.
  2. The matrix C from the center points, denoted in coded variables as (0,0,0,...,0), where there are k zeros.
  3. A matrix E from the axial points, with 2k rows. Each factor is sequentially placed at ±α and all other factors are at zero. The value of α is determined by the designer; while arbitrary, some values may give the design desirable properties. This part would look like:
E=[α000α0000α000α000000α0000α].

Then d is the vertical concatenation:

d=[FCE].

The design matrix X used in linear regression is the horizontal concatenation of a column of 1s (intercept), d, and all elementwise products of a pair of columns of d:

X=[1dd(1)×d(2)d(1)×d(3)d(k1)×d(k)d(1)2d(2)2d(k)2],

where d(i) represents the ith column in d.

Choosing α

There are many different methods to select a useful value of α. Let F be the number of points due to the factorial design and T = 2k + n, the number of additional points, where n is the number of central points in the design. Common values are as follows (Myers, 1971):

  1. Orthogonal design:: α=(Q×F/4)1/4, where Q=(F+TF)2;
  2. Rotatable design: α = F1/4 (the design implemented by MATLAB’s ccdesign function).

References

Myers, Raymond H. Response Surface Methodology. Boston: Allyn and Bacon, Inc., 1971 Template:Experimental design Template:Statistics