org.matsim.contrib.socnetgen.sna.snowball.analysis.EstimatorTask.java Source code

Introduction

Here is the source code for org.matsim.contrib.socnetgen.sna.snowball.analysis.EstimatorTask.java

Source

/* *********************************************************************** *
 * project: org.matsim.*
 * EstimatorTask.java
 *                                                                         *
 * *********************************************************************** *
 *                                                                         *
 * copyright       : (C) 2010 by the members listed in the COPYING,        *
 *                   LICENSE and WARRANTY file.                            *
 * email           : info at matsim dot org                                *
 *                                                                         *
 * *********************************************************************** *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *   See also COPYING, LICENSE and WARRANTY file                           *
 *                                                                         *
 * *********************************************************************** */
package org.matsim.contrib.socnetgen.sna.snowball.analysis;

import gnu.trove.iterator.TIntDoubleIterator;
import gnu.trove.map.hash.TIntDoubleHashMap;
import gnu.trove.map.hash.TIntIntHashMap;

import java.io.BufferedWriter;
import java.io.FileWriter;
import java.io.IOException;
import java.util.Arrays;
import java.util.Map;

import org.apache.commons.math.stat.descriptive.DescriptiveStatistics;
import org.matsim.contrib.socnetgen.sna.graph.Edge;
import org.matsim.contrib.socnetgen.sna.graph.Graph;
import org.matsim.contrib.socnetgen.sna.graph.Vertex;
import org.matsim.contrib.socnetgen.sna.graph.analysis.AnalyzerTask;
import org.matsim.contrib.socnetgen.sna.snowball.SampledVertex;

/**
 * @author illenberger
 * 
 */
public class EstimatorTask extends AnalyzerTask {

    private final PiEstimator estimator;

    public EstimatorTask(PiEstimator estimator) {
        this.estimator = estimator;
    }

    @Override
    public void analyze(Graph graph, Map<String, DescriptiveStatistics> stats) {
        TIntDoubleHashMap probas = new TIntDoubleHashMap(graph.getVertices().size());
        TIntIntHashMap counts = new TIntIntHashMap(graph.getVertices().size());

        double N_estim = 0;
        for (Vertex vertex : graph.getVertices()) {
            if (((SampledVertex) vertex).isSampled()) {
                int k = vertex.getNeighbours().size();
                double p = estimator.probability((SampledVertex) vertex);

                probas.adjustOrPutValue(k, p, p);
                counts.adjustOrPutValue(k, 1, 1);

                if (p > 0)
                    N_estim += 1 / p;
            }
        }
        DescriptiveStatistics ds = new DescriptiveStatistics();
        ds.addValue(N_estim);
        stats.put("N_estim", ds);

        double M_estim = 0;
        for (Edge edge : graph.getEdges()) {
            SampledVertex v_i = (SampledVertex) edge.getVertices().getFirst();
            SampledVertex v_j = (SampledVertex) edge.getVertices().getSecond();
            if (v_i.isSampled() && v_j.isSampled()) {
                double p_i = estimator.probability(v_i);
                double p_j = estimator.probability(v_j);
                if (p_i > 0 && p_j > 0) {
                    M_estim += 1 / ((p_i + p_j) - (p_i * p_j));
                    //               M_estim += 1/(p_i * p_j);
                }
            }
        }
        ds = new DescriptiveStatistics();
        ds.addValue(M_estim);
        stats.put("M_estim", ds);

        TIntDoubleIterator it = probas.iterator();
        for (int i = 0; i < probas.size(); i++) {
            it.advance();
            it.setValue(it.value() / counts.get(it.key()));
        }

        writeValues(probas, "proba");
    }

    private void writeValues(TIntDoubleHashMap values, String valName) {
        int[] keys = values.keys();
        Arrays.sort(keys);

        BufferedWriter writer;
        try {
            writer = new BufferedWriter(
                    new FileWriter(String.format("%1$s/%2$s.txt", getOutputDirectory(), valName)));
            writer.write("k\t");
            writer.write(valName);
            writer.newLine();
            for (int key : keys) {
                writer.write(String.valueOf(key));
                writer.write("\t");
                writer.write(String.valueOf(values.get(key)));
                writer.newLine();
            }
            writer.close();
        } catch (IOException e) {
            e.printStackTrace();
        }
    }
}