From couteau at codelutin.com  Mon Nov 30 16:50:10 2009
From: couteau at codelutin.com (couteau at codelutin.com)
Date: Mon, 30 Nov 2009 16:50:10 +0100 (CET)
Subject: [Rd] Bug in R evaluating a huge instruction (PR#14096)
Message-ID: <20091130155010.4EFF1282EFE3@mail.pubhealth.ku.dk>

Hello,

I do not know if i am facing a bug or a mistake (and if it is a bug i do =

not know in which category I can put it on the bug tracker), so I post=20
my message here.

I am trying to use the Morris function from the sensitivity package (but =

that is not the problem). I created a (big) instruction that should=20
generate a morris object in which there is a data.frame i am interested=20
in. My instruction is :

a<-morris(model=3DNULL,factors=3Dc("PecheAmateurs=20
d=E9cembre.proportionMetier","PechePetitMetierMer=20
mai.proportionMetier","PecheAmateurs=20
mai.proportionMetier","PecheChalutiere=20
janvier.minInactivityDays","PechePetitMetierEtang=20
octobre.minInactivityDays","PechePetitMetierEtang=20
d=E9cembre.minInactivityDays","FiletMer.standardisationFactor","PecheAmat=
eurs=20
juin.minInactivityDays","PechePetitMetierEtang=20
novembre.minInactivityDays","PechePetitMetierMer=20
septembre.minInactivityDays","PechePetitMetierMer=20
avril.minInactivityDays","PecheAmateurs=20
ao=FBt.proportionMetier","PecheChalutiere=20
novembre.minInactivityDays","PechePetitMetierMer=20
f=E9vrier.proportionMetier","PecheChalutiere=20
juillet.minInactivityDays","PechePetitMetierMer=20
octobre.proportionMetier","CSar Groupe=20
5.reproductionRate","PechePetitMetierMer=20
ao=FBt.proportionMetier","PechePetitMetierMer=20
janvier.minInactivityDays","PechePetitMetierEtang=20
mai.minInactivityDays","PecheAmateurs=20
octobre.minInactivityDays","PecheAmateurs=20
septembre.minInactivityDays","PecheChalutiere=20
juin.minInactivityDays","PecheChalutiere=20
avril.minInactivityDays","PecheAmateurs=20
juillet.proportionMetier","PecheChalutiere=20
mai.minInactivityDays","PechePetitMetierEtang=20
septembre.minInactivityDays","PechePetitMetierMer=20
janvier.proportionMetier","PecheAmateurs=20
septembre.proportionMetier","PechePetitMetierEtang=20
f=E9vrier.minInactivityDays","PechePetitMetierEtang=20
mars.minInactivityDays","PechePetitMetierEtang=20
juillet.minInactivityDays","PechePetitMetierMer=20
septembre.proportionMetier","CSar.capturability","PechePetitMetierMer=20
juin.proportionMetier","PecheChalutiere=20
mars.minInactivityDays","PecheChalutiere=20
septembre.minInactivityDays","PecheAmateurs=20
janvier.minInactivityDays","PechePetitMetierMer=20
ao=FBt.minInactivityDays","PechePetitMetierMer=20
juillet.proportionMetier","PecheChalutiere=20
d=E9cembre.minInactivityDays","PecheAmateurs=20
avril.minInactivityDays","PecheAmateurs=20
mars.proportionMetier","PechePetitMetierMer=20
novembre.minInactivityDays","PechePetitMetierEtang=20
ao=FBt.minInactivityDays","CSar.recruitmentDistribution","PecheChalutiere=
=20
f=E9vrier.minInactivityDays","CSar.naturalDeathRate","PecheAmateurs=20
octobre.proportionMetier","PechePetitMetierMer=20
mai.minInactivityDays","PechePetitMetierMer=20
mars.proportionMetier","PecheAmateurs=20
ao=FBt.minInactivityDays","FiletEtang.standardisationFactor","PecheChalut=
iere=20
ao=FBt.minInactivityDays","PecheAmateurs=20
janvier.proportionMetier","PecheAmateurs=20
d=E9cembre.minInactivityDays","PechePetitMetierEtang=20
avril.minInactivityDays","PecheAmateurs=20
mars.minInactivityDays","PechePetitMetierMer=20
d=E9cembre.minInactivityDays","PecheAmateurs=20
mai.minInactivityDays","PechePetitMetierEtang=20
juin.minInactivityDays","Harpon.standardisationFactor","PecheAmateurs=20
novembre.minInactivityDays","PecheAmateurs=20
novembre.proportionMetier","PechePetitMetierMer=20
mars.minInactivityDays","Ligneemb.standardisationFactor","PecheAmateurs=20
juin.proportionMetier","PecheAmateurs=20
avril.proportionMetier","PechePetitMetierMer=20
juillet.minInactivityDays","PechePetitMetierMer=20
avril.proportionMetier","PechePetitMetierMer=20
octobre.minInactivityDays","PechePetitMetierMer=20
f=E9vrier.minInactivityDays","PecheAmateurs=20
juillet.minInactivityDays","PechePetitMetierMer=20
d=E9cembre.proportionMetier","PecheAmateurs=20
f=E9vrier.proportionMetier","PecheChalutiere=20
octobre.minInactivityDays","CSar Groupe 4.reproductionRate","CSar Groupe =

7.reproductionRate","Lignecote.standardisationFactor","CSar.meanWeight","=
CSar=20
Groupe 6.reproductionRate","PechePetitMetierMer=20
juin.minInactivityDays","PechePetitMetierMer=20
novembre.proportionMetier","PechePetitMetierEtang=20
janvier.minInactivityDays","PecheAmateurs=20
f=E9vrier.minInactivityDays","CSar.growth"),r=3D10,design=3Dlist(type=3D"=
oat",levels=3Dc(5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5=
,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,=
5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5),grid.jump=3Dc(2,2,2,2,2,2,2,2,=
2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2=
,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,=
2,2,2,2,2)),binf=3Dc(0.0,0.0,0.0,11.0,10.0,12.0,0.06,10.0,11.0,0.0,0.0,0.=
0,10.0,0.0,11.0,0.0,1120.0,0.0,10.0,22.0,12.0,10.0,10.0,11.0,0.0,11.0,15.=
0,0.0,0.0,23.0,16.0,22.0,0.0,0.0,0.0,11.0,10.0,18.0,0.0,0.0,11.0,11.0,0.0=
,3.0,22.0,0.0,8.0,0.0,0.0,0.0,0.0,8.0,0.1456,11.0,0.0,18.0,17.0,11.0,7.0,=
11.0,20.0,0.0464,16.0,0.0,1.0,0.1515,0.0,0.0,0.0,0.0,1.0,11.0,9.0,0.0,0.0=
,11.0,1120.0,1120.0,0.0749,0.0,1120.0,0.0,0.0,13.0,16.0,0.0),bsup=3Dc(1.0=
,1.0,1.0,15.0,14.0,16.0,0.09,14.0,15.0,2.0,2.0,1.0,14.0,1.0,15.0,1.0,1680=
=2E0,1.0,14.0,26.0,16.0,14.0,14.0,14.0,1.0,15.0,19.0,1.0,1.0,27.0,20.0,26=
=2E0,1.0,1.0,1.0,15.0,14.0,22.0,2.0,1.0,15.0,15.0,1.0,7.0,26.0,1.0,12.0,1=
=2E0,1.0,2.0,1.0,12.0,0.2184,15.0,1.0,22.0,21.0,15.0,11.0,15.0,24.0,0.069=
6,20.0,1.0,5.0,0.22,1.0,1.0,2.0,1.0,5.0,15.0,13.0,1.0,1.0,15.0,1680.0,168=
0.0,0.11244,1.0,1680.0,2.0,1.0,17.0,20.0,1.0))

I checked several times and this intruction seams to have the right=20
syntax but R asks me to finish my instruction (I need to add two=20
parenthesis at the end so that it can be evaluated, but do not give me=20
the right results).

I replaced the factors vector by its size (the morris method allows=20
that) and the evaluation is ok (with right results).

I tryed also replacing the binf and bsup vectors by a single value=20
(authorized by the morris function also), leaving the factors vector,=20
and the evaluation is also ok (so it is not the factors vector that is=20
creating the problem).

This seems really odd to me and I cannot see where is the problem so I=20
am thinking about a R bug (even if the sensitivity package is not=20
loaded, i cannot enter the instruction, thus it is not a sensitivity=20
problem).

Did somebody already encountered such a bug (if it is a bug) ?

I am using R 2.9.2 under linux/Xubuntu 9.10 and this has also been tried =

on R 2.9.0 on Windows 7.

I can provide more info if needed.

Jean Couteau
Code Lutin



[Rd] Bug in R evaluating a huge instruction (PR#14096)
couteau at codelutin.com