Author: jcouteau Date: 2009-04-29 08:39:07 +0000 (Wed, 29 Apr 2009) New Revision: 126 Modified: trunk/sensitivity/SensitivityCalculatorROptimumLHS.java trunk/sensitivity/SensitivityCalculatorRRandomLHS.java Log: Adaptation des m?\195?\169thodes Hypercube latin au nouveau CDC. Possibilit?\195?\169 de modifier les instructions R. Pop-up d'erreur si facteurs non corrects. Modified: trunk/sensitivity/SensitivityCalculatorROptimumLHS.java =================================================================== --- trunk/sensitivity/SensitivityCalculatorROptimumLHS.java 2009-04-28 13:43:11 UTC (rev 125) +++ trunk/sensitivity/SensitivityCalculatorROptimumLHS.java 2009-04-29 08:39:07 UTC (rev 126) @@ -19,9 +19,14 @@ package sensitivity; import java.io.File; -import java.io.Serializable; import java.util.List; +import javax.swing.Box; +import javax.swing.JLabel; +import javax.swing.JOptionPane; +import javax.swing.JScrollPane; +import javax.swing.JTextPane; + import org.apache.commons.logging.Log; import org.apache.commons.logging.LogFactory; import org.codelutin.j2r.REngine; @@ -31,38 +36,42 @@ import org.codelutin.util.FileUtil; import fr.ifremer.isisfish.datastore.SimulationStorage; +import fr.ifremer.isisfish.simulator.sensitivity.AbstractSensitivityCalculator; import fr.ifremer.isisfish.simulator.sensitivity.DesignPlan; +import fr.ifremer.isisfish.simulator.sensitivity.Domain; import fr.ifremer.isisfish.simulator.sensitivity.Factor; import fr.ifremer.isisfish.simulator.sensitivity.Scenario; -import fr.ifremer.isisfish.simulator.sensitivity.SensitivityCalculator; -import fr.ifremer.isisfish.simulator.sensitivity.AbstractSensitivityCalculator; import fr.ifremer.isisfish.simulator.sensitivity.SensitivityException; import fr.ifremer.isisfish.simulator.sensitivity.SensitivityScenarios; import fr.ifremer.isisfish.simulator.sensitivity.domain.ContinuousDomain; +import fr.ifremer.isisfish.simulator.sensitivity.domain.DiscreteDomain; import fr.ifremer.isisfish.simulator.sensitivity.domain.EquationContinuousDomain; import fr.ifremer.isisfish.simulator.sensitivity.domain.MatrixContinuousDomain; import fr.ifremer.isisfish.util.Doc; -/**save.image() - * Implementation of Optimum Latin Hypercube method using R. +/** + * save.image() Implementation of Optimum Latin Hypercube method using R. * * @author jcouteau * @version $Revision: 1.0 $ * */ @Doc("Optimum Latin Hypercube") -public class SensitivityCalculatorROptimumLHS extends AbstractSensitivityCalculator { +public class SensitivityCalculatorROptimumLHS extends + AbstractSensitivityCalculator { /** to use log facility, just put in your code: log.info("..."); */ static private Log log = LogFactory .getLog(SensitivityCalculatorROptimumLHS.class); @Doc("Number of simulations (default=10)") - public int param_simulationNumber; + public int param_simulationNumber = 10; @Doc("The maximum number of times the Columnwise Pairwise algorithm is applied to all the columns(default=2).") - public int param_MaxSweeps; + public int param_MaxSweeps = 2; @Doc("The optimal stopping criterion (between 0 and 1) (default=0.1).") - public double param_eps; + public double param_eps = 0.1; + @Doc("True to be able to modify the code sent to R") + public boolean param_modifR = false; /** * Retourne vrai si le calculateur sait gerer la cardinalité des facteurs @@ -77,12 +86,28 @@ public SensitivityScenarios compute(DesignPlan plan, File outputdirectory) throws SensitivityException { + setIsisFactorsR(plan, outputdirectory); + int factornumber = plan.getFactors().size(); List<Factor> factors = plan.getFactors(); double[] dataframe = new double[0]; SensitivityScenarios thisExperiment = new SensitivityScenarios(); List<Scenario> thisExperimentScenarios = thisExperiment.getScenarios(); + //Test all factors, if one is discrete, return null + for (int i = 0; i < factornumber; i++) { + if (factors.get(i).getDomain() instanceof DiscreteDomain) { + JOptionPane + .showMessageDialog( + null, + "Error", + factors.get(i).getName() + + " has a discrete domain, this is not acceptable for this method.", + JOptionPane.ERROR_MESSAGE); + return null; + } + } + REngine engine = new RProxy(); if ((Integer) param_simulationNumber == 0) { @@ -98,14 +123,46 @@ try { + //Set working directory to get Isis R session + log.info("setwd(\"" + outputdirectory.getParent() + "\")"); + engine.voidEval("setwd(\"" + outputdirectory.getParent() + "\")"); + + //Get Isis R session + log.info("load(\".RData\")"); + engine.voidEval("load(\".RData\")"); + //Load the lhs library engine.voidEval("library(lhs)"); log.info("Message sent to R : " + "library(lhs)"); - String rInstruction = "x<-optimumLHS(n=" + param_simulationNumber - + ",k=" + factornumber + ",maxSweeps=" + param_MaxSweeps - + ",eps=" + param_eps + ")"; + //Remove all objects that will be used + engine.voidEval("remove(\"isis.factor.distribution\")"); + engine.voidEval("remove(\"isis.simule\")"); + engine.voidEval("remove(\"isis.MethodExp\")"); + engine.voidEval("remove(\"isis.methodAnalyse\")"); + String rInstruction = "x<-optimumLHS(n=%s,k=%s,maxSweeps=%s,eps=%s)"; + + String rCall = String.format(rInstruction, param_simulationNumber, + factornumber, param_MaxSweeps, param_eps); + + if (param_modifR) { + JLabel label = new JLabel( + "Modifier le code R envoyé si vous le souhaitez"); + JTextPane text = new JTextPane(); + text.setText(rCall); + text.setSize(400, 400); + text.setPreferredSize(text.getSize()); + + Box box = Box.createVerticalBox(); + box.add(label); + box.add(new JScrollPane(text)); + + JOptionPane.showMessageDialog(null, box, "R modif", + JOptionPane.QUESTION_MESSAGE); + rCall = text.getText(); + } + //Create the scenarios engine.voidEval(rInstruction); log.info("Message sent to R : " + rInstruction); @@ -122,8 +179,7 @@ for (int j = 0; j < param_simulationNumber; j++) { Scenario experimentScenario = new Scenario(); for (int i = 0; i < factornumber; i++) { - Factor factor = plan.getFactors() - .get(i); + Factor factor = plan.getFactors().get(i); if ((factor.getDomain() instanceof MatrixContinuousDomain) || (factor.getDomain() instanceof EquationContinuousDomain)) { factor @@ -145,13 +201,78 @@ thisExperiment.setScenarios(thisExperimentScenarios); } - //Create the factors vectors - for (int j = 0; j < factornumber; j++) { - Factor factor = thisExperimentScenarios - .get(0).getFactors().get(j); + String isisFactorDistribution = "isis.factor.distribution<-data.frame(NomFacteur=c(%s),NomDistribution=c(%s),ParametreDistribution=c(%s))"; + + String distribution = ""; + String parameters = ""; + String factorNames = ""; + + for (int i = 0; i < factors.size(); i++) { + Domain domain = factors.get(i).getDomain(); + if (i != 0) { + distribution += ","; + factorNames += ","; + parameters += ","; + } + + factorNames += "\"" + factors.get(i).getName() + "\""; + distribution += "\"qunif\""; + parameters += "\"[" + ((ContinuousDomain) domain).getMinBound() + + ";" + ((ContinuousDomain) domain).getMaxBound() + + "]\""; + } + + log.info("Message sent to R : " + + String.format(isisFactorDistribution, factorNames, + distribution, parameters)); + engine.voidEval(String.format(isisFactorDistribution, factorNames, + distribution, parameters)); + + log.info("Message sent to R : " + "call<-" + "\"x<-optimumLHS(n=" + + param_simulationNumber + ",k=" + factornumber + + ",maxSweeps=" + param_MaxSweeps + ",eps=" + param_eps + + ")\""); + engine.voidEval("call<-" + "\"x<-optimumLHS(n=" + + param_simulationNumber + ",k=" + factornumber + + ",maxSweeps=" + param_MaxSweeps + ",eps=" + param_eps + + ")\""); + + log + .info("Message sent to R : " + + "isis.MethodExp<-list(\"isis.factors\"=isis.factors,\"isis.factor.distribution\"=isis.factor.distribution,\"call\"=call)"); + engine + .voidEval("isis.MethodExp<-list(\"isis.factors\"=isis.factors,\"isis.factor.distribution\"=isis.factor.distribution,\"call\"=call)"); + + log + .info("Message sent to R : " + + "attr(isis.MethodExp,\"nomModel\")<-\"isis-fish-externe-R\""); + engine + .voidEval("attr(isis.MethodExp,\"nomModel\")<-\"isis-fish-externe-R\""); + + log.info("Message sent to R : " + "isis.simule<-data.frame(x)"); + engine.voidEval("isis.simule<-data.frame(x)"); + + log + .info("Message sent to R : " + + "attr(isis.simule,\"nomModel\")<-\"isis-fish-externe-R\""); + engine + .voidEval("attr(isis.simule,\"nomModel\")<-\"isis-fish-externe-R\""); + + log.info("Message sent to R : " + + "names(isis.simule)<-isis.factors[[1]]"); + engine.voidEval("names(isis.simule)<-isis.factors[[1]]"); + + String data = "data<-data.frame("; + + //Create the factors vectors and the dataframe instruction + for (int j = 0; j < factors.size(); j++) { + Factor factor = thisExperimentScenarios.get(0).getFactors() + .get(j); if (factor.getDomain() instanceof EquationContinuousDomain) { - String vector = "factor" + j + "<-c("; + //String vector = "factor" + j + "<-c("; + String vector = factor.getName().replaceAll(" ", "") + + "<-c("; for (int i = 0; i < param_simulationNumber; i++) { if (i < (param_simulationNumber - 1)) { vector = vector @@ -170,7 +291,9 @@ engine.voidEval(vector); log.info("Message sent to R : " + vector); } else if (factor.getDomain() instanceof MatrixContinuousDomain) { - String vector = "factor" + j + "<-c("; + //String vector = "factor" + j + "<-c("; + String vector = factor.getName().replaceAll(" ", "") + + "<-c("; for (int i = 0; i < param_simulationNumber; i++) { if (i < (param_simulationNumber - 1)) { vector = vector @@ -189,7 +312,9 @@ engine.voidEval(vector); log.info("Message sent to R : " + vector); } else { - String vector = "factor" + j + "<-c("; + //String vector = "factor" + j + "<-c("; + String vector = factor.getName().replaceAll(" ", "") + + "<-c("; for (int i = 0; i < param_simulationNumber; i++) { if (i < (param_simulationNumber - 1)) { vector = vector @@ -207,15 +332,13 @@ engine.voidEval(vector); log.info("Message sent to R : " + vector); } - } - //Create the data data.frame from the factors - String data = "data<-data.frame("; - for (int j = 0; j < factornumber; j++) { - if (j < factornumber - 1) { - data = data + "factor" + j + "=factor(factor" + j + "),"; + if (j < factors.size() - 1) { + data += factor.getName().replaceAll(" ", "") + "=factor(" + + factor.getName().replaceAll(" ", "") + "),"; } else { - data = data + "factor" + j + "=factor(factor" + j + "))"; + data += factor.getName().replaceAll(" ", "") + "=factor(" + + factor.getName().replaceAll(" ", "") + "))"; } } @@ -237,7 +360,7 @@ rInstruction += ")"; engine.voidEval(rInstruction); - log.info("Message sent to R : " + rInstruction); + log.info("Message sent to R" + rInstruction); // Set output directory engine.voidEval("setwd(\"" + outputdirectory.getAbsolutePath() @@ -245,7 +368,7 @@ log.info("Message sent to R : " + "setwd(\"" + outputdirectory.getAbsolutePath() + "\")"); - // Export the scenario matrix for the second run in a .randomlhs.csv file + // Export the scenario matrix for the second run in a .optimumlhs.csv file engine.voidEval("write.csv(data,file=\".optimumlhs.csv\")"); log.info("Message sent to R : " + "write.csv(data,file=\".optimumlhs.csv\")"); @@ -255,6 +378,14 @@ log.info("Message sent to R" + "dput(factornames,file=\".factornames\")"); + //Set working directory to save Isis R session + log.info("setwd(\"" + outputdirectory.getParent() + "\")"); + engine.voidEval("setwd(\"" + outputdirectory.getParent() + "\")"); + + // Save Isis R session + log.info("save.image()"); + engine.voidEval("save.image()"); + engine.terminate(); } catch (Exception e) { @@ -282,36 +413,48 @@ REngine engine = new RProxy(); try { - int sensitivityNumber = simulationStorages.get(0).getParameter() - .getSensitivityExport().size(); + //Set working directory to get Isis R session + log.info("setwd(\"" + outputdirectory.getParent() + "\")"); + engine.voidEval("setwd(\"" + outputdirectory.getParent() + "\")"); - for (int k = 0; k < sensitivityNumber; k++) { + //Get Isis R session + log.info("load(\".RData\")"); + engine.voidEval("load(\".RData\")"); - // Set output directory - engine.voidEval("setwd(\"" + outputdirectory.getAbsolutePath() - + "\")"); - log.info("Message sent to R : setwd(\"" - + outputdirectory.getAbsolutePath() + "\")"); + // Set output directory + engine.voidEval("setwd(\"" + outputdirectory.getAbsolutePath() + + "\")"); + log.info("Message sent to R : setwd(\"" + + outputdirectory.getAbsolutePath() + "\")"); - //Get back the scenarios - //engine.voidEval("factors<-read.csv(\".expandgrid.csv\",row.names=1,col.names=1)"); - engine.voidEval("factors<-read.csv(\".optimumlhs.csv\")"); - log - .info("Message sent to R : factors<-read.csv(\".optimumlhs.csv\")"); + //Get back the scenarios + engine.voidEval("factors<-read.csv(\".optimumlhs.csv\")"); + log + .info("Message sent to R : factors<-read.csv(\".optimumlhs.csv\")"); - // Get the factornames object in the .factornames file - engine.voidEval("factornames<-dget(\".factornames\")"); - log.info("Message sent to R : " - + "factornames<-dget(\".factornames\")"); + // Get the factornames object in the .factornames file + engine.voidEval("factornames<-dget(\".factornames\")"); + log.info("Message sent to R : " + + "factornames<-dget(\".factornames\")"); - //Get back the factors number - int factorNumber = ((Double) engine.eval("length(factors)-1")) - .intValue(); - //factorNumber=factorNumber-1; + //Get back the factors number + int factorNumber = ((Double) engine.eval("length(factors)-1")) + .intValue(); - //Create the results vectors - String result = "result<-c("; - for (int l = 0; l < simulationStorages.size(); l++) { + //Get back the simulation number + param_simulationNumber = (Integer) (engine + .eval("length(factors[,1])")); + + int sensitivityNumber = simulationStorages.get(0).getParameter() + .getSensitivityExport().size(); + + for (int k = 0; k < sensitivityNumber; k++) { + + // Creates the R expression to import results in R + String name = simulationStorages.get(0).getParameter() + .getSensitivityExport().get(k).getExportFilename(); + String rInstruction = name + "<-c("; + for (int l = 0; l < param_simulationNumber; l++) { File importFile = new File(simulationStorages.get(l) .getDirectory().toString() + File.separator @@ -324,28 +467,65 @@ .getExtensionFilename()); String simulResult = FileUtil.readAsString(importFile); double simulationResult = Double.valueOf(simulResult); - if (l < simulationStorages.size() - 1) { - result = result + simulationResult + ","; + if (l < param_simulationNumber - 1) { + rInstruction = rInstruction + simulationResult + ","; } else { - result = result + simulationResult; + rInstruction = rInstruction + simulationResult; } } - result = result + ")"; - engine.voidEval(result); - log.info("Message sent to R : " + result); + rInstruction = rInstruction + ")"; + log.info("Message sent to R : " + rInstruction); + + // Send the simulation results + engine.voidEval(rInstruction); + + //Put results in isis.simule + engine.voidEval("isis.simule<-data.frame(isis.simule," + name + + ")"); + log.info("Message sent to R : " + + "isis.simule<-data.frame(isis.simule," + name + ")"); + } + + //adding attribute to isis.Simule + log + .info("Message sent to R : " + + "attr(isis.simule,\"nomModel\")<-\"isis-fish-externe-R\""); + engine + .voidEval("attr(isis.simule,\"nomModel\")<-\"isis-fish-externe-R\""); + + log.info("Message sent to R : " + + "attr(isis.simule,\"call\")<-isis.MethodExp$call"); + engine.voidEval("attr(isis.simule,\"call\")<-isis.MethodExp$call"); + + for (int k = 0; k < sensitivityNumber; k++) { + + // Creates the R expression to import results in R + String name = simulationStorages.get(0).getParameter() + .getSensitivityExport().get(k).getExportFilename(); + //Create the dataforaov data.frame - String dataframe = "dataforaov<-data.frame(factors,result=result)"; + String dataframe = "dataforaov<-data.frame(factors," + name + + "=" + name + ")"; engine.voidEval(dataframe); log.info("Message sent to R : " + dataframe); //Call aov() - String aovCall = "aovresult<-aov(result~"; + String aovCall = "aovresult<-aov(" + name + "~"; for (int j = 0; j < factorNumber; j++) { + log.info("Message sent to R : " + "names(factors)[" + + (j + 2) + "]"); + if (j < (factorNumber - 1)) { - aovCall = aovCall + "factor" + j + "+"; + aovCall = aovCall + + engine + .eval("names(factors)[" + (j + 2) + "]") + + "+"; } else { - aovCall = aovCall + "factor" + j + ",data=dataforaov)"; + aovCall = aovCall + + engine + .eval("names(factors)[" + (j + 2) + "]") + + ",data=dataforaov)"; } } engine.voidEval(aovCall); @@ -388,14 +568,6 @@ log.info("Message sent to R : " + "row.names(exportsensitivity)<-factornames"); - /*Set the export directory - *Export directory is the first simulation export directory. - */ - engine.voidEval("setwd(\"" + outputdirectory.getAbsolutePath() - + "\")"); - log.info("Message sent to R : setwd(\"" - + outputdirectory.getAbsolutePath() + "\")"); - //Set dataforaov names engine .voidEval("resultsnames<-c(\"Simulation\",factornames,\"Result\")"); @@ -406,6 +578,14 @@ log.info("Message sent to R : " + "names(dataforaov)<-resultsnames"); + /*Set the export directory + *Export directory is the first simulation export directory. + */ + engine.voidEval("setwd(\"" + outputdirectory.getAbsolutePath() + + "\")"); + log.info("Message sent to R : setwd(\"" + + outputdirectory.getAbsolutePath() + "\")"); + //Save the results with the scenarios. engine.voidEval("write.csv(dataforaov,\"" + simulationStorages.get(0).getParameter() @@ -428,10 +608,71 @@ .getExportFilename() + "_SensitivityIndices.csv\")"); //FIXME export through java to enable export when using Rserve (when distant Rserve). - engine.terminate(); + + //creating isis.methodAnalyse + log + .info("Message sent to R : " + + "isis.methodAnalyse<-list(\"isis.factors\"=isis.factors,\"isis.factor.distribution\"=isis.factor.distribution,\"isis.simule\"=isis.simule,call_method=\"" + + aovCall + "\"" + ",aovresult)"); + engine + .voidEval("isis.methodAnalyse<-list(\"isis.factors\"=isis.factors,\"isis.factor.distribution\"=isis.factor.distribution,\"isis.simule\"=isis.simule,call_method=\"" + + aovCall + "\"" + ",aovresult)"); + + log + .info("Message sent to R : " + + "attr(isis.methodAnalyse,\"nomModel\")<-\"isis-fish-externe-R\""); + engine + .voidEval("attr(isis.methodAnalyse,\"nomModel\")<-\"isis-fish-externe-R\""); + + String renameIsisMethodAnalyse = "%s.isis.methodAnalyse<-isis.methodAnalyse"; + String simulationName = simulationStorages.get(0).getName() + .replaceAll("-", ""); + log.info("Message sent to R : " + + String.format(renameIsisMethodAnalyse, simulationName + + "." + name)); + engine.voidEval(String.format(renameIsisMethodAnalyse, + simulationName + "." + name)); + } - log.info("end"); + //Rename R objects for saving purpose + + String renameIsisSimule = "%s.isis.simule<-isis.simule"; + String renameIsisFactorDistribution = "%s.isis.factor.distribution<-isis.factor.distribution"; + String renameIsisFactor = "%s.isis.factor<-isis.factors"; + String renameIsisMethodExp = "%s.isis.methodExp<-isis.MethodExp"; + + String simulationName = simulationStorages.get(0).getName() + .replaceAll("-", ""); + + log.info("Message sent to R : " + + String.format(renameIsisSimule, simulationName)); + engine.voidEval(String.format(renameIsisSimule, simulationName)); + + log.info("Message sent to R : " + + String.format(renameIsisFactorDistribution, + simulationName)); + engine.voidEval(String.format(renameIsisFactorDistribution, + simulationName)); + + log.info("Message sent to R : " + + String.format(renameIsisFactor, simulationName)); + engine.voidEval(String.format(renameIsisFactor, simulationName)); + + log.info("Message sent to R : " + + String.format(renameIsisMethodExp, simulationName)); + engine.voidEval(String.format(renameIsisMethodExp, simulationName)); + + //Set working directory to save Isis R session + log.info("setwd(\"" + outputdirectory.getParent() + "\")"); + engine.voidEval("setwd(\"" + outputdirectory.getParent() + "\")"); + + // Save Isis R session + log.info("save.image()"); + engine.voidEval("save.image()"); + + engine.terminate(); + } catch (Exception e) { e.printStackTrace(); // Error while processing Modified: trunk/sensitivity/SensitivityCalculatorRRandomLHS.java =================================================================== --- trunk/sensitivity/SensitivityCalculatorRRandomLHS.java 2009-04-28 13:43:11 UTC (rev 125) +++ trunk/sensitivity/SensitivityCalculatorRRandomLHS.java 2009-04-29 08:39:07 UTC (rev 126) @@ -21,7 +21,11 @@ import java.io.File; import java.util.List; +import javax.swing.Box; +import javax.swing.JLabel; import javax.swing.JOptionPane; +import javax.swing.JScrollPane; +import javax.swing.JTextPane; import org.apache.commons.logging.Log; import org.apache.commons.logging.LogFactory; @@ -62,6 +66,8 @@ @Doc("Simulation number (default=10)") public int param_simulationNumber = 10; + @Doc("True to be able to modify the code sent to R") + public boolean param_modifR = false; /** * Retourne vrai si le calculateur sait gerer la cardinalité des facteurs @@ -112,17 +118,39 @@ //Get Isis R session log.info("load(\".RData\")"); engine.voidEval("load(\".RData\")"); + + //Set working directory for other results + log.info("setwd(\"" + outputdirectory + "\")"); + engine.voidEval("setwd(\"" + outputdirectory + "\")"); //Load the lhs library engine.voidEval("library(lhs)"); log.info("Message sent to R : " + "library(lhs)"); //Create the scenarios - engine.voidEval("x<-randomLHS(" + param_simulationNumber + "," - + factornumber + ")"); - log.info("Message sent to R : " + "x<-randomLHS(" - + param_simulationNumber + "," + factornumber + ")"); + String rInstruction = "x<-randomLHS(%s,%s)"; + String rCall = String.format(rInstruction,param_simulationNumber,factornumber); + + if (param_modifR) { + JLabel label = new JLabel( + "Modifier le code R envoyé si vous le souhaitez"); + JTextPane text = new JTextPane(); + text.setText(rCall); + text.setSize(400, 400); + text.setPreferredSize(text.getSize()); + Box box = Box.createVerticalBox(); + box.add(label); + box.add(new JScrollPane(text)); + + JOptionPane.showMessageDialog(null, box, "R modif", + JOptionPane.QUESTION_MESSAGE); + rCall = text.getText(); + } + + engine.voidEval(rCall); + log.info("Message sent to R : " + rCall); + // Get back experiment plan dataframe = (double[]) engine.eval("x"); log.info("Message sent to R" + "x"); @@ -139,9 +167,7 @@ if ((factor.getDomain() instanceof MatrixContinuousDomain) || (factor.getDomain() instanceof EquationContinuousDomain)) { factor - .setValueForIdentifier(Double.valueOf( - morris.getValue(new int[] { i, j })) - .toString()); + .setValueForIdentifier(morris.getValue(new int[] { i, j })); } else { Double value = (Double) ((ContinuousDomain) factor .getDomain()).getMinBound() @@ -385,13 +411,16 @@ //Get back the factors number int factorNumber = ((Double) engine.eval("length(factors)-1")) .intValue(); + + //Get back the simulation number + param_simulationNumber = (Integer)(engine.eval("length(factors[,1])")); - param_simulationNumber = simulationStorages.get(0).getParameter() + int sensitivityNumber = simulationStorages.get(0).getParameter() .getSensitivityExport().size(); - for (int k = 0; k < param_simulationNumber; k++) { + for (int k = 0; k < sensitivityNumber; k++) { - // Creates the R expression to import results in R + // Creates the R expression to import results in R String name = simulationStorages.get(0).getParameter() .getSensitivityExport().get(k).getExportFilename(); String rInstruction = name + "<-c("; @@ -408,7 +437,7 @@ .getExtensionFilename()); String simulResult = FileUtil.readAsString(importFile); double simulationResult = Double.valueOf(simulResult); - if (l < simulationStorages.size() - 1) { + if (l < param_simulationNumber - 1) { rInstruction = rInstruction + simulationResult + ","; } else { rInstruction = rInstruction + simulationResult; @@ -439,7 +468,7 @@ + "attr(isis.simule,\"call\")<-isis.MethodExp$call"); engine.voidEval("attr(isis.simule,\"call\")<-isis.MethodExp$call"); - for (int k = 0; k < param_simulationNumber; k++) { + for (int k = 0; k < sensitivityNumber; k++) { // Creates the R expression to import results in R String name = simulationStorages.get(0).getParameter()