From izmirlig at mail.nih.gov  Fri Nov  4 01:32:37 2005
From: izmirlig at mail.nih.gov (izmirlig@mail.nih.gov)
Date: Fri,  4 Nov 2005 01:32:37 +0100 (CET)
Subject: [Rd] small bug in gl1ce, package lasso2 (PR#8280)
Message-ID: <20051104003237.A00B8238E5@slim.kubism.ku.dk>

Full_Name: Grant Izmirlian
Version: 2.2.0
OS: SuSe Linux version 9.2
Submission from: (NULL) (156.40.34.177)


Sorry about the last submission, my bug-fix had an error in it because ifelse 
doesn't vectorize.  I'll repost with the correct bug-fix.
-------------------------------------------------------------------------------
The option exists to include all parameters, including the intercept, in the L-1
constraint, by specifying the argument, sweep.out=NULL, explicitly in the
call. However, upon doing this, the function stops with an error report 
"Matrix build from transformed variables has a constant column"

I was able to fix it as follows. In the following line, taken from the same
block of code producing the stop error, modify the existing line from

      X.to.C.stds <- sqrt(apply(X.to.C.w, 2, var))

to 


      X.to.C.stds <- (!is.null(sweep.out)) *
                     c(1,sqrt(apply(X.to.C.w[,-is.null(sweep.out)], 2, var))) +

                     (is.null(sweep.out)) * sqrt(apply(X.to.C.w,2,var)))

That should work.



[Rd] small bug in gl1ce, package lasso2 (PR#8280)
izmirlig@mail.nih.gov