wtaoptbtib {Aoptbdtvc}

Weighted A-optimal balanced treatment incomplete block designs

Description

This function generates weighted A-optimal balanced treatment incomplete block design for test vs control comparisons with specified parameters

Usage

wtaoptbtib(v,b,k,alpha,rho=0,ntrial=5)

Arguments

v

number of test treatments

b

number of blocks

k

block size

alpha

Weight for test versus test comparisons. Should be between 0 to 1

rho

rho=0

ntrial

number of trials, default is 5

Value

It either returns a text message or a design. If a design is found, it returns a list with following components

parameters

parameters of the design

design

generated weighted A-optmal BTIB design

N

incidence matrix of the generated weighted A-optmal BTIB design

NNP

concurrence matrix of the generated design

Note

The function is useful to construct weighted A-optimal BTIB designs upto 30 test treatments and up to block size 10. May not be very useful beyond 30 test treatments. For k<=3, designs with larger number of test treatments may be obtained.

Author(s)

Baidya Nath Mandal <[e-mail]>

References

Gupta, V., Ramana, D. and Parsad, R. (1999). Weighted A-efficiency of block designs for making treatment-control and treatment-treatment comparisons, Journal of statistical planning and inference 77(2), 301-319.

Mandal, B. N., Parsad, R. and Dash, S. (2017). A-optimal block designs for comparing test treatments with control treatment(s) - an algorithmic approach, upcoming project report, ICAR-Indian Agricultural Statistics Research Institute, New Delhi, India.

Examples

##construct a weighted A-optimal BTIB design with 4 test treatments in 6 blocks each of size 4 
##with weights to test vs test treatments comparisons as 0.6
wtaoptbtib(v=4,b=6,k=4,alpha=0.6,rho=0)
##construct an A-optimal BTIB design with 9 test treatments in 12 blocks each of size 4 
##with weights to test vs test treatments comparisons as 0
wtaoptbtib(v=9,b=12,k=4,alpha=0,rho=0)
##design not found
## Not run: wtaoptbtib(v=3,b=6,k=5,alpha=0.2,rho=0)
##BTIB design does not exist for these parameters
#Not run
wtaoptbtib(3,4,3,0.2,0)