Examples with RunProfileLikelihoodGeneral org.vcell.vmicro.workflow.task.RunProfileLikelihoodGeneral used on opensource projects

Example 1 with RunProfileLikelihoodGeneral

use of org.vcell.vmicro.workflow.task.RunProfileLikelihoodGeneral in project vcell by virtualcell.

the class OptModelParamPanel method computeProfileLikelihood.

public void computeProfileLikelihood() {
    final Repository repository = new MemoryRepository();
    Workflow workflow = new Workflow("profileLikelihoodWorkflow");
    final TaskContext context = new TaskContext(workflow, repository, localWorkspace);
    System.err.println("OptModelParamPanel.showParameterEvaluation(): how do we pass in the initial guess to ProfileLikelihood code??? should be an independent input to ProfileLikelihood so that it is explicit ... and OptContext is immutable.????");
    final RunProfileLikelihoodGeneral runProfileLikelihoodGeneral = new RunProfileLikelihoodGeneral("internal");
    WorkflowParameter<OptContext> optContextParam = workflow.addParameter(OptContext.class, "optContext", repository, optContext);
    workflow.connectParameter(optContextParam, runProfileLikelihoodGeneral.optContext);
    workflow.addTask(runProfileLikelihoodGeneral);
    final DisplayProfileLikelihoodPlots displayProfileLikelihoodPlots = new DisplayProfileLikelihoodPlots("displayProfileLikihood");
    workflow.connect2(runProfileLikelihoodGeneral.profileData, displayProfileLikelihoodPlots.profileData);
    WorkflowParameter<String> titleParam = workflow.addParameter(String.class, "title", repository, "profile likelihood");
    workflow.connectParameter(titleParam, displayProfileLikelihoodPlots.title);
    workflow.addTask(displayProfileLikelihoodPlots);
    AsynchClientTask evaluateTask = new AsynchClientTask("Prepare to evaluate parameters ...", AsynchClientTask.TASKTYPE_NONSWING_BLOCKING) {

        public void run(Hashtable<String, Object> hashTable) throws Exception {
            runProfileLikelihoodGeneral.compute(context, getClientTaskStatusSupport());
        }
    };
    AsynchClientTask showResultTask = new AsynchClientTask("Showing profile likelihood and confidence intervals ...", AsynchClientTask.TASKTYPE_SWING_BLOCKING) {

        public void run(Hashtable<String, Object> hashTable) throws Exception {
            displayProfileLikelihoodPlots.compute(context, getClientTaskStatusSupport());
        }
    };
    // dispatch
    ClientTaskDispatcher.dispatch(OptModelParamPanel.this, new Hashtable<String, Object>(), new AsynchClientTask[] { evaluateTask, showResultTask }, false, true, null, true);
}

Example 2 with RunProfileLikelihoodGeneral

use of org.vcell.vmicro.workflow.task.RunProfileLikelihoodGeneral in project vcell by virtualcell.

the class WorkflowTest method getVFrapExample.

public static Workflow getVFrapExample(Repository repository, File workingDirectory, File vfrapFile) {
    // 
    // construct the dataflow graph
    // 
    LocalWorkspace localWorkspace = new LocalWorkspace(workingDirectory);
    Workflow workflow = new Workflow("main");
    TaskContext context = new TaskContext(workflow, repository, localWorkspace);
    // 
    // workflow parameters
    // 
    WorkflowParameter<Double> maxIntensity = workflow.addParameter(Double.class, "maxIntensity");
    context.setParameterValue(maxIntensity, 10000.0);
    WorkflowParameter<Boolean> bNoise = workflow.addParameter(Boolean.class, "bNoise");
    context.setParameterValue(bNoise, true);
    WorkflowParameter<File> vfrapFileParam = workflow.addParameter(File.class, "vfrapFile");
    context.setParameterValue(vfrapFileParam, vfrapFile);
    WorkflowParameter<Double> bleachThreshold = workflow.addParameter(Double.class, "bleachThreshold");
    context.setParameterValue(bleachThreshold, 0.80);
    WorkflowParameter<Double> cellThreshold = workflow.addParameter(Double.class, "cellThreshold");
    context.setParameterValue(cellThreshold, 0.5);
    WorkflowParameter<String> modelTypeOne = workflow.addParameter(String.class, "modelType");
    context.setParameterValue(modelTypeOne, GenerateRefSimOptModel.ModelType.DiffOne.toString());
    WorkflowParameter<String> modelTypeTwo = workflow.addParameter(String.class, "modelType2");
    context.setParameterValue(modelTypeTwo, GenerateRefSimOptModel.ModelType.DiffTwoWithPenalty.toString());
    WorkflowParameter<String> displayROITitle = workflow.addParameter(String.class, "displayROITitle");
    context.setParameterValue(displayROITitle, "rois");
    // input data from VFrap file
    ImportRawTimeSeriesFromVFrap timeSeriesFromVCell = new ImportRawTimeSeriesFromVFrap("frapImport");
    workflow.connectParameter(vfrapFileParam, timeSeriesFromVCell.vfrapFile);
    workflow.addTask(timeSeriesFromVCell);
    // // input data from HDF5 postprocessing (includes PSF)
    // ImportRawTimeSeriesFromHdf5Fluor timeSeriesFromVCell = new ImportRawTimeSeriesFromHdf5Fluor("t1");
    // timeSeriesFromVCell.bNoise.setSource(bNoise);
    // timeSeriesFromVCell.maxIntensity.setSource(maxIntensity);
    // timeSeriesFromVCell.fluorDataName.setSource(hdf5FluorName);
    // timeSeriesFromVCell.vcellHdf5File.setSource(hdf5File);
    // timeSeriesFromVCell.zSliceIndex.setSource(hdf5ZSliceIndex);
    // workflow.addTask(timeSeriesFromVCell);
    GenerateCellROIsFromRawTimeSeries generateCellROIs = new GenerateCellROIsFromRawTimeSeries("generateCellROIs");
    workflow.connectParameter(cellThreshold, generateCellROIs.cellThreshold);
    workflow.connect2(timeSeriesFromVCell.rawTimeSeriesImages, generateCellROIs.rawTimeSeriesImages);
    workflow.addTask(generateCellROIs);
    GenerateNormalizedFrapData generateNormalizedFrapData = new GenerateNormalizedFrapData("generateNormalizedFrapData");
    workflow.connect2(generateCellROIs.backgroundROI_2D, generateNormalizedFrapData.backgroundROI_2D);
    workflow.connect2(generateCellROIs.indexOfFirstPostbleach, generateNormalizedFrapData.indexOfFirstPostbleach);
    workflow.connect2(timeSeriesFromVCell.rawTimeSeriesImages, generateNormalizedFrapData.rawImageTimeSeries);
    workflow.addTask(generateNormalizedFrapData);
    GenerateBleachROI generateROIs = new GenerateBleachROI("generateROIs");
    workflow.connectParameter(bleachThreshold, generateROIs.bleachThreshold);
    workflow.connect2(generateCellROIs.cellROI_2D, generateROIs.cellROI_2D);
    workflow.connect2(generateNormalizedFrapData.normalizedFrapData, generateROIs.normalizedTimeSeries);
    workflow.addTask(generateROIs);
    GenerateDependentImageROIs generateDependentROIs = new GenerateDependentImageROIs("generateDependentROIs");
    workflow.connect2(generateCellROIs.cellROI_2D, generateDependentROIs.cellROI_2D);
    workflow.connect2(generateROIs.bleachedROI_2D, generateDependentROIs.bleachedROI_2D);
    workflow.addTask(generateDependentROIs);
    DisplayDependentROIs displayDependentROIs = new DisplayDependentROIs("displayDependentROIs");
    workflow.connect2(generateDependentROIs.imageDataROIs, displayDependentROIs.imageROIs);
    workflow.connect2(generateCellROIs.cellROI_2D, displayDependentROIs.cellROI);
    workflow.connectParameter(displayROITitle, displayDependentROIs.title);
    workflow.addTask(displayDependentROIs);
    GenerateReducedData generateReducedNormalizedData = new GenerateReducedData("generateReducedNormalizedData");
    workflow.connect2(generateNormalizedFrapData.normalizedFrapData, generateReducedNormalizedData.imageTimeSeries);
    workflow.connect2(generateDependentROIs.imageDataROIs, generateReducedNormalizedData.imageDataROIs);
    workflow.addTask(generateReducedNormalizedData);
    ComputeMeasurementError computeMeasurementError = new ComputeMeasurementError("computeMeasurementError");
    workflow.connect2(generateDependentROIs.imageDataROIs, computeMeasurementError.imageDataROIs);
    workflow.connect2(generateCellROIs.indexOfFirstPostbleach, computeMeasurementError.indexFirstPostbleach);
    workflow.connect2(generateNormalizedFrapData.prebleachAverage, computeMeasurementError.prebleachAverage);
    workflow.connect2(timeSeriesFromVCell.rawTimeSeriesImages, computeMeasurementError.rawImageTimeSeries);
    workflow.addTask(computeMeasurementError);
    // 
    // SLOW WAY
    // 
    // RunRefSimulation runRefSimulationFull = new RunRefSimulation("runRefSimulationFull");
    // runRefSimulationFull.cellROI_2D.setSource(generateCellROIs.cellROI_2D);
    // runRefSimulationFull.normalizedTimeSeries.setSource(generateNormalizedFrapData.normalizedFrapData);
    // workflow.addTask(runRefSimulationFull);
    // GenerateReducedRefData generateReducedRefSimData = new GenerateReducedRefData("generateReducedRefSimData");
    // generateReducedRefSimData.imageTimeSeries.setSource(runRefSimulationFull.refSimTimeSeries);
    // generateReducedRefSimData.imageDataROIs.setSource(generateDependentROIs.imageDataROIs);
    // workflow.addTask(generateReducedRefSimData);
    // Generate2DOptContext generate2DOptContextFull = new Generate2DOptContext("generate2DOptContextFull");
    // generate2DOptContextFull.normalizedMeasurementErrors.setSource(computeMeasurementError.normalizedMeasurementError);
    // generate2DOptContextFull.modelType.setSource(modelType);
    // generate2DOptContextFull.normExpData.setSource(generateReducedNormalizedData.reducedROIData);
    // generate2DOptContextFull.refSimData.setSource(generateReducedRefSimData.reducedROIData);
    // generate2DOptContextFull.refSimDiffusionRate.setSource(runRefSimulationFull.refSimDiffusionRate);
    // workflow.addTask(generate2DOptContextFull);
    // 
    // FAST WAY
    // 
    GenerateTrivial2DPsf psf_2D = new GenerateTrivial2DPsf("psf_2D");
    workflow.addTask(psf_2D);
    RunRefSimulationFast runRefSimulationFast = new RunRefSimulationFast("runRefSimulationFast");
    workflow.connect2(generateCellROIs.cellROI_2D, runRefSimulationFast.cellROI_2D);
    workflow.connect2(generateNormalizedFrapData.normalizedFrapData, runRefSimulationFast.normalizedTimeSeries);
    workflow.connect2(generateDependentROIs.imageDataROIs, runRefSimulationFast.imageDataROIs);
    workflow.connect2(psf_2D.psf_2D, runRefSimulationFast.psf);
    workflow.addTask(runRefSimulationFast);
    // 
    // model with One mobile fraction
    // 
    {
        GenerateRefSimOptModel generateRefSimOptModelOne = new GenerateRefSimOptModel("generateRefSimOptModel");
        workflow.connectParameter(modelTypeOne, generateRefSimOptModelOne.modelType);
        workflow.connect2(runRefSimulationFast.reducedROIData, generateRefSimOptModelOne.refSimData);
        workflow.connect2(runRefSimulationFast.refSimDiffusionRate, generateRefSimOptModelOne.refSimDiffusionRate);
        workflow.addTask(generateRefSimOptModelOne);
        Generate2DOptContext generate2DOptContextOne = new Generate2DOptContext("generate2DOptContextOne");
        workflow.connect2(computeMeasurementError.normalizedMeasurementError, generate2DOptContextOne.normalizedMeasurementErrors);
        workflow.connect2(generateReducedNormalizedData.reducedROIData, generate2DOptContextOne.normExpData);
        workflow.connect2(generateRefSimOptModelOne.optModel, generate2DOptContextOne.optModel);
        workflow.addTask(generate2DOptContextOne);
        RunProfileLikelihoodGeneral runProfileLikelihoodOne = new RunProfileLikelihoodGeneral("runProfileLikelihoodOne");
        workflow.connect2(generate2DOptContextOne.optContext, runProfileLikelihoodOne.optContext);
        workflow.addTask(runProfileLikelihoodOne);
    }
    // 
    // model with One mobile fraction
    // 
    GenerateRefSimOptModel generateRefSimOptModelTwo = new GenerateRefSimOptModel("generateRefSimOptModelTwo");
    workflow.connectParameter(modelTypeTwo, generateRefSimOptModelTwo.modelType);
    workflow.connect2(runRefSimulationFast.reducedROIData, generateRefSimOptModelTwo.refSimData);
    workflow.connect2(runRefSimulationFast.refSimDiffusionRate, generateRefSimOptModelTwo.refSimDiffusionRate);
    workflow.addTask(generateRefSimOptModelTwo);
    Generate2DOptContext generate2DOptContextTwo = new Generate2DOptContext("generate2DOptContextTwo");
    workflow.connect2(computeMeasurementError.normalizedMeasurementError, generate2DOptContextTwo.normalizedMeasurementErrors);
    workflow.connect2(generateRefSimOptModelTwo.optModel, generate2DOptContextTwo.optModel);
    workflow.connect2(generateReducedNormalizedData.reducedROIData, generate2DOptContextTwo.normExpData);
    workflow.addTask(generate2DOptContextTwo);
    RunProfileLikelihoodGeneral runProfileLikelihoodTwo = new RunProfileLikelihoodGeneral("runProfileLikelihoodTwo");
    workflow.connect2(generate2DOptContextTwo.optContext, runProfileLikelihoodTwo.optContext);
    workflow.addTask(runProfileLikelihoodTwo);
    return workflow;
}