METHOD AND SYSTEM FOR INDIVIDUAL DEMAND FORECASTING

1 . A method for individual forecasting of a future event for a subject using historical data, the historical data comprising a plurality of historical events associated with the subject, the method executed on at least one processing unit, the method comprising:

receiving the historical data associated with the subject;

determining a random variable representing a remaining time until the future event;

predicting a time to the future event using a distribution function that is determined using a recurrent neural network, the distribution function comprising a learned density with peaks that approximate the times of the historical events in the historical data;

determining a log-likelihood function based on a probability that the random variable exceeds an amount of time remaining until a next historical event in the historical data and parameterized by the distribution function; and

outputting a forecast of a time to the future event as the log-likelihood function.

2 . The method of claim 1 , wherein a loss function for the recurrent neural network comprises a negative of the log-likelihood function.

3 . The method of claim 1 , wherein the random variable is conditioned based on inter-arrival times of the historical events in the historical data.

4 . The method of claim 1 , wherein the random variable is conditioned based on excess times since arrival of preceding historical events in the historical data.

5 . The method of claim 1 , wherein the log-likelihood function at each time is the log of the probability that the random variable is in the set of time until the next historical event when the next historical event has been observed, and the log of the survival function otherwise.

6 . The method of claim 5 , wherein the distribution function follows a Weibull distribution.

7 . The method of claim 6 , wherein the distribution function is determined as (k/λ)((s+t)/λ) k−1 S W (t), where k is the shape of the Weibull distribution, λ is the scale of the Weibull distribution, t is the time-step, and S W (t) is the survival function.

8 . The method of claim 1 , wherein outputting the forecast of the time to the future event as the log-likelihood function comprises determining a sum of log-likelihoods at each time-step.

9 . The method of claim 8 , further comprising transforming the sum of log-likelihoods as a function of recurrent neural network parameters and historical data, and determining a minimizer of an overall observed loss of the recurrent neural network using such function.

10 . The method of claim 1 , further comprising outputting derivative values of the log-likelihood function.

11 . A system for individual forecasting of a future event for a subject using historical data, the historical data comprising a plurality of historical events associated with the subject, the system comprising one or more processors in communication with a data storage, the one or more processors configurable to execute:

a data acquisition module to receive the historical data associated with the subject;

a conditional excess module to determine a random variable representing a remaining time until the future event;

a machine learning module 120 to predict a time to the future event using a distribution function that is determined using a recurrent neural network, the distribution function comprising a learned density with peaks that approximate the times of the historical events in the historical data; and

a forecasting module to determine a log-likelihood function based on a probability that the random variable exceeds an amount of time remaining until a next historical event in the historical data and parameterized by the distribution function, and to output a forecast of a time to the future event as the log-likelihood function.

12 . The system of claim 11 , wherein a loss function for the recurrent neural network comprises a negative of the log-likelihood function.

13 . The system of claim 11 , wherein the random variable is conditioned based on inter-arrival times of the historical events in the historical data.

14 . The system of claim 11 , wherein the random variable is conditioned based on excess times since arrival of preceding historical events in the historical data.

15 . The system of claim 11 , wherein the log-likelihood function at each time is the log of the probability that the random variable is in the set of time until the next historical event when the next historical event has been observed, and the log of the survival function otherwise.

16 . The system of claim 15 , wherein the distribution function follows a Weibull distribution.

17 . The system of claim 16 , wherein the distribution function is determined as (k/λ)((s+t)/λ) k−1 S W (t), where k is the shape of the Weibull distribution, λ is the scale of the Weibull distribution, t is the time-step, and S W (t) is the survival function.

18 . The system of claim 11 , wherein outputting the forecast of the time to the future event as the log-likelihood function comprises determining a sum of log-likelihoods at each time-step.

19 . The system of claim 18 , wherein the forecasting module further transforms the sum of log-likelihoods as a function of recurrent neural network parameters and historical data, and determining a minimizer of an overall observed loss of the recurrent neural network using such function.

20 . The system of claim 11 , wherein the forecasting module further outputs derivative values of the log-likelihood function.