Return-Path: X-Original-To: apmail-commons-dev-archive@www.apache.org Delivered-To: apmail-commons-dev-archive@www.apache.org Received: from mail.apache.org (hermes.apache.org [140.211.11.3]) by minotaur.apache.org (Postfix) with SMTP id B428865F5 for ; Thu, 21 Jul 2011 05:41:39 +0000 (UTC) Received: (qmail 868 invoked by uid 500); 21 Jul 2011 05:41:36 -0000 Delivered-To: apmail-commons-dev-archive@commons.apache.org Received: (qmail 240 invoked by uid 500); 21 Jul 2011 05:41:23 -0000 Mailing-List: contact dev-help@commons.apache.org; run by ezmlm Precedence: bulk List-Help: List-Unsubscribe: List-Post: List-Id: Reply-To: "Commons Developers List" Delivered-To: mailing list dev@commons.apache.org Received: (qmail 231 invoked by uid 99); 21 Jul 2011 05:41:16 -0000 Received: from nike.apache.org (HELO nike.apache.org) (192.87.106.230) by apache.org (qpsmtpd/0.29) with ESMTP; Thu, 21 Jul 2011 05:41:16 +0000 X-ASF-Spam-Status: No, hits=-0.7 required=5.0 tests=FREEMAIL_FROM,RCVD_IN_DNSWL_LOW,SPF_PASS X-Spam-Check-By: apache.org Received-SPF: pass (nike.apache.org: domain of phil.steitz@gmail.com designates 209.85.210.41 as permitted sender) Received: from [209.85.210.41] (HELO mail-pz0-f41.google.com) (209.85.210.41) by apache.org (qpsmtpd/0.29) with ESMTP; Thu, 21 Jul 2011 05:41:06 +0000 Received: by pzk4 with SMTP id 4so1388051pzk.0 for ; Wed, 20 Jul 2011 22:40:44 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=message-id:date:from:user-agent:mime-version:to:subject:references :in-reply-to:content-type:content-transfer-encoding; bh=byxzo1EtfI94i82+UUOrn4myAZ2oMjS81LhmdFiNE4s=; b=etFiVPHcM7cy1Db+md1qOIIQDFUL+v8YY1EoX/pwJHnmzcyJDV5hmn6jHsK1rgET9q Dgke8fuFaEPKLLhrb4nlpmEgS1hk9HUYqLn4q4xhQJKQCDKYDaLiTrXv7DV0JFnPKUhI lMhaqRu7lo9oYvBmO5uwsPh6vlM3SWMsSwlD0= Received: by 10.142.65.24 with SMTP id n24mr4294168wfa.281.1311226844621; Wed, 20 Jul 2011 22:40:44 -0700 (PDT) Received: from a.local (71-223-74-208.phnx.qwest.net [71.223.74.208]) by mx.google.com with ESMTPS id e15sm472836wfd.3.2011.07.20.22.40.42 (version=SSLv3 cipher=OTHER); Wed, 20 Jul 2011 22:40:43 -0700 (PDT) Message-ID: <4E27BBD9.9060805@gmail.com> Date: Wed, 20 Jul 2011 22:40:41 -0700 From: Phil Steitz User-Agent: Mozilla/5.0 (Macintosh; Intel Mac OS X 10.6; rv:5.0) Gecko/20110624 Thunderbird/5.0 MIME-Version: 1.0 To: Commons Developers List Subject: Re: [math] Pivoting in QR decomposition References: In-Reply-To: Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: 7bit X-Virus-Checked: Checked by ClamAV on apache.org On 7/20/11 9:27 PM, Greg Sterijevski wrote: > Hello, > > In attempting to build constrained ols into the OLSMultipleLinearRegression > class, I rediscovered the fact that to keep the regression coefficients in > their canonical order, I need to know how the QR decomposition pivoted. In > other words, while I get the correct constrained parameters, their order is > changed. Is there any way to get this info out (pivot order)? Can you explain a little more exactly what you are trying to do? How exactly is the parameter vector getting permuted? The QRDecomp solver certainly should not do that. See the QRDecomp sources for what is exposed in the public API. We can talk about adding more, but we should be careful about exposing more implementation details. You can now get the Householder reflection vectors and obviously the decomp itself. As you can see from the sources, the current OLS impl just uses the solver. Phil > > Thanks, > > -Greg > --------------------------------------------------------------------- To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org For additional commands, e-mail: dev-help@commons.apache.org