-
Notifications
You must be signed in to change notification settings - Fork 32
/
pair-fmi.gplot
73 lines (61 loc) · 2.21 KB
/
pair-fmi.gplot
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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
# set term png small color
#################################################
#
set term png
set out 'pair-fmi.png'
# set term postscript eps enhanced color 20
# set out 'pair-fmi.eps'
set key right top
set title "MI Distribution of word pairs"
set xlabel "MI"
set ylabel "Frequency"
set xrange [-12:27]
set logscale y
set format y "10^{%L}"
set yrange [4e-8:5e-2]
g(x,m,s) = (1/sqrt(2*3.14159*s)) * exp(- (x-m)**2 / (2*s**2))
norm = 1.0 / 28184319
plot "pair-fmi.dat" using 2:(norm * $3) with lines lw 2 title "data_{ }", \
"pair-fmi.dat" using 2:(0.017*2**(-0.35*$2)) with lines lw 2 title "2^{-0.35 MI}", \
"pair-fmi.dat" using 2:(0.2*2**(2.0*$2)) with lines lw 2 title "4^{MI}"
# "pair-fmi.dat" using 2:(0.04*g($2,0.4,2)) with lines lw 2 title "Gauss(0.4,2)", \
#################################################
#
set term png
set out 'pair-fmi-signal.png'
# set term postscript eps enhanced color 20
# set out 'pair-fmi-signal.eps'
set key right top
set title "MI Distribution of Word Pairs"
set xlabel "MI"
set ylabel "Frequency"
set xrange [-12:27]
set logscale y
set format y "10^{%L}"
set yrange [4e-8:5e-2]
g(x,m,s) = (1/sqrt(2*3.14159*s)) * exp(- (x-m)**2 / (2*s**2))
lnm(x,m,s) = (1/(x*s*sqrt(2*3.14159))) * exp(- (log(x)-m)**2 / (2*s**2))
norm = 1.0 / 28184319
plot "pair-fmi.dat" using 2:(norm * $3) with lines lw 2 title "data_{ }", \
"pair-fmi.dat" using 2:(0.033*g($2,0.32,2)) with lines lw 2 title "Gauss(0.4,2)", \
"pair-fmi.dat" using 2:(0.028*lnm($2,1.85,0.54)) with lines lw 2 title "Log Normal(1.85,0.54)", \
#################################################
#
set term png
set out 'pair-rmi.png'
#
set term postscript eps enhanced color 20
set out 'pair-rmi.eps'
set key right top
set title "Ranked-MI Distribution of word pairs"
set xlabel "Ranked-MI"
set ylabel "Frequency"
set xrange [-30:30]
set logscale y
set format y "10^{%L}"
set yrange [4e-8:2e-2]
norm = 1.0 / 28184319
plot "pair-rmi.dat" using 2:(norm * $3) with lines lw 2 title "Ranked-MI data_{ }", \
"pair-rmi.dat" using 2:(0.0033*2**(-0.24*$2)) with lines lw 2 title "2^{-0.24 MI}", \
"pair-fmi.dat" using (1.5*$2-10):(norm * $3) with lines lw 2 title "MI_{shifted}"
#################################################