J'ai retrouvé ce bout de code dans un vieux script R. Le dataframe sur lequel tu dois travailler est isis.simule dans l'objet isis.methodAnalyse. dataforaov doit etre une variable intermediaire detruite à la la fin de l'analyse. Jean?? temp= as_SOBOL1500_200912071656_0.SensitivityGenitorBiomassRelativeY5.isis.methodAnalyse$isis.simule fmla=as.formula(paste("temp[,10]~(",paste(names(temp)[1:5],collapse="+"),")^2")) test=aov(fmla,data=temp) plot(summary(test)$coefficient[1:5,2]) a+ steph Le 21/05/2012 15:29, Loic GASCHE a écrit :
Bonjour,
En voulant faire d'autres anova que celles faites par ISIS à la suite de mes simulations (pour par exemple tester les interactions sans avoir à tout refaire tourner), j'ai observé dans le RData qu'elles sont appelées de la manière suivante par ISIS :
aovresult<-aov(SensitivityBiomassRelativeY6_Plaice7D~X1+X2+X3+X4+X5+X6+X7+X8+X9+X10+X11+X12+X13+X14+X15+X16+X17+X18+
X19+X20+X21+X22+X23+X24+X25+X26+X27+X28+X29+X30+X31+X32+X33+X34,data=dataforaov)"
La question est d'où sort ce dataforaov ? Et existe-t-il encore dans le RData final ou est-il supprimé à un moment (et remplacé par isis.simule par exemple) ?
Loïc _______________________________________________ Isis-fish-devel mailing list Isis-fish-devel@list.isis-fish.org http://list.isis-fish.org/cgi-bin/mailman/listinfo/isis-fish-devel
-- ...................................................................... Stephanie MAHEVAS (Stephanie.Mahevas@ifremer.fr) IFREMER/EMH (Ecologie et Modèles pour l'Halieutique) Tel: (33) 2 40 37 41 81 Fax: (33) 2 40 37 40 75 o \ o / _ o __| \ / |__ o _ \ o / o /|\ | /\ ___\o \o | o/ o/__ /\ | /|\ / \ / \ | \ /) | ( \ /o\ / ) | (\ / | / \ / \ ......................................................................