Need to figure out SIS parameters for comparison, remove matrix inverse (use factorization instead) for faster run, and then ready to run. Still to add analysis steps (Kendall tau, etc) but that can be run separately.