Fix R2 calculation

1 files changed, 8 insertions, 3 deletions

diff --git a/src/statistics.c b/src/statistics.c
index 50bfd513..84453ffa 100644
--- a/src/statistics.c
+++ b/src/statistics.c
@@ -59,13 +59,18 @@ double stat_r2(double *obs, double *calc, unsigned int *c, int size,
 	for ( i=1; i<size; i++ ) {
 
 		if ( c[i] > 0 ) {
+
 			double obsi;
+
 			obsi = obs[i] / (double)c[i];
-			top += pow(fabs(obsi - scale*calc[i]), 2.0);
+			obsi = obsi / scale;
+
+			top += pow(fabs(obsi - calc[i]), 2.0);
 			bot += pow(obsi, 2.0);
-		}
 
-	} /* else reflection not measured so don't worry about it */
+		} /* else reflection not measured so don't worry about it */
+
+	}
 
 	return sqrt(top/bot);
 }