Fix percentile calculation

Author: ishii-norimi

lib/model/kernel_density_estimator.js

@@ -55,7 +55,7 @@ export default class KernelDensityEstimator {
 			const np = n * p
 			const np_l = Math.floor(np)
 			const np_h = Math.ceil(np)
-			return k[np_l] + (np_h - np_l) * (k[np_h] - k[np_l])
+			return k[np_l] + (np - np_l) * (k[np_h] - k[np_l])
 		}
 		const sgm = Math.min(std, (q(0.75) - q(0.25)) / 1.34)
 

lib/model/nadaraya_watson.js

@@ -38,7 +38,7 @@ export default class NadarayaWatson {
 				const np = n * p
 				const np_l = Math.floor(np)
 				const np_h = Math.ceil(np)
-				return k[np_l] + (np_h - np_l) * (k[np_h] - k[np_l])
+				return k[np_l] + (np - np_l) * (k[np_h] - k[np_l])
 			}
 			const sgm = Math.min(std, (q(0.75) - q(0.25)) / 1.34)
 

lib/model/priestley_chao.js

@@ -37,7 +37,7 @@ export default class PriestleyChao {
 				const np = n * p
 				const np_l = Math.floor(np)
 				const np_h = Math.ceil(np)
-				return k[np_l] + (np_h - np_l) * (k[np_h] - k[np_l])
+				return k[np_l] + (np - np_l) * (k[np_h] - k[np_l])
 			}
 			const sgm = Math.min(std, (q(0.75) - q(0.25)) / 1.34)
 

lib/model/tukeys_fences.js

@@ -26,7 +26,7 @@ export default class TukeysFences {
 				const np = n * p
 				const np_l = Math.floor(np)
 				const np_h = Math.ceil(np)
-				return x[np_l] + (np_h - np_l) * (x[np_h] - x[np_l])
+				return x[np_l] + (np - np_l) * (x[np_h] - x[np_l])
 			}
 
 			const q1 = q(0.25)