Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

[REF] neurosynth decode #892

Open
wants to merge 1 commit into
base: main
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
32 changes: 17 additions & 15 deletions nimare/decode/base.py
Original file line number Diff line number Diff line change
Expand Up @@ -21,6 +21,14 @@ class Decoder(NiMAREBase):

__id_cols = ["id", "study_id", "contrast_id"]

def __init__(self):
self._required_inputs = None
self.inputs_ = {}
self.feature_group = None
self.features = None
self.features_ = []
self.frequency_threshold = 0.001

def _collect_inputs(self, dataset, drop_invalid=True):
"""Search for, and validate, required inputs as necessary."""
if not hasattr(dataset, "slice"):
Expand Down Expand Up @@ -53,34 +61,27 @@ def _preprocess_input(self, dataset):
This also takes into account which features have at least one study in the
Dataset with the feature.
"""
# Reduce feature list as desired
if self.feature_group is not None:
if not self.feature_group.endswith("__"):
self.feature_group += "__"
feature_names = self.inputs_["annotations"].columns.values
feature_names = [f for f in feature_names if f.startswith(self.feature_group)]
if self.features is not None:
features = [f.split("__")[-1] for f in feature_names if f in self.features]
else:
features = feature_names
features = [
f.split("__")[-1] for f in feature_names if self.features and f in self.features
]
else:
if self.features is None:
features = self.inputs_["annotations"].columns.values
else:
features = self.features
features = self.features or self.inputs_["annotations"].columns.values

features = [f for f in features if f not in self.__id_cols]
n_features_orig = len(features)

# At least one study in the dataset much have each label
counts = (self.inputs_["annotations"][features] > self.frequency_threshold).sum(0)
features = counts[counts > 0].index.tolist()
if not len(features):
self.features_ = counts[counts > 0].index.tolist()
if not self.features_:
raise Exception("No features identified in Dataset!")
elif len(features) < n_features_orig:
LGR.info(f"Retaining {len(features)}/{n_features_orig} features.")

self.features_ = features
elif len(self.features_) < n_features_orig:
LGR.info(f"Retaining {len(self.features_)}/{n_features_orig} features.")

def fit(self, dataset, drop_invalid=True):
"""Fit Decoder to Dataset.
Expand Down Expand Up @@ -113,3 +114,4 @@ def _fit(self, dataset):

Must return a DataFrame, with one row for each feature.
"""
pass
2 changes: 1 addition & 1 deletion nimare/decode/discrete.py
Original file line number Diff line number Diff line change
Expand Up @@ -449,6 +449,7 @@ def __init__(
u=0.05,
correction="fdr_bh",
):
super().__init__()
self.feature_group = feature_group
self.features = features
self.frequency_threshold = frequency_threshold
Expand Down Expand Up @@ -499,7 +500,6 @@ def neurosynth_decode(
annotations,
ids,
ids2=None,
feature_group=None,
features=None,
frequency_threshold=0.001,
prior=0.5,
Expand Down
Loading