/* file: physionet2012.c	G. Moody	22 March 2012

Sample entry for the PhysioNet/CinC Challenge 2012

This sample requires saps_score.c.  Compile it using gcc by
    gcc -o physionet2012 physionet2012.c saps_score.c
*/

#include <stdio.h>
#include <string.h>
#define PLEN 16	/* max length of data type string, including trailing zero */

struct data {
    int time;   	/* minutes since admission */
    char p[PLEN];	/* type of data */
    double val; 	/* parameter value */
};

/* Risk of in-hospital death as a function of calculated SAPS score [0-56].
   These values are *examples*;  they are *not* good estimates of risk! */
double risk[] = {			/*  SAPS */
    .1, .1, .1, .1, .1, .1, .1, .1,	/*  0- 7 */
    .1, .1, .1, .1, .1, .1, .1, .1,	/*  8-15 */
    .1, .1, .2, .2, .3, .3, .3, .4,	/* 16-23 */
    .4, .4, .5, .5, .5, .6, .6, .6,	/* 24-31 */
    .7, .7, .7, .8, .8, .9, .9, .9,	/* 32-39 */
    .9, .9, .9, .9, .9, .9, .9, .9,	/* 40-47 */
    .9, .9, .9, .9, .9, .9, .9, .9, .9	/* 48-56 */
};

main(int argc, char **argv)
{
    char buf[128], RecordID[20], pname[PLEN], *s, *t;
    double value;
    int h, m, n = 0, saps, survival;
    static struct data Data[2000], *d = Data;

    /* Look for the RecordID. */
    while (fgets(buf, sizeof(buf), stdin)) {
	if (strncmp(buf, "00:00,RecordID,", 15) == 0) {
	    strncpy(RecordID, buf+15, sizeof(RecordID));
	    RecordID[strlen(RecordID)-1] = '\0';  /* discard trailing newline */
	    break;
	}
    }

    /* Read the remainder of the input, parse it, store into Data */
    while (fgets(buf, sizeof(buf), stdin)) {
	for (s = buf; *s && *s != ','; s++)
	    ;
	if (*s == '\0') continue; /* skip this line if no comma found */
	*s++ = '\0';
	if (sscanf(buf, "%d:%d", &h, &m) == 2) d->time = 60*h + m;
	else continue;		  /* skip if first field not in hh:mm form */
	for (t = s; *t && *t != ','; t++)
	    ;
	if (*t == '\0') continue; /* skip if second comma is missing */
	*t++ = '\0';
	strncpy(d->p, s, PLEN);
	if (sscanf(t, "%lf", &value) == 1) {
	    d->val = value;
	    d++;
	    n++;
	}
    }

    /* Analyse the data.  This sample entry calculates SAPS-I from the data. */
    saps = saps_score(&Data, n);
    
    /* The saps_score function calculates a score even if one or more SAPS
       parameters are missing, but returns a negative score in this case.
       Here, we treat a negative (truncated) score as if it were positive. */
   if (saps < 0) saps = -saps;
 
   /* This sample entry predicts death (survival = 1) if SAPS is 20 or more. */
   survival = (saps >= 20);

   /* It uses risk[saps] to estimate risk for event 2. */
   printf("%s,%d,%g\n", RecordID, survival, risk[saps]);
}