[Rd] Bug in R evaluating a huge instruction (PR#14096)

From: <couteau_at_codelutin.com>
Date: Mon, 30 Nov 2009 16:50:10 +0100 (CET)


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



R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Mon 30 Nov 2009 - 22:12:15 GMT

This archive was generated by hypermail 2.2.0 : Mon 30 Nov 2009 - 23:00:53 GMT