Re: [bug #13259] full analysis-high Te values -- March 26, 2009 - 15:01

There's a system test for this, so I'll play around with this.  We
need to work out if removing this statement will cause problems
elsewhere.  Note I've made changes to the 1.3 line already which
modifies the behaviour of model selection and model elimination in
model-free analysis a little.

Regards,

Edward


On Thu, Mar 26, 2009 at 2:36 PM, Sébastien Morin
<sebastien.morin.1@xxxxxxxxx> wrote:
> Hi,
>
> I commented the if statement for the 'final' data pipe name:
>
>        # Model elimination.
>        #if modsel_pipe != 'final':
>        #    eliminate()
>        eliminate()
>
> This as the same effect as changing the modsel_pipe argument from
> 'final' to 'aic' at line 303: model 'm8' is eliminated and model 'm1' is
> chosen.
>
> What is the best option ?
> Remove the if statement at line 602 or change the data pipe name at line
> 303 ?
>
> Regards,
>
>
> Séb  :)
>
>
>
>
>
> Edward d'Auvergne wrote:
> > I've been trying to debug it, and the script may need to be modified.
> > If you look at the self.model_selection() method in full_analysis.py,
> > there is an 'if' statement for the 'final' data pipe name.  If you
> > remove this statement, does the bug go away?
> >
> > Cheers,
> >
> > Edward
> >
> >
> > On Thu, Mar 26, 2009 at 4:13 AM, Sébastien Morin
> > <sebastien.morin.1@xxxxxxxxx> wrote:
> >
> > > Sébastien Morin wrote:
> > >
> > > > Hi Ed,
> > > >
> > > > I just tried something using the files in the bug report, as well as with
> > > > the up-to-date 1.3 repository (r9001)...
> > > >
> > > > I modified line 303 of the script from:
> > > >
> > > >                    self.model_selection(modsel_pipe='final',
> > > > dir=self.base_dir + 'aic')
> > > >
> > > > to:
> > > >
> > > >                    self.model_selection(modsel_pipe='aic',
> > > > dir=self.base_dir + 'aic')
> > > >
> > > > With this change, model 'm8' is eliminated and model 'm1' is selected
> > > > instead.
> > > >
> > > > Is it possible that model-free model elimination was not performed during
> > > > normal rounds of minimisation (with modsel_pipe='aic') ?
> > > Sorry, correct the sentence above with:
> > >
> > > ---
> > > Is it possible that model-free model elimination was not performed during
> > > normal rounds of minimisation (with  line 303 including the argument:
> > > modsel_pipe='final') ?
> > > ---
> > >
> > > Séb
> > >
> > > > I was sure I had tried this change yesterday with Pierre-Yves (on his 
> > > > full
> > > > datasets). However, probably something was wrong or maybe the changes you
> > > > made today were also needed. Anyway...
> > > >
> > > > What do you think ?
> > > >
> > > > Regards,
> > > >
> > > >
> > > > Séb  :)
> > > >
> > > >
> > > >
> > > >
> > > > Edward d Auvergne wrote:
> > > >
> > > >
> > > > > Update of bug #13259 (project relax):
> > > > >
> > > > >                  Status:                    None => Confirmed
> > > > >               Assigned to:                    None => bugman
> > > > >        Operating System:               GNU/Linux => All systems
> > > > >
> > > > >    _______________________________________________________
> > > > >
> > > > > Reply to this item at:
> > > > >
> > > > >  <http://gna.org/bugs/?13259>
> > > > >
> > > > > _______________________________________________
> > > > >  Message sent via/by Gna!
> > > > >  http://gna.org/
> > > > >
> > > > >
> > > > > _______________________________________________
> > > > > relax (http://nmr-relax.com)
> > > > >
> > > > > This is the relax-devel mailing list
> > > > > relax-devel@xxxxxxx
> > > > >
> > > > > To unsubscribe from this list, get a password
> > > > > reminder, or change your subscription options,
> > > > > visit the list information page at
> > > > > https://mail.gna.org/listinfo/relax-devel
> > > --
> > > Sébastien Morin
> > > PhD Student
> > > S. Gagné NMR Laboratory
> > > Université Laval & PROTEO
> > > Québec, Canada
>
>
> --
> Sébastien Morin
> PhD Student
> S. Gagné NMR Laboratory
> Université Laval & PROTEO
> Québec, Canada
>
>
>
> _______________________________________________
> relax (http://nmr-relax.com)
>
> This is the relax-devel mailing list
> relax-devel@xxxxxxx
>
> To unsubscribe from this list, get a password
> reminder, or change your subscription options,
> visit the list information page at
> https://mail.gna.org/listinfo/relax-devel