bugGNU Octave - Bugs: bug #43699, [optim] Error in curvefit_stat...

 
 

bug #43699: [optim] Error in curvefit_stat when requesting covp with multiple correlated parameters

Submitted by:  Markus Appel <mappel>
Submitted on:  Tue 25 Nov 2014 04:14:29 PM UTC  
 
Category: Octave Forge PackageSeverity: 3 - Normal
Priority: 5 - NormalItem Group: Segfault, Bus Error, etc.
Status: FixedAssigned to: None
Originator Name: Open/Closed: Closed
Release: 3.8.0Operating System: GNU/Linux

Add a New Comment(Rich Markup)
   

You are not logged in

Please log in, so followups can be emailed to you.

 

Fri 28 Nov 2014 01:47:55 PM UTC, comment #2:

Sorry, I've forgotten to commit the last change before pushing. It is now in a changeset with some additional changes addressing the remaining issues I mentioned.

Olaf Till <i7tiol>
Project Member
Wed 26 Nov 2014 02:29:15 PM UTC, comment #1:

Thanks for the report. The triggering error was just usage of a vector in the wrong orientation, but there turned out to be an additional problem with usage of error bounds for basis vectors of more than one-dimensional null spaces. A fix for the bug is now pushed; tell me if you need a summarizing patch for the time before a new package release.

Though this bug should be fixed, there seem to be still some issues with the procedure computing covd and corp, but I'll need a bit time for this.

Olaf Till <i7tiol>
Project Member
Tue 25 Nov 2014 04:14:29 PM UTC, original submission:

When running the attached minimal test script, I get the following error from the curvefit_stat call:
with optim 1.4.0 in octave 3.8.0 on Ubuntu 12.04.
In the script, I have added dummy parameters on which the model function does not depend to provoke the error. In reality I encountered this error when multiple parameters of a model function became fully correlated.

From playing around it seems that _null_optim_ cannot handle an input matrix with more than one singular value being zero (within tol). Commenting out line 20: "PIN.dummy2 = 1;" in the demonstration script makes it work and results in the expected correlation matrix INFO.covp with one row/column being NA.

Markus Appel <mappel>

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach File(s):
   
   
Comment:
   

Attached Files
file #32542:  curvefittest.m added by mappel (742B - text/x-objcsrc)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by i7tiol (Posted a comment)
  • -unavailable- added by mappel (Submitted the item)
  • -unavailable- added by mappel
  •  

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only project members can vote.

     

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 4 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Fri 28 Nov 2014 02:36:44 PM UTCcarandraugStatusNone=>Fixed
      Open/ClosedOpen=>Closed
    Tue 25 Nov 2014 04:14:29 PM UTCmappelAttached File-=>Added curvefittest.m, #32542
      Carbon-Copy-=>Added -unavailable-

    Back to the top


    Powered by Savane 3.1-cleanup1