cmyriam | Bonjour,
je travail sous linux et plus particulierement avec un logiciel appele ms qui permet de simuler des milliers de donnee biologique. Meme si vous ne le connaissez pas, je pense que vous pouvez quand meme m'aider avec mon pb.
Voila dans mon shell, je rentre cette commande pour lancer une simulation:
Code :
- > ./ms 68 10000 -t 45.94 -I 3 23 8 37 -n 2 1.14 -g 2 25098 -n 3 3.09 -g 3 1304876 -es 0.000005298 3 0.85 -ej 0.000005298 3 2 -ej 0.000005298 4 1 -ej 0.00028 2 1 -en 0.019171 1 0.0000163 -en 0.019175 1 0.4134 -r 300 2000 | ./statsMS_2pop 23 8 3 > Output.csv
|
Avec ./statsMS_2pop 23 8 3 qui est un script qui me permet de calculer des stats a partir des donnees obrenu grace a la simulation.
Mon probleme est que je voudrais relancer a la fin de chaque simulation une autre automatiquement tout en changeant certaines donnees. Je veux changer a chaque fois des valeurs (celle entre "" ) de "68", -t "45.95", -I 3 "23 8 37", -r "300" 2000.
Voila le script que j'ai ecrit:
Code :
- t <- c(45.94,54.6,28.45);
- nSample <- c(68,57,34);
- nsubPop <- matrix(c(3,23,8,37,12,14,3,10,19,23,14,16),nrow=3,ncol=4,byrow=TRUE)
- r <- c(300,157,278)
- for(iCount in 1:3)
- {
- system(paste("./ms",nSample[iCount],"10000 -t",t[iCount],"-I",nsubPop[iCount,],"-n 2 1.14 -g 2 25098 -n 3 3.09 -g 3 1304876 -es 0.000005298 3 0.85 -ej 0.000005298 3 2 -ej 0.000005298 4 1 -ej 0.00028 2 1 -en 0.019171 1 0.0000163 -en 0.019175 1 0.4134 -r",
- r[iCount],"| ./statsMS_2pop",nsubPop[iCount,2],nsubPop[iCount,3],"3 >>Results.csv" ))
- }
|
Helas il ne fonctionne pas, de plus n'y aurait-il pas un moyen plus facile de faire ca comme peut etre recuperer les valeurs a partir d'un autre fichier.
Enfin derniere question, devrais-je faire le fichier en bash?
Merci par avance pour l'aide,
Myriam
|