[R] R: fractional factorial design in R

From: <statistical.model_at_googlemail.com>
Date: Wed 25 Jan 2006 - 08:03:00 EST


sorry, some small mistakes in the previuos syntax. This works!

design.test <- gen.orthogonal.design(c(2,4,3),numCards=16) design.test

gen.orthogonal.design <- function(listFactors,numCards){

	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=FALSE,nTrials=numCards)

	design<-desPB$design#[,2:(numFactors+1)]
	cat("Number of trials: ", fill=T, length(design[,1]), append=T)
	print(cor(design))
	return(design)

}

However, it is necessary to run the function and guess numCards until the correlation matrix is diagonal and all levels are selected for the final design.
Any idea how to solve this problem without an iterative function?

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!



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jan 25 10:22:21 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:10 EST