1 2R version 3.6.0 (2019-04-26) -- "Planting of a Tree" 3Copyright (C) 2019 The R Foundation for Statistical Computing 4Platform: x86_64-w64-mingw32/x64 (64-bit) 5 6R is free software and comes with ABSOLUTELY NO WARRANTY. 7You are welcome to redistribute it under certain conditions. 8Type 'license()' or 'licence()' for distribution details. 9 10R is a collaborative project with many contributors. 11Type 'contributors()' for more information and 12'citation()' on how to cite R or R packages in publications. 13 14Type 'demo()' for some demos, 'help()' for on-line help, or 15'help.start()' for an HTML browser interface to help. 16Type 'q()' to quit R. 17 18> ## test column order options of oa.design 19> require(DoE.base) 20Loading required package: DoE.base 21Loading required package: grid 22Loading required package: conf.design 23Registered S3 method overwritten by 'DoE.base': 24 method from 25 factorize.factor conf.design 26 27Attaching package: 'DoE.base' 28 29The following objects are masked from 'package:stats': 30 31 aov, lm 32 33The following object is masked from 'package:graphics': 34 35 plot.design 36 37The following object is masked from 'package:base': 38 39 lengths 40 41> 42> P3.3(oa.design(L18, nlevels=c(2,3,3,3),columns="order"), detailed=TRUE) 43The columns of the array have been used in order of appearance. 44For designs with relatively few columns, 45the properties can sometimes be substantially improved 46using option columns with min3 or even min34. 47 48 length3 frequency 49 0.0000 2 50 0.5000 1 51 0.6667 1 52attr(,"A3") 53[1] 1.166667 54attr(,"detail") 55 1:2:3 1:2:4 1:3:4 2:3:4 560.0000 0.0000 0.6667 0.5000 57> P3.3(oa.design(L18, nlevels=c(2,3,3,3),columns="min3"), detailed=TRUE) 58 length3 frequency 59 0.0 3 60 0.5 1 61attr(,"A3") 62[1] 0.5 63attr(,"detail") 641:2:3 1:2:4 1:3:4 2:3:4 65 0.0 0.0 0.0 0.5 66> P3.3(oa.design(L18, nlevels=c(2,3,3,3),columns="min34"), detailed=TRUE) 67 length3 frequency 68 0.0 3 69 0.5 1 70attr(,"A3") 71[1] 0.5 72attr(,"detail") 731:2:3 1:2:4 1:3:4 2:3:4 74 0.0 0.0 0.0 0.5 75> P3.3(oa.design(L18, nlevels=c(2,3,3,3),columns="minRPFT"), detailed=TRUE) 76 length3 frequency 77 0.0 3 78 0.5 1 79attr(,"A3") 80[1] 0.5 81attr(,"detail") 821:2:3 1:2:4 1:3:4 2:3:4 83 0.0 0.0 0.0 0.5 84> P3.3(oa.design(L18, nlevels=c(2,3,3,3),columns="minRelProjAberr"), detailed=TRUE) 85 length3 frequency 86 0.0 3 87 0.5 1 88attr(,"A3") 89[1] 0.5 90attr(,"detail") 911:2:3 1:2:4 1:3:4 2:3:4 92 0.0 0.0 0.0 0.5 93> 94> 95> P3.3(oa.design(L18, nlevels=c(2,3,3,3,3),columns="min3"), parft=TRUE) 96 length3.parft frequency 97 0.00 6 98 0.25 1 99 0.50 3 100attr(,"A3") 101[1] 3.5 102attr(,"sumPARFT3") 103[1] 1.75 104> P3.3(oa.design(L18, nlevels=c(2,3,3,3,3),columns="min34"), parft=TRUE) 105 length3.parft frequency 106 0.00 6 107 0.25 1 108 0.50 3 109attr(,"A3") 110[1] 3.5 111attr(,"sumPARFT3") 112[1] 1.75 113> P3.3(oa.design(L18, nlevels=c(2,3,3,3,3),columns="minRPFT"), parft=TRUE) 114 length3.parft frequency 115 0.00 6 116 0.25 1 117 0.50 3 118attr(,"A3") 119[1] 3.5 120attr(,"sumPARFT3") 121[1] 1.75 122> P3.3(oa.design(L18, nlevels=c(2,3,3,3,3),columns="minRelProjAberr"), parft=TRUE) 123 length3.parft frequency 124 0.00 6 125 0.25 1 126 0.50 3 127attr(,"A3") 128[1] 3.5 129attr(,"sumPARFT3") 130[1] 1.75 131> 132> P3.3(oa.design(L18, nlevels=c(3,3,2,3,3),columns="min3"), parftdf=TRUE) 133 length3.parftdf frequency 134 0.00 6 135 0.25 1 136 0.50 3 137attr(,"A3") 138[1] 3.5 139attr(,"sumPARFTdf3") 140[1] 1.75 141> P3.3(oa.design(L18, nlevels=c(3,3,2,3,3),columns="min34"), parftdf=TRUE) 142 length3.parftdf frequency 143 0.00 6 144 0.25 1 145 0.50 3 146attr(,"A3") 147[1] 3.5 148attr(,"sumPARFTdf3") 149[1] 1.75 150> P3.3(oa.design(L18, nlevels=c(3,3,2,3,3),columns="minRPFT"), parftdf=TRUE) 151 length3.parftdf frequency 152 0.00 6 153 0.25 1 154 0.50 3 155attr(,"A3") 156[1] 3.5 157attr(,"sumPARFTdf3") 158[1] 1.75 159> P3.3(oa.design(L18, nlevels=c(3,3,2,3,3),columns="minRelProjAberr"), parftdf=TRUE) 160 length3.parftdf frequency 161 0.00 6 162 0.25 1 163 0.50 3 164attr(,"A3") 165[1] 3.5 166attr(,"sumPARFTdf3") 167[1] 1.75 168> 169> P3.3(oa.design(L18, nlevels=c(2,3,3,3,3,3),columns="min3"), rela=TRUE) 170 length3.rela frequency 171 0.0000 7 172 0.2500 7 173 0.5000 3 174 0.6667 3 175attr(,"A3") 176[1] 8.5 177attr(,"rA3") 178[1] 5.25 179attr(,"GR") 180[1] 3.1835 181> P3.3(oa.design(L18, nlevels=c(2,3,3,3,3,3),columns="min34"), rela=TRUE) 182 length3.rela frequency 183 0.0000 7 184 0.2500 7 185 0.5000 3 186 0.6667 3 187attr(,"A3") 188[1] 8.5 189attr(,"rA3") 190[1] 5.25 191attr(,"GR") 192[1] 3.1835 193> P3.3(oa.design(L18, nlevels=c(2,3,3,3,3,3),columns="minRPFT"), rela=TRUE) 194 length3.rela frequency 195 0.0000 7 196 0.2500 7 197 0.5000 3 198 0.6667 3 199attr(,"A3") 200[1] 8.5 201attr(,"rA3") 202[1] 5.25 203attr(,"GR") 204[1] 3.1835 205> P3.3(oa.design(L18, nlevels=c(2,3,3,3,3,3),columns="minRelProjAberr"), rela=TRUE) 206 length3.rela frequency 207 0.0000 7 208 0.2500 7 209 0.5000 3 210 0.6667 3 211attr(,"A3") 212[1] 8.5 213attr(,"rA3") 214[1] 5.25 215attr(,"GR") 216[1] 3.1835 217> 218> GRind(oa.design(L18, nlevels=c(2,3,3,3,3,3,3))) 219The columns of the array have been used in order of appearance. 220For designs with relatively few columns, 221the properties can sometimes be substantially improved 222using option columns with min3 or even min34. 223 224$GRs 225 GR GRind 226 3 3 227 228$GR.i 229 A B C D E F G 230GRtot.i 3.184 3 3.293 3 3 3.293 3.293 231GRind.i 3.184 3 3.293 3 3 3.293 3.293 232 233$ARFT 234 aveR2 frequency 235 0.000 27 236 0.250 48 237 0.333 12 238 0.500 9 239 0.667 6 240 1.000 3 241 242$SCFT 243 SC frequency 244 0.000 48 245 0.250 96 246 0.333 24 247 0.500 12 248 0.667 6 249 1.000 9 250 251attr(,"class") 252[1] "GRind" "list" 253> GRind(oa.design(L18, nlevels=c(2,3,3,3,3,3,3)), arft=FALSE) 254The columns of the array have been used in order of appearance. 255For designs with relatively few columns, 256the properties can sometimes be substantially improved 257using option columns with min3 or even min34. 258 259$GRs 260 GR GRind 261 3 3 262 263$GR.i 264 A B C D E F G 265GRtot.i 3.184 3 3.293 3 3 3.293 3.293 266GRind.i 3.184 3 3.293 3 3 3.293 3.293 267 268$SCFT 269 SC frequency 270 0.000 48 271 0.250 96 272 0.333 24 273 0.500 12 274 0.667 6 275 1.000 9 276 277attr(,"class") 278[1] "GRind" "list" 279> GRind(oa.design(L18, nlevels=c(2,3,3,3,3,3,3)), scft=FALSE) 280The columns of the array have been used in order of appearance. 281For designs with relatively few columns, 282the properties can sometimes be substantially improved 283using option columns with min3 or even min34. 284 285$GRs 286 GR GRind 287 3 3 288 289$GR.i 290 A B C D E F G 291GRtot.i 3.184 3 3.293 3 3 3.293 3.293 292GRind.i 3.184 3 3.293 3 3 3.293 3.293 293 294$ARFT 295 aveR2 frequency 296 0.000 27 297 0.250 48 298 0.333 12 299 0.500 9 300 0.667 6 301 1.000 3 302 303attr(,"class") 304[1] "GRind" "list" 305> GRind(oa.design(L18, nlevels=c(2,3,3,3,3,3,3)), arft=FALSE, scft=FALSE) 306The columns of the array have been used in order of appearance. 307For designs with relatively few columns, 308the properties can sometimes be substantially improved 309using option columns with min3 or even min34. 310 311$GRs 312 GR GRind 313 3 3 314 315$GR.i 316 A B C D E F G 317GRtot.i 3.184 3 3.293 3 3 3.293 3.293 318GRind.i 3.184 3 3.293 3 3 3.293 3.293 319 320attr(,"class") 321[1] "GRind" "list" 322> GRind(oa.design(L18, nlevels=c(2,3,3,3,3,3,3)), arft=FALSE, scft=FALSE, cancor=TRUE) 323The columns of the array have been used in order of appearance. 324For designs with relatively few columns, 325the properties can sometimes be substantially improved 326using option columns with min3 or even min34. 327 328$GRs 329 GR GRind 330 3 3 331 332$GR.i 333 A B C D E F G 334GRtot.i 3.184 3 3.293 3 3 3.293 3.293 335GRind.i 3.184 3 3.293 3 3 3.293 3.293 336 337$cancors 338, , = 1 339 340 others 341factor 1:2 1:3 1:4 1:5 1:6 2:3 2:4 2:5 2:6 3:4 3:5 3:6 4:5 342 A 0 0.000 0.000 0.000 0.000 0.667 0.667 0.00 0.00 0.00 0.667 0.667 0.667 343 B 0 0.000 0.000 0.000 0.000 0.250 0.250 1.00 1.00 1.00 0.250 0.250 0.250 344 C 0 0.333 0.333 0.000 0.000 0.250 0.250 0.50 0.50 0.25 0.250 0.250 0.250 345 D 0 0.333 0.000 0.333 0.333 0.250 1.000 0.25 0.25 0.25 0.250 0.250 0.250 346 E 0 0.333 0.000 0.333 0.333 0.250 1.000 0.25 0.25 0.25 0.250 0.250 0.250 347 F 0 0.000 0.333 0.333 0.000 0.500 0.250 0.25 0.50 0.25 0.250 0.250 0.250 348 G 0 0.000 0.333 0.333 0.000 0.500 0.250 0.25 0.50 0.25 0.250 0.250 0.250 349 others 350factor 4:6 5:6 351 A 0.667 0.00 352 B 0.250 1.00 353 C 0.250 0.25 354 D 0.250 0.25 355 E 0.250 0.25 356 F 0.250 0.25 357 G 0.250 0.25 358 359, , = 2 360 361 others 362factor 1:2 1:3 1:4 1:5 1:6 2:3 2:4 2:5 2:6 3:4 3:5 3:6 4:5 4:6 363 A NA NA NA NA NA NA NA NA NA NA NA NA NA NA 364 B 0 0.000 0.000 0.000 0.000 0.25 0.25 0.00 0.00 1.00 0.25 0.25 0.25 0.25 365 C 0 0.333 0.333 0.000 0.000 0.25 0.25 0.50 0.50 0.25 0.25 0.25 0.25 0.25 366 D 0 0.333 0.000 0.333 0.333 0.25 1.00 0.25 0.25 0.25 0.25 0.25 0.25 0.25 367 E 0 0.333 0.000 0.333 0.333 0.25 1.00 0.25 0.25 0.25 0.25 0.25 0.25 0.25 368 F 0 0.000 0.333 0.333 0.000 0.50 0.25 0.25 0.50 0.25 0.25 0.25 0.25 0.25 369 G 0 0.000 0.333 0.333 0.000 0.50 0.25 0.25 0.50 0.25 0.25 0.25 0.25 0.25 370 others 371factor 5:6 372 A NA 373 B 0.00 374 C 0.25 375 D 0.25 376 E 0.25 377 F 0.25 378 G 0.25 379 380 381attr(,"class") 382[1] "GRind" "list" 383> GRind(oa.design(L18, nlevels=c(2,3,3,3,3,3,3)), cancor=TRUE) 384The columns of the array have been used in order of appearance. 385For designs with relatively few columns, 386the properties can sometimes be substantially improved 387using option columns with min3 or even min34. 388 389$GRs 390 GR GRind 391 3 3 392 393$GR.i 394 A B C D E F G 395GRtot.i 3.184 3 3.293 3 3 3.293 3.293 396GRind.i 3.184 3 3.293 3 3 3.293 3.293 397 398$ARFT 399 aveR2 frequency 400 0.000 27 401 0.250 48 402 0.333 12 403 0.500 9 404 0.667 6 405 1.000 3 406 407$SCFT 408 SC frequency 409 0.000 48 410 0.250 96 411 0.333 24 412 0.500 12 413 0.667 6 414 1.000 9 415 416$cancors 417, , = 1 418 419 others 420factor 1:2 1:3 1:4 1:5 1:6 2:3 2:4 2:5 2:6 3:4 3:5 3:6 4:5 421 A 0 0.000 0.000 0.000 0.000 0.667 0.667 0.00 0.00 0.00 0.667 0.667 0.667 422 B 0 0.000 0.000 0.000 0.000 0.250 0.250 1.00 1.00 1.00 0.250 0.250 0.250 423 C 0 0.333 0.333 0.000 0.000 0.250 0.250 0.50 0.50 0.25 0.250 0.250 0.250 424 D 0 0.333 0.000 0.333 0.333 0.250 1.000 0.25 0.25 0.25 0.250 0.250 0.250 425 E 0 0.333 0.000 0.333 0.333 0.250 1.000 0.25 0.25 0.25 0.250 0.250 0.250 426 F 0 0.000 0.333 0.333 0.000 0.500 0.250 0.25 0.50 0.25 0.250 0.250 0.250 427 G 0 0.000 0.333 0.333 0.000 0.500 0.250 0.25 0.50 0.25 0.250 0.250 0.250 428 others 429factor 4:6 5:6 430 A 0.667 0.00 431 B 0.250 1.00 432 C 0.250 0.25 433 D 0.250 0.25 434 E 0.250 0.25 435 F 0.250 0.25 436 G 0.250 0.25 437 438, , = 2 439 440 others 441factor 1:2 1:3 1:4 1:5 1:6 2:3 2:4 2:5 2:6 3:4 3:5 3:6 4:5 4:6 442 A NA NA NA NA NA NA NA NA NA NA NA NA NA NA 443 B 0 0.000 0.000 0.000 0.000 0.25 0.25 0.00 0.00 1.00 0.25 0.25 0.25 0.25 444 C 0 0.333 0.333 0.000 0.000 0.25 0.25 0.50 0.50 0.25 0.25 0.25 0.25 0.25 445 D 0 0.333 0.000 0.333 0.333 0.25 1.00 0.25 0.25 0.25 0.25 0.25 0.25 0.25 446 E 0 0.333 0.000 0.333 0.333 0.25 1.00 0.25 0.25 0.25 0.25 0.25 0.25 0.25 447 F 0 0.000 0.333 0.333 0.000 0.50 0.25 0.25 0.50 0.25 0.25 0.25 0.25 0.25 448 G 0 0.000 0.333 0.333 0.000 0.50 0.25 0.25 0.50 0.25 0.25 0.25 0.25 0.25 449 others 450factor 5:6 451 A NA 452 B 0.00 453 C 0.25 454 D 0.25 455 E 0.25 456 F 0.25 457 G 0.25 458 459 460attr(,"class") 461[1] "GRind" "list" 462> 463> ## interesting, but take too long 464> #P3.3(oa.design(L36.2.11.3.12, nlevels=c(2,2,2,3,3,3),columns="min3")) 465> #P3.3(oa.design(L36.2.11.3.12, nlevels=c(2,2,2,3,3,3),columns="min3.rela"), rela=TRUE) 466> #P3.3(oa.design(L36.2.11.3.12, nlevels=c(2,2,2,3,3,3),columns="min34")) 467> #P3.3(oa.design(L36.2.11.3.12, nlevels=c(2,2,2,3,3,3),columns="min34.rela"), rela=TRUE) 468> #P3.3(oa.design(L36.2.11.3.12, nlevels=c(2,2,2,3,3,3),columns="minRPFT"), rela=TRUE) 469> #P3.3(oa.design(L36.2.11.3.12, nlevels=c(2,2,2,3,3,3),columns="minRelProjAberr"), rela=TRUE) 470> 471> ## also interesting but also take too long 472> #P3.3(oa.design(L32.2.10.4.7, nlevels=c(2,2,2,4,4,4,4,4), columns="order"), rela=TRUE) 473> #P3.3(oa.design(L32.2.10.4.7, nlevels=c(2,2,2,4,4,4,4,4), columns="min34"), rela=TRUE) 474> #P3.3(oa.design(L32.2.10.4.7, nlevels=c(2,2,2,4,4,4,4,4), columns="min34.rela"), rela=TRUE) 475> #P3.3(oa.design(L32.2.10.4.7, nlevels=c(2,2,2,4,4,4,4,4), columns="minRPFT"), rela=TRUE) 476> #P3.3(oa.design(L32.2.10.4.7, nlevels=c(2,2,2,4,4,4,4,4), columns="minRelProjAberr"), rela=TRUE) 477> 478> proc.time() 479 user system elapsed 480 2.26 0.03 2.29 481