PBIBD (version 1.3)

apbibd: Calculates the efficiencies of PBIB designs with any number of associate classes.

Description

This function calculates the different kinds of efficiencies and the overall efficiency factor of Partially Balanced Incomplete Block Designs with any number of associate classes. The total number of treatments i.e. v, replications i.e. r, block size i.e. k, vector l of lambda's (lambda i being the ith element of vector l), vector n of number of associates (n i, i.e. number of ith associates, being the ith element of vector n), a list P of P-matrices of the association scheme of the design (Pi being the ith matrix of the list P) are to be supplied as input to the function.

Usage

apbibd(v, r, k, l, n, P)

Arguments

v

Total number of treatments of the design

r

Replication of the treatments in the design

k

Block size of the design

l

A vector containing lambda 1, lambda 2, lambda 3,..., lambda m as its first, second, third,..., mth elements

n

A vector containing n1, n 2, n 3, ..., n m as its first, second, third,..., mth elements

P

A list containing P-matrices of the association scheme of the design such that P1 is its first element, P2 is second element, P3 is third element,..., Pm is the mth element

Value

Returns a list with (m+1) components:

E1

Efficiency E1 of the design

E2

Efficiency E2 of the design and so on …

Em

Efficiency Em of the design

E

Overall efficiency factor of the design

Examples

# NOT RUN {
v<-25
r<-9
k<-9
l<-c(5,2,5,2,5)
n<-c(2,8,2,8,4)
P1<-matrix(c(0,0,1,0,0,0,0,0,4,4,1,0,1,0,0,0,4,0,4,0,0,4,0,0,0),nrow=5,ncol=5)
P2<-matrix(c(0,0,0,1,1,0,0,1,3,3,0,1,0,1,0,1,3,1,3,0,1,3,0,0,0),nrow=5,ncol=5)
P3<-matrix(c(1,0,1,0,0,0,4,0,4,0,1,0,0,0,0,0,4,0,0,4,0,0,0,4,0),nrow=5,ncol=5)
P4<-matrix(c(0,1,0,1,0,1,3,1,3,0,0,1,0,0,1,1,3,0,0,3,0,0,1,3,0),nrow=5,ncol=5)
P5<-matrix(c(0,2,0,0,0,2,6,0,0,0,0,0,0,2,0,0,0,2,6,0,0,0,0,0,3),nrow=5,ncol=5)
P<-list(P1,P2,P3,P4,P5)
apbibd(v,r,k,l,n,P)
# }