using System;

using MathNet.Numerics;

using System.Collections.Generic;

using System.Linq;

public class Program

{

    public static void Main()

    {

        Console.WriteLine("Calculating linear regression");

        List<(double, DateTime)> collatedEgfrMeasurements = new List<(double, DateTime)>()

        {

            /*

                PSEUDOCODE REQUIREMENTS

                * At least 3 egfr measurements in the last 12 months

                * All measurements in past 24 months will be included

            */

            //Doesn't matter which order the measurements are specified but order from oldest first to most recent last for clarity

            //ADD DATES HERE

            (100, DateTime.Parse("2023-12-13T11:52:18")),

            (85, DateTime.Parse("2024-09-12T11:52:18")),

            (65, DateTime.Parse("2024-12-11T11:52:18"))     

        };

        //collatedEgfrMeasurements will have all egfr measurements in the last 2 years from DateTime.Now

        double[] egfrMeasurements = collatedEgfrMeasurements.Select(x => x.Item1).ToArray();

        List<DateTime> egfrDates = collatedEgfrMeasurements.Select(x => x.Item2).ToList();

​

        //Gets the number of days since two years ago (from DateTime.Now) and converts to a decimal year

        DateTime twentyFourMonthsAgo = DateTime.Now.AddYears(-2);

        double[] elapsedDatesInYears = egfrDates.Select(x => Convert.ToDouble(((x - twentyFourMonthsAgo).TotalDays)/365)).ToArray();

​

        //Do a linear regression to calculate the slope

        (double A, double B) regression = Fit.Line(elapsedDatesInYears, egfrMeasurements);

        double intercept = regression.A;

        double slope = regression.B;

        Console.WriteLine("Intercept is: " + intercept);

        Console.WriteLine("Slope is: " + slope + " ml/min/year");

    }

}