Are there 'higher order moments' in physics? This may be a rather noob question but please let me clarify: I'm struggling to understand the use of the word 'moments' w.r.t., probability distributions. It seems after some research and poking around it seems to have been derived from physics when trying to solve/prove something about/related to binomial distribution and the method was called method of moments. I've asked the corresponding question here: https://stats.stackexchange.com/q/17595/4426
Now 'Pearson' (one of the very famous statisticians) comments:

We shall now proceed to find the first four moments of the system of 
  rectangles round GN. If the inertia of each rectangle might be
  considered  as concentrated along its mid vertical, we should have for
  the sth moment round NG, writing d = c(1 + nq).

Here are some of the details of the proof (as in the above post):

Now Pearson talks about calculating the 'rth' moment and uses a derivative function to do so:

Question: I'm not aware of such a function from my knowledge of elementary physics. What kind of moments are being calculated here? How do you calculate 'higher order moments'? Is there any such thing? 
Basically looking to clarify something in statistics but was historically alluded to physics and hence just want to get it ironed out :)
UPDATE: Intent of question: What I want to know is does the above derivation have anything at all to do with the concept of moments in physics and how is it related? Since the 'word' moment (and its intent) seem to be borrowed from physics when the author is making the derivation. I personally want to know if something like this does exist in the field of physics and how are these two derivations (and 'moments' related)
 A: I have to start by saying I don't know anything about the derivative method shown in this excerpt. I tried some calculations but it doesn't even seem to give the same result as the standard definition, so I'm guessing he is calculating something different from what we call "moments" in modern physics. Anyway, by way of explanation:

The word "moment" is used for several different purposes in physics, so it can be kind of a confusing term because you have to know what is meant by the context. But all the various meanings of moment stem from its definition in math.
In math, a moment is a way of characterizing some distribution. It could be a probability distribution, a mass distribution, a charge distribution, or anything similar; all you need is some function $f(x)$ which defines the density of the quantity (mass/charge/probability) in question. In other words, $\int_a^b f(x)\;\mathrm{d}x$ is the amount of "stuff" between $a$ and $b$.
The $n$th mathematical moment of a distribution with density function $f(x)$ around a point $c$ is computed by a very simple formula:
$$I^{(n)}(c) = \int (x - c)^n f(x)\ \mathrm{d}x$$
(as a slight abuse of notation, when $n = 0$ the moment is independent of $c$ so I'll write it as $I^{(0)}$) This generalizes to higher-dimensional spaces, but then the moment becomes an $n$-index tensor:
$$I_{i_1\cdots i_n}^{(n)}(\mathbf{c}) = \idotsint \prod_{j=1}^{d}(r_{i_j} - c_{i_j}) f(\mathbf{r})\ \mathrm{d}^d\mathbf{r}$$
In physical applications, the definitions used are a little different, but in general an $n$th moment involves the integral of some $n$th power of position multiplied by the distribution function $f(\mathbf{r})$. (The aforementioned differences show up in how you use the various components of $\mathbf{r}$ to compute that $n$th power.)
Many typical measures used to describe physical systems or mathematical distributions can be represented as moments. For example:

*

*If $f(x)$ is a 1D probability distribution:

*

*The normalization constant (which is 1) is $I^{(0)}$

*The mean value is $\langle x\rangle = I^{(1)}(0)$

*The variance is $I^{(2)}(\langle x\rangle)$



*If $f(\mathbf{r})$ is a mass distribution:

*

*The total mass is $I^{(0)}$

*The center of mass is $I^{(1)}(0)$ (from which comes the term "weighted average")

*The moment of inertia around any point $\mathbf{c}$ is a second moment



*If $f(\mathbf{r})$ is a charge distribution:

*

*The total charge, or monopole moment, is $I^{(0)}$

*The dipole moment is $I^{(1)}(0)$

*The quadrupole moment is a second moment

*and so on



For charge distributions, the quantities $I^{(n)}(0),\ n=0,1,2,\ldots$ (as modified with the required extra terms) are called the electric multipole moments $Q^{(n)}$. These quantities are of particular interest because you can expand the electric potential of an arbitrary charge distribution in terms involving successive moments:
$$\Phi(\mathbf{r}) = \sum_{n=0}^{\infty} \sum_{\{i,j\}}\frac{C_n Q_{i_1\cdots i_n}^{(n)}x_{i_1}\cdots x_{i_n}}{r^{2n+1}} \sim \sum_n \frac{C_n Q^{(n)}}{r^{n+1}}$$
In many situations, $r$ is relatively large so it's sufficient to use only the first nonzero term of this series in a calculation. In a sense, higher moments incorporate more detailed features of the charge distribution, which "blur out" and thus have little effect at large distances.

For the example you're looking at here, it sounds like Pearson is calculating the moments of area in the $x$ dimension around the origin - in other words, the density function $f(x)$ is the function that would trace along the tops of the rectangles.
$$f(x) = a\binom{n}{k}p^{n-k}q^k,\quad \tfrac{(2k + 1)c}{2} \le x < \tfrac{(2k + 3)c}{2}$$
(you could think of this as calculating the moments of mass of a cardboard cutout of the binomial distribution, assuming the cardboard is uniform density).
You can plug this into the integral definition of a moment, although the resulting expression is rather complicated, and as I said, it doesn't seem to give the same results as the derivative method Pearson is using. So I believe he's calculating something different.
A: The $n$th moment of a distribution is
$$ \sum_i w_i x_i^n $$
or
$$ \int \mathrm{d}x \, \rho(x) x^n$$
where $w$ is the "weight" of each discrete point (or $\rho$ is the continuous density) with "distance" $x$.
What physical values should be used for $x$ and $w$ ($\rho$) depends on what moment you are calculating. In the above calculation the weight appears to be a measured probability and the distance the position of the bin.
Note that the zeroth moment is just the sum of the weights (integral of the density), the first moment is the kind of calculation you see in the "moment of inertia" if we let $x$ be "the distance from the center".
A: Here's an example of higher moments being useful:
In heavy-ion collisions, we count the number of protons and antiprotons detected during each collision. We also record the net proton number, which is simply the difference of the former two numbers. We assemble this data into an approximation of the probability distribution of seeing a given net proton number.
We calculate the first- through eighth-order central moments of this distribution, and use these to calculate the first- through eighth-order cumulants. The ratios of these cumulants are an indicator of the presence of quark-gluon plasma, in that they would be significantly different if nuclei were melting. For a more detailed explanation, see https://arxiv.org/abs/1607.06602.
A: We can say that higher order moments represent the shape of the distribution: the mean the second moment, represent how much is the scattering, where this can be seen by a point scatter. the third represents the skewness /asymmetry, the fourth how much is the flatness of the curve. Further moments continue this way but it is as far as I know inexpressible in language or you can consider the first moments of the square of the variable, the cube...
In general, these higher orders represent the details of the variable: Going to higher order, you represent the data in more and more detail view. This observation has practical implications: if you have a sample of finite amount of data, you can calculate the moments up to an order: computing the higher order is useless and error because you don't have enough data. Or  better stated, the variance of the moment considered in large.
I f you are  interested, my book(PhD thesis) "Construction of Random Signals from their Higher Order moments", Ismail Chamseddine, 1997.
