Skip to content

HMM Model

Hidden Markov Model implementation using hmmlearn.

pitch_sequencing.models.hmm_model

HMM model wrapper using hmmlearn CategoricalHMM.

HMMModel

Bases: BaseModel

Hidden Markov Model for pitch sequence prediction.

Source code in src/pitch_sequencing/models/hmm_model.py 
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
class HMMModel(BaseModel):
    """Hidden Markov Model for pitch sequence prediction."""

    def __init__(self, config=None):
        config = config or {}
        self.min_components = config.get("min_components", 1)
        self.max_components = config.get("max_components", 8)
        self.n_iter = config.get("n_iter", 100)
        self._model = None
        self._best_n = None

    @property
    def name(self) -> str:
        return "HMM"

    @property
    def model_type(self) -> str:
        return "sequence"

    def fit(self, X_train, y_train, X_val=None, y_val=None, **kwargs):
        """Train HMM by sweeping n_components and picking best by validation accuracy.

        For HMM, X_train is expected to be a flat 2D array of shape (n_samples, 1)
        with encoded pitch types (the HMM uses its own flat encoding, not windowed).
        y_train is the same flat array (self-supervised next-token prediction).
        """
        from hmmlearn import hmm as hmmlearn_hmm

        best_accuracy = 0
        best_model = None

        for n_components in range(self.min_components, self.max_components + 1):
            model = hmmlearn_hmm.CategoricalHMM(
                n_components=n_components,
                n_iter=self.n_iter,
                random_state=42,
            )
            model.fit(X_train)

            if X_val is not None:
                predicted = model.predict(X_val)
                actual = X_val.flatten()
                accuracy = np.mean(predicted == actual)
            else:
                predicted = model.predict(X_train)
                actual = X_train.flatten()
                accuracy = np.mean(predicted == actual)

            if accuracy > best_accuracy:
                best_accuracy = accuracy
                best_model = model
                self._best_n = n_components

        self._model = best_model

    def predict(self, X) -> np.ndarray:
        return self._model.predict(X)

    def predict_proba(self, X) -> np.ndarray:
        """Return emission probabilities for each sample given predicted state."""
        states = self._model.predict(X)
        emission = self._model.emissionprob_
        return emission[states]

    def get_params(self) -> dict:
        return {
            "n_components": self._best_n,
            "min_components": self.min_components,
            "max_components": self.max_components,
            "n_iter": self.n_iter,
        }

fit(X_train, y_train, X_val=None, y_val=None, **kwargs)

Train HMM by sweeping n_components and picking best by validation accuracy.

For HMM, X_train is expected to be a flat 2D array of shape (n_samples, 1) with encoded pitch types (the HMM uses its own flat encoding, not windowed). y_train is the same flat array (self-supervised next-token prediction).

Source code in src/pitch_sequencing/models/hmm_model.py 
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
def fit(self, X_train, y_train, X_val=None, y_val=None, **kwargs):
    """Train HMM by sweeping n_components and picking best by validation accuracy.

    For HMM, X_train is expected to be a flat 2D array of shape (n_samples, 1)
    with encoded pitch types (the HMM uses its own flat encoding, not windowed).
    y_train is the same flat array (self-supervised next-token prediction).
    """
    from hmmlearn import hmm as hmmlearn_hmm

    best_accuracy = 0
    best_model = None

    for n_components in range(self.min_components, self.max_components + 1):
        model = hmmlearn_hmm.CategoricalHMM(
            n_components=n_components,
            n_iter=self.n_iter,
            random_state=42,
        )
        model.fit(X_train)

        if X_val is not None:
            predicted = model.predict(X_val)
            actual = X_val.flatten()
            accuracy = np.mean(predicted == actual)
        else:
            predicted = model.predict(X_train)
            actual = X_train.flatten()
            accuracy = np.mean(predicted == actual)

        if accuracy > best_accuracy:
            best_accuracy = accuracy
            best_model = model
            self._best_n = n_components

    self._model = best_model

predict_proba(X)

Return emission probabilities for each sample given predicted state.

Source code in src/pitch_sequencing/models/hmm_model.py 
65
66
67
68
69
def predict_proba(self, X) -> np.ndarray:
    """Return emission probabilities for each sample given predicted state."""
    states = self._model.predict(X)
    emission = self._model.emissionprob_
    return emission[states]