cmatbd.mae {optbdmaeAT}

Computes the treatment information matrix

Description

Computes the information matrix (C-matrix) for treatment effects under either the linear fixed effects model or the linear mixed effects model setting for a given block design of size 2.

Usage

cmatbd.mae(trt.N, blk.N, theta, des)

Arguments

trt.N

integer, specifying number of treatments, v.

blk.N

integer, specifying number of arrays, b.

theta

numeric, representing a function of the ratio of random array variance and random error variance. It takes any value between 0 and 1, inclusive.

des

matrix, a 2 x b block design with b blocks of size k = 2 and v treatments.

Value

Returns a v x v treatment information matrix (C-matrix).

Author(s)

Dibaba Bayisa Gemechu, Legesse Kassa Debusho, and Linda Haines

References

Debusho, L. K., Gemechu, D. B., and Haines, L. M. (2016). Algorithmic construction of optimal block designs for two-colour cDNA microarray experiments using the linear mixed model. Under review.

Gemechu D. B., Debusho L. K. and Haines L. M. (2014). A-optimal designs for two-colour cDNA microarray experiments using the linear mixed effects model. Peer-reviewed Proceedings of the Annual Conference of the South African Statistical Association for 2014 (SASA 2014), Rhodes University, Grahamstown, South Africa. pp 33-40, ISBN: 978-1-86822-659-7.

See Also

optbdmaeAT, fixparbd.mae, intcbd.mae

Examples

##Information matrix

     trt.N <- 3 
     blk.N <- 3 
     theta <- 0.2 
     dsgn <- intcbd.mae(trt.N = 3, blk.N = 3)

     cmatbd.mae(trt.N = 3, blk.N = 3, theta = 0.2, des = dsgn)