
Your message dated Sun, 22 Oct 2006 12:00:53 +0200
with messageid <[email protected]>
has caused the Debian Bug report #393495,
regarding octaveforge: help of leasqr need improvements
to be marked as having been forwarded to the upstream software
author(s) [email protected]
(NB: If you are a system administrator and have no idea what I am
talking about this indicates a serious mail system misconfiguration
somewhere. Please contact me immediately.)
Debian bug tracking system administrator
(administrator, Debian Bugs database)
 Begin Message 
Subject: 
: Bug#393495: octaveforge: help of leasqr need improvements] 
From: 
Rafael Laboissiere 
Date: 
Sun, 22 Oct 2006 12:00:53 +0200 
package octaveforge
tags 393495 upstream
thanks
I am forwarding below a bug report filed against the octaveforge Debian
package, regarding the documentation for the leasqr function. Please, tell
me if you think this patch is okay. Since I am a member of the octave
project at SF, I can apply the patch to CVS if you wish.
For your convenience, I am attaching below a file containing the patch.
[BTW, Francisco, the next time please attach a standalone file with the
patch. It makes things easier.]
Rafael
 Forwarded message from Francesco Potorti` <[email protected]> 
From: Francesco Potorti` <[email protected]>
Subject: [Pkgoctavedevel] Bug#393495: octaveforge: help of leasqr need
improvements
Date: Mon, 16 Oct 2006 17:49:45 +0200
To: Debian Bug Tracking System <[email protected]>
ReplyTo: Francesco Potorti` <[email protected]>,
[email protected]
Organization: Debian GNU/Linux
Messageid: <[email protected]>
Package: octaveforge
Version: 2006.07.09+dfsg12
Severity: minor
I suggest these improvements to the leasqr help test, which make the
text more consistent and readable and add some information.
 leasqr.m~ 20061016 17:28:25.000000000 +0200
+++ leasqr.m 20061016 17:48:43.000000000 +0200
@@ 22,24 +22,30 @@ function [f,p,kvg,iter,corp,covp,covr,st
% LevenbergMarquardt nonlinear regression of f(x,p) to y(x).
%
% Version 3.beta
% {}= optional parameters
% x=vec or mat of indep variables, 1 row/observation: x=[x0 x1....xm]
% y=vec of obs values, same no. of rows as x.
% wt=vec(dim=length(x)) of statistical weights. These should be set
% to be proportional to (sqrt of var(y))^1; (That is, the covariance
% matrix of the data is assumed to be proportional to diagonal with diagonal
% equal to (wt.^2)^1. The constant of proportionality will be estimated.),
% default=ones(length(y),1).
% pin=vector of initial parameters to be adjusted by leasqr.
% dp=fractional incr of p for numerical partials,default= .001*ones(size(pin))
% dp(j)>0 means central differences.
% dp(j)<0 means onesided differences.
% Note: dp(j)=0 holds p(j) fixed i.e. leasqr wont change initial guess: pin(j)
% F=name of function in quotes,of the form y=f(x,p)
% dFdp=name of partials Mfile in quotes default is prt=dfdp(x,f,p,dp,F)
% stol=scalar tolerances on fractional improvement in ss,default stol=.0001
% niter=scalar max no. of iterations, default = 20
% options=matrix of n rows (same number of rows as pin) containing
+% Optional parameters are in braces {}.
+% x = column vector or matrix of independent variables, 1 row per
+% observation: x = [x0 x1....xm].
+% y = column vector of observed values, same number of rows as x.
+% wt = column vector (dim=length(x)) of statistical weights. These
+% should be set to be proportional to (sqrt of var(y))^1; (That is,
+% the covariance matrix of the data is assumed to be proportional to
+% diagonal with diagonal equal to (wt.^2)^1. The constant of
+% proportionality will be estimated.); default = ones(length(y),1).
+% pin = column vec of initial parameters to be adjusted by leasqr.
+% dp = fractional increment of p for numerical partial derivatives;
+% default = .001*ones(size(pin))
+% dp(j) > 0 means central differences on jth parameter p(j).
+% dp(j) < 0 means onesided differences on jth parameter p(j).
+% dp(j) = 0 holds p(j) fixed i.e. leasqr wont change initial guess: pin(j)
+% F = name of function in quotes; the function shall be of the form y=f(x,p),
+% with y, x, p of the form y, x, pin as described above.
+% dFdp = name of partial derivative function in quotes; default is "dfdp", a
+% slow but general partial derivatives function; the function shall be
+% of the form prt=dfdp(x,f,p,dp,F) (see dfdp.m).
+% stol = scalar tolerance on fractional improvement in scalar sum of
+% squares = sum((wt.*(yf))^2); default stol = .0001;
+% niter = scalar maximum number of iterations; default = 20;
+% options = matrix of n rows (same number of rows as pin) containing
% column 1: desired fractional precision in parameter estimates.
% Iterations are terminated if change in parameter vector (chg) on two
% consecutive iterations is less than their corresponding elements
@@ 52,16 +58,16 @@ function [f,p,kvg,iter,corp,covp,covr,st
% default = Inf*ones().
%
% OUTPUT VARIABLES
% f=vec function values computed in function func.
% p=vec trial or final parameters. i.e, the solution.
% kvg=scalar: =1 if convergence, =0 otherwise.
% iter=scalar no. of interations used.
% corp= correlation matrix for parameters
% covp= covariance matrix of the parameters
% covr = diag(covariance matrix of the residuals)
% stdresid= standardized residuals
% Z= matrix that defines confidence region
% r2= coefficient of multiple determination
+% f = column vector of values computed: f = F(x,p).
+% p = column vector trial or final parameters. i.e, the solution.
+% kvg = scalar: = 1 if convergence, = 0 otherwise.
+% iter = scalar number of iterations used.
+% corp = correlation matrix for parameters.
+% covp = covariance matrix of the parameters.
+% covr = diag(covariance matrix of the residuals).
+% stdresid = standardized residuals.
+% Z = matrix that defines confidence region (see comments in the source).
+% r2 = coefficient of multiple determination.
%
% All Zero guesses not acceptable
@@ 113,7 +119,7 @@ function [f,p,kvg,iter,corp,covp,covr,st
% 8) Change to more efficent algorithm of Bard for selecting epsL.
% 9) Tighten up memory usage by making use of sparse matrices (if
% MATLAB version >= 4.0) in computation of covp, corp, stdresid.
% Modified by Francesco Potorti
+% Modified by Francesco Potortì
% for use in Octave
%
% References:
 System Information:
Debian Release: testing/unstable
APT prefers testing
APT policy: (990, 'testing'), (500, 'testingproposedupdates')
Architecture: amd64 (x86_64)
Shell: /bin/sh linked to /bin/bash
Kernel: Linux 2.6.151amd64k8smp
Locale: LANG=C, [email protected] (charmap=ANSI_X3.41968) (ignored: LC_ALL
set to C)
Versions of packages octaveforge depends on:
ii octave2.1forge 2006.03.17+dfsg11 Contributed functions from the GNU
octaveforge recommends no packages.
 no debconf information
_______________________________________________
Pkgoctavedevel mailing list
[email protected]
http://lists.alioth.debian.org/mailman/listinfo/pkgoctavedevel
 End forwarded message 

Rafael
leasqr.mpatch
Description: Text document
 End Message 

