import numpy as np
import pandas as pd
from sklearn.ensemble import ExtraTreesClassifier, GradientBoostingClassifier
from sklearn.feature_selection import RFE
from sklearn.model_selection import train_test_split
from sklearn.pipeline import make_pipeline
from tpot.export_utils import set_param_recursive

# NOTE: Make sure that the outcome column is labeled 'target' in the data file
tpot_data = pd.read_csv("PATH/TO/DATA/FILE",
                        sep="COLUMN_SEPARATOR", dtype=np.float64)
features = tpot_data.drop("target", axis=1)
training_features, testing_features, training_target, testing_target = train_test_split(
    features, tpot_data["target"], random_state=42
)

# Average CV score on the training set was: 0.8490676876315174
exported_pipeline = make_pipeline(
    RFE(
        estimator=ExtraTreesClassifier(
            criterion="entropy", max_features=0.2, n_estimators=100
        ),
        step=0.9000000000000001,
    ),
    GradientBoostingClassifier(
        learning_rate=0.01,
        max_depth=3,
        max_features=0.2,
        min_samples_leaf=1,
        min_samples_split=4,
        n_estimators=100,
        subsample=0.4,
    ),
)
# Fix random state for all the steps in exported pipeline
set_param_recursive(exported_pipeline.steps, "random_state", 42)

exported_pipeline.fit(training_features, training_target)
results = exported_pipeline.predict(testing_features)
