Statistical Methods in Linguistics
University of Rochester
August 31, September 2 and 9, 2026
Probability and statistics is all about describing what happens precisely.
That means stating what can happen and how often it happens.
August 31
I am an associate professor in Linguistics and Computer Science.
My office is 511A Lattimore Hall. My website is aaronstevenwhite.io.
I study how words, constructions, and discourses encode events, attitudes, and inferences.
That includes work on how verb meanings relate to syntactic selection, how event structure can be inferred from semantic annotations, and how uncertain inferences update a discourse.
My data include acceptability and inference judgments, corpus distributions, and decompositional semantic annotations.
I use probability and statistics to connect patterns in those observations to hypotheses about the linguistic representations that could have produced them.
Probability and statistics is all about describing what happens precisely.
That means stating what can happen and how often it happens.
What can happen? A speaker can produce many tongue and lip trajectories while articulating a vowel.
How often? We can describe how those trajectories are distributed across vowels, speakers, and speaking contexts.
What can happen? A vowel token can have different durations, formant values, and formant trajectories.
How often? We can describe how those measurements are distributed across vowels, speakers, and speaking contexts.
What can happen? A comprehender can spend different amounts of time reading a word or sentence region.
How often? We can describe how reading times vary with word frequency, syntactic structure, readers, and passages.
What can happen? A producer can choose the double object or prepositional dative construction to describe a transfer.
How often? We can describe how the probability of each choice varies with the theme, recipient, verb, and discourse context.
What can happen? A participant can choose the left endpoint, the right endpoint, or an interior slider value.
How often? We can describe the probability of each endpoint and the distribution of values in the interior.
What can happen? Each pairing of a participant and an item can yield a different response.
How often? We can describe how the response distribution varies across participants and items as well as across the conditions we manipulate.
What can happen? A design determines which conditions each participant and item can appear in and which observations can enter the analysis.
How often? Assignment probabilities and replication determine how often the comparisons supported by the design occur.
What can happen? An acoustic token can arise from one of several vowel categories that we do not observe directly.
How often? We can describe how frequently each category generates a token and how measurements are distributed within each category.
What can happen? A verb-frame combination can receive a range of acceptability judgments.
How often? We can describe judgments across verb-frame combinations and ask which recurring patterns support a structured acceptability lexicon.
What can happen? A language-feature combination can be observed or missing, and an observed combination can take different values.
How often? We can describe how missingness varies with sampling and documentation and how recorded values are distributed.
What can happen? An earlier response can determine which question an annotator sees next, and the annotator can then choose among the available responses.
How often? We can describe the probability of each path through the task and the responses given by different annotators to different items.
The thing we want to describe and the data we collect to learn about it need not be identical.
We may observe the verb in that frame in a corpus.
We may also observe a participant judge a sentence containing that verb and frame.
Neither observation is identical to the verb’s acceptability in that frame.
A corpus token tells us that a producer used the form in one context.
A judgment tells us how one participant responded to one sentence under one task.
Both may inform an account of acceptability, but through different processes.
The course describes both the data we collect and the processes that could have produced those data.
Fall 2026
Class meets from 12:30 PM to 1:45 PM in Lattimore 513.
The first meeting is Monday, August 31. The final class meeting is Monday, December 14.
Monday, October 26 is the review class.
Wednesday, October 28 is the midterm and the LING 414 proposal deadline.
We meet for two 75 minute sessions each week.
Plan for at least 480 minutes of work outside class each week for notes, readings, problem sets, and project work.
LING 110 with a grade of C minus or better is the prerequisite.
Each chapter introduces one principal concept and links back to the chapters it requires.
Readings appear in the chapter where they become useful, together with an explanation of what to take from them.
The course also uses linguistics papers when their data or arguments bear directly on the method.
Install R and either RStudio or Visual Studio Code.
Installation instructions and troubleshooting belong on the course Zulip.
| component | weight |
|---|---|
| five problem sets | 65 percent |
| midterm on October 28 | 25 percent |
| two office hour meetings | 10 percent |
| component | weight |
|---|---|
| five problem sets | 40 percent |
| midterm on October 28 | 20 percent |
| final project | 40 percent |
Each problem set develops an analysis that the notes prepare you to understand.
Working through the code and interpreting its output are part of learning the next portion of the course.
Students are selected at random to present and will not know in advance whether they are presenting.
All students must attend and be prepared to present on any assessment day following a due date.
The assessment concerns what the code does, why the analysis answers the question, and what the result means.
The final project rubric states the requirements for each option.
Office hours are by appointment through the scheduler on my website.
The two meetings may concern course material, a problem set, a possible project, or a statistical question from other work.
Do not use email for course communication.
Post general questions in the appropriate Zulip channel so that everyone can use the answer.
Use a direct message for grades, absences, accommodations, or other private matters.
I generally reply within one to two business days.
I do not monitor Zulip after 5 PM or on weekends.
A problem set may be submitted at most three days late.
Project components are not accepted late without prior approval. Ask for an extension at least 48 hours before the deadline, except in an emergency.
You may discuss problem sets with classmates, but each student must write and submit their own solution.
Shared discussion should help you understand the analysis well enough to explain it independently.
You may use generative AI for coding assistance on problem sets and projects.
You must understand the resulting code and cite the tool and the specific prompts in code comments. Generative AI may not be used for written assignments.
| A | A minus | B plus | B | B minus | C plus | C | C minus | D plus | D | E |
|---|---|---|---|---|---|---|---|---|---|---|
| 93 | 90 | 87 | 83 | 80 | 77 | 73 | 70 | 67 | 60 | below 60 |
Students seeking accommodations should contact the University Office of Disability Resources.
Faculty are mandatory reporters for Title IX matters. The syllabus links to the relevant confidential and nonconfidential resources.
Piantadosi 2014 and Zipf’s law
It relates each word type’s frequency to its frequency rank.
A corpus records a large collection of usages produced by many speakers and writers in many contexts.
We can ask for a precise description of that entire collection.
George Washington’s 1793 inaugural address is included with quanteda.
quanteda, data_corpus_inaugural
| word type | token count | frequency rank |
|---|---|---|
| the | 13 | 1 |
| of | 11 | 2 |
| i | 6 | 3 |
| to | 5 | 4 |
| in | 3 | 5 |
| shall | 3 | 6 |
Equal counts are ordered alphabetically here.
If r is a word’s frequency rank and f(r) is its frequency, then
f(r)\propto\frac{1}{r^\alpha},
where \alpha is often near 1 for word frequencies.
Piantadosi 2014, Equation 1
In the Washington address, f(1)=13. If \alpha=1, then
f(2)=6.5 \qquad\text{and}\qquad f(3)\approx4.33.
The observed counts are 11 and 6. The inverse relation describes the direction of the pattern, but it does not reproduce these counts.
Piantadosi treats the Zipf and Mandelbrot form as the working description:
f(r)\propto\frac{1}{(r+\beta)^\alpha}.
The additional parameter \beta changes the curve among the highest frequency words.
Piantadosi 2014, Equation 2
Both rank and frequency range over several orders of magnitude.
A logarithmic axis allocates comparable space to comparable ratios, such as 10 to 100 and 100 to 1,000.
The red curve is the Zipf and Mandelbrot description. The gray curve follows the local average more closely.
Piantadosi 2014, Figure 1a
If rank and frequency are estimated from the same corpus, their measurement errors are related.
Even equally probable word types will receive different observed ranks after chance differences in their counts.
The paper randomly splits the word tokens into two parts.
One part estimates each word’s rank. The other estimates its frequency.
A chance fluctuation in one half cannot affect both estimated rank and estimated frequency.
The vertical axis shows the difference between observed log frequency and the curve’s prediction.
Piantadosi 2014, Figure 1b
The differences form long runs above and below zero, including a large scoop among lower frequency words.
Even the highest frequency words show systematic local structure that the simple curve does not describe.
The full word frequency distribution is more complex than the large scale relation between frequency and rank.
Deriving a power law is not enough because many incompatible processes can yield a similar curve.
Piantadosi thus examines additional properties of word frequencies that distinguish possible explanations.
The common rank ordering is estimated across languages. Similar meanings tend to occupy similar portions of the frequency distribution.
Piantadosi 2014, Figure 2
The horizontal axis is cardinality rather than an independently assigned rank. Smaller numbers receive higher frequencies across English, Russian, and Italian.
Piantadosi 2014, Figure 3
Taboo words with roughly shared referential content still show highly unequal frequencies.
Social register and other aspects of use may thus matter alongside reference.
Piantadosi 2014, Figure 4
Months, planets, and elements give a language relatively little freedom to choose the relevant referents.
Piantadosi 2014, Figure 5
Syntactic category frequencies in the Penn Treebank also form a strongly unequal distribution.
Piantadosi 2014, Figure 6
Determiners, prepositions, modals, nouns, and two verb categories differ in their fitted curves and in their systematic departures from those curves.
Piantadosi 2014, Figure 7
The probability of using Dallas differs between a discussion of Lyndon Johnson and a discussion of Carl Sagan.
Topics, social groups, technologies, and historical events all change word frequencies.
The simple rank frequency curve describes an aggregate whose components need not follow one unchanging distribution.
Twenty-five participants wrote stories of at least 2,000 words using eight novel alien names. The average across participants was near Zipfian.
Piantadosi 2014, Figure 8
Piantadosi averages across participants after ordering each participant’s names by frequency.
Larger studies would be needed to establish the pattern reliably within individual participants.
Piantadosi reviews near Zipfian frequency distributions in music, computer programs, and internet systems.
Near-Zipfian distributions across unrelated domains could reflect a general process, or different processes could produce similar aggregate descriptions.
A process that emits characters independently and occasionally emits a space can produce a near Zipfian distribution of resulting strings.
Humans do not produce words by emitting independent characters until a space happens to occur.
Treating the letter e as a boundary in the American National Corpus still produces a near Zipfian curve.
Piantadosi 2014, Figure 9
A model can describe the aggregate pattern while misdescribing how speakers produce words.
If a word becomes more likely to recur after it has already occurred, a strongly unequal frequency distribution can develop.
But discourse topic may explain both the earlier use and the later reuse. Reuse itself need not be the underlying cause.
Meaning strongly predicts frequency, so an account of word frequency must represent meaning.
But meanings constrained by nature and newly introduced names also produce near Zipfian patterns. Semantic organization alone is thus incomplete.
Optimization models can yield Zipfian frequencies by balancing assumptions about speaker and listener costs.
The paper argues that the psychological assumptions and the required parameter values must be tested independently of the resulting curve.
General arguments from information, computation, or entropy may explain why power laws occur in many systems.
Without new predictions, the observed curve cannot distinguish these accounts from more specific psychological processes.
The novel name experiment suggests that memory may contribute to unequal reuse even without an established lexicon.
Piantadosi presents this as a possibility that requires a fuller model and new evidence.
What can happen? Each speaker or writer can choose a word type in a particular context.
How often? Zipf’s law describes the aggregate frequencies of those choices across an entire corpus.
Describing what can happen and how often it happens does not by itself identify the process that produced the pattern.
September 2
On August 31, we chose what one observation represents.
On September 2, we state what can happen to one represented observation and which distinctions receive probabilities.
Its interpretation depends on the range of possibilities under consideration.
Suppose one observation is one represented English personal pronoun token.
Which forms can the model produce as outcomes?
This module assumes the basic set theory reviewed in Languages as formal objects.
You should already be comfortable with:
\Omega=\{\textit{I},\textit{me},\textit{you},\textit{he},\textit{him}, \textit{she},\textit{her},\textit{we},\textit{us},\textit{they},\textit{them}\}
The sample space is the set of possible outcomes under the declared representation.
\textit{herself}\notin\Omega says that the current model omits herself.
It does not say that speakers cannot produce herself.
| represented unit | one possible outcome |
|---|---|
| vowel category | /i/ |
| vowel token | one production of /i/ |
| acoustic measurement | an F1 value |
| acoustic trajectory | a sequence of F1 values |
A vowel study asks whether formant trajectories differ by dialect.
Should one outcome be a vowel category, one F1 value, or a complete trajectory?
Which outcomes answer the question in the same way?
Several outcomes give the same answer.
The forms he, her, and them are distinct outcomes, but each is third person.
An event groups outcomes that make a proposition true.
T=\{\textit{he},\textit{him},\textit{she},\textit{her},\textit{they},\textit{them}\}
Let A be the event that the pronoun is accusative.
The outcome them belongs to both T and A.
It is still one outcome.
If T is third person and A is accusative, then
T\cap A=\{\textit{him},\textit{her},\textit{them}\}.
The entire sample space \Omega is the event that some represented outcome occurs.
The empty set \varnothing is the event that no represented outcome occurs.
Let P be plural and A be accusative.
Write P\cap A, P\cup A, and P^c for the four outcome pronoun space \{\textit{he},\textit{him},\textit{they},\textit{them}\}.
Sigma-algebras
Suppose the model records plural number but not case.
It can distinguish plural from nonplural forms without distinguishing they from them.
A sigma-algebra is the collection of events to which a probability model assigns probabilities.
We will write this collection as \mathcal F.
\Omega\in\mathcal F
The model must be able to measure the event that some represented outcome occurs.
If B\in\mathcal F, then
B^c\in\mathcal F.
Measuring “plural” requires measuring “not plural.”
If B_1,B_2,\ldots\in\mathcal F, then
\bigcup_i B_i\in\mathcal F.
Measurable alternatives remain measurable when combined.
For plural P in \Omega_4=\{\textit{he},\textit{him},\textit{they},\textit{them}\},
\mathcal F_P=\{\varnothing,P,P^c,\Omega_4\}.
The collection \{\Omega_4,P\} is not a sigma-algebra.
It omits P^c and \varnothing.
\Omega states which outcomes exist in the representation.
\mathcal F states which groupings of those outcomes can receive probabilities.
If the event “accusative” is not in \mathcal F_P, can the model assign a probability to accusative form?
What additional distinction must the event space retain?
Begin from the distinctions the question needs
Let P be plural and A be accusative.
We want the smallest sigma-algebra that contains both events.
The collection \{P,A\} is a generating set.
\sigma(P,A) is the smallest sigma-algebra containing both events.
Each outcome must fall into one combination of plural or nonplural and accusative or nonaccusative.
These smallest nonempty intersections are the atoms.
P\cap A=\{\textit{them}\},\qquad P\cap A^c=\{\textit{they}\},
P^c\cap A=\{\textit{him}\},\qquad P^c\cap A^c=\{\textit{he}\}.
Every atom is a singleton in this four outcome space.
Every subset is thus a union of atoms, so
\sigma(P,A)=2^{\Omega_4}.
Generating from P alone cannot distinguish they from them.
More observations do not repair a distinction omitted by the representation.
Probability measures
| outcome | probability |
|---|---|
| he | .20 |
| him | .30 |
| they | .10 |
| them | .40 |
A probability measure maps measurable events to numbers.
\mathbb P:\mathcal F\rightarrow[0,1]
The triple
\langle\Omega,\mathcal F,\mathbb P\rangle
is a probability space.
For every B\in\mathcal F,
\mathbb P(B)\geq0.
\mathbb P(\Omega)=1.
The represented possibilities exhaust the probability mass.
For pairwise disjoint events B_1,B_2,\ldots,
\mathbb P\left(\bigcup_i B_i\right)=\sum_i\mathbb P(B_i).
For A=\{\textit{him},\textit{them}\},
\mathbb P(A)=.30+.40=.70.
Four exhaustive outcomes receive probabilities .40, .35, .20, and .10.
Can these values define a probability measure? Which requirement fails?
September 9
| outcome | number | case | probability |
|---|---|---|---|
| he | singular | nominative | .20 |
| him | singular | accusative | .30 |
| they | plural | nominative | .10 |
| them | plural | accusative | .40 |
Mutual exclusivity
One pronoun outcome cannot be both nominative and accusative in this representation.
The two events have no outcomes in common.
Events B and C are mutually exclusive when
B\cap C=\varnothing.
If B\cap C=\varnothing, then
\mathbb P(B\cup C)=\mathbb P(B)+\mathbb P(C).
For events that can occur together,
\mathbb P(B\cup C) =\mathbb P(B)+\mathbb P(C)-\mathbb P(B\cap C).
“Plural” and “accusative” are different labels, but them belongs to both events.
The represented outcomes determine whether the intersection is empty.
Plurality and case together
The joint probability of P and A is
\mathbb P(P\cap A).
The event P\cap A contains only them.
Thus
\mathbb P(P\cap A)=.40.
| accusative A | nominative A^c | total | |
|---|---|---|---|
| plural P | .40 | .10 | .50 |
| singular P^c | .30 | .20 | .50 |
| total | .70 | .30 | 1 |
Marginal probabilities appear in the table margins.
\mathbb P(P)=.40+.10=.50
Conditional probability
Among plural pronouns, how much of the probability belongs to accusative forms?
The phrase “among plural pronouns” makes P the reference event.
| plural outcome | probability in \Omega_4 | probability within P |
|---|---|---|
| they | .10 | .10/.50=.20 |
| them | .40 | .40/.50=.80 |
The conditional probability of A given P is
\mathbb P(A\mid P)=\frac{\mathbb P(A\cap P)}{\mathbb P(P)}.
\begin{aligned} \mathbb P(A\mid P) &=\frac{.40}{.50}\\ &=.80. \end{aligned}
\mathbb P(P\mid A)=\frac{.40}{.70}\approx.57.
The overlap stays the same. The reference event changes.
For \mathbb P(A\mid P),
Most passive clauses may have animate subjects.
It does not follow that most clauses with animate subjects are passive.
The product rule
Starting from
\mathbb P(A\mid P)=\frac{\mathbb P(A\cap P)}{\mathbb P(P)},
multiply both sides by \mathbb P(P).
The product rule gives
\mathbb P(A\cap P)=\mathbb P(A\mid P)\mathbb P(P).
\begin{aligned} \mathbb P(A\cap P) &=.80\times.50\\ &=.40. \end{aligned}
The same intersection can be written
\mathbb P(A\cap P)=\mathbb P(P\mid A)\mathbb P(A).
Both products recover .40.
The product rule contains a conditional probability. It allows the events to be associated.
Reversing a conditional probability
\mathbb P(B\mid C)\mathbb P(C) =\mathbb P(C\mid B)\mathbb P(B).
Both sides equal \mathbb P(B\cap C).
Bayes’ rule solves the equality for one conditional probability.
\mathbb P(B\mid C) =\frac{\mathbb P(C\mid B)\mathbb P(B)}{\mathbb P(C)}.
| feature F | no feature F^c | total | |
|---|---|---|---|
| dialect group D | 90 | 10 | 100 |
| other group D^c | 180 | 720 | 900 |
| total | 270 | 730 | 1,000 |
\mathbb P(F\mid D)=\frac{90}{100}=.90.
The denominator is the 100 members of the dialect group.
\mathbb P(D\mid F)=\frac{90}{270}=\frac13.
The denominator is the 270 tokens with the feature.
\mathbb P(D\mid F) =\frac{\mathbb P(F\mid D)\mathbb P(D)}{\mathbb P(F)} =\frac{.90\times.10}{.27} =\frac13.
Does conditioning change probability?
Events B and C are independent when
\mathbb P(B\cap C)=\mathbb P(B)\mathbb P(C).
\mathbb P(P)\mathbb P(A)=.50\times.70=.35
but
\mathbb P(P\cap A)=.40.
When \mathbb P(C)>0, independence implies
\mathbb P(B\mid C)=\mathbb P(B).
Learning that C occurred does not change the probability of B.
| accusative A | nominative A^c | total | |
|---|---|---|---|
| plural P | .35 | .15 | .50 |
| singular P^c | .35 | .15 | .50 |
| total | .70 | .30 | 1 |
The event labels are unchanged across the two tables.
The arrangement of probability mass determines whether the events are independent.
Mutually exclusive events with positive probability cannot be independent.
Observing one event rules out the other.
Explain why
From records to probability spaces
Those distinctions must be retained by the data representation and declared by its probability space.