Description Usage Arguments Details Value Author(s) References See Also Examples
View source: R/userfcnnonlinearHVDM.R
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
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | #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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.