Skip to content

Implementation of MoStress: a Sequence Model for Stress Classification

License

Notifications You must be signed in to change notification settings

ramos-ai/MoStress

Repository files navigation

README

1. Introduction

This is the implementation of MoStress: a sequence model for stress classification. For more information about the model itself and all the processes and experiments around it, please see the articles on Publications section . For more questions and discussions, please contact any of the authors of this repository.

The resume of MoStress is shown on the image and the diagram below below:

    flowchart LR

        %%{init: {'flowchart' : {'curve' : 'linear'}}}%%
        %%{init: {'theme':'neutral'}}%%

        classDef circles fill:#DAE8FC
        classDef circlesGreen fill:#D5E8D4
        classDef subgraphBackground fill:#FFF2CC
        classDef roundRectangles fill:#FFE6CC
        classDef ifs fill:#F5F5F5

        inputData((Input Data)):::circles
        preprocessedData((Preprocessed Data)):::circles
        weightsArray((Weights Array)):::circlesGreen
        if{" "}:::ifs

        Preprocessing:::subgraphBackground

        subgraph Preprocessing
            direction LR  

            fourier(Fourier Analysis):::roundRectangles
            normalization(Windows Normalization):::roundRectangles
            label(Windows Labelling):::roundRectangles
            weights(Weights Calculation):::roundRectangles

            fourier --> normalization --> label --> weights

        end


        inputData --> Preprocessing
        Preprocessing -->|Generate| if
        if --> preprocessedData
        if --> weightsArray
Loading
    flowchart LR

    %%{init: {'flowchart' : {'curve' : 'linear'}}}%%
    %%{init: {'theme':'neutral'}}%%

    classDef circles fill:#DAE8FC
    classDef circlesGreen fill:#D5E8D4
    classDef rectangle fill:#E1D5E7
    classDef subgraphBackground fill:#FFF2CC
    classDef roundRectangles fill:#F8CECC
    classDef ifs fill:#F5F5F5

    preprocessedData1((Preprocessed Data)):::circles
    weightsArray((Weights Array)):::circlesGreen
    if{" "}:::ifs
    if2{" "}:::ifs
    prediction[Predictions]:::rectangle
    metrics[Learning Metrics]:::rectangle

    preprocessedData1 --> if
    weightsArray --> if

    NeuralNetwork:::subgraphBackground

    subgraph NeuralNetwork

            rnn("Recurrent Neural Network"):::roundRectangles
            esn(Echo State Network):::roundRectangles
            nbeats(NBetas Feature Extractor):::roundRectangles

        end

    if -->|Input| rnn
    if -->|Input| esn
    if -->|Input| nbeats

    rnn -->|Generate| if2
    esn -->|Generate| if2
    nbeats -->|Generate| if2

    if2 --> prediction
    if2 --> metrics
Loading

The code with the results can be found on folder main, and the folder experiments we have testing of new features.

Also it is important to mention that this code it doesn't contain exploratory functions or methods, here is just the execution of MoStress.

Although , contributions to add those exploratory features or any other contribution are really welcome, please, create your fork, branch and pull request to add it!

2. Prerequisites

Python Version Pip Version Tensorflow Version Sklearn Version Pandas Version Numpy Version ReservoirPy Version KerasBeats Version Matplotlib Version Seaborn Version

At this point in time, MoStress only support WESAD dataset, therefore, in order to run the code successfully, it is important to obtain the WESAD data in advance.

3. Datasets

Right now, there is only the implementation of the physiologic data collected from chest sensor of the WESAD. Thus, to help deal with the possible configuration of different datasets, we use the json configs/wesadDatasetOptions.json to add custom configurations and if you want to add new datasets, we suggest to do the same.

If you need to add a new dataset, just create the class which implement it and add the class call on the datasets/DatasetFactory.py and also make your class inherit the abstract class datasets/Dataset.py and don't forget to implement the _getData() method.

4. MoStress

4.1 Building your setup

MoStress has a builder (MoSB - MoStress Builder), which can be used for setup preparation. Make sure you have the properly versions for Python, Pip, and Python-Env. With MoStress Builder, your environment is prepared automatically, installing the necessary libraries and creating directories. To run the builder, move to the root directory of the downloaded project and execute one of the following commands:

./builder
sh builder
bash builder

After preparing the environment, just run the environment according to the your operational system. Do not forget to make sure the environment created (.mostress-env) is being used by Jupyter or VSCode.

Linux/MacOS

source $PWD.mostress-env/bin/activate

4.2 Preprocessing Step

4.2.1 Config File

On the folder configFiles, there is a configs/wesadDatasetOptions.json, where we set all the parameters needed on the preprocessing.

4.2.2 Implemented Steps

On moStress/preprocessing we have the main class MoStressPreprocessing.py, which has the call for all the steps that has to be done.

Each step were implemented as classes on moStress/preprocessing/implementedSteps folder, and they all extends the abstract class: moStress/preprocessing/implementedSteps/Steps.py, so if you want to implement a new step, please, inherit this class also.

4.3 Neural Network Architecture

The class MoStressNeuralNetwork basically is a wrapper which takes different models architecture, therefore, to add a new model, create a class which implements your architecture on models/architectures, make your new architecture inherits of the abstract class models/architectures/BaseArchitecture.py, and add the architecture call on models/architectures/ModelFactory.py.

We also advice that you pass a instance of the class MoStressNeuralNetwork as the constructor of your new model, since this class may already have utils properties to your model.

Currently we have the follow architectures implemented:

  • REGULARIZER-GRU,
  • REGULARIZER-LSTM,
  • BASELINE-GRU,
  • BASELINE-LSTM,
  • BASELINE-RESERVOIR
  • NBEATS-FEATURE-EXTRACTOR

The regularizer prefix means that we add a Ridge Regularizer and a Max Norm Bias on the output layer of the sequential model.

Also, after we train the model, we save it on models/saved

NOTE1: we don't save the reservoir model yet.

NOTE2: it is possible that the code of NBEATS-FEATURE-EXTRACTOR throw some error.

5 Improvements

  1. Implement a resume report after the data preprocessing;
  2. Implement the k-fold validation to improve the results;
  3. Check if it is possible to lower the number of for loops;
  4. Implement unit tests to ensure secure changes.

6. Future Work

  1. Change the RNN for reservoir computing or spiking neural networks;
  2. Use N-BEATS, with or without the preprocessing step.
  3. Use MoStress with new datasets

7. Publications

Here is our publications:

  1. A. de Souza, M. B. Melchiades, S. J. Rigo, G. d. O. Ramos, Mostress: a sequence model for stress classification, in: 2022 International Joint Conference on Neural Networks (IJCNN), IEEE, Padova, Italy, 2022, pp. 1–8. doi:10.1109/IJCNN55064.2022.9892953. Link

And here is how you can cite this research:

@INPROCEEDINGS{deSouza2022ijcnn,
  author={de Souza, Arturo and Melchiades, Mateus B. and Rigo, Sandro J. and Ramos, Gabriel de O.},
  title = {MoStress: a Sequence Model for Stress Classification},
  booktitle = {2022 International Joint Conference on Neural Networks (IJCNN)},
  year = {2022},
  address = {Padova, Italy},
  month = {July},
  publisher = {IEEE},
  pages={1-8},
  doi={10.1109/IJCNN55064.2022.9892953}
}

Appendix

A. Code Architecture

classDiagram
  class DatasetFactory {
    +make(datasetName~string~, datasetPath~string~, datasetOptions~string~)
  }
  class  Dataset{
    <<abstract>>
    loadData(dataPath~string~, dataEncoding~string~ = 'latin1')$
    saveData(dataPath~string~, data~any~)$
    #_adjustUnnecessaryLabelCode(labelData~any~, oldIndex~int~, newIndex~int~)
    #_getData()*
  }
  class WesadPhysioChest {
    +String dataPath
    +Dictionary dataOpts
    +Any data

    -_getLabel(labelCode~int|string~) String
    -_getData() Any
  }
  class Steps {
    <<abstract>>
    +execute()*
  }
  class Normalization {
    +MoStressPreprocessing moStressPreprocessing
    -Boolean _hasNormalizationFinished

    -_rollingZScore(df~pd.DataFrame~)
    +execute()
  }
  class WindowsLabelling {
    +MoStressPreprocessing moStressPreprocessing
    -Boolean _hasWindowsLabellingFinished

    -_getElementArrayFrequency(npArray~np.Array~)
    -_labellingWindows(df~pd.DataFrame~, labels~List~)
    +execute()
  }
  class WeightsCalculation {
    +MoStressPreprocessing moStressPreprocessing
    -Boolean _hasWeightsCalculationFinished

    -_getValidationData()
    -_getTrainingData()
    -_getWeights()
    +execute()
  }
  class MoStressPreprocessing {
    +List features
    +List targets
    +List featuresModelValidation
    +List targetsModelValidation
    +Dict quantityOfSets
    +Dict discardedWindowsCounter

    -_applyNormalization()
    -_applyWindowsLabelling()
    -_applyWeightsCalculation()
    +execute()
  }

  Normalization ..> Steps: inherits
  WindowsLabelling ..> Steps: inherits
  WeightsCalculation ..> Steps: inherits

  MoStressPreprocessing *-- Normalization: compose
  MoStressPreprocessing *-- WindowsLabelling: compose
  MoStressPreprocessing *-- WeightsCalculation: compose

  DatasetFactory ..|> MoStressPreprocessing: input

  WesadPhysioChest ..> Dataset: inherits
  DatasetFactory --|> WesadPhysioChest: make

  class ModelFactory {
    +make(architectureName~string~, architectureOptions~dictionary~)
  }
  class BaseArchitecture {
    <<abstract>>     
    -_compileModel()*
    -_setModelCallbacks()*
    -_fitModel()*
    -_makePredictions()*
    -_saveModel()*
  }
  class SequentialModels {
    +MoStressNeuralNetwork moStressNeuralNetwork

    +gruRegularizerMoStress()
    +lstmRegularizerMoStress()
    +gruBaselineMostress()
    +lstmBaselineMostress()
  }
  class ReservoirModels {
    +MoStressNeuralNetwork moStressNeuralNetwork

    +baseline()
  }
  class EvaluateModel {
    +Any model
    +String modelName
    +List featuresValidation
    +List targetsValidation
    +List classes

    -_makePredictions()
    +getClassesPredicted()
    -_getConfusionMatrix()
    +getConfusionMatrix()
    +printConfusionMatrix()
    +executeEvaluation()
  }
  class MoStressNeuralNetwork {
    +Any modelOpts
    +Dict weights
    -Tensor _xTrain
    -Tensor _xTest
    -Tensor _yTrain
    -Tensor _xTest
    -Int _winSize
    -Int _numFeatures
    -Int _numClasses
    -String _modelName
    -String _optimizerName
    -List _callbacks
    +Any model
    +Dict history
    +String modelFullName
    -Int _reservoirRandomSeed
    -Int _reservoirVerbosityState
    -List _allTrainFeatures
    -List _allTrainTargets

    +execute()
  }

ModelFactory --|> SequentialModels:make
ModelFactory --|> ReservoirModels:make
MoStressNeuralNetwork *-- ModelFactory: compose
MoStressNeuralNetwork *-- EvaluateModel: evaluate
MoStressPreprocessing --|> MoStressNeuralNetwork: input
SequentialModels ..> BaseArchitecture: inherits
ReservoirModels ..> BaseArchitecture: inherits

  class GenericBlock {
    +Number lookback
    +Number horizon
    +Number num_neurons
    +Number block_layers

    +call(inputs)
  }
  class TrendBlock {
    +Number lookback
    +Number horizon
    +Number num_neurons
    +Number block_layers
    +Number polynomial_term

    +call(inputs)
  }
  class SeasonalBlock {
    +Number lookback
    +Number horizon
    +Number num_neurons
    +Number block_layers
    +Number num_harmonics

    +call(inputs)
  }
  class NBeatsBlock {
    +String model_type          
    +Number lookback            
    +Number horizon             
    +Number num_generic_neurons 
    +Number num_generic_stacks  
    +Number num_generic_layers  
    +Number num_trend_neurons   
    +Number num_trend_stacks    
    +Number num_trend_layers    
    +Number num_seasonal_neurons
    +Number num_seasonal_stacks 
    +Number num_seasonal_layers 
    +Number num_harmonics       
    +Number polynomial_term     

    +get_config()
    +call(inputs)
  }
  class NBeatsModel {
    +String model_type          
    +Number lookback            
    +Number horizon             
    +Number num_generic_neurons 
    +Number num_generic_stacks  
    +Number num_generic_layers  
    +Number num_trend_neurons   
    +Number num_trend_stacks    
    +Number num_trend_layers    
    +Number num_seasonal_neurons
    +Number num_seasonal_stacks 
    +Number num_seasonal_layers 
    +Number num_harmonics       
    +Number polynomial_term    
    +String loss
    +Number learning_rate
    +Number batch_size

    +get_residual_model()
    +get_config()
    +build_layer()
    +build_model()
    +compile_model()
    +fit()
    +predict()
    +evaluate()
  }
  class NBeatsFeatureExtractor {
    +MoStressNeuralNetwork moStressNeuralNetwork
    +NBeatsModel nBeats
    +String residualsFolderPath
    +String residualsTrainingFolderPath
    +String residualsValidationFolderPath
    +String nBeatsSavedModelBasePath
    -List callbacks
    +Boolean hasToCollectResiduals
    +Boolean hasToFitBasicClassifier

    +classificationModel()
    -_setModelCallbacksPath()
    +prepareResidualsDataset(epochs, verbose)
    -_trainNbeatsModel(epochs, verbose)
    -_generateResiduals()
    #collectValidationResiduals()
  }

  NBeatsBlock *-- GenericBlock: compose
  NBeatsBlock *-- TrendBlock: compose
  NBeatsBlock *-- SeasonalBlock: compose
  NBeatsModel *-- NBeatsBlock: compose
  NBeatsFeatureExtractor *-- NBeatsModel: compose
  NBeatsFeatureExtractor ..> BaseArchitecture: inherits
  ModelFactory --|> NBeatsFeatureExtractor:make


Loading

License

This project uses the following license: MIT.

Releases

No releases published

Packages

No packages published