tmpdir<-tempdir()
setwd(tmpdir)
#load data
data(CNOlistToy,package="CellNOptR")
data(ToyModel,package="CellNOptR")
#pre-process model (partial)
indicesToy<-indexFinder(CNOlistToy,ToyModel,verbose=TRUE)
ToyNCNOcutComp<-compressModel(ToyModel,indicesToy)
indicesToyNCNOcutComp<-indexFinder(CNOlistToy,ToyNCNOcutComp)
ToyNCNOcutCompExp<-expandGates(ToyNCNOcutComp)
#optimise
ToyFields4Sim<-prep4sim(ToyNCNOcutCompExp)
initBstring<-rep(1,length(ToyNCNOcutCompExp$reacID))
ToyT1opt<-gaBinaryT1(
CNOlist=CNOlistToy,
model=ToyNCNOcutCompExp,
initBstring=initBstring,
maxGens=2,
popSize=5,
verbose=TRUE)
#write report
namesFilesToy<-list(
dataPlot=NA,
evolFitT1=NA,
evolFitT2=NA,
simResultsT1=NA,
simResultsT2=NA,
scaffold=NA,
scaffoldDot=NA,
tscaffold=NA,
wscaffold=NA,
PKN=NA,
PKNdot=NA,
wPKN=NA,
nPKN=NA)
writeReport(
modelOriginal=ToyModel,
modelOpt=ToyNCNOcutCompExp,
optimResT1=ToyT1opt,
optimResT2=NA,
CNOlist=CNOlistToy,
directory="testToy",
namesFiles=namesFilesToy,
namesData=list(CNOlist="Toy",model="ToyModel"),
resE=NA)
Run the code above in your browser using DataLab