#data<-gen.factorial(levels=3,nVars=3,varNames=c("A","B","C"))
#d<-matrix(c("SpecialOps","DieselEngines","Screws","1","Reconn","GasTurbines","Waterjets","2","NA","NA","NA","3","NA","NA","NA","4","NA","NA","NA","5","NA","NA","NA","6","NA","NA","NA","7","NA","NA","NA","8"), nrow=8, ncol=4, byrow=TRUE, dimnames=list(c("row1", "row2", "row3", "row4", "row5", "row6", "row7", "row8"),c("Payload", "PropulsionPlant", "PropellerType", "VesselDetectability")))

rm(list=ls())
d<-gen.factorial(levels=c(2,2,2,8),nVars=4,factors="all",varNames=c("Payload","PropulsionPlant","PropellerType","VesselDetectability"))
optFederov(~.,d,approximate=FALSE,criterion="I",augment=FALSE,args=TRUE)