BIBsize: Balanced incomplete blocksize

In daewr: Design and Analysis of Experiments with R

Description

This function computes the number of blocks, treatment frequency and lambda for a potential BIB design

Usage

BIBsize(t,k)

Arguments

t	input - number of levels of the treatment factor
k	input - blocksize or number of experimental units per block

Value

a list containing the b=number of blocks, r=number of treatment replicates and lambda for a potential BIB design with t levels of treatment factor and blocksize k.

Author(s)

John Lawson

Examples

BIBsize(6,3)
## The function is currently defined as
BIBsize<-function(t,k)
{
  b<-t
  r<-0
  lambda<-0
  check<-0
  while (check==0) {
   while (r==0) {
     #cat("r=",r)
     testr<-(b*k)/t
     #cat("testr=",testr,"b=",b)
     if (testr==floor(testr)) {
       r<-testr
       } else {
       b<-b+1
       }
     }
      #cat("b=",b, "r=",r)
      testl<-(r*(k-1))/(t-1)
      #cat("testl=",testl,"b=",b)
      if (testl==floor(testl)) {
       lambda<-testl
       check=1
       } else {
       r<-0
       b<-b+1
      #cat("b=",b, "r=",r)
        }

    #cat("lambda=",lambda)
    }
  cat("Posible BIB design with b=",b," and r=",r," lambda=",lambda,"\n") 
}
       

daewr documentation built on March 13, 2021, 3:01 a.m.