training.nl: Performs the HVDM training step and returns a list containing...
In rHVDM: Hidden Variable Dynamic Modeling

Description Usage Arguments Details Value Author(s) References See Also Examples

This method performs the nonlinear training step of the HVDM algorithm. It returns a list that will then be used in the subsequent screening step.

1	training.nl(inputHVDM,transforms,constraints,forcetransforms,genemodels,firstguess)

`inputHVDM`	an HVDM training object
`transforms`	a vector containing the kinetic parameter identifiers that have to be transformed during optimisation (optional)
`constraints`	"known" values for the activator signal
`forcetransforms`	Boolean, whether the transformation in argument transforms have to be applied
`genemodels`	the type of model that has to be used for each gene
`firstguess`	a first guess for all the parameters

Contrary to the linear training function (without .nl suffix), this function takes as main input another training object (either a linear training object or a non-linear one, where all the genes have been fitted witha MIchelis-Menten model ("MM").

An exponential transform is set by default for the basal (Bj) and degradation (Dj) rates, as well as for the kinetic parameters (Vj and Kj) in the production function (through the transforms argument). This forces the values for these parameters to be positive For the exponent Nj, the expp1 function ensures it is greater than one, when the hill formulation is used. To turn this off, set the forcetransforms switch to FALSE. Even in this case the degradation rate will not be allowed to take non-positive values as it causes problems with the differential operator used internally. The value in the vector indicates the parameter to be transformed: "Bj": basal rate of transcription, "Sj": sensitivity, "Dj": degrdation rate, etc.. The entry label indicates the transform to be applied.

The constraints argument is used to specify values for know values of some of the parameters, typically the strength of the activator for some time points. See example below for the syntax.

The genemodels argument is compulsory and used to specify the model used for each individual gene. MM means Michaelis-Menten model whereas hill means a hill function is used. The general form of the production function is Bj+Vj*f(t)^Nj/(Kj^Nj+f(t)^Nj) where Nj=1 for the MM model and Nj>1 in the hill case. See below for the syntax of the input.

The firstguess argument is not in use yet and will offer the possibility to enter a first guess for the fitting.

a list containing the results.

Martino Barenco

M. Barenco, D. Tomescu, D. Brewer, R. Callard, J. Stark, M. Hubank (2006) Ranked predictions of p53 targets using Hidden Variable Dynamic Modelling. Genome Biology, V7(3), R25.

HVDMcheck,screening.nl,fitgene,HVDMreport,training

#load data and fit a linear model
data(HVDMexample)
rm(fiveGyMAS5)
data(HVDMexample2)
tp532<-training(eset=twodosesMAS5,genes=p53traingenes,degrate=0.8)
#formulate constraints
CONSTRAINTS<-c(350,35)
names(CONSTRAINTS)<-c("trfact1.5Gy.1.4","trfact1.hGy.1.4")
#specify individual gene models
GENEMODELS<-rep("MM",5)
GENEMODELS[c(1,3)]<-"hill"
names(GENEMODELS)<-p53traingenes
#run the model
tp53hyb<-training.nl(inputHVDM=tp532,constraints=CONSTRAINTS,genemodels=GENEMODELS)