Learn R Programming

MXM (version 0.8.7)

model: Regression model(s) obtained from SES

Description

One or more regression models obtained from SES are returned.

Usage

model(target, dataset, sesObject, nsignat = 1, test = NULL)

Arguments

target
A numeric vector containing the values of the target variable. If the values are proportions or percentages, i.e. strictly within 0 and 1 they are mapped into R using log( target/(1 - target) ). It can also discrete data, binary data (as factor), nomi
dataset
A numeric matrix containing the variables. Rows are samples and columns are features. If you have categorical variables, this should be a data frame.
sesObject
An object with the results of a SES run.
nsignat
How many signatures to use. If nsignat = 1 (default value) the first set of variables will be used for the model. If you want more, then specify the nubmer of signatures you want. If you want the models based on all signatures, specify "all". If you p
test
If you know the test used in SES put it here, otherwise leave it NULL. It will take this information from the SEs object. If you used a robust version of a test (wherever possible), robust model(s) will be created.

Value

  • A list including:
  • modDepending on the number of signatures requested, one or models will be returned.
  • signatureA matrix (or just one vector if one signature only) with the variables of each signature, along with the BIC of the corresponding regression model.

Details

This command is useful if you want to see all models and check for example their fitting ability, MSE in linear models for exmaple.

References

Aitchison J. (1986). The Statistical Analysis of Compositional Data, Chapman & Hall; reprinted in 2003, with additional material, by The Blackburn Press. Cox D.R. (1972). Regression models and life-tables. J. R. Stat. Soc., 34, 187-220. Draper, N.R. and Smith H. (1988). Applied regression analysis. New York, Wiley, 3rd edition. Ferrari S.L.P. and Cribari-Neto F. (2004). Beta Regression for Modelling Rates and Proportions. Journal of Applied Statistics, 31(7): 799-815. Gutenbrunner C., Jureckova J., Koenker R. and Portnoy S. (1993). Tests of Linear Hypothesis based on Regression Rank Scores, Journal of NonParametric Statistics 2, 307-331. Joseph M.H. (2011). Negative Binomial Regression. Cambridge University Press, 2nd edition. Koenker R.W. (2005). Quantile Regression, Cambridge University Press. Lagani V., Kortas G. and Tsamardinos I. (2013). Biomarker signature identification in "omics" with multiclass outcome. Computational and Structural Biotechnology Journal, 6(7): 1-7. Lagani V. and Tsamardinos I. (2010). Structure-based variable selection for survival data. Bioinformatics Journal 16(15): 1887-1894. Lambert, Diane (1992). Zero-inflated Poisson regression, with an application to defects in manufacturing. Technometrics 34(1)1: 1-14. Mardia K.V., Kent J.T. and Bibby J.M. (1979). Multivariate Analysis, Academic Press, New York, USA. Maronna R.D. Yohai M.V. (2006). Robust Statistics, Theory and Methods. Wiley. McCullagh P., and Nelder J.A. (1989). Generalized linear models. CRC press, USA, 2nd edition.

See Also

SES, cv.ses

Examples

Run this code
# simulate a dataset with continuous data
dataset <- matrix(runif(1000 * 101, 1, 100), nrow = 1000 ) 
#the target feature is the last column of the dataset as a vector
target <- dataset[, 100]
dataset <- dataset[, -100]
sesObject <- SES(target , dataset , max_k=3 , threshold = 0.05);
model(target, dataset, sesObject, nsignat = 1, test = NULL) 
model(target, dataset, sesObject, nsignat = 40, test = NULL)

Run the code above in your browser using DataLab