prototorch_models/examples/gmlvq_spiral.py

95 lines
2.2 KiB
Python
Raw Normal View History

2021-11-17 14:01:38 +00:00
"""GMLVQ example using the spiral dataset."""
2021-05-21 15:55:55 +00:00
import argparse
import warnings
2021-05-21 15:55:55 +00:00
2021-06-04 20:21:28 +00:00
import prototorch as pt
import pytorch_lightning as pl
2021-05-21 16:54:47 +00:00
import torch
from prototorch.models import (
GMLVQ,
PruneLoserPrototypes,
VisGLVQ2D,
)
from pytorch_lightning.callbacks import EarlyStopping
from pytorch_lightning.utilities.seed import seed_everything
from pytorch_lightning.utilities.warnings import PossibleUserWarning
from torch.utils.data import DataLoader
warnings.filterwarnings("ignore", category=PossibleUserWarning)
warnings.filterwarnings("ignore", category=UserWarning)
if __name__ == "__main__":
# Reproducibility
seed_everything(seed=4)
2021-05-21 15:55:55 +00:00
# Command-line arguments
parser = argparse.ArgumentParser()
parser = pl.Trainer.add_argparse_args(parser)
args = parser.parse_args()
# Dataset
2021-06-04 20:21:28 +00:00
train_ds = pt.datasets.Spiral(num_samples=500, noise=0.5)
# Dataloaders
train_loader = DataLoader(train_ds, batch_size=256)
# Hyperparameters
2021-05-25 13:41:10 +00:00
num_classes = 2
2021-06-04 20:21:28 +00:00
prototypes_per_class = 10
hparams = dict(
2021-05-25 13:41:10 +00:00
distribution=(num_classes, prototypes_per_class),
2021-06-04 20:21:28 +00:00
transfer_function="swish_beta",
2021-05-04 18:56:16 +00:00
transfer_beta=10.0,
2021-06-04 20:21:28 +00:00
proto_lr=0.1,
bb_lr=0.1,
input_dim=2,
latent_dim=2,
)
# Initialize the model
model = GMLVQ(
2021-06-04 20:21:28 +00:00
hparams,
optimizer=torch.optim.Adam,
prototypes_initializer=pt.initializers.SSCI(train_ds, noise=1e-2),
2021-06-04 20:21:28 +00:00
)
# Callbacks
vis = VisGLVQ2D(
2021-06-04 20:21:28 +00:00
train_ds,
show_last_only=False,
block=False,
)
pruning = PruneLoserPrototypes(
threshold=0.01,
2021-06-04 20:21:28 +00:00
idle_epochs=10,
prune_quota_per_epoch=5,
frequency=5,
2021-06-04 20:21:28 +00:00
replace=True,
prototypes_initializer=pt.initializers.SSCI(train_ds, noise=1e-1),
2021-06-04 20:21:28 +00:00
verbose=True,
)
es = EarlyStopping(
2021-06-04 20:21:28 +00:00
monitor="train_loss",
min_delta=1.0,
patience=5,
mode="min",
check_on_train_epoch_end=True,
)
# Setup trainer
2021-05-21 15:55:55 +00:00
trainer = pl.Trainer.from_argparse_args(
args,
2021-06-04 20:21:28 +00:00
callbacks=[
vis,
es,
2021-06-04 20:21:28 +00:00
pruning,
],
max_epochs=1000,
log_every_n_steps=1,
detect_anomaly=True,
)
# Training loop
trainer.fit(model, train_loader)