| Type: | Package |
| Title: | Tool for Construction of Two-Phase Experimental Designs |
| Version: | 1.0.0 |
| Maintainer: | Akhilesh Jha <jha.akhilesh09@gmail.com> |
| Description: | Provides functions to construct two-phase design layouts, compute treatment- and block-incidence matrices, derive C-matrices for residual, direct, and interaction effects, and calculate the efficiency factor for two-phase experimental designs with factorial treatment structure. |
| License: | GPL-3 |
| Encoding: | UTF-8 |
| RoxygenNote: | 7.3.2 |
| Imports: | MASS, Matrix |
| Depends: | R (≥ 4.0) |
| Suggests: | testthat (≥ 3.0.0) |
| Config/testthat/edition: | 3 |
| NeedsCompilation: | no |
| Packaged: | 2025-10-03 20:20:37 UTC; Akhilesh Jha |
| Author: | Akhilesh Jha [aut, cre], Cini Varghese [aut], Seema Jaggi [aut], Eldho Varghese [aut], Mohd Harun [aut] |
| Repository: | CRAN |
| Date/Publication: | 2025-10-08 19:40:17 UTC |
Two-Phase Experimental Design Development
Description
'TPDDev()' constructs a two-phase experimental design layout, computes the information matrices for residual, direct, and interaction effects, and calculates the efficiency factor.
Usage
TPDDev(d1, base_d2)
Arguments
d1 |
Numeric matrix. Phase-I layout (rows = blocks, columns = plots per block). |
base_d2 |
Numeric matrix. Base Phase-II layout (rows = blocks, columns = plots per block). |
Details
Steps performed:
Checks that 'ncol(d1)' equals 'nrow(base_d2)'.
Expands the Phase-II layout for each Phase-I block.
Builds the combined layout and treatment-incidence matrices.
Computes block-incidence matrix and overall mean matrix.
Derives overall C-matrix for factorial treatment effects.
Projects C-matrix to obtain residual, direct, and interaction effect matrices.
Calculates the efficiency factor of the design.
Value
A list containing design parameters, combined two-phase design layout, information matrices, and efficiency factor.
Examples
d1 <- matrix(c(
1,2,3,4,5,6,7,
2,3,4,5,6,7,8,
3,4,5,6,7,8,1,
4,5,6,7,8,1,2,
5,6,7,8,1,2,3,
6,7,8,1,2,3,4,
7,8,1,2,3,4,5,
8,1,2,3,4,5,6
), nrow = 8, byrow = TRUE)
base_d2 <- matrix(c(
1,2,4,
2,3,5,
3,4,6,
4,5,7,
5,6,1,
6,7,2,
7,1,3
), ncol = 3, byrow = TRUE)
result <- TPDDev(d1, base_d2)
result$efficiency_factor
result$C_dir