Boltzmann used a ρ ρ ln ⁡ ⁡ ρ ρ formula as early as 1866. He interpreted ρ as a density in phase space—without mentioning probability—but since this satisfies the axiomatic definition of a probability measure we can retrospectively interpret it as a probability anyway. Gibbs gave an explicitly probabilistic interpretation in 1878.