Article 4VKZY Set theory at the command line

Set theory at the command line

by
John
from John D. Cook on (#4VKZY)

Often you have two lists, and you want to know what items belong to both lists, or which things belong to one list but not the other. The command line utility comm was written for this task.

Given two files, A and B, comm returns its output in three columns:

  1. A - B
  2. B - A
  3. A a(C) B.

Here the minus sign means set difference, i.e. A - B is the set of things in A but not in B.

comm_venn.png
The numbering above will be correspond to command line options for comm. More on that shortly.

Difference and intersection

Here's an example. Suppose we have a file states.txt containing a list of US states

 Alabama California Georgia Idaho Virginia

and a file names.txt containing female names.

 Charlotte Della Frances Georgia Margaret Marilyn Virginia

Then the command

 comm states.txt names.txt

returns the following.

 Alabama California Charlotte Della Frances Georgia Idaho Margaret Marilyn Virginia

The first column is states which are not female names. The second column is female names which are not states. The third column is states which are also female names.

Filtering output

The output of comm is tab-separated. So you could pull out one of the columns by piping the output through cut, but you probably don't want to do that for a couple reasons. First, you might get unwanted spaces. For example,

 comm states.txt names.txt | cut -f 1

returns

AlabamaCaliforniaIdaho

Worse, if you ask cut to return the second column, you won't get what you expect.

AlabamaCaliforniaCharlotteDellaFrancesIdahoMargaretMarilyn

This is because although comm uses tabs, it doesn't produce a typical tab-separated file.

The way to filter comm output is to tell comm which columns to suppress.

If you only want the first column, states that are not names, use the option -23, i.e. to select column 1, tell comm not to print columns 2 or 3. So

 comm -23 states.txt names.txt

returns

 Alabama California Idaho

with no blank lines. Similarly, if you just want column 2, use comm -13. And if you want the intersection, column 3, use comm -12.

Sorting

The comm utility assumes its input files are sorted. If they're not, it will warn you.

If your files are not already sorted, you could sort them and send them to comm in a one-liner [1].

 comm <(sort states.txt) <(sort names.txt)
Multisets

If your input files are truly sets, i.e. no item appears twice, comm will act as you expect. If you actually have a multiset, i.e. items may appear more than once, the output may surprise you at first, though on further reflection you may agree it does what you'd hope.

Suppose you have a file places.txt of places

 Cincinnati Cleveland Washington Washington

and a file of US Presidents.

 Adams Adams Cleveland Cleveland Jefferson Washington

Then

 comm places.txt presidents.txt

produces the following output.

 Adams Adams Cincinnati Cleveland Cleveland Jefferson Washington Washington

Cleveland was on the list of presidents twice. (Grover Cleveland was the 22nd and the 24th president, the only US president to serve non-consecutive terms.) The output of comm lists Cleveland once in the intersection of places and presidents (say 22nd president and a city in Ohio), but also lists him once as a president not corresponding to a place. That is because there are two Clevelands on our list of presidents, there is only one in our list of places.

Suppose we had included five places named Cleveland (There are cities named Cleveland in several states). Then comm would list Cleveland 3 times as a place not a president, and 2 times as a place and a president.

In general, comm utility follows the mathematical conventions for multisets. Suppose an item x appears m times in multiset A, and n times in a multiset B. Then x appears max(m-n, 0) times in A-B, max(n-m, 0) times in B-A, and min(m, n) times in A a(C) B.

Union

To form the union of two files as multisets of lines, just combine them into one file, with duplicates. You can join file1 and file2 with cat (short for "concatenate").

 cat file1 file2 > multiset_union_file

To find the union of two files as sets, first find the union as multisets, then remove duplicates.

 cat file1 file2 | sort | uniq > set_union_file

Note that uniq, like comm, assumes files are sorted.

Related

[1] The commutility is native to Unix-like systems, but has been ported to Windows. The examples in this post will work on Windows with ports of the necessary utilities, except for where we sort a file before sending it on to comm with <(sort file). That's not a feature of sort or comm but of the shell, bash in my case. The Windows command line doesn't support this syntax. (But bash ported to Windows would.)

e-LiAPd26RA
External Content
Source RSS or Atom Feed
Feed Location http://feeds.feedburner.com/TheEndeavour?format=xml
Feed Title John D. Cook
Feed Link https://www.johndcook.com/blog
Reply 0 comments