graphoptBbd function

Creates the graphical layout of resultant optimal design

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

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.

Returns

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.

Author(s)

Dibaba Bayisa Gemechu, Legesse Kassa Debusho, and Linda Haines

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.

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")