File: stats.ml

package info (click to toggle)
ocamlgsl 0.3.5-3
  • links: PTS
  • area: main
  • in suites: sarge
  • size: 3,444 kB
  • ctags: 2,901
  • sloc: ml: 7,956; ansic: 6,796; makefile: 303; sh: 87; awk: 13
file content (28 lines) | stat: -rw-r--r-- 1,195 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28

let _ = 
  let data = [| 17.2; 18.1; 16.5; 18.3; 12.6 |] in
  let mean = Gsl_stats.mean data in
  let variance = Gsl_stats.variance data in
  let largest  = Gsl_stats.max data in
  let smallest = Gsl_stats.min data in
  Printf.printf "The dataset is %g, %g, %g, %g, %g\n"
    data.(0) data.(1) data.(2) data.(3) data.(4) ;
  Printf.printf "The sample mean is %g\n" mean ;
  Printf.printf "The estimated variance is %g\n"  variance ;
  Printf.printf "The largest value is %g\n"  largest ;
  Printf.printf "The smallest value is %g\n"  smallest

let _ = 
  let data = [| 17.2; 18.1; 16.5; 18.3; 12.6 |] in
  Printf.printf "Original dataset is %g, %g, %g, %g, %g\n"
    data.(0) data.(1) data.(2) data.(3) data.(4) ;
  Array.sort compare data ;
  Printf.printf "Sorted dataset is %g, %g, %g, %g, %g\n"
    data.(0) data.(1) data.(2) data.(3) data.(4) ;
  let median = Gsl_stats.quantile_from_sorted_data data 0.5 in
  let upperq = Gsl_stats.quantile_from_sorted_data data 0.75 in
  let lowerq = Gsl_stats.quantile_from_sorted_data data 0.25 in
  Printf.printf "The median is %g\n" median ;
  Printf.printf "The upper quartile is %g\n" upperq ;
  Printf.printf "The lower quartile is %g\n" lowerq