edu.cmu.tetrad.algcomparison.algorithm.pairwise.Skew Maven / Gradle / Ivy
package edu.cmu.tetrad.algcomparison.algorithm.pairwise;
import edu.cmu.tetrad.algcomparison.algorithm.Algorithm;
import edu.cmu.tetrad.algcomparison.utils.TakesExternalGraph;
import edu.cmu.tetrad.annotation.AlgType;
import edu.cmu.tetrad.annotation.Bootstrapping;
import edu.cmu.tetrad.data.DataModel;
import edu.cmu.tetrad.data.DataSet;
import edu.cmu.tetrad.data.DataType;
import edu.cmu.tetrad.data.SimpleDataLoader;
import edu.cmu.tetrad.graph.EdgeListGraph;
import edu.cmu.tetrad.graph.Graph;
import edu.cmu.tetrad.search.Lofs;
import edu.cmu.tetrad.util.Parameters;
import edu.cmu.tetrad.util.Params;
import edu.pitt.dbmi.algo.resampling.GeneralResamplingTest;
import java.util.ArrayList;
import java.util.List;
/**
* Skew.
*
* @author josephramsey
*/
@edu.cmu.tetrad.annotation.Algorithm(
name = "Skew",
command = "skew",
algoType = AlgType.orient_pairwise,
dataType = DataType.Continuous
)
@Bootstrapping
public class Skew implements Algorithm, TakesExternalGraph {
private static final long serialVersionUID = 23L;
private Algorithm algorithm;
private Graph externalGraph;
public Skew() {
}
public Skew(Algorithm algorithm) {
this.algorithm = algorithm;
}
@Override
public Graph search(DataModel dataSet, Parameters parameters) {
if (parameters.getInt(Params.NUMBER_RESAMPLING) < 1) {
Graph graph = this.algorithm.search(dataSet, parameters);
if (graph != null) {
this.externalGraph = graph;
} else {
throw new IllegalArgumentException("This Skew algorithm needs both data and a graph source as inputs; it \n"
+ "will orient the edges in the input graph using the data");
}
List dataSets = new ArrayList<>();
dataSets.add(SimpleDataLoader.getContinuousDataSet(dataSet));
Lofs lofs = new Lofs(this.externalGraph, dataSets);
lofs.setRule(Lofs.Rule.Skew);
return lofs.orient();
} else {
Skew skew = new Skew(this.algorithm);
if (this.externalGraph != null) {
skew.setExternalGraph(this.algorithm);
}
DataSet data = (DataSet) dataSet;
GeneralResamplingTest search = new GeneralResamplingTest(data, skew, parameters.getInt(Params.NUMBER_RESAMPLING), parameters.getDouble(Params.PERCENT_RESAMPLE_SIZE), parameters.getBoolean(Params.RESAMPLING_WITH_REPLACEMENT), parameters.getInt(Params.RESAMPLING_ENSEMBLE), parameters.getBoolean(Params.ADD_ORIGINAL_DATASET));
search.setParameters(parameters);
search.setVerbose(parameters.getBoolean(Params.VERBOSE));
return search.search();
}
}
@Override
public Graph getComparisonGraph(Graph graph) {
return new EdgeListGraph(graph);
}
@Override
public String getDescription() {
return "Skew" + (this.algorithm != null ? " with initial graph from "
+ this.algorithm.getDescription() : "");
}
@Override
public DataType getDataType() {
return DataType.Continuous;
}
@Override
public List getParameters() {
List parameters = new ArrayList<>();
if (this.algorithm != null && !this.algorithm.getParameters().isEmpty()) {
parameters.addAll(this.algorithm.getParameters());
}
parameters.add(Params.VERBOSE);
return parameters;
}
@Override
public void setExternalGraph(Algorithm algorithm) {
if (algorithm == null) {
throw new IllegalArgumentException("This Skew algorithm needs both data and a graph source as inputs; it \n"
+ "will orient the edges in the input graph using the data.");
}
this.algorithm = algorithm;
}
}