[R] Linear Models with positive coefficients?

Douglas Bates bates at stat.wisc.edu
Thu Jul 8 15:14:32 CEST 1999


John Logsdon <j.logsdon at lancaster.ac.uk> writes:

> Just a thought.  Otherwise a constrained optimisation algorithm, which is
> not available in R.  Or parameterise so that the log of the parameter is
> used, as in the gnlr suite for the shape parameter.

There are simpler ways of going about the task because it is
essentially an integer programming problem.  You are just searching
for the subset of the coefficients that should be set equal to zero.

Dennis Wolf and I wrote about it several years ago including some
Fortran code.

@article{bate:wolf:1984,
    Author = {Bates, Douglas M. and Wolf, Dennis A.},
    Title = {Non-negative Regression By {G}ivens Rotations},
    Year = 1984,
    Journal = {Communications in Statistics, Part B -- Simulation and
              Computation},
    Volume = 13,
    Pages = {841--850}
}

Given the number of things I am supposed to be working on, I will not
be able to make an R implementation of this soon.  If anyone wants to
try I would be happy to offer what advice I can.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list