1 use plotters::prelude::*;
2
3 use rand::SeedableRng;
4 use rand_distr::{Distribution, Normal};
5 use rand_xorshift::XorShiftRng;
6
main() -> Result<(), Box<dyn std::error::Error>>7 fn main() -> Result<(), Box<dyn std::error::Error>> {
8 let root =
9 BitMapBackend::new("plotters-doc-data/normal-dist.png", (1024, 768)).into_drawing_area();
10
11 root.fill(&WHITE)?;
12
13 let sd = 0.13;
14
15 let random_points: Vec<(f64, f64)> = {
16 let norm_dist = Normal::new(0.5, sd).unwrap();
17 let mut x_rand = XorShiftRng::from_seed(*b"MyFragileSeed123");
18 let mut y_rand = XorShiftRng::from_seed(*b"MyFragileSeed321");
19 let x_iter = norm_dist.sample_iter(&mut x_rand);
20 let y_iter = norm_dist.sample_iter(&mut y_rand);
21 x_iter.zip(y_iter).take(5000).collect()
22 };
23
24 let areas = root.split_by_breakpoints([944], [80]);
25
26 let mut x_hist_ctx = ChartBuilder::on(&areas[0])
27 .y_label_area_size(40)
28 .build_cartesian_2d((0.0..1.0).step(0.01).use_round().into_segmented(), 0..250)?;
29 let mut y_hist_ctx = ChartBuilder::on(&areas[3])
30 .x_label_area_size(40)
31 .build_cartesian_2d(0..250, (0.0..1.0).step(0.01).use_round())?;
32 let mut scatter_ctx = ChartBuilder::on(&areas[2])
33 .x_label_area_size(40)
34 .y_label_area_size(40)
35 .build_cartesian_2d(0f64..1f64, 0f64..1f64)?;
36 scatter_ctx
37 .configure_mesh()
38 .disable_x_mesh()
39 .disable_y_mesh()
40 .draw()?;
41 scatter_ctx.draw_series(
42 random_points
43 .iter()
44 .map(|(x, y)| Circle::new((*x, *y), 2, GREEN.filled())),
45 )?;
46 let x_hist = Histogram::vertical(&x_hist_ctx)
47 .style(GREEN.filled())
48 .margin(0)
49 .data(random_points.iter().map(|(x, _)| (*x, 1)));
50 let y_hist = Histogram::horizontal(&y_hist_ctx)
51 .style(GREEN.filled())
52 .margin(0)
53 .data(random_points.iter().map(|(_, y)| (*y, 1)));
54 x_hist_ctx.draw_series(x_hist)?;
55 y_hist_ctx.draw_series(y_hist)?;
56
57 Ok(())
58 }
59 #[test]
entry_point()60 fn entry_point() {
61 main().unwrap()
62 }
63
