Source code for pygod.detector.radar

# -*- coding: utf-8 -*-
""" Residual Analysis for Anomaly Detection in Attributed Networks
# Author: Kay Liu <>
# License: BSD 2 clause

import time
import torch
import warnings
from torch_geometric.utils import to_dense_adj

from . import Detector
from ..utils import logger, validate_device

[docs] class Radar(Detector): """ Residual Analysis for Anomaly Detection in Attributed Networks Radar is an anomaly detector with residual analysis. .. note:: This detector is transductive only. Using ``predict`` with unseen data will train the detector from scratch. See :cite:`li2017radar` for details. Parameters ---------- gamma : float, optional Weight of node feature w.r.t. structure. Default: ``1``. weight_decay : float, optional Weight decay (L2 penalty). Default: ``0.``. lr : float, optional Learning rate. Default: ``0.004``. epoch : int, optional Maximum number of training epoch. Default: ``100``. gpu : int GPU Index, -1 for using CPU. Default: ``-1``. contamination : float, optional Valid in (0., 0.5). The proportion of outliers in the data set. Used when fitting to define the threshold on the decision function. Default: ``0.1``. verbose : int, optional Verbosity mode. Range in [0, 3]. Larger value for printing out more log information. Default: ``0``. """ def __init__(self, gamma=1., weight_decay=0., lr=0.004, epoch=100, gpu=-1, contamination=0.1, verbose=0): super(Radar, self).__init__(contamination=contamination, verbose=verbose) # model param self.gamma = gamma self.weight_decay = weight_decay # training param = lr self.epoch = epoch self.device = validate_device(gpu) self.model = None
[docs] def fit(self, data, label=None): data.s = to_dense_adj(data.edge_index)[0] x, s, l, w_init, r_init = self.process_graph(data) self.model = RadarBase(w_init, r_init) optimizer = torch.optim.Adam(self.model.parameters(),, weight_decay=self.weight_decay) for epoch in range(self.epoch): start_time = time.time() x_, r = self.model(x) loss = self._loss(x, x_, r, l) optimizer.zero_grad() loss.backward() optimizer.step() decision_score = torch.sum(torch.pow(r, 2), dim=1).detach().cpu() logger(epoch=epoch, loss=loss.item(), score=self.decision_score_, target=label, time=time.time() - start_time, verbose=self.verbose, train=True) self.decision_score_ = decision_score self._process_decision_score() return self
def decision_function(self, data, label=None): if data is not None: warnings.warn("This detector is transductive only. " "Training from scratch with the input data."), label) return self.decision_score_ def process_graph(self, data): x = s = s = torch.max(s, s.T) laplacian = torch.diag(torch.sum(s, dim=1)) - s w_init = torch.eye(x.shape[0]).to(self.device) r_init = torch.inverse((1 + self.weight_decay) * torch.eye(x.shape[0]).to(self.device) + self.gamma * laplacian) @ x return x, s, laplacian, w_init, r_init def _loss(self, x, x_, r, l): return torch.norm(x - x_ - r, 2) + \ self.gamma * torch.trace(r.T @ l @ r)
class RadarBase(torch.nn.Module): def __init__(self, w, r): super(RadarBase, self).__init__() self.w = torch.nn.Parameter(w) self.r = torch.nn.Parameter(r) def forward(self, x): return self.w @ x, self.r