Example #1
0
        public double Classify(object dataSource)
        {
            var enumerable = ClassificationModel
                             .Properties
                             .Select(p => p.Scaler[Convert.ToDouble(PatientPropertyProvider.GetPropertyValue(dataSource, p.Name))] * p.Coefficient)
                             .ToList();

            var intermediateResult = enumerable.Sum() + ClassificationModel.FreeCoefficient;

            return(Math.Exp(intermediateResult) / (1 + Math.Exp(intermediateResult)));
        }
        public double CalculateOptimalCutOff(ClassificationModel model, List <Patient> patients)
        {
            //ToDo: filter patients which could be used for current model
            var applicablePatients = patients.Where(patient =>
            {
                var dataSource = new { Patient = patient, PatientVisit = patient.LastVisit };
                return(model.Properties.All(p => PatientPropertyProvider.GetPropertyValue(dataSource, p.Name) != null));
            })
                                     .ToList();

            foreach (var p in applicablePatients)
            {
                var lastVisitData = p.LastVisit;
                Console.Write($@"{p.Surname} bmi {lastVisitData.ObesityBmi} waist {lastVisitData.ObesityWaistCircumference} gene {p.AgtAgtr2}  gender {p.Gender} phiz {lastVisitData.PhysicalActivity} maleHered {p.MaleHeredity}; \n");
            }

            var classificator = new PatientClassificator(model);

            var healthyCorrect   = 0;
            var illCorrect       = applicablePatients.Count(p => p.HypertensionStage != HypertensionStage.Healthy);
            var healthyIncorrect = applicablePatients.Count(p => p.HypertensionStage == HypertensionStage.Healthy);
            var illIncorrect     = 0;

            var basePatientStageScoresGroupped =
                applicablePatients.Select(patient => new
            {
                Patient = patient,
                Score   = classificator.Classify(new { Patient = patient, PatientVisit = patient.LastVisit })
            })
                .GroupBy(patientScore => patientScore.Score)
                .OrderBy(scoreGroup => scoreGroup.Key)
                .ToList();

            var cutoffs = new List <CutOffValues>
            {
                new CutOffValues
                {
                    HealthyCorrect   = healthyCorrect,
                    HealthyIncorrect = healthyIncorrect,
                    IllCorrect       = illCorrect,
                    IllIncorrect     = illIncorrect,
                    Score            = 0
                }
            };


            foreach (var scoreGroup in basePatientStageScoresGroupped)
            {
                foreach (var patientScore in scoreGroup)
                {
                    if (patientScore.Patient.HypertensionStage == HypertensionStage.Healthy)
                    {
                        healthyCorrect++;
                        healthyIncorrect--;
                    }
                    else
                    {
                        illCorrect--;
                        illIncorrect++;
                    }
                }

                cutoffs.Add(new CutOffValues
                {
                    HealthyCorrect   = healthyCorrect,
                    HealthyIncorrect = healthyIncorrect,
                    IllCorrect       = illCorrect,
                    IllIncorrect     = illIncorrect,
                    Score            = scoreGroup.Key
                }
                            );
            }
            var cutoff = cutoffs.OrderByDescending(c => c.SumPercent).ThenBy(c => c.PercentDifference);
//            foreach (var coutoffTest in cutoffs)
//            {
//                Console.Write(coutoffTest.SumPercent + " " + coutoffTest.HealthPercent + " " + coutoffTest.IllPercent + "; \n");
//            }
            var cutoff2 = cutoffs.OrderBy(c => c.PercentDifference).ThenByDescending(c => c.SumPercent);

            return(cutoff.FirstOrDefault().Score);
        }