<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Hello,<br>
<br>
You can use the numderivative macro, with your residual function
as first argument.<br>
<br>
hth<br>
<br>
S.<br>
<br>
<br>
Le 17/09/2015 16:18, Pablo Caron a écrit :<br>
</div>
<blockquote
cite="mid:CAF_hC28eS1R9yrVEu5p250A9VSb0H3uLxkXyiampL-fQ67f+zQ@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_quote">
<div dir="ltr">Dear user group,
<div><br>
</div>
<div>I need to extract the Jacobian from the lsqrsolve
function. The fortran routine lmdif (<a
moz-do-not-send="true"
href="http://www.netlib.org/minpack/lmdif.f"
target="_blank"><a class="moz-txt-link-freetext" href="http://www.netlib.org/minpack/lmdif.f">http://www.netlib.org/minpack/lmdif.f</a></a>)
provide this as the output fjac, but lsqrsolve does not
have this option.</div>
<div><br>
</div>
<div>I need to evaluate the goodness of a regression I made
using lsqrsolve and I found that this is done through the
correlation matrix. I need either the correlation matrix
or the jacobian.<br>
</div>
<div><br>
</div>
<div>I would like to get an output similar to the following
<div><font face="monospace, monospace"><br>
</font></div>
<div><font face="monospace, monospace">Final set of
parameters Asymptotic Standard Error</font></div>
<div><font face="monospace, monospace">=======================
==========================</font></div>
<div><font face="monospace, monospace">m = 0.174844
+/- 0.0002964 (0.1695%)</font></div>
<div><font face="monospace, monospace">b = 1506.34
+/- 5.756 (0.3821%)</font></div>
<div><font face="monospace, monospace"><br>
</font></div>
<div><font face="monospace, monospace">correlation matrix
of the fit parameters:</font></div>
<div><font face="monospace, monospace"> m b</font></div>
<div><font face="monospace, monospace">m 1.000</font></div>
<div><font face="monospace, monospace">b -0.448 1.000</font></div>
</div>
<div><font face="monospace, monospace"><br>
</font></div>
Regards<span class="HOEnZb"><font color="#888888"><br>
<br>
Pablo</font></span></div>
</div>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:users@lists.scilab.org">users@lists.scilab.org</a>
<a class="moz-txt-link-freetext" href="http://lists.scilab.org/mailman/listinfo/users">http://lists.scilab.org/mailman/listinfo/users</a>
</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
Département de Génie Informatique
EA 4297 Transformations Intégrées de la Matière Renouvelable
Université de Technologie de Compiègne - CS 60319
60203 Compiègne cedex</pre>
</body>
</html>