similar to: glmpath and coxpath variables

Hi, I am new to model selection by coefficient shrinkage method such as lasso. And I became particularly interested in variable selection in Cox regression by lasso. I became aware of the coxpath() in R package glmpath does lasso on Cox model. I have tried the sample script on the help page of coxpath(), but I have difficult time understanding the output. Therefore, I would greatly appreciate if

We have uploaded to CRAN a new version of glmpath, a package which fits the L1 regularization path for generalized linear models. The revision includes: - coxpath, a function for fitting the L1-regularization path for the Cox ph model; - bootstrap functions for analyzing sparse solutions; - the ability to mix in L2 regularization along with L1 (elasticnet). We have also completed a report that

We have uploaded to CRAN a new version of glmpath, a package which fits the L1 regularization path for generalized linear models. The revision includes: - coxpath, a function for fitting the L1-regularization path for the Cox ph model; - bootstrap functions for analyzing sparse solutions; - the ability to mix in L2 regularization along with L1 (elasticnet). We have also completed a report that

Hi all, I'm relatively new to using R, and have been trying to fit an L1 regularization path using coxpath from the glmpath library. I'm interested in using a cross validation framework, where I crossvalidate on a training set to select the lambda that achieves the lowest error, then use that value of lambda on the entire training set, before applying to a test set. This seems to entail

The problem first appeared in R 2.6.1 and is still there in R 2.6.2 On Windows running R CMD check command for glmpath package fails. The reason seems to be that when R is running the examples file (glmpath-Ex.R), it skips about 50 lines and as a result gives a syntax error. I'm working with a modified version of the CRAN glmpath 0.94. My version happens to give a more clear example of a

Hello R helpers: *My first message didn't pass trough filter so here it's again* I would like to obtain probability of an event for one single patient as a function of time (from survfit.coxph) object, as I want to find what is the probability of an event say at 1 month and what is the probability of an event at 80 months and compare. So I tried the following but it fails miserably. I

Hello: I would like to obtain probability of an event for one single patient as a function of time (from survfit.coxph) object, as I want to find what is the probability of an event say at 1 month and what is the probability of an event at 80 months and compare. So I tried the following but it fails miserably. I looked at some old posts but could not figure out the solution. Here's what I did

All, I was having trouble trying to create a new class of data and pass it on to splom (in the lattice library). I mentioned this to Martin Morgan after a talk he gave. Following is not so much a question, but rather an answer from Morgan that might be useful to others. Here is the edited part of an email conversation with him: On Thursday, August 26, 2010 1:36 PM, Martin Morgan wrote:

Dear colleagues, I hope someone can help me with my problem. I have fitted a cox model with the following syntax: # cox01def <-cph(Surv(TEVENT,EVENT) ~ ifelse(AGE>50, (AGE-50)^2,0) + BMI + # HDL+DIABETES +HISTCAR2 + log(CREAT)+ as.factor(ALBUMIN)+STENOSIS+IMT,data # = XC, x=T, y=T, surv=T) *1 Furthermore I have estimated my beta's also with a Lasso method - Coxpath ( from

Dear R users, Please can you help me with a relatively straightforward problem that I am struggling with? I am simply trying to plot a baseline survivor and hazard function for a simple data set of lung cancer survival where `futime' is follow up time in months and status is 1=dead and 0=alive. Using the survival package: lung.wbs <- survreg( Surv(futime, status)~ 1, data=lung,

Dear List, I have a theoretical question related to epidemiological data analysis: If the treatment status (tx = 0,1) changes over time for the patients in a non-randomized cohort, is there a way to estimate the treatment effect? (i.e., after joining the study, some patients may have to wait for a period of time before receiving the treatment, i.e., the situation of patient with id == 2 for the

Hello everyone, These are some questions about the 'pec' function in R. These questions deal with prediction error curves and their derivation. Prediction error curves are documented in, for example, "Efron-type measures of prediction error for survival analysis" by Gerds and Schumacher. I have detailed some syntax that I have used at the bottom of this email. The associated

Hi, perhaps you can help me to find out, how to find the best Lambda in a LASSO-model. I have a feature selection problem with 150 proteins potentially predicting Cancer or Noncancer. With a lasso model fit.glm <- glmpath(x=as.matrix(X), y=target, family="binomial") (target is 0, 1 <- Cancer non cancer, X the proteins, numerical in expression), I get following path (PICTURE

I'm trying to do LASSO in R with the package glmpath. However, I'm not sure if I am using the accompanying prediction function *predict.glmpath()* correctly. Suppose I fit some regularized binomial regression model like so: library(glmpath);load(heart.data);attach(heart.data); fit <- glmpath(x, y, family=binomial) Then I can use predict.glmpath() to estimate the value of the

We have uploaded to CRAN the first version of glmpath, which fits the L1 regularization path for generalized linear models. The lars package fits the entire piecewise-linear L1 regularization path for the lasso. The coefficient paths for L1 regularized glms, however, are not piecewise linear. glmpath uses convex optimization - in particular predictor-corrector methods- to fit the

We have uploaded to CRAN the first version of glmpath, which fits the L1 regularization path for generalized linear models. The lars package fits the entire piecewise-linear L1 regularization path for the lasso. The coefficient paths for L1 regularized glms, however, are not piecewise linear. glmpath uses convex optimization - in particular predictor-corrector methods- to fit the

Hi R-users! I am trying to learn how to use the glmpath package. I have a dataframe like this > dim(data) [1] 605 109 and selected the following > response <- data[,1] > features<-as.matrix(data[,3:109]) > mymodel <- glmpath(features,response, family = binomial) Error in if (lambda <= min.lambda) { : missing value where TRUE/FALSE expected Reading the glmpath pdf, I

Hi all, I am using glampath package for L1 regularized logistic regression. I have read the article " L1 regularization path algorithm for GLM" by park and Hastie (2006). One thing I can't understand that how to find best lambda for my prediction. I want to use that lambda for the prediction not the entire set. thanks. -- View this message in context:

First, here is your message as it appears on R-help. On 10/14/2012 05:00 AM, r-help-request@r-project.org wrote: > I?m trying to set up proportional hazard model that is stratified with > respect to covariate 1 and has an interaction between covariate 1 and > another variable, covariate 2. Both variables are categorical. In the > following, I try to illustrate the two problems that

Hi Claire, I'm replying and CC-ing to the R-help list to get more eyes on your question since others will likely have more/better advice, and perhaps someone else in the future will have a similar question, and might find this thread handy. I've removed your specific research aim since that might be private information, but you can include that later if others find it necessary to know