Re: [R] OptFederov and

From: Michael Griffiths <>
Date: Thu, 10 Mar 2011 11:23:35 +0000

Dear R users,

I have used the AlgDesign package to construct a D-optimal exp. design from a file containing a set of allowed runs. The code for the optFederov call is below:

#Using optFederov ------------------------------



file<- read.table("C:/Upstream/Designs/Vodacom Tanzania phase1 grammar/JMP_allowed.txt",sep="\t",header=TRUE,strip.white=TRUE)

formula<-as.formula("~structure+intro + intro/intro1 + concept + concept/concept1 + prize")

design<-optFederov(formula, data=file, nTrials=48, evaluateI=T, criterion="D", maxIteration=10000, nRepeats=300,



I also use the Rcmdr interface at time for ease, when using the Rcmdr interface is used the same data file, formula and desired number of runs. I also set the seed to that of the optFederov call above. I was expected the same design to be generated. However, when I examine the balance of the generated design, and optFederov yield different designs. Given that is a DoE wrapper for optFerderov I am confused as to why this is happening.

#Using R-commander

Dataset.Dopt <-, data= file, formula,

   nRepeat= 300 ,randomize= TRUE, seed=403)

Apologies if this has a simple answer, I have consulted the R-help archives, but have not found an answer as yet. If the forum could offer reasons as to why this is happening I would be most grateful.


M. Griffiths

        [[alternative HTML version deleted]] mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Thu 10 Mar 2011 - 11:51:32 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 10 Mar 2011 - 12:00:19 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive