summaryrefslogtreecommitdiff
path: root/research/lesci/analyze.go
blob: 5fc3b2de0102d7e9f187b2c22baf77073acb1b43 (plain)
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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
package lesci

import (
	"context"
	"database/sql"
	"fmt"
)

var (
	sqlPairs = `
		CREATE TEMP TABLE pairs AS
		SELECT uid, pos, token, logprob,
		       LEAD(token, 1) OVER w AS next_tok,
		       LAG(token, 1)  OVER w AS prev_tok
		FROM context
		WINDOW w AS (PARTITION BY uid ORDER BY pos)`

	sqlMatched = `
		CREATE TEMP TABLE matched AS
		SELECT p.*, r1.c AS merged_as_first, r2.c AS merged_as_second
		FROM pairs p
		LEFT JOIN oov_rules r1 ON p.token = r1.a AND p.next_tok = r1.b
		LEFT JOIN oov_rules r2 ON p.prev_tok = r2.a AND p.token = r2.b`

	sqlSummedOOV = `
		CREATE TEMP TABLE summed_oov AS
		WITH processed AS (
		    SELECT
		        uid,
		        logprob,
		        COALESCE(merged_as_first, merged_as_second, token) AS tok,
		        CASE
		            WHEN merged_as_first  IS NOT NULL THEN pos
		            WHEN merged_as_second IS NOT NULL THEN pos - 1
		            ELSE pos
		        END AS adjusted_pos
		    FROM matched
		)
		SELECT
		    uid,
		    tok,
		    adjusted_pos AS pos,
		    SUM(logprob) AS token_logprob
		FROM processed
		WHERE tok >= ? AND tok < ?
		GROUP BY uid, tok, adjusted_pos`

	sqlResults = `
		CREATE TABLE lesci_results AS
		SELECT tok,
		       QUANTILE_CONT(token_logprob, 0.5)  AS median,
		       QUANTILE_CONT(token_logprob, 0.75) AS q75,
		       QUANTILE_CONT(token_logprob, 0.25) AS q25,
		       QUANTILE_CONT(token_logprob, 0.75)
		         - QUANTILE_CONT(token_logprob, 0.25) AS iqr,
		       AVG(token_logprob)        AS mean,
		       STDDEV_SAMP(token_logprob) AS std,
		       COUNT(*)                  AS num,
		       tok < ?                   AS treat
		FROM summed_oov
		GROUP BY tok`
)

func (e *Experiment) Analyze(db *sql.DB) error {
	ctx := context.Background()

	var tx *sql.Tx

	if t, err := db.BeginTx(ctx, nil); err != nil {
		return err
	} else {
		tx = t
	}

	defer tx.Rollback()

	if _, err := tx.ExecContext(ctx, `DROP TABLE IF EXISTS lesci_results`); err != nil {
		return err
	}

	if _, err := tx.ExecContext(ctx, sqlPairs); err != nil {
		return fmt.Errorf("pairs: %w", err)
	}

	if _, err := tx.ExecContext(ctx, sqlMatched); err != nil {
		return fmt.Errorf("matched: %w", err)
	}

	low := e.cutoff - e.window
	high := e.cutoff + e.window

	if e.window == -1 {
		low = 0
		high = 10000000 // TODO use counterfactual
	}

	if _, err := tx.ExecContext(ctx, sqlSummedOOV, low, high); err != nil {
		return fmt.Errorf("summed_oov: %w", err)
	}

	if _, err := tx.ExecContext(ctx, sqlResults, e.cutoff); err != nil {
		return fmt.Errorf("lesci_results: %w", err)
	}

	if err := numSamples(ctx, tx); err != nil {
		return err
	}

	return tx.Commit()
}

func numSamples(ctx context.Context, tx *sql.Tx) error {
	var rows *sql.Rows

	if r, err := tx.QueryContext(ctx, `SELECT treat, COUNT(*), SUM(num) FROM lesci_results GROUP BY treat ORDER BY treat DESC`); err != nil {
		return err
	} else {
		rows = r
	}

	defer rows.Close()

	for rows.Next() {
		var treat bool
		var tokens, occurrences int

		if err := rows.Scan(&treat, &tokens, &occurrences); err != nil {
			return err
		}

		side := "right"

		if treat {
			side = "left"
		}

		fmt.Printf("%s: %d tokens, %d occurrences\n", side, tokens, occurrences)
	}

	return rows.Err()
}