 speed dist
set loc speed qua
set loc dist qua
read cases inline
1 4 2
2 4 10
3 7 4
4 7 22
5 8 16
6 9 10
7 10 18
8 10 26
9 10 34
10 11 17
11 11 28
12 12 14
13 12 20
14 12 24
15 12 28
16 13 26
17 13 34
18 13 34
19 13 46
20 14 26
21 14 36
22 14 60
23 14 80
24 15 20
25 15 26
26 15 54
27 16 32
28 16 40
29 17 32
30 17 40
31 17 50
32 18 42
33 18 56
34 18 76
35 18 84
36 19 36
37 19 46
38 19 68
39 20 32
40 20 48
41 20 52
42 20 56
43 20 64
44 22 66
45 23 54
46 24 70
47 24 92
48 24 93
49 24 120
50 25 85
;;;;
run
eval (require "plot")
eval (define speed (data "speed"))
eval (define dist (data "dist"))
eval (define ys (lowess speed dist 'frac (/ 2 3)))
eval (define ys2 (lowess speed dist 'frac 0.2))
eval (define fig1 (create-plot speed dist 500 "speed" "dist"))
eval (add-points fig1 speed dist "red")
eval (add-lines fig1 speed ys "blue" 0)
eval (add-lines fig1 speed ys2 "blue" 1)
eval (plot-title fig1 "Lowess for cars dataset")
eval (add-lines fig1 '(2 5) '(130 130) "blue" 0)
eval (add-lines fig1 '(2 6) '(120 120) "blue" 1)
eval (add-text fig1 6 130 "f=2/3" 14)
eval (add-text fig1 6 120 "f=1/5" 14)
