Central composite design

From Wikipedia, the free encyclopedia

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 (perhaps 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 −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 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:

${\mathbf 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}}.$

Then d is the vertical concatenation:

${\mathbf d}={\begin{bmatrix}{\mathbf F}\\{\mathbf C}\\{\mathbf E}\end{bmatrix}}.$

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:

${\mathbf X}={\begin{bmatrix}{\mathbf 1}&{\mathbf d}&{\mathbf d}(1)\times {\mathbf d}(2)&{\mathbf d}(1)\times {\mathbf d}(3)&\cdots &{\mathbf d}(k-1)\times {\mathbf d}(k)&{\mathbf d}(1)^{2}&{\mathbf d}(2)^{2}&\cdots &{\mathbf d}(k)^{2}\end{bmatrix}},$

where d(i) represents the i^th 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):

Orthogonal design:: $\alpha =(Q\times F/4)^{{1/4}}\,\!$ , where $Q=({\sqrt {F+T}}-{\sqrt {F}})^{2}$ ;
Rotatable design: α = F^1/4 (the design implemented by MATLAB’s ccdesign function).

References

Myers, Raymond H. Response Surface Methodology. Boston: Allyn and Bacon, Inc., 1971

v t e Design of experiments

Scientific Method	Scientific experiment Statistical design Control Internal & external validity Experimental unit Blinding Optimal design: Bayesian Random assignment Randomization Restricted randomization Replication versus subsampling Sample size

Treatment & Blocking	Treatment Effect size Contrast Interaction Confounding Orthogonality Blocking Covariate Nuisance variable

Models & Inference	Linear regression Ordinary least squares Bayesian Random effect Mixed model Hierarchical model: Bayesian Analysis of variance (Anova) Cochran's theorem Manova (multivariate) Ancova (covariance) Compare means Multiple comparison

Designs: Completely Randomized	Factorial Fractional factorial Plackett-Burman Taguchi Response surface methodology Polynomial & rational modeling Box-Behnken Central composite Block Generalized randomized block design (GRBD) Latin square Graeco-Latin square Orthogonal array Latin hypercube Repeated measures design Crossover study Randomized controlled trial Sequential analysis Sequential probability ratio test

Glossary Category Statistics portal Statistical outline Statistical topics

Statistics

Descriptive statistics

Continuous data

Location	Mean (Arithmetic, Geometric, Harmonic) Median Mode

Dispersion	Range Standard deviation Coefficient of variation Percentile Interquartile range

Shape	Variance Skewness Kurtosis Moments L-moments

Count data

Index of dispersion

Summary tables

Dependence

Statistical graphics

Data collection

Designing studies	Effect size Standard error Statistical power Sample size determination

Survey methodology	Sampling Stratified sampling Cluster sampling Opinion poll Questionnaire

Controlled experiment	Design of experiments Randomized experiment Random assignment Replication Blocking Factorial experiment Optimal design

Uncontrolled studies	Natural experiment Quasi-experiment Observational study

Statistical inference

Statistical theory	Sampling distribution Order statistic Scan statistic Record value Sufficiency Completeness Exponential family Permutation test (Randomization test) Empirical distribution Bootstrap U statistic Efficiency Asymptotics Robustness

Frequentist inference	Unbiased estimator (Mean unbiased minimum variance, Median unbiased) Biased estimators (Maximum likelihood, Method of moments, Minimum distance, Density estimation) Confidence interval Testing hypotheses Power Parametric tests (Likelihood-ratio, Wald, Score)

Specific tests	Z (normal) Student's t-test F Goodness of fit (Chi-squared, G, Sample source, sample normality, Skewness & kurtosis Normality, Model comparison, Model quality) Signed-rank (1-sample, 2-sample, 1-way anova) Shapiro–Wilk Kolmogorov–Smirnov

Bayesian inference	Bayesian probability Prior Posterior Credible interval Bayes factor Bayesian estimator Maximum posterior estimator

Correlation and regression analysis

Correlation	Pearson product–moment correlation Partial correlation Confounding variable Coefficient of determination

Regression analysis	Errors and residuals Regression model validation Mixed effects models Simultaneous equations models MARS

Linear regression	Simple linear regression Ordinary least squares General linear model Bayesian regression

Non-standard predictors	Nonlinear regression Nonparametric Semiparametric Isotonic Robust Heteroscedasticity Homoscedasticity

Generalized linear model	Exponential families Logistic (Bernoulli) Binomial Poisson

Partition of variance	Analysis of variance (ANOVA) Analysis of covariance Multivariate ANOVA Degrees of freedom

Categorical, multivariate, time-series, or survival analysis

Categorical data

Multivariate statistics

Time series analysis

General	Decomposition Trend Stationarity Seasonal adjustment Exponential smoothing Cointegration

Specific tests	Granger causality Q-Statistic Durbin–Watson

Time domain	ACF PACF XCF ARMA model ARIMA model ARCH Vector autoregression

Frequency domain	Spectral density estimation Fourier analysis

Survival analysis

Applications

Biostatistics	Bioinformatics Clinical trials & studies Epidemiology Medical statistics

Engineering statistics	Chemometrics Methods engineering Probabilistic design Process & Quality control Reliability System identification

Social statistics	Actuarial science Census Crime statistics Demography Econometrics National accounts Official statistics Population Psychometrics

Spatial statistics	Cartography Environmental statistics Geographic information system Geostatistics Kriging

Category
Portal
Outline
Index

This article is issued from Wikipedia. The text is available under the Creative Commons Attribution/Share Alike; additional terms may apply for the media files.