plot.clogitL1 {clogitL1} | R Documentation |
Takes a clogitL1
object and plots the parameter profile associated with it.
## S3 method for class 'clogitL1' plot(x, logX=T, add.legend=F, add.labels=T, lty=1:ncol(x$beta), col=1:ncol(x$beta), ...)
x |
an object of type |
logX |
should the horizontal axis be on log scale? |
add.legend |
set to TRUE if legend should be printed in top right hand corner. Legend will contain names of variables in data.frame, if specified, otherwise will be numbered from 1 to p in order encountered in original input matrix x |
add.labels |
set to TRUE if labels are to be added to curves at leftmost side. If variable names are available, these are plotted, otherwise, curves are numbered from 1 to p in order encountered in original input matrix x |
lty |
usual 'lty' plotting parameter. |
col |
usual 'col' plotting parameter. |
... |
additional arguments to |
http://www.jstatsoft.org/v58/i12/
set.seed(145) # data parameters K = 10 # number of strata n = 5 # number in strata m = 2 # cases per stratum p = 20 # predictors # generate data y = rep(c(rep(1, m), rep(0, n-m)), K) X = matrix (rnorm(K*n*p, 0, 1), ncol = p) # pure noise strata = sort(rep(1:K, n)) par(mfrow = c(1,2)) # fit the conditional logistic model clObj = clogitL1(y=y, x=X, strata) plot(clObj, logX=TRUE) # cross validation clcvObj = cv.clogitL1(clObj) plot(clcvObj)