>If an orthogonal main effect plan exists for the number of trials you
specify, optFederov() in AlgDesign will more than likely find it for
you, since such a design should be an optimal design.
thanks very much to you and the others!!
I have written a little function, and now the syntax to obtain an orthogonal
design such as 2x4x3x2 is:
design.test <- gen.orthogonal.design(c(2,4,3,2))
gen.orthogonal.design <- function(listFactors){
library(AlgDesign)
FactorsNames<-c("A","B","C","D","E","F","G","H","J","K","L")
numFactors<-length(listFactors)
dat<-gen.factorial(listFactors,center=FALSE,varNames=FactorsNames[1:numFacto
rs])
desPB<-optFederov(~.,dat,nRepeats=20,approximate=TRUE)
design<-desPB$design[,2:numFactors]
cat("Minimum number of trials: ", fill=T, length(design[,1]),
append=T)
#cor(design)
return(design)
}
Roberto Furlan
University of Turin, Italy
----------------------------------------
La mia Cartella di Posta in Arrivo รจ protetta con SPAMfighter
188 messaggi contenenti spam sono stati bloccati con successo.
Scarica gratuitamente SPAMfighter!
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
