Si tu construit une variable ctemp <-c("PecheAmateurs décembre.proportionMetier", "PechePetitMetierMer mai.proportionMetier","PecheAmateurs mai.proportionMetier", "PecheChalutiere janvier.minInactivityDays","PechePetitMetierEtang octobre.minInactivityDays", "PechePetitMetierEtang décembre.minInactivityDays","FiletMer.standardisationFactor", "PecheAmateurs juin.minInactivityDays","PechePetitMetierEtang novembre.minInactivityDays", "PechePetitMetierMer septembre.minInactivityDays","PechePetitMetierMer avril.minInactivityDays", "PecheAmateurs août.proportionMetier","PecheChalutiere novembre.minInactivityDays", "PechePetitMetierMer février.proportionMetier","PecheChalutiere juillet.minInactivityDays", "PechePetitMetierMer octobre.proportionMetier","CSar Groupe 5.reproductionRate", "PechePetitMetierMer août.proportionMetier","PechePetitMetierMer janvier.minInactivityDays", "PechePetitMetierEtang mai.minInactivityDays","PecheAmateurs octobre.minInactivityDays", "PecheAmateurs septembre.minInactivityDays","PecheChalutiere juin.minInactivityDays", "PecheChalutiere avril.minInactivityDays","PecheAmateurs juillet.proportionMetier", "PecheChalutiere mai.minInactivityDays","PechePetitMetierEtang septembre.minInactivityDays", "PechePetitMetierMer janvier.proportionMetier","PecheAmateurs septembre.proportionMetier", "PechePetitMetierEtang février.minInactivityDays","PechePetitMetierEtang mars.minInactivityDays", "PechePetitMetierEtang juillet.minInactivityDays","PechePetitMetierMer septembre.proportionMetier", "CSar.capturability","PechePetitMetierMer juin.proportionMetier","PecheChalutiere mars.minInactivityDays", "PecheChalutiere septembre.minInactivityDays","PecheAmateurs janvier.minInactivityDays", "PechePetitMetierMer août.minInactivityDays","PechePetitMetierMer juillet.proportionMetier", "PecheChalutiere décembre.minInactivityDays","PecheAmateurs avril.minInactivityDays", "PecheAmateurs mars.proportionMetier","PechePetitMetierMer novembre.minInactivityDays", "PechePetitMetierEtang août.minInactivityDays","CSar.recruitmentDistribution","PecheChalutiere février.minInactivityDays", "CSar.naturalDeathRate","PecheAmateurs octobre.proportionMetier","PechePetitMetierMer mai.minInactivityDays", "PechePetitMetierMer mars.proportionMetier","PecheAmateurs août.minInactivityDays", "FiletEtang.standardisationFactor","PecheChalutiere août.minInactivityDays","PecheAmateurs janvier.proportionMetier", "PecheAmateurs décembre.minInactivityDays","PechePetitMetierEtang avril.minInactivityDays", "PecheAmateurs mars.minInactivityDays","PechePetitMetierMer décembre.minInactivityDays", "PecheAmateurs mai.minInactivityDays","PechePetitMetierEtang juin.minInactivityDays", "Harpon.standardisationFactor","PecheAmateurs novembre.minInactivityDays", "PecheAmateurs novembre.proportionMetier","PechePetitMetierMer mars.minInactivityDays", "Ligneemb.standardisationFactor","PecheAmateurs juin.proportionMetier","PecheAmateurs avril.proportionMetier", "PechePetitMetierMer juillet.minInactivityDays", "PechePetitMetierMer avril.proportionMetier","PechePetitMetierMer octobre.minInactivityDays", "PechePetitMetierMer février.minInactivityDays","PecheAmateurs juillet.minInactivityDays", "PechePetitMetierMer décembre.proportionMetier","PecheAmateurs février.proportionMetier", "PecheChalutiere octobre.minInactivityDays","CSar Groupe 4.reproductionRate","CSar Groupe 7.reproductionRate", "Lignecote.standardisationFactor","CSar.meanWeight","CSar Groupe 6.reproductionRate", "PechePetitMetierMer juin.minInactivityDays","PechePetitMetierMer novembre.proportionMetier", "PechePetitMetierEtang janvier.minInactivityDays","PecheAmateurs février.minInactivityDays","CSar.growth") et que tu lances a<-morris(model=NULL,factors=ctemp, r=10,design=list(type="oat",levels=c(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=c(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=c(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=c(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.0,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.0,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.0,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.0696,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,1680.0,0.11244,1.0,1680.0,2.0,1.0,17.0,20.0,1.0) ) ca semble ne plus poser de probleme steph Jean Couteau <couteau@codelutin.com> a écrit :
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 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 generate a morris object in which there is a data.frame i am interested in. My instruction is :
a<-morris(model=NULL,factors=c("PecheAmateurs décembre.proportionMetier","PechePetitMetierMer mai.proportionMetier","PecheAmateurs mai.proportionMetier","PecheChalutiere janvier.minInactivityDays","PechePetitMetierEtang octobre.minInactivityDays","PechePetitMetierEtang décembre.minInactivityDays","FiletMer.standardisationFactor","PecheAmateurs juin.minInactivityDays","PechePetitMetierEtang novembre.minInactivityDays","PechePetitMetierMer septembre.minInactivityDays","PechePetitMetierMer avril.minInactivityDays","PecheAmateurs août.proportionMetier","PecheChalutiere novembre.minInactivityDays","PechePetitMetierMer février.proportionMetier","PecheChalutiere juillet.minInactivityDays","PechePetitMetierMer octobre.proportionMetier","CSar Groupe 5.reproductionRate","PechePetitMetierMer août.proportionMetier","PechePetitMetierMer janvier.minInactivityDays","PechePetitMetierEtang mai.minInactivityDays","PecheAmateurs octobre.minInactivityDays","PecheAmateurs septembre.minInactivityDays","PecheChalutiere juin.minInactivityDays","PecheChalutiere avril.minInactivityDays","PecheAmateurs juillet.proportionMetier","PecheChalutiere mai.minInactivityDays","PechePetitMetierEtang septembre.minInactivityDays","PechePetitMetierMer janvier.proportionMetier","PecheAmateurs septembre.proportionMetier","PechePetitMetierEtang février.minInactivityDays","PechePetitMetierEtang mars.minInactivityDays","PechePetitMetierEtang juillet.minInactivityDays","PechePetitMetierMer septembre.proportionMetier","CSar.capturability","PechePetitMetierMer juin.proportionMetier","PecheChalutiere mars.minInactivityDays","PecheChalutiere septembre.minInactivityDays","PecheAmateurs janvier.minInactivityDays","PechePetitMetierMer août.minInactivityDays","PechePetitMetierMer juillet.proportionMetier","PecheChalutiere décembre.minInactivityDays","PecheAmateurs avril.minInactivityDays","PecheAmateurs mars.proportionMetier","PechePetitMetierMer novembre.minInactivityDays","PechePetitMetierEtang août.minInactivityDays","CSar.recruitmentDistribution","PecheChalutiere février.minInactivityDays","CSar.naturalDeathRate","PecheAmateurs octobre.proportionMetier","PechePetitMetierMer mai.minInactivityDays","PechePetitMetierMer mars.proportionMetier","PecheAmateurs août.minInactivityDays","FiletEtang.standardisationFactor","PecheChalutiere août.minInactivityDays","PecheAmateurs janvier.proportionMetier","PecheAmateurs décembre.minInactivityDays","PechePetitMetierEtang avril.minInactivityDays","PecheAmateurs mars.minInactivityDays","PechePetitMetierMer décembre.minInactivityDays","PecheAmateurs mai.minInactivityDays","PechePetitMetierEtang juin.minInactivityDays","Harpon.standardisationFactor","PecheAmateurs novembre.minInactivityDays","PecheAmateurs novembre.proportionMetier","PechePetitMetierMer mars.minInactivityDays","Ligneemb.standardisationFactor","PecheAmateurs juin.proportionMetier","PecheAmateurs avril.proportionMetier","PechePetitMetierMer juillet.minInactivityDays","PechePetitMetierMer avril.proportionMetier","PechePetitMetierMer octobre.minInactivityDays","PechePetitMetierMer février.minInactivityDays","PecheAmateurs juillet.minInactivityDays","PechePetitMetierMer décembre.proportionMetier","PecheAmateurs février.proportionMetier","PecheChalutiere octobre.minInactivityDays","CSar Groupe 4.reproductionRate","CSar Groupe 7.reproductionRate","Lignecote.standardisationFactor","CSar.meanWeight","CSar Groupe 6.reproductionRate","PechePetitMetierMer juin.minInactivityDays","PechePetitMetierMer novembre.proportionMetier","PechePetitMetierEtang janvier.minInactivityDays","PecheAmateurs février.minInactivityDays","CSar.growth"),r=10,design=list(type="oat",levels=c(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=c(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=c(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=c(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.0,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.0,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.0,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.0696,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,1680.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 syntax but R asks me to finish my instruction (I need to add two parenthesis at the end so that it can be evaluated, but do not give me the right results).
I replaced the factors vector by its size (the morris method allows that) and the evaluation is ok (with right results).
I tryed also replacing the binf and bsup vectors by a single value (authorized by the morris function also), leaving the factors vector, and the evaluation is also ok (so it is not the factors vector that is creating the problem).
This seems really odd to me and I cannot see where is the problem so I am thinking about a R bug (even if the sensitivity package is not loaded, i cannot enter the instruction, thus it is not a sensitivity 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
_______________________________________________ Isis-fish-devel mailing list Isis-fish-devel@list.isis-fish.org http://list.isis-fish.org/cgi-bin/mailman/listinfo/isis-fish-devel