estimate_risk_dream6: Expected risk estimation (comparison with litterature).

In pauwels2014: Bayesian Experimental Design for Systems Biology.

Description

Implements the strategy of Dream6 challenge winning team.

Usage

estimate_risk_dream6(thetas, knobj, experiment_fun)

Arguments

thetas	A sample from the posterior associated to the knoweldge list knobjs.
knobj	A knowledge list. See knobjs.
experiment_fun	A function that represents the molecular perturbation to be performed. See experiment_list1.

Details

The global variable observables should be defined.

Value

A dataframe with the following columns

Measurement	Factor representing possible measurements. See observables.
Risk	The risk associated to this measurement.
Cost	The cost associated to this measurement.

Author(s)

Edouard Pauwels

Examples

data(experiment_list1)
data(observables)
data(knobjs)
sapply(	
	1:length(knobjs),
	function(k){
		assign(names(knobjs)[k], knobjs[[k]], envir = .GlobalEnv)
	}
)

knobjActMult1$global_parameters$n_simu_weights <- 3

estimate_risk_dream6(knobjActMult1$datas[[1]]$thetas[1:10,], 
	knobjActMult1, experiment_list1$nothing)

pauwels2014 documentation built on May 1, 2019, 6:29 p.m.