graphoptBbd {Boptbd}

Creates the graphical layout of resultant optimal design

Description

The function graphoptBbd creates the graphical layout of resultant Bayesian A-, or D-optimal block design on a separate pop-up GUI tcltk window.

Usage

graphoptBbd(trt.N, blk.N, alpha, beta, OptdesF, Optcrit)

Arguments

trt.N

integer, specifying number of treatments, v.

blk.N

integer, specifying number of blocks, b.

alpha

numeric, representing parameter of beta distribution.

beta

numeric, representing parameter of beta distribution.

OptdesF

matrix, a 2 x b obtained Bayesian optimal block design.

Optcrit

character specifying the optimality criteria to be used. Thus, Optcrit takes the letter "A" and "D" for Bayesian A- and D-optimal block designs, respectively.

Details

Detail discussions concerning the constructions of a graphs can be found in igraph R package.

Value

Returns the graphical layout of the resultant Bayesian optimal block design 'OptdesF' on a separate pop-up window. Furthermore, the function saves the graphical layout of the resultant Bayesain optimal block design in .pdf format in a temporary directory.

Note

When closing a pop-up window for graphical layout of the resultant designs (Graph plot), if the window is closed by clicking on the red button with "X" sign (top-right), the warning message "Warning message: In rm(list = cmd, envir = .tkplot.env): object 'tkp ...' not found" will occur in R-console irrespective of what command is executed next. To resolve this warning message, click on "close" menu that is located at the top-left of the graph plot pop-up window when closing this window.

Author(s)

Dibaba Bayisa Gemechu, Legesse Kassa Debusho, and Linda Haines

See Also

igraph

Examples

    ##To create the graphical layout of the D-optimal block design
    ##obtained using the treatment exchange algorithm for 
    trt.N <- 10  #Number of treatments
    blk.N <- 10  #Number of blocks
    alpha <- 0.1 #alpha value 
    beta  <- 0.1 #beta value 
    OptdesF <- rbind(1:10, c(2:10,1)) #Bayesian A-optimal block design (loop design)

    graphoptBbd(trt.N = 10, blk.N = 10, alpha = 0.1, beta = 0.1, OptdesF, Optcrit = "A")