stripchart 0.0.3

Sign up to get free protection for your applications and to get access to all the features.
Files changed (504) hide show
  1. checksums.yaml +7 -0
  2. data/LICENSE +21 -0
  3. data/README.md +21 -0
  4. data/exe/stripchart +5 -0
  5. data/lib/stripchart/app.rb +71 -0
  6. data/lib/stripchart/public/app.css +37 -0
  7. data/lib/stripchart/public/app.js +167 -0
  8. data/lib/stripchart/public/components/d3/LICENSE +26 -0
  9. data/lib/stripchart/public/components/d3/Makefile +45 -0
  10. data/lib/stripchart/public/components/d3/README.md +7 -0
  11. data/lib/stripchart/public/components/d3/bower.json +25 -0
  12. data/lib/stripchart/public/components/d3/component.json +30 -0
  13. data/lib/stripchart/public/components/d3/d3.js +8798 -0
  14. data/lib/stripchart/public/components/d3/d3.min.js +5 -0
  15. data/lib/stripchart/public/components/d3/index-browserify.js +3 -0
  16. data/lib/stripchart/public/components/d3/index.js +15 -0
  17. data/lib/stripchart/public/components/d3/lib/colorbrewer/LICENSE +38 -0
  18. data/lib/stripchart/public/components/d3/lib/colorbrewer/colorbrewer.css +1690 -0
  19. data/lib/stripchart/public/components/d3/lib/colorbrewer/colorbrewer.js +302 -0
  20. data/lib/stripchart/public/components/d3/lib/geographiclib/LICENSE +23 -0
  21. data/lib/stripchart/public/components/d3/lib/jit/LICENSE +27 -0
  22. data/lib/stripchart/public/components/d3/lib/penner/LICENSE +31 -0
  23. data/lib/stripchart/public/components/d3/lib/polymaps/LICENSE +28 -0
  24. data/lib/stripchart/public/components/d3/lib/protovis/LICENSE +27 -0
  25. data/lib/stripchart/public/components/d3/lib/science/LICENSE +26 -0
  26. data/lib/stripchart/public/components/d3/lib/science/science.js +225 -0
  27. data/lib/stripchart/public/components/d3/lib/science/science.lin.js +27 -0
  28. data/lib/stripchart/public/components/d3/lib/science/science.lin.min.js +1 -0
  29. data/lib/stripchart/public/components/d3/lib/science/science.min.js +1 -0
  30. data/lib/stripchart/public/components/d3/lib/science/science.stats.js +720 -0
  31. data/lib/stripchart/public/components/d3/lib/science/science.stats.min.js +1 -0
  32. data/lib/stripchart/public/components/d3/package.json +53 -0
  33. data/lib/stripchart/public/components/d3/src/arrays/ascending.js +3 -0
  34. data/lib/stripchart/public/components/d3/src/arrays/bisect.js +28 -0
  35. data/lib/stripchart/public/components/d3/src/arrays/descending.js +3 -0
  36. data/lib/stripchart/public/components/d3/src/arrays/entries.js +5 -0
  37. data/lib/stripchart/public/components/d3/src/arrays/extent.js +21 -0
  38. data/lib/stripchart/public/components/d3/src/arrays/index.js +22 -0
  39. data/lib/stripchart/public/components/d3/src/arrays/keys.js +5 -0
  40. data/lib/stripchart/public/components/d3/src/arrays/map.js +50 -0
  41. data/lib/stripchart/public/components/d3/src/arrays/max.js +14 -0
  42. data/lib/stripchart/public/components/d3/src/arrays/mean.js +15 -0
  43. data/lib/stripchart/public/components/d3/src/arrays/median.js +9 -0
  44. data/lib/stripchart/public/components/d3/src/arrays/merge.js +3 -0
  45. data/lib/stripchart/public/components/d3/src/arrays/min.js +14 -0
  46. data/lib/stripchart/public/components/d3/src/arrays/nest.js +97 -0
  47. data/lib/stripchart/public/components/d3/src/arrays/permute.js +7 -0
  48. data/lib/stripchart/public/components/d3/src/arrays/quantile.js +8 -0
  49. data/lib/stripchart/public/components/d3/src/arrays/range.js +24 -0
  50. data/lib/stripchart/public/components/d3/src/arrays/set.js +38 -0
  51. data/lib/stripchart/public/components/d3/src/arrays/shuffle.js +8 -0
  52. data/lib/stripchart/public/components/d3/src/arrays/sum.js +14 -0
  53. data/lib/stripchart/public/components/d3/src/arrays/transpose.js +5 -0
  54. data/lib/stripchart/public/components/d3/src/arrays/values.js +5 -0
  55. data/lib/stripchart/public/components/d3/src/arrays/zip.js +15 -0
  56. data/lib/stripchart/public/components/d3/src/behavior/behavior.js +1 -0
  57. data/lib/stripchart/public/components/d3/src/behavior/drag.js +79 -0
  58. data/lib/stripchart/public/components/d3/src/behavior/index.js +3 -0
  59. data/lib/stripchart/public/components/d3/src/behavior/zoom.js +209 -0
  60. data/lib/stripchart/public/components/d3/src/color/color.js +5 -0
  61. data/lib/stripchart/public/components/d3/src/color/hcl.js +42 -0
  62. data/lib/stripchart/public/components/d3/src/color/hsl.js +64 -0
  63. data/lib/stripchart/public/components/d3/src/color/index.js +6 -0
  64. data/lib/stripchart/public/components/d3/src/color/lab.js +68 -0
  65. data/lib/stripchart/public/components/d3/src/color/rgb.js +307 -0
  66. data/lib/stripchart/public/components/d3/src/color/xyz.js +7 -0
  67. data/lib/stripchart/public/components/d3/src/compat/date.js +3 -0
  68. data/lib/stripchart/public/components/d3/src/compat/index.js +2 -0
  69. data/lib/stripchart/public/components/d3/src/compat/style.js +20 -0
  70. data/lib/stripchart/public/components/d3/src/core/array.js +19 -0
  71. data/lib/stripchart/public/components/d3/src/core/class.js +12 -0
  72. data/lib/stripchart/public/components/d3/src/core/document.js +3 -0
  73. data/lib/stripchart/public/components/d3/src/core/functor.js +5 -0
  74. data/lib/stripchart/public/components/d3/src/core/identity.js +3 -0
  75. data/lib/stripchart/public/components/d3/src/core/index.js +3 -0
  76. data/lib/stripchart/public/components/d3/src/core/noop.js +1 -0
  77. data/lib/stripchart/public/components/d3/src/core/ns.js +22 -0
  78. data/lib/stripchart/public/components/d3/src/core/rebind.js +16 -0
  79. data/lib/stripchart/public/components/d3/src/core/source.js +3 -0
  80. data/lib/stripchart/public/components/d3/src/core/subclass.js +11 -0
  81. data/lib/stripchart/public/components/d3/src/core/target.js +3 -0
  82. data/lib/stripchart/public/components/d3/src/core/true.js +3 -0
  83. data/lib/stripchart/public/components/d3/src/core/vendor.js +12 -0
  84. data/lib/stripchart/public/components/d3/src/d3.js +23 -0
  85. data/lib/stripchart/public/components/d3/src/dsv/csv.js +3 -0
  86. data/lib/stripchart/public/components/d3/src/dsv/dsv.js +136 -0
  87. data/lib/stripchart/public/components/d3/src/dsv/index.js +3 -0
  88. data/lib/stripchart/public/components/d3/src/dsv/tsv.js +3 -0
  89. data/lib/stripchart/public/components/d3/src/end.js +2 -0
  90. data/lib/stripchart/public/components/d3/src/event/dispatch.js +69 -0
  91. data/lib/stripchart/public/components/d3/src/event/drag.js +27 -0
  92. data/lib/stripchart/public/components/d3/src/event/event.js +50 -0
  93. data/lib/stripchart/public/components/d3/src/event/index.js +5 -0
  94. data/lib/stripchart/public/components/d3/src/event/mouse.js +39 -0
  95. data/lib/stripchart/public/components/d3/src/event/timer.js +82 -0
  96. data/lib/stripchart/public/components/d3/src/event/touches.js +12 -0
  97. data/lib/stripchart/public/components/d3/src/format/collapse.js +3 -0
  98. data/lib/stripchart/public/components/d3/src/format/format-locale.js +4 -0
  99. data/lib/stripchart/public/components/d3/src/format/format-localized.js +5 -0
  100. data/lib/stripchart/public/components/d3/src/format/format.js +144 -0
  101. data/lib/stripchart/public/components/d3/src/format/formatPrefix.js +22 -0
  102. data/lib/stripchart/public/components/d3/src/format/index.js +4 -0
  103. data/lib/stripchart/public/components/d3/src/format/requote.js +5 -0
  104. data/lib/stripchart/public/components/d3/src/format/round.js +5 -0
  105. data/lib/stripchart/public/components/d3/src/geo/albers-usa.js +129 -0
  106. data/lib/stripchart/public/components/d3/src/geo/albers.js +11 -0
  107. data/lib/stripchart/public/components/d3/src/geo/area.js +67 -0
  108. data/lib/stripchart/public/components/d3/src/geo/azimuthal-equal-area.js +12 -0
  109. data/lib/stripchart/public/components/d3/src/geo/azimuthal-equidistant.js +13 -0
  110. data/lib/stripchart/public/components/d3/src/geo/azimuthal.js +25 -0
  111. data/lib/stripchart/public/components/d3/src/geo/bounds.js +169 -0
  112. data/lib/stripchart/public/components/d3/src/geo/cartesian.js +47 -0
  113. data/lib/stripchart/public/components/d3/src/geo/centroid.js +149 -0
  114. data/lib/stripchart/public/components/d3/src/geo/circle.js +81 -0
  115. data/lib/stripchart/public/components/d3/src/geo/clip-antimeridian.js +101 -0
  116. data/lib/stripchart/public/components/d3/src/geo/clip-circle.js +182 -0
  117. data/lib/stripchart/public/components/d3/src/geo/clip-polygon.js +94 -0
  118. data/lib/stripchart/public/components/d3/src/geo/clip-view.js +212 -0
  119. data/lib/stripchart/public/components/d3/src/geo/clip.js +135 -0
  120. data/lib/stripchart/public/components/d3/src/geo/compose.js +12 -0
  121. data/lib/stripchart/public/components/d3/src/geo/conic-conformal.js +36 -0
  122. data/lib/stripchart/public/components/d3/src/geo/conic-equal-area.js +33 -0
  123. data/lib/stripchart/public/components/d3/src/geo/conic-equidistant.js +35 -0
  124. data/lib/stripchart/public/components/d3/src/geo/conic.js +16 -0
  125. data/lib/stripchart/public/components/d3/src/geo/distance.js +13 -0
  126. data/lib/stripchart/public/components/d3/src/geo/equirectangular.js +10 -0
  127. data/lib/stripchart/public/components/d3/src/geo/geo.js +1 -0
  128. data/lib/stripchart/public/components/d3/src/geo/gnomonic.js +12 -0
  129. data/lib/stripchart/public/components/d3/src/geo/graticule.js +102 -0
  130. data/lib/stripchart/public/components/d3/src/geo/greatArc.js +42 -0
  131. data/lib/stripchart/public/components/d3/src/geo/index.js +32 -0
  132. data/lib/stripchart/public/components/d3/src/geo/interpolate.js +38 -0
  133. data/lib/stripchart/public/components/d3/src/geo/length.js +43 -0
  134. data/lib/stripchart/public/components/d3/src/geo/mercator.js +48 -0
  135. data/lib/stripchart/public/components/d3/src/geo/orthographic.js +12 -0
  136. data/lib/stripchart/public/components/d3/src/geo/path-area.js +40 -0
  137. data/lib/stripchart/public/components/d3/src/geo/path-bounds.js +21 -0
  138. data/lib/stripchart/public/components/d3/src/geo/path-buffer.js +59 -0
  139. data/lib/stripchart/public/components/d3/src/geo/path-centroid.js +78 -0
  140. data/lib/stripchart/public/components/d3/src/geo/path-context.js +49 -0
  141. data/lib/stripchart/public/components/d3/src/geo/path.js +97 -0
  142. data/lib/stripchart/public/components/d3/src/geo/point-in-polygon.js +71 -0
  143. data/lib/stripchart/public/components/d3/src/geo/projection.js +132 -0
  144. data/lib/stripchart/public/components/d3/src/geo/resample.js +98 -0
  145. data/lib/stripchart/public/components/d3/src/geo/rotation.js +72 -0
  146. data/lib/stripchart/public/components/d3/src/geo/spherical.js +12 -0
  147. data/lib/stripchart/public/components/d3/src/geo/stereographic.js +12 -0
  148. data/lib/stripchart/public/components/d3/src/geo/stream.js +71 -0
  149. data/lib/stripchart/public/components/d3/src/geo/transverse-mercator.js +23 -0
  150. data/lib/stripchart/public/components/d3/src/geom/delaunay.js +31 -0
  151. data/lib/stripchart/public/components/d3/src/geom/geom.js +1 -0
  152. data/lib/stripchart/public/components/d3/src/geom/hull.js +120 -0
  153. data/lib/stripchart/public/components/d3/src/geom/index.js +6 -0
  154. data/lib/stripchart/public/components/d3/src/geom/polygon.js +105 -0
  155. data/lib/stripchart/public/components/d3/src/geom/quadtree.js +190 -0
  156. data/lib/stripchart/public/components/d3/src/geom/voronoi.js +545 -0
  157. data/lib/stripchart/public/components/d3/src/interpolate/array.js +19 -0
  158. data/lib/stripchart/public/components/d3/src/interpolate/ease.js +110 -0
  159. data/lib/stripchart/public/components/d3/src/interpolate/hcl.js +23 -0
  160. data/lib/stripchart/public/components/d3/src/interpolate/hsl.js +23 -0
  161. data/lib/stripchart/public/components/d3/src/interpolate/index.js +13 -0
  162. data/lib/stripchart/public/components/d3/src/interpolate/interpolate.js +25 -0
  163. data/lib/stripchart/public/components/d3/src/interpolate/lab.js +20 -0
  164. data/lib/stripchart/public/components/d3/src/interpolate/number.js +6 -0
  165. data/lib/stripchart/public/components/d3/src/interpolate/object.js +25 -0
  166. data/lib/stripchart/public/components/d3/src/interpolate/rgb.js +20 -0
  167. data/lib/stripchart/public/components/d3/src/interpolate/round.js +6 -0
  168. data/lib/stripchart/public/components/d3/src/interpolate/string.js +88 -0
  169. data/lib/stripchart/public/components/d3/src/interpolate/transform.js +56 -0
  170. data/lib/stripchart/public/components/d3/src/interpolate/uninterpolate.js +9 -0
  171. data/lib/stripchart/public/components/d3/src/layout/bundle.js +59 -0
  172. data/lib/stripchart/public/components/d3/src/layout/chord.js +157 -0
  173. data/lib/stripchart/public/components/d3/src/layout/cluster.js +92 -0
  174. data/lib/stripchart/public/components/d3/src/layout/force.js +359 -0
  175. data/lib/stripchart/public/components/d3/src/layout/hierarchy.js +116 -0
  176. data/lib/stripchart/public/components/d3/src/layout/histogram.js +110 -0
  177. data/lib/stripchart/public/components/d3/src/layout/index.js +13 -0
  178. data/lib/stripchart/public/components/d3/src/layout/layout.js +1 -0
  179. data/lib/stripchart/public/components/d3/src/layout/pack.js +212 -0
  180. data/lib/stripchart/public/components/d3/src/layout/partition.js +51 -0
  181. data/lib/stripchart/public/components/d3/src/layout/pie.js +99 -0
  182. data/lib/stripchart/public/components/d3/src/layout/stack.js +245 -0
  183. data/lib/stripchart/public/components/d3/src/layout/tree.js +251 -0
  184. data/lib/stripchart/public/components/d3/src/layout/treemap.js +230 -0
  185. data/lib/stripchart/public/components/d3/src/math/adder.js +34 -0
  186. data/lib/stripchart/public/components/d3/src/math/index.js +2 -0
  187. data/lib/stripchart/public/components/d3/src/math/number.js +3 -0
  188. data/lib/stripchart/public/components/d3/src/math/random.js +28 -0
  189. data/lib/stripchart/public/components/d3/src/math/transform.js +64 -0
  190. data/lib/stripchart/public/components/d3/src/math/trigonometry.js +29 -0
  191. data/lib/stripchart/public/components/d3/src/scale/bilinear.js +7 -0
  192. data/lib/stripchart/public/components/d3/src/scale/category.js +57 -0
  193. data/lib/stripchart/public/components/d3/src/scale/identity.js +33 -0
  194. data/lib/stripchart/public/components/d3/src/scale/index.js +11 -0
  195. data/lib/stripchart/public/components/d3/src/scale/linear.js +128 -0
  196. data/lib/stripchart/public/components/d3/src/scale/log.js +92 -0
  197. data/lib/stripchart/public/components/d3/src/scale/nice.js +28 -0
  198. data/lib/stripchart/public/components/d3/src/scale/ordinal.js +92 -0
  199. data/lib/stripchart/public/components/d3/src/scale/polylinear.js +24 -0
  200. data/lib/stripchart/public/components/d3/src/scale/pow.js +58 -0
  201. data/lib/stripchart/public/components/d3/src/scale/quantile.js +54 -0
  202. data/lib/stripchart/public/components/d3/src/scale/quantize.js +44 -0
  203. data/lib/stripchart/public/components/d3/src/scale/scale.js +10 -0
  204. data/lib/stripchart/public/components/d3/src/scale/sqrt.js +6 -0
  205. data/lib/stripchart/public/components/d3/src/scale/threshold.js +36 -0
  206. data/lib/stripchart/public/components/d3/src/selection/append.js +16 -0
  207. data/lib/stripchart/public/components/d3/src/selection/attr.js +62 -0
  208. data/lib/stripchart/public/components/d3/src/selection/call.js +8 -0
  209. data/lib/stripchart/public/components/d3/src/selection/classed.js +72 -0
  210. data/lib/stripchart/public/components/d3/src/selection/data.js +118 -0
  211. data/lib/stripchart/public/components/d3/src/selection/datum.js +7 -0
  212. data/lib/stripchart/public/components/d3/src/selection/each.js +16 -0
  213. data/lib/stripchart/public/components/d3/src/selection/empty.js +5 -0
  214. data/lib/stripchart/public/components/d3/src/selection/enter-insert.js +20 -0
  215. data/lib/stripchart/public/components/d3/src/selection/enter-select.js +27 -0
  216. data/lib/stripchart/public/components/d3/src/selection/enter.js +21 -0
  217. data/lib/stripchart/public/components/d3/src/selection/filter.js +28 -0
  218. data/lib/stripchart/public/components/d3/src/selection/html.js +10 -0
  219. data/lib/stripchart/public/components/d3/src/selection/index.js +1 -0
  220. data/lib/stripchart/public/components/d3/src/selection/insert.js +9 -0
  221. data/lib/stripchart/public/components/d3/src/selection/node.js +11 -0
  222. data/lib/stripchart/public/components/d3/src/selection/on.js +104 -0
  223. data/lib/stripchart/public/components/d3/src/selection/order.js +13 -0
  224. data/lib/stripchart/public/components/d3/src/selection/property.js +43 -0
  225. data/lib/stripchart/public/components/d3/src/selection/remove.js +11 -0
  226. data/lib/stripchart/public/components/d3/src/selection/select.js +32 -0
  227. data/lib/stripchart/public/components/d3/src/selection/selectAll.js +27 -0
  228. data/lib/stripchart/public/components/d3/src/selection/selection.js +67 -0
  229. data/lib/stripchart/public/components/d3/src/selection/size.js +7 -0
  230. data/lib/stripchart/public/components/d3/src/selection/sort.js +15 -0
  231. data/lib/stripchart/public/components/d3/src/selection/style.js +56 -0
  232. data/lib/stripchart/public/components/d3/src/selection/text.js +10 -0
  233. data/lib/stripchart/public/components/d3/src/selection/transition.js +19 -0
  234. data/lib/stripchart/public/components/d3/src/start.js +2 -0
  235. data/lib/stripchart/public/components/d3/src/svg/arc.js +99 -0
  236. data/lib/stripchart/public/components/d3/src/svg/area-radial.js +14 -0
  237. data/lib/stripchart/public/components/d3/src/svg/area.js +120 -0
  238. data/lib/stripchart/public/components/d3/src/svg/axis.js +222 -0
  239. data/lib/stripchart/public/components/d3/src/svg/brush.js +379 -0
  240. data/lib/stripchart/public/components/d3/src/svg/chord.js +90 -0
  241. data/lib/stripchart/public/components/d3/src/svg/diagonal-radial.js +26 -0
  242. data/lib/stripchart/public/components/d3/src/svg/diagonal.js +43 -0
  243. data/lib/stripchart/public/components/d3/src/svg/index.js +12 -0
  244. data/lib/stripchart/public/components/d3/src/svg/line-radial.js +26 -0
  245. data/lib/stripchart/public/components/d3/src/svg/line.js +439 -0
  246. data/lib/stripchart/public/components/d3/src/svg/svg.js +1 -0
  247. data/lib/stripchart/public/components/d3/src/svg/symbol.js +105 -0
  248. data/lib/stripchart/public/components/d3/src/time/day.js +21 -0
  249. data/lib/stripchart/public/components/d3/src/time/format-iso.js +20 -0
  250. data/lib/stripchart/public/components/d3/src/time/format-locale.js +10 -0
  251. data/lib/stripchart/public/components/d3/src/time/format-localized.js +11 -0
  252. data/lib/stripchart/public/components/d3/src/time/format-utc.js +31 -0
  253. data/lib/stripchart/public/components/d3/src/time/format.js +303 -0
  254. data/lib/stripchart/public/components/d3/src/time/hour.js +14 -0
  255. data/lib/stripchart/public/components/d3/src/time/index.js +15 -0
  256. data/lib/stripchart/public/components/d3/src/time/interval.js +71 -0
  257. data/lib/stripchart/public/components/d3/src/time/minute.js +13 -0
  258. data/lib/stripchart/public/components/d3/src/time/month.js +16 -0
  259. data/lib/stripchart/public/components/d3/src/time/scale-utc.js +43 -0
  260. data/lib/stripchart/public/components/d3/src/time/scale.js +150 -0
  261. data/lib/stripchart/public/components/d3/src/time/second.js +13 -0
  262. data/lib/stripchart/public/components/d3/src/time/time.js +35 -0
  263. data/lib/stripchart/public/components/d3/src/time/week.js +32 -0
  264. data/lib/stripchart/public/components/d3/src/time/year.js +16 -0
  265. data/lib/stripchart/public/components/d3/src/transition/attr.js +58 -0
  266. data/lib/stripchart/public/components/d3/src/transition/delay.js +9 -0
  267. data/lib/stripchart/public/components/d3/src/transition/duration.js +9 -0
  268. data/lib/stripchart/public/components/d3/src/transition/each.js +23 -0
  269. data/lib/stripchart/public/components/d3/src/transition/ease.js +10 -0
  270. data/lib/stripchart/public/components/d3/src/transition/filter.js +22 -0
  271. data/lib/stripchart/public/components/d3/src/transition/index.js +1 -0
  272. data/lib/stripchart/public/components/d3/src/transition/remove.js +8 -0
  273. data/lib/stripchart/public/components/d3/src/transition/select.js +27 -0
  274. data/lib/stripchart/public/components/d3/src/transition/selectAll.js +30 -0
  275. data/lib/stripchart/public/components/d3/src/transition/style.js +53 -0
  276. data/lib/stripchart/public/components/d3/src/transition/subtransition.js +25 -0
  277. data/lib/stripchart/public/components/d3/src/transition/text.js +11 -0
  278. data/lib/stripchart/public/components/d3/src/transition/transition.js +118 -0
  279. data/lib/stripchart/public/components/d3/src/transition/tween.js +17 -0
  280. data/lib/stripchart/public/components/d3/src/xhr/html.js +12 -0
  281. data/lib/stripchart/public/components/d3/src/xhr/index.js +5 -0
  282. data/lib/stripchart/public/components/d3/src/xhr/json.js +9 -0
  283. data/lib/stripchart/public/components/d3/src/xhr/text.js +5 -0
  284. data/lib/stripchart/public/components/d3/src/xhr/xhr.js +118 -0
  285. data/lib/stripchart/public/components/d3/src/xhr/xml.js +5 -0
  286. data/lib/stripchart/public/components/d3/test/XMLHttpRequest.js +55 -0
  287. data/lib/stripchart/public/components/d3/test/arrays/ascending-test.js +45 -0
  288. data/lib/stripchart/public/components/d3/test/arrays/bisect-test.js +257 -0
  289. data/lib/stripchart/public/components/d3/test/arrays/descending-test.js +45 -0
  290. data/lib/stripchart/public/components/d3/test/arrays/entries-test.js +39 -0
  291. data/lib/stripchart/public/components/d3/test/arrays/extent-test.js +50 -0
  292. data/lib/stripchart/public/components/d3/test/arrays/keys-test.js +27 -0
  293. data/lib/stripchart/public/components/d3/test/arrays/map-test.js +214 -0
  294. data/lib/stripchart/public/components/d3/test/arrays/max-test.js +51 -0
  295. data/lib/stripchart/public/components/d3/test/arrays/mean-test.js +39 -0
  296. data/lib/stripchart/public/components/d3/test/arrays/median-test.js +39 -0
  297. data/lib/stripchart/public/components/d3/test/arrays/merge-test.js +26 -0
  298. data/lib/stripchart/public/components/d3/test/arrays/min-test.js +51 -0
  299. data/lib/stripchart/public/components/d3/test/arrays/nest-test.js +256 -0
  300. data/lib/stripchart/public/components/d3/test/arrays/permute-test.js +49 -0
  301. data/lib/stripchart/public/components/d3/test/arrays/quantile-test.js +56 -0
  302. data/lib/stripchart/public/components/d3/test/arrays/range-test.js +99 -0
  303. data/lib/stripchart/public/components/d3/test/arrays/set-test.js +136 -0
  304. data/lib/stripchart/public/components/d3/test/arrays/sum-test.js +43 -0
  305. data/lib/stripchart/public/components/d3/test/arrays/transpose-test.js +28 -0
  306. data/lib/stripchart/public/components/d3/test/arrays/values-test.js +31 -0
  307. data/lib/stripchart/public/components/d3/test/arrays/zip-test.js +28 -0
  308. data/lib/stripchart/public/components/d3/test/assert.js +153 -0
  309. data/lib/stripchart/public/components/d3/test/color/hcl-test.js +112 -0
  310. data/lib/stripchart/public/components/d3/test/color/hsl-test.js +118 -0
  311. data/lib/stripchart/public/components/d3/test/color/lab-test.js +81 -0
  312. data/lib/stripchart/public/components/d3/test/color/rgb-test.js +91 -0
  313. data/lib/stripchart/public/components/d3/test/core/functor-test.js +25 -0
  314. data/lib/stripchart/public/components/d3/test/core/ns-test.js +57 -0
  315. data/lib/stripchart/public/components/d3/test/core/rebind-test.js +50 -0
  316. data/lib/stripchart/public/components/d3/test/core/version-test.js +16 -0
  317. data/lib/stripchart/public/components/d3/test/data/sample.csv +2 -0
  318. data/lib/stripchart/public/components/d3/test/data/sample.html +5 -0
  319. data/lib/stripchart/public/components/d3/test/data/sample.json +1 -0
  320. data/lib/stripchart/public/components/d3/test/data/sample.tsv +2 -0
  321. data/lib/stripchart/public/components/d3/test/data/sample.txt +1 -0
  322. data/lib/stripchart/public/components/d3/test/data/sample.xml +4 -0
  323. data/lib/stripchart/public/components/d3/test/data/us-counties.json +3186 -0
  324. data/lib/stripchart/public/components/d3/test/dsv/csv-test.js +214 -0
  325. data/lib/stripchart/public/components/d3/test/dsv/tsv-test.js +185 -0
  326. data/lib/stripchart/public/components/d3/test/event/dispatch-test.js +163 -0
  327. data/lib/stripchart/public/components/d3/test/event/mouse-test.html +22 -0
  328. data/lib/stripchart/public/components/d3/test/event/timer-test.js +74 -0
  329. data/lib/stripchart/public/components/d3/test/event/touch-test.html +61 -0
  330. data/lib/stripchart/public/components/d3/test/format/format-test.js +360 -0
  331. data/lib/stripchart/public/components/d3/test/format/formatPrefix-test.js +112 -0
  332. data/lib/stripchart/public/components/d3/test/format/requote-test.js +49 -0
  333. data/lib/stripchart/public/components/d3/test/format/round-test.js +66 -0
  334. data/lib/stripchart/public/components/d3/test/geo/albers-test.js +40 -0
  335. data/lib/stripchart/public/components/d3/test/geo/albers-usa-test.js +33 -0
  336. data/lib/stripchart/public/components/d3/test/geo/area-benchmark.js +12 -0
  337. data/lib/stripchart/public/components/d3/test/geo/area-test.js +185 -0
  338. data/lib/stripchart/public/components/d3/test/geo/azimuthal-equal-area-test.js +38 -0
  339. data/lib/stripchart/public/components/d3/test/geo/azimuthal-equidistant-test.js +38 -0
  340. data/lib/stripchart/public/components/d3/test/geo/benchmark.js +83 -0
  341. data/lib/stripchart/public/components/d3/test/geo/bounds-test.js +265 -0
  342. data/lib/stripchart/public/components/d3/test/geo/centroid-test.js +187 -0
  343. data/lib/stripchart/public/components/d3/test/geo/circle-test.js +29 -0
  344. data/lib/stripchart/public/components/d3/test/geo/clip-view-mock.js +4 -0
  345. data/lib/stripchart/public/components/d3/test/geo/clip-view-test.js +44 -0
  346. data/lib/stripchart/public/components/d3/test/geo/conic-conformal-test.js +38 -0
  347. data/lib/stripchart/public/components/d3/test/geo/conic-equal-area-test.js +38 -0
  348. data/lib/stripchart/public/components/d3/test/geo/conic-equidistant-test.js +38 -0
  349. data/lib/stripchart/public/components/d3/test/geo/distance-test.js +20 -0
  350. data/lib/stripchart/public/components/d3/test/geo/equirectangular-test.js +38 -0
  351. data/lib/stripchart/public/components/d3/test/geo/gnomonic-test.js +38 -0
  352. data/lib/stripchart/public/components/d3/test/geo/graticule-test.js +188 -0
  353. data/lib/stripchart/public/components/d3/test/geo/greatArc-test.js +35 -0
  354. data/lib/stripchart/public/components/d3/test/geo/interpolate-test.js +22 -0
  355. data/lib/stripchart/public/components/d3/test/geo/length-test.js +44 -0
  356. data/lib/stripchart/public/components/d3/test/geo/mercator-test.js +38 -0
  357. data/lib/stripchart/public/components/d3/test/geo/orthographic-test.js +34 -0
  358. data/lib/stripchart/public/components/d3/test/geo/path-test.js +832 -0
  359. data/lib/stripchart/public/components/d3/test/geo/point-in-polygon-mock.js +18 -0
  360. data/lib/stripchart/public/components/d3/test/geo/point-in-polygon-test.js +352 -0
  361. data/lib/stripchart/public/components/d3/test/geo/projection-test-suite.js +29 -0
  362. data/lib/stripchart/public/components/d3/test/geo/projection-test.js +58 -0
  363. data/lib/stripchart/public/components/d3/test/geo/rotation-test.js +35 -0
  364. data/lib/stripchart/public/components/d3/test/geo/stereographic-test.js +38 -0
  365. data/lib/stripchart/public/components/d3/test/geo/stream-test.js +264 -0
  366. data/lib/stripchart/public/components/d3/test/geo/transverse-mercator-test.js +38 -0
  367. data/lib/stripchart/public/components/d3/test/geom/hull-test.js +68 -0
  368. data/lib/stripchart/public/components/d3/test/geom/polygon-test.js +151 -0
  369. data/lib/stripchart/public/components/d3/test/geom/quadtree-test.js +139 -0
  370. data/lib/stripchart/public/components/d3/test/geom/voronoi-test.js +230 -0
  371. data/lib/stripchart/public/components/d3/test/interpolate/array-test.js +24 -0
  372. data/lib/stripchart/public/components/d3/test/interpolate/ease-test.js +135 -0
  373. data/lib/stripchart/public/components/d3/test/interpolate/hcl-test.js +47 -0
  374. data/lib/stripchart/public/components/d3/test/interpolate/hsl-test.js +47 -0
  375. data/lib/stripchart/public/components/d3/test/interpolate/interpolate-test.js +157 -0
  376. data/lib/stripchart/public/components/d3/test/interpolate/interpolate-transform-test.js +60 -0
  377. data/lib/stripchart/public/components/d3/test/interpolate/lab-test.js +31 -0
  378. data/lib/stripchart/public/components/d3/test/interpolate/number-test.js +20 -0
  379. data/lib/stripchart/public/components/d3/test/interpolate/object-test.js +35 -0
  380. data/lib/stripchart/public/components/d3/test/interpolate/rgb-test.js +31 -0
  381. data/lib/stripchart/public/components/d3/test/interpolate/round-test.js +17 -0
  382. data/lib/stripchart/public/components/d3/test/interpolate/string-test.js +43 -0
  383. data/lib/stripchart/public/components/d3/test/layout/chord-test.js +164 -0
  384. data/lib/stripchart/public/components/d3/test/layout/cluster-test.js +44 -0
  385. data/lib/stripchart/public/components/d3/test/layout/force-test.js +119 -0
  386. data/lib/stripchart/public/components/d3/test/layout/hierarchy-test.js +27 -0
  387. data/lib/stripchart/public/components/d3/test/layout/histogram-test.js +92 -0
  388. data/lib/stripchart/public/components/d3/test/layout/pack-test.js +129 -0
  389. data/lib/stripchart/public/components/d3/test/layout/partition-test.js +45 -0
  390. data/lib/stripchart/public/components/d3/test/layout/pie-test.js +24 -0
  391. data/lib/stripchart/public/components/d3/test/layout/tree-test.js +45 -0
  392. data/lib/stripchart/public/components/d3/test/layout/treemap-test.js +181 -0
  393. data/lib/stripchart/public/components/d3/test/load.js +58 -0
  394. data/lib/stripchart/public/components/d3/test/math/random-test.js +37 -0
  395. data/lib/stripchart/public/components/d3/test/math/transform-null-matrix-test.html +26 -0
  396. data/lib/stripchart/public/components/d3/test/math/transform-null-test.html +26 -0
  397. data/lib/stripchart/public/components/d3/test/math/transform-rotate-origin-test.html +40 -0
  398. data/lib/stripchart/public/components/d3/test/math/transform-rotate-test.html +138 -0
  399. data/lib/stripchart/public/components/d3/test/math/transform-test.html +88 -0
  400. data/lib/stripchart/public/components/d3/test/scale/category-test.js +76 -0
  401. data/lib/stripchart/public/components/d3/test/scale/identity-test.js +163 -0
  402. data/lib/stripchart/public/components/d3/test/scale/linear-test.js +281 -0
  403. data/lib/stripchart/public/components/d3/test/scale/log-test.js +378 -0
  404. data/lib/stripchart/public/components/d3/test/scale/ordinal-test.js +272 -0
  405. data/lib/stripchart/public/components/d3/test/scale/pow-test.js +267 -0
  406. data/lib/stripchart/public/components/d3/test/scale/quantile-test.js +91 -0
  407. data/lib/stripchart/public/components/d3/test/scale/quantize-test.js +92 -0
  408. data/lib/stripchart/public/components/d3/test/scale/sqrt-test.js +260 -0
  409. data/lib/stripchart/public/components/d3/test/scale/threshold-test.js +64 -0
  410. data/lib/stripchart/public/components/d3/test/selection/append-test.js +153 -0
  411. data/lib/stripchart/public/components/d3/test/selection/attr-test.js +182 -0
  412. data/lib/stripchart/public/components/d3/test/selection/call-test.js +75 -0
  413. data/lib/stripchart/public/components/d3/test/selection/classed-test.js +269 -0
  414. data/lib/stripchart/public/components/d3/test/selection/data-test.js +270 -0
  415. data/lib/stripchart/public/components/d3/test/selection/datum-test.js +87 -0
  416. data/lib/stripchart/public/components/d3/test/selection/each-test.js +87 -0
  417. data/lib/stripchart/public/components/d3/test/selection/empty-test.js +53 -0
  418. data/lib/stripchart/public/components/d3/test/selection/enter-test.js +24 -0
  419. data/lib/stripchart/public/components/d3/test/selection/filter-test.js +72 -0
  420. data/lib/stripchart/public/components/d3/test/selection/html-test.js +134 -0
  421. data/lib/stripchart/public/components/d3/test/selection/insert-test.js +169 -0
  422. data/lib/stripchart/public/components/d3/test/selection/node-test.js +55 -0
  423. data/lib/stripchart/public/components/d3/test/selection/on-test.js +135 -0
  424. data/lib/stripchart/public/components/d3/test/selection/order-test.js +32 -0
  425. data/lib/stripchart/public/components/d3/test/selection/property-test.js +118 -0
  426. data/lib/stripchart/public/components/d3/test/selection/remove-test.js +38 -0
  427. data/lib/stripchart/public/components/d3/test/selection/select-test.js +69 -0
  428. data/lib/stripchart/public/components/d3/test/selection/selectAll-test.js +52 -0
  429. data/lib/stripchart/public/components/d3/test/selection/selection-select-test.js +128 -0
  430. data/lib/stripchart/public/components/d3/test/selection/selection-selectAll-test.js +131 -0
  431. data/lib/stripchart/public/components/d3/test/selection/selection-test.js +38 -0
  432. data/lib/stripchart/public/components/d3/test/selection/size-test.js +52 -0
  433. data/lib/stripchart/public/components/d3/test/selection/sort-test.js +81 -0
  434. data/lib/stripchart/public/components/d3/test/selection/style-test.js +129 -0
  435. data/lib/stripchart/public/components/d3/test/selection/text-test.js +115 -0
  436. data/lib/stripchart/public/components/d3/test/svg/arc-test.js +143 -0
  437. data/lib/stripchart/public/components/d3/test/svg/area-radial-test.js +199 -0
  438. data/lib/stripchart/public/components/d3/test/svg/area-test.js +192 -0
  439. data/lib/stripchart/public/components/d3/test/svg/axis-test.js +447 -0
  440. data/lib/stripchart/public/components/d3/test/svg/brush-test.js +96 -0
  441. data/lib/stripchart/public/components/d3/test/svg/diagonal-test.js +62 -0
  442. data/lib/stripchart/public/components/d3/test/svg/line-radial-test.js +120 -0
  443. data/lib/stripchart/public/components/d3/test/svg/line-test.js +207 -0
  444. data/lib/stripchart/public/components/d3/test/svg/symbol-test.js +91 -0
  445. data/lib/stripchart/public/components/d3/test/time/day-test.js +178 -0
  446. data/lib/stripchart/public/components/d3/test/time/dayOfYear-test.js +16 -0
  447. data/lib/stripchart/public/components/d3/test/time/days-test.js +95 -0
  448. data/lib/stripchart/public/components/d3/test/time/format-iso-test.js +33 -0
  449. data/lib/stripchart/public/components/d3/test/time/format-test.js +300 -0
  450. data/lib/stripchart/public/components/d3/test/time/format-utc-test.js +219 -0
  451. data/lib/stripchart/public/components/d3/test/time/hour-test.js +215 -0
  452. data/lib/stripchart/public/components/d3/test/time/hours-test.js +109 -0
  453. data/lib/stripchart/public/components/d3/test/time/minute-test.js +127 -0
  454. data/lib/stripchart/public/components/d3/test/time/minutes-test.js +91 -0
  455. data/lib/stripchart/public/components/d3/test/time/month-test.js +150 -0
  456. data/lib/stripchart/public/components/d3/test/time/months-test.js +95 -0
  457. data/lib/stripchart/public/components/d3/test/time/scale-test.js +585 -0
  458. data/lib/stripchart/public/components/d3/test/time/second-test.js +143 -0
  459. data/lib/stripchart/public/components/d3/test/time/seconds-test.js +91 -0
  460. data/lib/stripchart/public/components/d3/test/time/time.js +30 -0
  461. data/lib/stripchart/public/components/d3/test/time/week-test.js +162 -0
  462. data/lib/stripchart/public/components/d3/test/time/weeks-test.js +95 -0
  463. data/lib/stripchart/public/components/d3/test/time/year-test.js +137 -0
  464. data/lib/stripchart/public/components/d3/test/time/years-test.js +74 -0
  465. data/lib/stripchart/public/components/d3/test/transition/transition-test-attr.js +56 -0
  466. data/lib/stripchart/public/components/d3/test/transition/transition-test-attrTween.js +59 -0
  467. data/lib/stripchart/public/components/d3/test/transition/transition-test-call.js +30 -0
  468. data/lib/stripchart/public/components/d3/test/transition/transition-test-delay.js +56 -0
  469. data/lib/stripchart/public/components/d3/test/transition/transition-test-duration.js +60 -0
  470. data/lib/stripchart/public/components/d3/test/transition/transition-test-each.js +163 -0
  471. data/lib/stripchart/public/components/d3/test/transition/transition-test-filter.js +84 -0
  472. data/lib/stripchart/public/components/d3/test/transition/transition-test-id.js +22 -0
  473. data/lib/stripchart/public/components/d3/test/transition/transition-test-node.js +39 -0
  474. data/lib/stripchart/public/components/d3/test/transition/transition-test-remove.js +45 -0
  475. data/lib/stripchart/public/components/d3/test/transition/transition-test-select.js +60 -0
  476. data/lib/stripchart/public/components/d3/test/transition/transition-test-selectAll.js +56 -0
  477. data/lib/stripchart/public/components/d3/test/transition/transition-test-size.js +37 -0
  478. data/lib/stripchart/public/components/d3/test/transition/transition-test-style.js +51 -0
  479. data/lib/stripchart/public/components/d3/test/transition/transition-test-styleTween.js +65 -0
  480. data/lib/stripchart/public/components/d3/test/transition/transition-test-text.js +27 -0
  481. data/lib/stripchart/public/components/d3/test/transition/transition-test-time.js +38 -0
  482. data/lib/stripchart/public/components/d3/test/transition/transition-test-transition.js +59 -0
  483. data/lib/stripchart/public/components/d3/test/transition/transition-test-tween.js +62 -0
  484. data/lib/stripchart/public/components/d3/test/transition/transition-test.js +89 -0
  485. data/lib/stripchart/public/components/d3/test/xhr/html-test.js +37 -0
  486. data/lib/stripchart/public/components/d3/test/xhr/json-test.js +51 -0
  487. data/lib/stripchart/public/components/d3/test/xhr/text-test.js +46 -0
  488. data/lib/stripchart/public/components/d3/test/xhr/xhr-test.js +54 -0
  489. data/lib/stripchart/public/components/d3/test/xhr/xml-test.js +46 -0
  490. data/lib/stripchart/public/components/jquery/README.md +11 -0
  491. data/lib/stripchart/public/components/jquery/bower.json +11 -0
  492. data/lib/stripchart/public/components/jquery/component.json +23 -0
  493. data/lib/stripchart/public/components/jquery/composer.json +35 -0
  494. data/lib/stripchart/public/components/jquery/jquery-migrate.js +511 -0
  495. data/lib/stripchart/public/components/jquery/jquery-migrate.min.js +3 -0
  496. data/lib/stripchart/public/components/jquery/jquery.js +8829 -0
  497. data/lib/stripchart/public/components/jquery/jquery.min.js +6 -0
  498. data/lib/stripchart/public/components/jquery/jquery.min.map +1 -0
  499. data/lib/stripchart/public/components/jquery/package.json +7 -0
  500. data/lib/stripchart/version.rb +3 -0
  501. data/lib/stripchart/views/index.erb +19 -0
  502. data/lib/stripchart/web.rb +22 -0
  503. data/lib/stripchart/websocket.rb +22 -0
  504. metadata +602 -0
@@ -0,0 +1,225 @@
1
+ (function(){science = {version: "1.7.0"}; // semver
2
+ science.ascending = function(a, b) {
3
+ return a - b;
4
+ };
5
+ // Euler's constant.
6
+ science.EULER = .5772156649015329;
7
+ // Compute exp(x) - 1 accurately for small x.
8
+ science.expm1 = function(x) {
9
+ return (x < 1e-5 && x > -1e-5) ? x + .5 * x * x : Math.exp(x) - 1;
10
+ };
11
+ science.functor = function(v) {
12
+ return typeof v === "function" ? v : function() { return v; };
13
+ };
14
+ // Based on:
15
+ // http://www.johndcook.com/blog/2010/06/02/whats-so-hard-about-finding-a-hypotenuse/
16
+ science.hypot = function(x, y) {
17
+ x = Math.abs(x);
18
+ y = Math.abs(y);
19
+ var max,
20
+ min;
21
+ if (x > y) { max = x; min = y; }
22
+ else { max = y; min = x; }
23
+ var r = min / max;
24
+ return max * Math.sqrt(1 + r * r);
25
+ };
26
+ science.quadratic = function() {
27
+ var complex = false;
28
+
29
+ function quadratic(a, b, c) {
30
+ var d = b * b - 4 * a * c;
31
+ if (d > 0) {
32
+ d = Math.sqrt(d) / (2 * a);
33
+ return complex
34
+ ? [{r: -b - d, i: 0}, {r: -b + d, i: 0}]
35
+ : [-b - d, -b + d];
36
+ } else if (d === 0) {
37
+ d = -b / (2 * a);
38
+ return complex ? [{r: d, i: 0}] : [d];
39
+ } else {
40
+ if (complex) {
41
+ d = Math.sqrt(-d) / (2 * a);
42
+ return [
43
+ {r: -b, i: -d},
44
+ {r: -b, i: d}
45
+ ];
46
+ }
47
+ return [];
48
+ }
49
+ }
50
+
51
+ quadratic.complex = function(x) {
52
+ if (!arguments.length) return complex;
53
+ complex = x;
54
+ return quadratic;
55
+ };
56
+
57
+ return quadratic;
58
+ };
59
+ // Constructs a multi-dimensional array filled with zeroes.
60
+ science.zeroes = function(n) {
61
+ var i = -1,
62
+ a = [];
63
+ if (arguments.length === 1)
64
+ while (++i < n)
65
+ a[i] = 0;
66
+ else
67
+ while (++i < n)
68
+ a[i] = science.zeroes.apply(
69
+ this, Array.prototype.slice.call(arguments, 1));
70
+ return a;
71
+ };
72
+ science.vector = {};
73
+ science.vector.cross = function(a, b) {
74
+ // TODO how to handle non-3D vectors?
75
+ // TODO handle 7D vectors?
76
+ return [
77
+ a[1] * b[2] - a[2] * b[1],
78
+ a[2] * b[0] - a[0] * b[2],
79
+ a[0] * b[1] - a[1] * b[0]
80
+ ];
81
+ };
82
+ science.vector.dot = function(a, b) {
83
+ var s = 0,
84
+ i = -1,
85
+ n = Math.min(a.length, b.length);
86
+ while (++i < n) s += a[i] * b[i];
87
+ return s;
88
+ };
89
+ science.vector.length = function(p) {
90
+ return Math.sqrt(science.vector.dot(p, p));
91
+ };
92
+ science.vector.normalize = function(p) {
93
+ var length = science.vector.length(p);
94
+ return p.map(function(d) { return d / length; });
95
+ };
96
+ // 4x4 matrix determinant.
97
+ science.vector.determinant = function(matrix) {
98
+ var m = matrix[0].concat(matrix[1]).concat(matrix[2]).concat(matrix[3]);
99
+ return (
100
+ m[12] * m[9] * m[6] * m[3] - m[8] * m[13] * m[6] * m[3] -
101
+ m[12] * m[5] * m[10] * m[3] + m[4] * m[13] * m[10] * m[3] +
102
+ m[8] * m[5] * m[14] * m[3] - m[4] * m[9] * m[14] * m[3] -
103
+ m[12] * m[9] * m[2] * m[7] + m[8] * m[13] * m[2] * m[7] +
104
+ m[12] * m[1] * m[10] * m[7] - m[0] * m[13] * m[10] * m[7] -
105
+ m[8] * m[1] * m[14] * m[7] + m[0] * m[9] * m[14] * m[7] +
106
+ m[12] * m[5] * m[2] * m[11] - m[4] * m[13] * m[2] * m[11] -
107
+ m[12] * m[1] * m[6] * m[11] + m[0] * m[13] * m[6] * m[11] +
108
+ m[4] * m[1] * m[14] * m[11] - m[0] * m[5] * m[14] * m[11] -
109
+ m[8] * m[5] * m[2] * m[15] + m[4] * m[9] * m[2] * m[15] +
110
+ m[8] * m[1] * m[6] * m[15] - m[0] * m[9] * m[6] * m[15] -
111
+ m[4] * m[1] * m[10] * m[15] + m[0] * m[5] * m[10] * m[15]);
112
+ };
113
+ // Performs in-place Gauss-Jordan elimination.
114
+ //
115
+ // Based on Jarno Elonen's Python version (public domain):
116
+ // http://elonen.iki.fi/code/misc-notes/python-gaussj/index.html
117
+ science.vector.gaussjordan = function(m, eps) {
118
+ if (!eps) eps = 1e-10;
119
+
120
+ var h = m.length,
121
+ w = m[0].length,
122
+ y = -1,
123
+ y2,
124
+ x;
125
+
126
+ while (++y < h) {
127
+ var maxrow = y;
128
+
129
+ // Find max pivot.
130
+ y2 = y; while (++y2 < h) {
131
+ if (Math.abs(m[y2][y]) > Math.abs(m[maxrow][y]))
132
+ maxrow = y2;
133
+ }
134
+
135
+ // Swap.
136
+ var tmp = m[y];
137
+ m[y] = m[maxrow];
138
+ m[maxrow] = tmp;
139
+
140
+ // Singular?
141
+ if (Math.abs(m[y][y]) <= eps) return false;
142
+
143
+ // Eliminate column y.
144
+ y2 = y; while (++y2 < h) {
145
+ var c = m[y2][y] / m[y][y];
146
+ x = y - 1; while (++x < w) {
147
+ m[y2][x] -= m[y][x] * c;
148
+ }
149
+ }
150
+ }
151
+
152
+ // Backsubstitute.
153
+ y = h; while (--y >= 0) {
154
+ var c = m[y][y];
155
+ y2 = -1; while (++y2 < y) {
156
+ x = w; while (--x >= y) {
157
+ m[y2][x] -= m[y][x] * m[y2][y] / c;
158
+ }
159
+ }
160
+ m[y][y] /= c;
161
+ // Normalize row y.
162
+ x = h - 1; while (++x < w) {
163
+ m[y][x] /= c;
164
+ }
165
+ }
166
+ return true;
167
+ };
168
+ // Find matrix inverse using Gauss-Jordan.
169
+ science.vector.inverse = function(m) {
170
+ var n = m.length
171
+ i = -1;
172
+
173
+ // Check if the matrix is square.
174
+ if (n !== m[0].length) return;
175
+
176
+ // Augment with identity matrix I to get AI.
177
+ m = m.map(function(row, i) {
178
+ var identity = new Array(n),
179
+ j = -1;
180
+ while (++j < n) identity[j] = i === j ? 1 : 0;
181
+ return row.concat(identity);
182
+ });
183
+
184
+ // Compute IA^-1.
185
+ science.vector.gaussjordan(m);
186
+
187
+ // Remove identity matrix I to get A^-1.
188
+ while (++i < n) {
189
+ m[i] = m[i].slice(n);
190
+ }
191
+
192
+ return m;
193
+ };
194
+ science.vector.multiply = function(a, b) {
195
+ var m = a.length,
196
+ n = b[0].length,
197
+ p = b.length,
198
+ i = -1,
199
+ j,
200
+ k;
201
+ if (p !== a[0].length) throw {"error": "columns(a) != rows(b); " + a[0].length + " != " + p};
202
+ var ab = new Array(m);
203
+ while (++i < m) {
204
+ ab[i] = new Array(n);
205
+ j = -1; while(++j < n) {
206
+ var s = 0;
207
+ k = -1; while (++k < p) s += a[i][k] * b[k][j];
208
+ ab[i][j] = s;
209
+ }
210
+ }
211
+ return ab;
212
+ };
213
+ science.vector.transpose = function(a) {
214
+ var m = a.length,
215
+ n = a[0].length,
216
+ i = -1,
217
+ j,
218
+ b = new Array(n);
219
+ while (++i < n) {
220
+ b[i] = new Array(m);
221
+ j = -1; while (++j < m) b[i][j] = a[j][i];
222
+ }
223
+ return b;
224
+ };
225
+ })()
@@ -0,0 +1,27 @@
1
+ (function(){science.lin = {};
2
+ /**
3
+ * Solves tridiagonal systems of linear equations.
4
+ *
5
+ * Source: http://en.wikipedia.org/wiki/Tridiagonal_matrix_algorithm
6
+ *
7
+ * @param {number[]} a
8
+ * @param {number[]} b
9
+ * @param {number[]} c
10
+ * @param {number[]} d
11
+ * @param {number[]} x
12
+ * @param {number} n
13
+ */
14
+ science.lin.tridag = function(a, b, c, d, x, n) {
15
+ var i,
16
+ m;
17
+ for (i = 1; i < n; i++) {
18
+ m = a[i] / b[i - 1];
19
+ b[i] -= m * c[i - 1];
20
+ d[i] -= m * d[i - 1];
21
+ }
22
+ x[n - 1] = d[n - 1] / b[n - 1];
23
+ for (i = n - 2; i >= 0; i--) {
24
+ x[i] = (d[i] - c[i] * x[i + 1]) / b[i];
25
+ }
26
+ };
27
+ })()
@@ -0,0 +1 @@
1
+ (function(){science.lin={},science.lin.tridag=function(a,b,c,d,e,f){var g,h;for(g=1;g<f;g++)h=a[g]/b[g-1],b[g]-=h*c[g-1],d[g]-=h*d[g-1];e[f-1]=d[f-1]/b[f-1];for(g=f-2;g>=0;g--)e[g]=(d[g]-c[g]*e[g+1])/b[g]}})()
@@ -0,0 +1 @@
1
+ (function(){science={version:"1.7.0"},science.ascending=function(a,b){return a-b},science.EULER=.5772156649015329,science.expm1=function(a){return a<1e-5&&a>-0.00001?a+.5*a*a:Math.exp(a)-1},science.functor=function(a){return typeof a=="function"?a:function(){return a}},science.hypot=function(a,b){a=Math.abs(a),b=Math.abs(b);var c,d;a>b?(c=a,d=b):(c=b,d=a);var e=d/c;return c*Math.sqrt(1+e*e)},science.quadratic=function(){function b(b,c,d){var e=c*c-4*b*d;if(e>0){e=Math.sqrt(e)/(2*b);return a?[{r:-c-e,i:0},{r:-c+e,i:0}]:[-c-e,-c+e]}if(e===0){e=-c/(2*b);return a?[{r:e,i:0}]:[e]}if(a){e=Math.sqrt(-e)/(2*b);return[{r:-c,i:-e},{r:-c,i:e}]}return[]}var a=!1;b.complex=function(c){if(!arguments.length)return a;a=c;return b};return b},science.zeroes=function(a){var b=-1,c=[];if(arguments.length===1)while(++b<a)c[b]=0;else while(++b<a)c[b]=science.zeroes.apply(this,Array.prototype.slice.call(arguments,1));return c},science.vector={},science.vector.cross=function(a,b){return[a[1]*b[2]-a[2]*b[1],a[2]*b[0]-a[0]*b[2],a[0]*b[1]-a[1]*b[0]]},science.vector.dot=function(a,b){var c=0,d=-1,e=Math.min(a.length,b.length);while(++d<e)c+=a[d]*b[d];return c},science.vector.length=function(a){return Math.sqrt(science.vector.dot(a,a))},science.vector.normalize=function(a){var b=science.vector.length(a);return a.map(function(a){return a/b})},science.vector.determinant=function(a){var b=a[0].concat(a[1]).concat(a[2]).concat(a[3]);return b[12]*b[9]*b[6]*b[3]-b[8]*b[13]*b[6]*b[3]-b[12]*b[5]*b[10]*b[3]+b[4]*b[13]*b[10]*b[3]+b[8]*b[5]*b[14]*b[3]-b[4]*b[9]*b[14]*b[3]-b[12]*b[9]*b[2]*b[7]+b[8]*b[13]*b[2]*b[7]+b[12]*b[1]*b[10]*b[7]-b[0]*b[13]*b[10]*b[7]-b[8]*b[1]*b[14]*b[7]+b[0]*b[9]*b[14]*b[7]+b[12]*b[5]*b[2]*b[11]-b[4]*b[13]*b[2]*b[11]-b[12]*b[1]*b[6]*b[11]+b[0]*b[13]*b[6]*b[11]+b[4]*b[1]*b[14]*b[11]-b[0]*b[5]*b[14]*b[11]-b[8]*b[5]*b[2]*b[15]+b[4]*b[9]*b[2]*b[15]+b[8]*b[1]*b[6]*b[15]-b[0]*b[9]*b[6]*b[15]-b[4]*b[1]*b[10]*b[15]+b[0]*b[5]*b[10]*b[15]},science.vector.gaussjordan=function(a,b){b||(b=1e-10);var c=a.length,d=a[0].length,e=-1,f,g;while(++e<c){var h=e;f=e;while(++f<c)Math.abs(a[f][e])>Math.abs(a[h][e])&&(h=f);var i=a[e];a[e]=a[h],a[h]=i;if(Math.abs(a[e][e])<=b)return!1;f=e;while(++f<c){var j=a[f][e]/a[e][e];g=e-1;while(++g<d)a[f][g]-=a[e][g]*j}}e=c;while(--e>=0){var j=a[e][e];f=-1;while(++f<e){g=d;while(--g>=e)a[f][g]-=a[e][g]*a[f][e]/j}a[e][e]/=j,g=c-1;while(++g<d)a[e][g]/=j}return!0},science.vector.inverse=function(a){var b=a.length;i=-1;if(b===a[0].length){a=a.map(function(a,c){var d=Array(b),e=-1;while(++e<b)d[e]=c===e?1:0;return a.concat(d)}),science.vector.gaussjordan(a);while(++i<b)a[i]=a[i].slice(b);return a}},science.vector.multiply=function(a,b){var c=a.length,d=b[0].length,e=b.length,f=-1,g,h;if(e!==a[0].length)throw{error:"columns(a) != rows(b); "+a[0].length+" != "+e};var i=Array(c);while(++f<c){i[f]=Array(d),g=-1;while(++g<d){var j=0;h=-1;while(++h<e)j+=a[f][h]*b[h][g];i[f][g]=j}}return i},science.vector.transpose=function(a){var b=a.length,c=a[0].length,d=-1,e,f=Array(c);while(++d<c){f[d]=Array(b),e=-1;while(++e<b)f[d][e]=a[e][d]}return f}})()
@@ -0,0 +1,720 @@
1
+ (function(){science.stats = {};
2
+ // Bandwidth selectors for Gaussian kernels.
3
+ // Based on R's implementations in `stats.bw`.
4
+ science.stats.bandwidth = {
5
+
6
+ // Silverman, B. W. (1986) Density Estimation. London: Chapman and Hall.
7
+ nrd0: function(x) {
8
+ var hi = Math.sqrt(science.stats.variance(x));
9
+ if (!(lo = Math.min(hi, science.stats.iqr(x) / 1.34)))
10
+ (lo = hi) || (lo = Math.abs(x[1])) || (lo = 1);
11
+ return .9 * lo * Math.pow(x.length, -.2);
12
+ },
13
+
14
+ // Scott, D. W. (1992) Multivariate Density Estimation: Theory, Practice, and
15
+ // Visualization. Wiley.
16
+ nrd: function(x) {
17
+ var h = science.stats.iqr(x) / 1.34;
18
+ return 1.06 * Math.min(Math.sqrt(science.stats.variance(x)), h)
19
+ * Math.pow(x.length, -1/5);
20
+ }
21
+ };
22
+ science.stats.distance = {
23
+ euclidean: function(a, b) {
24
+ var n = a.length,
25
+ i = -1,
26
+ s = 0,
27
+ x;
28
+ while (++i < n) {
29
+ x = a[i] - b[i];
30
+ s += x * x;
31
+ }
32
+ return Math.sqrt(s);
33
+ },
34
+ manhattan: function(a, b) {
35
+ var n = a.length,
36
+ i = -1,
37
+ s = 0;
38
+ while (++i < n) s += Math.abs(a[i] - b[i]);
39
+ return s;
40
+ },
41
+ minkowski: function(p) {
42
+ return function(a, b) {
43
+ var n = a.length,
44
+ i = -1,
45
+ s = 0;
46
+ while (++i < n) s += Math.pow(Math.abs(a[i] - b[i]), p);
47
+ return Math.pow(s, 1 / p);
48
+ };
49
+ },
50
+ chebyshev: function(a, b) {
51
+ var n = a.length,
52
+ i = -1,
53
+ max = 0,
54
+ x;
55
+ while (++i < n) {
56
+ x = Math.abs(a[i] - b[i]);
57
+ if (x > max) max = x;
58
+ }
59
+ return max;
60
+ },
61
+ hamming: function(a, b) {
62
+ var n = a.length,
63
+ i = -1,
64
+ d = 0;
65
+ while (++i < n) if (a[i] !== b[i]) d++;
66
+ return d;
67
+ },
68
+ jaccard: function(a, b) {
69
+ var n = a.length,
70
+ i = -1,
71
+ s = 0;
72
+ while (++i < n) if (a[i] === b[i]) s++;
73
+ return s / n;
74
+ },
75
+ braycurtis: function(a, b) {
76
+ var n = a.length,
77
+ i = -1,
78
+ s0 = 0,
79
+ s1 = 0,
80
+ ai,
81
+ bi;
82
+ while (++i < n) {
83
+ ai = a[i];
84
+ bi = b[i];
85
+ s0 += Math.abs(ai - bi);
86
+ s1 += Math.abs(ai + bi);
87
+ }
88
+ return s0 / s1;
89
+ }
90
+ };
91
+ // Based on implementation in http://picomath.org/.
92
+ science.stats.erf = function(x) {
93
+ var a1 = 0.254829592,
94
+ a2 = -0.284496736,
95
+ a3 = 1.421413741,
96
+ a4 = -1.453152027,
97
+ a5 = 1.061405429,
98
+ p = 0.3275911;
99
+
100
+ // Save the sign of x
101
+ var sign = x < 0 ? -1 : 1;
102
+ if (x < 0) {
103
+ sign = -1;
104
+ x = -x;
105
+ }
106
+
107
+ // A&S formula 7.1.26
108
+ var t = 1 / (1 + p * x);
109
+ return sign * (
110
+ 1 - (((((a5 * t + a4) * t) + a3) * t + a2) * t + a1)
111
+ * t * Math.exp(-x * x));
112
+ };
113
+ science.stats.phi = function(x) {
114
+ return .5 * (1 + science.stats.erf(x / Math.SQRT2));
115
+ };
116
+ // See <http://en.wikipedia.org/wiki/Kernel_(statistics)>.
117
+ science.stats.kernel = {
118
+ uniform: function(u) {
119
+ if (u <= 1 && u >= -1) return .5;
120
+ return 0;
121
+ },
122
+ triangular: function(u) {
123
+ if (u <= 1 && u >= -1) return 1 - Math.abs(u);
124
+ return 0;
125
+ },
126
+ epanechnikov: function(u) {
127
+ if (u <= 1 && u >= -1) return .75 * (1 - u * u);
128
+ return 0;
129
+ },
130
+ quartic: function(u) {
131
+ if (u <= 1 && u >= -1) {
132
+ var tmp = 1 - u * u;
133
+ return (15 / 16) * tmp * tmp;
134
+ }
135
+ return 0;
136
+ },
137
+ triweight: function(u) {
138
+ if (u <= 1 && u >= -1) {
139
+ var tmp = 1 - u * u;
140
+ return (35 / 32) * tmp * tmp * tmp;
141
+ }
142
+ return 0;
143
+ },
144
+ gaussian: function(u) {
145
+ return 1 / Math.sqrt(2 * Math.PI) * Math.exp(-.5 * u * u);
146
+ },
147
+ cosine: function(u) {
148
+ if (u <= 1 && u >= -1) return Math.PI / 4 * Math.cos(Math.PI / 2 * u);
149
+ return 0;
150
+ }
151
+ };
152
+ // http://exploringdata.net/den_trac.htm
153
+ science.stats.kde = function() {
154
+ var kernel = science.stats.kernel.gaussian,
155
+ sample = [],
156
+ bandwidth = science.stats.bandwidth.nrd;
157
+
158
+ function kde(points, i) {
159
+ var bw = bandwidth.call(this, sample);
160
+ return points.map(function(x) {
161
+ var i = -1,
162
+ y = 0,
163
+ n = sample.length;
164
+ while (++i < n) {
165
+ y += kernel((x - sample[i]) / bw);
166
+ }
167
+ return [x, y / bw / n];
168
+ });
169
+ }
170
+
171
+ kde.kernel = function(x) {
172
+ if (!arguments.length) return kernel;
173
+ kernel = x;
174
+ return kde;
175
+ };
176
+
177
+ kde.sample = function(x) {
178
+ if (!arguments.length) return sample;
179
+ sample = x;
180
+ return kde;
181
+ };
182
+
183
+ kde.bandwidth = function(x) {
184
+ if (!arguments.length) return bandwidth;
185
+ bandwidth = science.functor(x);
186
+ return kde;
187
+ };
188
+
189
+ return kde;
190
+ };
191
+ // Based on figue implementation by Jean-Yves Delort.
192
+ // http://code.google.com/p/figue/
193
+ science.stats.kmeans = function() {
194
+ var distance = science.stats.distance.euclidean,
195
+ maxIterations = 1000,
196
+ k = 1;
197
+
198
+ function kmeans(vectors) {
199
+ var n = vectors.length,
200
+ assignments = [],
201
+ clusterSizes = [],
202
+ repeat = 1,
203
+ iterations = 0,
204
+ centroids = science_stats_kmeansRandom(k, vectors),
205
+ newCentroids,
206
+ i,
207
+ j,
208
+ x,
209
+ d,
210
+ min,
211
+ best;
212
+
213
+ while (repeat && iterations < maxIterations) {
214
+ // Assignment step.
215
+ j = -1; while (++j < k) {
216
+ clusterSizes[j] = 0;
217
+ }
218
+
219
+ i = -1; while (++i < n) {
220
+ x = vectors[i];
221
+ min = Infinity;
222
+ j = -1; while (++j < k) {
223
+ d = distance.call(this, centroids[j], x);
224
+ if (d < min) {
225
+ min = d;
226
+ best = j;
227
+ }
228
+ }
229
+ clusterSizes[assignments[i] = best]++;
230
+ }
231
+
232
+ // Update centroids step.
233
+ newCentroids = [];
234
+ i = -1; while (++i < n) {
235
+ x = assignments[i];
236
+ d = newCentroids[x];
237
+ if (d == null) newCentroids[x] = vectors[i].slice();
238
+ else {
239
+ j = -1; while (++j < d.length) {
240
+ d[j] += vectors[i][j];
241
+ }
242
+ }
243
+ }
244
+ j = -1; while (++j < k) {
245
+ x = newCentroids[j];
246
+ d = 1 / clusterSizes[j];
247
+ i = -1; while (++i < x.length) x[i] *= d;
248
+ }
249
+
250
+ // Check convergence.
251
+ repeat = 0;
252
+ j = -1; while (++j < k) {
253
+ if (!science_stats_kmeansCompare(newCentroids[j], centroids[j])) {
254
+ repeat = 1;
255
+ break;
256
+ }
257
+ }
258
+ centroids = newCentroids;
259
+ iterations++;
260
+ }
261
+ return {assignments: assignments, centroids: centroids};
262
+ }
263
+
264
+ kmeans.k = function(x) {
265
+ if (!arguments.length) return k;
266
+ k = x;
267
+ return kmeans;
268
+ };
269
+
270
+ kmeans.distance = function(x) {
271
+ if (!arguments.length) return distance;
272
+ distance = x;
273
+ return kmeans;
274
+ };
275
+
276
+ return kmeans;
277
+ };
278
+
279
+ function science_stats_kmeansCompare(a, b) {
280
+ if (!a || !b || a.length !== b.length) return false;
281
+ var n = a.length,
282
+ i = -1;
283
+ while (++i < n) if (a[i] !== b[i]) return false;
284
+ return true;
285
+ }
286
+
287
+ // Returns an array of k distinct vectors randomly selected from the input
288
+ // array of vectors. Returns null if k > n or if there are less than k distinct
289
+ // objects in vectors.
290
+ function science_stats_kmeansRandom(k, vectors) {
291
+ var n = vectors.length;
292
+ if (k > n) return null;
293
+
294
+ var selected_vectors = [];
295
+ var selected_indices = [];
296
+ var tested_indices = {};
297
+ var tested = 0;
298
+ var selected = 0;
299
+ var i,
300
+ vector,
301
+ select;
302
+
303
+ while (selected < k) {
304
+ if (tested === n) return null;
305
+
306
+ var random_index = Math.floor(Math.random() * n);
307
+ if (random_index in tested_indices) continue;
308
+
309
+ tested_indices[random_index] = 1;
310
+ tested++;
311
+ vector = vectors[random_index];
312
+ select = true;
313
+ for (i = 0; i < selected; i++) {
314
+ if (science_stats_kmeansCompare(vector, selected_vectors[i])) {
315
+ select = false;
316
+ break;
317
+ }
318
+ }
319
+ if (select) {
320
+ selected_vectors[selected] = vector;
321
+ selected_indices[selected] = random_index;
322
+ selected++;
323
+ }
324
+ }
325
+ return selected_vectors;
326
+ }
327
+ science.stats.hcluster = function() {
328
+ var distance = science.stats.distance.euclidean,
329
+ linkage = "simple"; // simple, complete or average
330
+
331
+ function hcluster(vectors) {
332
+ var n = vectors.length,
333
+ dMin = [],
334
+ cSize = [],
335
+ distMatrix = [],
336
+ clusters = [],
337
+ c1,
338
+ c2,
339
+ c1Cluster,
340
+ c2Cluster,
341
+ p,
342
+ root,
343
+ i,
344
+ j;
345
+
346
+ // Initialise distance matrix and vector of closest clusters.
347
+ i = -1; while (++i < n) {
348
+ dMin[i] = 0;
349
+ distMatrix[i] = [];
350
+ j = -1; while (++j < n) {
351
+ distMatrix[i][j] = i === j ? Infinity : distance(vectors[i] , vectors[j]);
352
+ if (distMatrix[i][dMin[i]] > distMatrix[i][j]) dMin[i] = j;
353
+ }
354
+ }
355
+
356
+ // create leaves of the tree
357
+ i = -1; while (++i < n) {
358
+ clusters[i] = [];
359
+ clusters[i][0] = {
360
+ left: null,
361
+ right: null,
362
+ dist: 0,
363
+ centroid: vectors[i],
364
+ size: 1,
365
+ depth: 0
366
+ };
367
+ cSize[i] = 1;
368
+ }
369
+
370
+ // Main loop
371
+ for (p = 0; p < n-1; p++) {
372
+ // find the closest pair of clusters
373
+ c1 = 0;
374
+ for (i = 0; i < n; i++) {
375
+ if (distMatrix[i][dMin[i]] < distMatrix[c1][dMin[c1]]) c1 = i;
376
+ }
377
+ c2 = dMin[c1];
378
+
379
+ // create node to store cluster info
380
+ c1Cluster = clusters[c1][0];
381
+ c2Cluster = clusters[c2][0];
382
+
383
+ newCluster = {
384
+ left: c1Cluster,
385
+ right: c2Cluster,
386
+ dist: distMatrix[c1][c2],
387
+ centroid: calculateCentroid(c1Cluster.size, c1Cluster.centroid,
388
+ c2Cluster.size, c2Cluster.centroid),
389
+ size: c1Cluster.size + c2Cluster.size,
390
+ depth: 1 + Math.max(c1Cluster.depth, c2Cluster.depth)
391
+ };
392
+ clusters[c1].splice(0, 0, newCluster);
393
+ cSize[c1] += cSize[c2];
394
+
395
+ // overwrite row c1 with respect to the linkage type
396
+ for (j = 0; j < n; j++) {
397
+ switch (linkage) {
398
+ case "single":
399
+ if (distMatrix[c1][j] > distMatrix[c2][j])
400
+ distMatrix[j][c1] = distMatrix[c1][j] = distMatrix[c2][j];
401
+ break;
402
+ case "complete":
403
+ if (distMatrix[c1][j] < distMatrix[c2][j])
404
+ distMatrix[j][c1] = distMatrix[c1][j] = distMatrix[c2][j];
405
+ break;
406
+ case "average":
407
+ distMatrix[j][c1] = distMatrix[c1][j] = (cSize[c1] * distMatrix[c1][j] + cSize[c2] * distMatrix[c2][j]) / (cSize[c1] + cSize[j]);
408
+ break;
409
+ }
410
+ }
411
+ distMatrix[c1][c1] = Infinity;
412
+
413
+ // infinity ­out old row c2 and column c2
414
+ for (i = 0; i < n; i++)
415
+ distMatrix[i][c2] = distMatrix[c2][i] = Infinity;
416
+
417
+ // update dmin and replace ones that previous pointed to c2 to point to c1
418
+ for (j = 0; j < n; j++) {
419
+ if (dMin[j] == c2) dMin[j] = c1;
420
+ if (distMatrix[c1][j] < distMatrix[c1][dMin[c1]]) dMin[c1] = j;
421
+ }
422
+
423
+ // keep track of the last added cluster
424
+ root = newCluster;
425
+ }
426
+
427
+ return root;
428
+ }
429
+
430
+ hcluster.distance = function(x) {
431
+ if (!arguments.length) return distance;
432
+ distance = x;
433
+ return hcluster;
434
+ };
435
+
436
+ return hcluster;
437
+ };
438
+
439
+ function calculateCentroid(c1Size, c1Centroid, c2Size, c2Centroid) {
440
+ var newCentroid = [],
441
+ newSize = c1Size + c2Size,
442
+ n = c1Centroid.length,
443
+ i = -1;
444
+ while (++i < n) {
445
+ newCentroid[i] = (c1Size * c1Centroid[i] + c2Size * c2Centroid[i]) / newSize;
446
+ }
447
+ return newCentroid;
448
+ }
449
+ science.stats.iqr = function(x) {
450
+ var quartiles = science.stats.quantiles(x, [.25, .75]);
451
+ return quartiles[1] - quartiles[0];
452
+ };
453
+ // Based on org.apache.commons.math.analysis.interpolation.LoessInterpolator
454
+ // from http://commons.apache.org/math/
455
+ science.stats.loess = function() {
456
+ var bandwidth = .3,
457
+ robustnessIters = 2,
458
+ accuracy = 1e-12;
459
+
460
+ function smooth(xval, yval, weights) {
461
+ var n = xval.length,
462
+ i;
463
+
464
+ if (n !== yval.length) throw {error: "Mismatched array lengths"};
465
+ if (n == 0) throw {error: "At least one point required."};
466
+
467
+ if (arguments.length < 3) {
468
+ weights = [];
469
+ i = -1; while (++i < n) weights[i] = 1;
470
+ }
471
+
472
+ science_stats_loessFiniteReal(xval);
473
+ science_stats_loessFiniteReal(yval);
474
+ science_stats_loessFiniteReal(weights);
475
+ science_stats_loessStrictlyIncreasing(xval);
476
+
477
+ if (n == 1) return [yval[0]];
478
+ if (n == 2) return [yval[0], yval[1]];
479
+
480
+ var bandwidthInPoints = Math.floor(bandwidth * n);
481
+
482
+ if (bandwidthInPoints < 2) throw {error: "Bandwidth too small."};
483
+
484
+ var res = [],
485
+ residuals = [],
486
+ robustnessWeights = [];
487
+
488
+ // Do an initial fit and 'robustnessIters' robustness iterations.
489
+ // This is equivalent to doing 'robustnessIters+1' robustness iterations
490
+ // starting with all robustness weights set to 1.
491
+ i = -1; while (++i < n) {
492
+ res[i] = 0;
493
+ residuals[i] = 0;
494
+ robustnessWeights[i] = 1;
495
+ }
496
+
497
+ var iter = -1;
498
+ while (++iter <= robustnessIters) {
499
+ var bandwidthInterval = [0, bandwidthInPoints - 1];
500
+ // At each x, compute a local weighted linear regression
501
+ var x;
502
+ i = -1; while (++i < n) {
503
+ x = xval[i];
504
+
505
+ // Find out the interval of source points on which
506
+ // a regression is to be made.
507
+ if (i > 0) {
508
+ science_stats_loessUpdateBandwidthInterval(xval, weights, i, bandwidthInterval);
509
+ }
510
+
511
+ var ileft = bandwidthInterval[0],
512
+ iright = bandwidthInterval[1];
513
+
514
+ // Compute the point of the bandwidth interval that is
515
+ // farthest from x
516
+ var edge = (xval[i] - xval[ileft]) > (xval[iright] - xval[i]) ? ileft : iright;
517
+
518
+ // Compute a least-squares linear fit weighted by
519
+ // the product of robustness weights and the tricube
520
+ // weight function.
521
+ // See http://en.wikipedia.org/wiki/Linear_regression
522
+ // (section "Univariate linear case")
523
+ // and http://en.wikipedia.org/wiki/Weighted_least_squares
524
+ // (section "Weighted least squares")
525
+ var sumWeights = 0,
526
+ sumX = 0,
527
+ sumXSquared = 0,
528
+ sumY = 0,
529
+ sumXY = 0,
530
+ denom = Math.abs(1 / (xval[edge] - x));
531
+
532
+ for (var k = ileft; k <= iright; ++k) {
533
+ var xk = xval[k],
534
+ yk = yval[k],
535
+ dist = k < i ? x - xk : xk - x,
536
+ w = science_stats_loessTricube(dist * denom) * robustnessWeights[k] * weights[k],
537
+ xkw = xk * w;
538
+ sumWeights += w;
539
+ sumX += xkw;
540
+ sumXSquared += xk * xkw;
541
+ sumY += yk * w;
542
+ sumXY += yk * xkw;
543
+ }
544
+
545
+ var meanX = sumX / sumWeights,
546
+ meanY = sumY / sumWeights,
547
+ meanXY = sumXY / sumWeights,
548
+ meanXSquared = sumXSquared / sumWeights;
549
+
550
+ var beta = (Math.sqrt(Math.abs(meanXSquared - meanX * meanX)) < accuracy)
551
+ ? 0 : ((meanXY - meanX * meanY) / (meanXSquared - meanX * meanX));
552
+
553
+ var alpha = meanY - beta * meanX;
554
+
555
+ res[i] = beta * x + alpha;
556
+ residuals[i] = Math.abs(yval[i] - res[i]);
557
+ }
558
+
559
+ // No need to recompute the robustness weights at the last
560
+ // iteration, they won't be needed anymore
561
+ if (iter === robustnessIters) {
562
+ break;
563
+ }
564
+
565
+ // Recompute the robustness weights.
566
+
567
+ // Find the median residual.
568
+ var sortedResiduals = residuals.slice();
569
+ sortedResiduals.sort();
570
+ var medianResidual = sortedResiduals[Math.floor(n / 2)];
571
+
572
+ if (Math.abs(medianResidual) < accuracy)
573
+ break;
574
+
575
+ var arg,
576
+ w;
577
+ i = -1; while (++i < n) {
578
+ arg = residuals[i] / (6 * medianResidual);
579
+ robustnessWeights[i] = (arg >= 1) ? 0 : ((w = 1 - arg * arg) * w);
580
+ }
581
+ }
582
+
583
+ return res;
584
+ }
585
+
586
+ smooth.bandwidth = function(x) {
587
+ if (!arguments.length) return x;
588
+ bandwidth = x;
589
+ return smooth;
590
+ };
591
+
592
+ smooth.robustnessIterations = function(x) {
593
+ if (!arguments.length) return x;
594
+ robustnessIters = x;
595
+ return smooth;
596
+ };
597
+
598
+ smooth.accuracy = function(x) {
599
+ if (!arguments.length) return x;
600
+ accuracy = x;
601
+ return smooth;
602
+ };
603
+
604
+ return smooth;
605
+ };
606
+
607
+ function science_stats_loessFiniteReal(values) {
608
+ var n = values.length,
609
+ i = -1;
610
+
611
+ while (++i < n) if (!isFinite(values[i])) return false;
612
+
613
+ return true;
614
+ }
615
+
616
+ function science_stats_loessStrictlyIncreasing(xval) {
617
+ var n = xval.length,
618
+ i = 0;
619
+
620
+ while (++i < n) if (xval[i - 1] >= xval[i]) return false;
621
+
622
+ return true;
623
+ }
624
+
625
+ // Compute the tricube weight function.
626
+ // http://en.wikipedia.org/wiki/Local_regression#Weight_function
627
+ function science_stats_loessTricube(x) {
628
+ return (x = 1 - x * x * x) * x * x;
629
+ }
630
+
631
+ // Given an index interval into xval that embraces a certain number of
632
+ // points closest to xval[i-1], update the interval so that it embraces
633
+ // the same number of points closest to xval[i], ignoring zero weights.
634
+ function science_stats_loessUpdateBandwidthInterval(
635
+ xval, weights, i, bandwidthInterval) {
636
+
637
+ var left = bandwidthInterval[0],
638
+ right = bandwidthInterval[1];
639
+
640
+ // The right edge should be adjusted if the next point to the right
641
+ // is closer to xval[i] than the leftmost point of the current interval
642
+ var nextRight = science_stats_loessNextNonzero(weights, right);
643
+ if ((nextRight < xval.length) && (xval[nextRight] - xval[i]) < (xval[i] - xval[left])) {
644
+ var nextLeft = science_stats_loessNextNonzero(weights, left);
645
+ bandwidthInterval[0] = nextLeft;
646
+ bandwidthInterval[1] = nextRight;
647
+ }
648
+ }
649
+
650
+ function science_stats_loessNextNonzero(weights, i) {
651
+ var j = i + 1;
652
+ while (j < weights.length && weights[j] === 0) j++;
653
+ return j;
654
+ }
655
+ // Welford's algorithm.
656
+ science.stats.mean = function(x) {
657
+ var n = x.length;
658
+ if (n === 0) return NaN;
659
+ var m = 0,
660
+ i = -1;
661
+ while (++i < n) m += (x[i] - m) / (i + 1);
662
+ return m;
663
+ };
664
+ science.stats.median = function(x) {
665
+ return science.stats.quantiles(x, [.5])[0];
666
+ };
667
+ science.stats.mode = function(x) {
668
+ x = x.slice().sort(science.ascending);
669
+ var mode,
670
+ n = x.length,
671
+ i = -1,
672
+ l = i,
673
+ last = null,
674
+ max = 0,
675
+ tmp,
676
+ v;
677
+ while (++i < n) {
678
+ if ((v = x[i]) !== last) {
679
+ if ((tmp = i - l) > max) {
680
+ max = tmp;
681
+ mode = last;
682
+ }
683
+ last = v;
684
+ l = i;
685
+ }
686
+ }
687
+ return mode;
688
+ };
689
+ // Uses R's quantile algorithm type=7.
690
+ science.stats.quantiles = function(d, quantiles) {
691
+ d = d.slice().sort(science.ascending);
692
+ var n_1 = d.length - 1;
693
+ return quantiles.map(function(q) {
694
+ if (q === 0) return d[0];
695
+ else if (q === 1) return d[n_1];
696
+
697
+ var index = 1 + q * n_1,
698
+ lo = Math.floor(index),
699
+ h = index - lo,
700
+ a = d[lo - 1];
701
+
702
+ return h === 0 ? a : a + h * (d[lo] - a);
703
+ });
704
+ };
705
+ // Unbiased estimate of a sample's variance.
706
+ // Also known as the sample variance, where the denominator is n - 1.
707
+ science.stats.variance = function(x) {
708
+ var n = x.length;
709
+ if (n < 1) return NaN;
710
+ if (n === 1) return 0;
711
+ var mean = science.stats.mean(x),
712
+ i = -1,
713
+ s = 0;
714
+ while (++i < n) {
715
+ var v = x[i] - mean;
716
+ s += v * v;
717
+ }
718
+ return s / (n - 1);
719
+ };
720
+ })()