data(Selection)
data(network)
#A nv value can chosen using the cutoff function
nv=.11
EGR1<-which(match(Selection@name,"EGR1")==1)
P<-position(network,nv=nv)
#We predict gene expression modulations within the network if EGR1 is experimentaly knocked-out.
prediction_ko5<-predict(Selection,network,nv=nv,targets=EGR1)
#Then we plot the results. Here for example we see changes at time point t2:
plot(prediction_ko5,time=2,ini=P,label_v=Selection@name)
Run the code above in your browser using DataLab