''' Created on Mar 15, 2012 .. codeauthor:: jhkwakkel <j.h.kwakkel (at) tudelft (dot) nl> ''' from analysis import clusterer from util import ema_logging from core import ModelEnsemble from test.scarcity_example import ScarcityModel if __name__ == "__main__": ema_logging.log_to_stderr(ema_logging.INFO) model = ScarcityModel(r'..\..\src\test', "fluCase") ensemble = ModelEnsemble() ensemble.set_model_structure(model) ensemble.parallel = True results = ensemble.perform_experiments(200) clusterer.cluster(data=results, outcome='relative market price', distance='gonenc', cMethod='maxclust', cValue=5, plotDendrogram=False)
''' Created on Mar 15, 2012 @author: jhkwakkel ''' from analysis import clusterer from expWorkbench import ema_logging from expWorkbench import ModelEnsemble from test.scarcity_example import ScarcityModel if __name__ == "__main__": ema_logging.log_to_stderr(ema_logging.INFO) model = ScarcityModel(r'..\..\src\test', "fluCase") ensemble = ModelEnsemble() ensemble.set_model_structure(model) ensemble.parallel = True results = ensemble.perform_experiments(200) clusterer.cluster(data=results, outcome='relative market price', distance='gonenc', cMethod='maxclust', cValue=5, plotDendrogram=False)
#load the data data = load_results(r'..\gallery\data\100 flu cases no policy.cPickle') # specify the number of desired clusters # note: the meaning of cValue is tied to the value for cMethod cValue = 5 #perform cluster analysis dist, clusteraloc, runlog, z = cluster(data=data, outcome='deceased population region 1', distance='gonenc', interClusterDistance='complete', cMethod = 'maxclust', cValue = cValue, plotDendrogram=False, plotClusters=False, groupPlot=False, sisterCount=100, tHoldCurvature = 0.1, tHoldSlope = 0.1 ) #the plotting fig = plt.figure() #for each cluster for i in range(1, cValue+1): #get the data values = data[1]['deceased population region 1'] values = values[clusteraloc==i]
''' Created on Sep 8, 2011 @author: gonengyucel ''' from analysis.clusterer import cluster from expWorkbench import load_results from expWorkbench import EMAlogging EMAlogging.log_to_stderr(EMAlogging.INFO) data = load_results(r'..\analysis\1000 flu cases.cPickle') cluster(data=data, outcome='deceased population region 1', distance='gonenc', interClusterDistance='complete', plotDendrogram=True, plotClusters=False, groupPlot=False, sisterCount=100, )