matrix_factorize.py

import torch
import torch.nn as nn
from torch.nn import init
from torch.autograd import Variable
import pickle
import numpy as np
import time
import random
from collections import defaultdict
from UV_Encoders import UV_Encoder
from UV_Aggregators import UV_Aggregator
from Social_Encoders import Social_Encoder
from Social_Aggregators import Social_Aggregator
import torch.nn.functional as F
import torch.utils.data
from sklearn.metrics import mean_squared_error
from sklearn.metrics import mean_absolute_error
from math import sqrt
import datetime
import argparse
import os
import json
"""
GraphRec: Graph Neural Networks for Social Recommendation. 
Wenqi Fan, Yao Ma, Qing Li, Yuan He, Eric Zhao, Jiliang Tang, and Dawei Yin. 
In Proceedings of the 28th International Conference on World Wide Web (WWW), 2019. Preprint[https://arxiv.org/abs/1902.07243]

If you use this code, please cite our paper:
```
@inproceedings{fan2019graph,
  title={Graph Neural Networks for Social Recommendation},
  author={Fan, Wenqi and Ma, Yao and Li, Qing and He, Yuan and Zhao, Eric and Tang, Jiliang and Yin, Dawei},
  booktitle={WWW},
  year={2019}
}
```

"""

def fx(x):
    return x
class Matrix(nn.Module):

    def __init__(self, u2e, v2e, embed_dim=256):
        super(Matrix, self).__init__()
        self.u2e = u2e
        self.v2e = v2e
        self.embed_dim = embed_dim

        
        self.criterion = nn.MSELoss()
        self.dropout = nn.Dropout(p=0.5)

    def forward(self, nodes_u, nodes_v):
        scores = torch.sum(self.u2e.weight[nodes_u]*self.v2e.weight[nodes_v],dim=-1)
        return scores.squeeze()

    def loss(self, nodes_u, nodes_v, labels_list):
        scores = self.forward(nodes_u, nodes_v)
        return self.criterion(scores, labels_list)

import sys
def train(model, device, train_loader, optimizer, epoch, best_rmse, best_mae):
    model.train()
    running_loss = 0.0
    expected_rmse = 999
    steps = []
    for i, data in enumerate(train_loader, 0):
        batch_nodes_u, batch_nodes_v, labels_list = data
        optimizer.zero_grad()
        loss = model.loss(batch_nodes_u.to(device), batch_nodes_v.to(device), labels_list.to(device))
        loss.backward(retain_graph=True)
        optimizer.step()
        running_loss += loss.item()
        steps.append(loss.item())
        sys.stdout.write('\r>>'+'step ' +str(i%100)+', Loss '+str(np.mean(steps)) )
         
        sys.stdout.flush()
        
        if i % 100 == 0:
            steps = []
            print('[%d, %5d] loss: %.3f, The best rmse/mae: %.6f / %.6f' % (
                epoch, i, running_loss / 100, best_rmse, best_mae))
            running_loss = 0.0
            sys.stdout.write('\n')
    return 0


def test(model, device, test_loader):
    model.eval()
    tmp_pred = []
    target = []
    with torch.no_grad():
        for test_u, test_v, tmp_target in test_loader:
          #  print(test_u.size(), test_v.size(), tmp_target.size() )
            test_u, test_v, tmp_target = test_u.to(device), test_v.to(device), tmp_target.to(device)
            val_output = model.forward(test_u, test_v)
            #print(val_output.size())
            tmp_pred.append(list(val_output.data.cpu().numpy()))
            target.append(list(tmp_target.data.cpu().numpy()))
    tmp_pred = np.array(sum(tmp_pred, []))
    target = np.array(sum(target, []))
    print(tmp_pred,target)
    expected_rmse = sqrt(mean_squared_error(tmp_pred, target))
    mae = mean_absolute_error(tmp_pred, target)
    return expected_rmse, mae

def rebuild(a):
    b = {}
    for u,v in a.items():
        b[int(u)] = v
    return b
def rebuild_list(a):
    res = []
    for i in a:
        res.append(int(i))
    return res
def main():
    # Training settings
    parser = argparse.ArgumentParser(description='Social Recommendation: GraphRec model')
    parser.add_argument('--batch_size', type=int, default=10000, metavar='N', help='input batch size for training')
    parser.add_argument('--embed_dim', type=int, default=256, metavar='N', help='embedding size')
    parser.add_argument('--phase', type=str, default="test", metavar='N', help='test phase')
    parser.add_argument('--dataset', type=str, default="100k", metavar='N', help='80 for 80-10-10 split dataset, 60 for 60-20-20 split dataset')
    parser.add_argument('--lr', type=float, default=0.001, metavar='LR', help='learning rate')
    parser.add_argument('--test_batch_size', type=int, default=1000, metavar='N', help='input batch size for testing')
    parser.add_argument('--epochs', type=int, default=200, metavar='N', help='number of epochs to train')
    args = parser.parse_args()

    os.environ['CUDA_VISIBLE_DEVICES'] = '-1'
    use_cuda = False
    if torch.cuda.is_available():
        use_cuda = True
    device = torch.device("cuda" if use_cuda else "cpu")

    embed_dim = args.embed_dim
    if args.dataset == "100k":
        data = json.load(open("data/100k/train_100k.json","r"))
        train_u = rebuild_list(data["users"])
        train_v = rebuild_list(data["items"])
        train_r = rebuild_list(data["rating"])
        data = json.load(open("data/100k/test_100k.json","r"))
        test_u = rebuild_list(data["users"])
        test_v = rebuild_list(data["items"])
        test_r = rebuild_list(data["rating"])
        user_cluster = rebuild(json.load(open("data/100k/user_cluster_100k.json","r")))
        user_neighbor = rebuild(json.load(open("data/100k/user_neighbor_100k_2.json","r")))
        movie_cluster = rebuild(json.load(open("data/100k/movie_cluster_100k.json","r")))
        movie_neighbor = rebuild(json.load(open("data/100k/movie_neighbor_100k_2.json","r")))


    elif args.dataset == "1m":
        data = json.load(open("data/1m/train_1m.json","r"))
        train_u = rebuild_list(data["users"])
        train_v = rebuild_list(data["items"])
        train_r = rebuild_list(data["rating"])
        data = json.load(open("data/1m/test_1m.json","r"))
        test_u = rebuild_list(data["users"])
        test_v = rebuild_list(data["items"])
        test_r = rebuild_list(data["rating"])
        user_cluster = rebuild(json.load(open("data/1m/user_cluster_1m.json","r")))
        user_neighbor = rebuild(json.load(open("data/1m/user_neighbor_1m.json","r")))
        movie_cluster = rebuild(json.load(open("data/1m/movie_cluster_1m.json","r")))
        movie_neighbor = rebuild(json.load(open("data/1m/movie_neighbor_1m.json","r")))
    mean = np.mean(test_r)

    print("predict mean rmse", np.mean(np.abs(test_r -mean)**2)**0.5)
    print("predict mean mae", np.mean(np.abs(test_r -mean)))
    print(len(train_u),len(test_v),len(test_r))
    res = []
    
    #print("user avg ",np.mean(np.abs(np.array(test_r)- np.array(res))**2)**0.5)
    
    trainset = torch.utils.data.TensorDataset(torch.LongTensor(train_u), torch.LongTensor(train_v),
                                              torch.FloatTensor(train_r))
    testset = torch.utils.data.TensorDataset(torch.LongTensor(test_u), torch.LongTensor(test_v),
                                             torch.FloatTensor(test_r))
    train_loader = torch.utils.data.DataLoader(trainset, batch_size=args.batch_size, shuffle=True)
    test_loader = torch.utils.data.DataLoader(testset, batch_size=args.test_batch_size, shuffle=True)
    num_users = max(train_u)+1 #history_u_lists.__len__()
    num_items = max(train_v)+1# history_v_lists.__len__()
    
    print(num_items)
    u2e = nn.Embedding(num_users, embed_dim).to(device)
    u2e.weight.data.uniform_(-0.1, 0.1)
    v2e = nn.Embedding(num_items, embed_dim).to(device)
    v2e.weight.data.uniform_(-0.1, 0.1)
    u_cluster2e = nn.Embedding(100, embed_dim).to(device)
    v_cluster2e = nn.Embedding(100, embed_dim).to(device)
    

    # user feature
    # features: item * rating
    
    # neighobrs
    agg_u_social = Social_Aggregator(None, u2e, embed_dim, cuda=device)
    enc_u = Social_Encoder(u2e, embed_dim, user_neighbor, agg_u_social,
                           base_model=u_cluster2e, cuda=device,cluster=user_cluster)

    # item feature: user * rating
    agg_v_social = Social_Aggregator(None, v2e, embed_dim, cuda=device)
    enc_v = Social_Encoder(v2e, embed_dim, movie_neighbor, agg_v_social,
                           base_model=v_cluster2e, cuda=device,cluster=movie_cluster)

    # model
    graphrec = Matrix(u2e, v2e, embed_dim).to(device)
    optimizer = torch.optim.Adam(graphrec.parameters(), lr=args.lr)
    
    graphrec.load_state_dict(torch.load("checkpoint/checkpoint_256_matrix.pth"), strict=False) #RMSprop , alpha=0.9
    #graphrec.eval()
    expected_rmse, mae = test(graphrec, device, test_loader)
    print("rmse: %.4f, mae:%.4f " % (expected_rmse, mae))
    
    best_rmse = 9
    best_mae = 9
    endure_count = 0
    
    for epoch in range(1, args.epochs + 1):

        train(graphrec, device, train_loader, optimizer, epoch, best_rmse, best_mae)
        expected_rmse, mae = test(graphrec, device, test_loader)
        # please add the validation set to tune the hyper-parameters based on your datasets.

        # early stopping (no validation set in toy dataset)
        if best_rmse > expected_rmse:
            best_rmse = expected_rmse
            torch.save(graphrec.state_dict(), "checkpoint/checkpoint_256_matrix.pth")
            #best_mae = mae
            endure_count = 0
        if best_mae > mae:
            best_mae = mae
            torch.save(graphrec.state_dict(), "checkpoint/checkpoint_best_mae_256.pth")
        else:
            endure_count += 1
        print("rmse: %.4f, mae:%.4f " % (expected_rmse, mae))
        print("best rmse: %.4f, mae:%.4f " % (best_rmse, best_mae))
        torch.save(graphrec.state_dict(), "checkpoint.pth")
        #if endure_count > 5:
        #    break
    

if __name__ == "__main__":
    main()

## 100k

## matrix 256 0.9567
## matrix 128 0.9440
##1m

# rmse: 0.8714