By the application of Rodrigues parameters, crystal orientations are represented as points on the unit semi-hypersphere $4+ c R 4 or equivalently on the projective hyperplane H3c R 4. For the statistical analysis of orientation data, probability models on $4+-H 3 are required. Among different hyperspherical analogs of the normal distribution in Euclidean space, corresponding to various of its characterizations, the Bingham model distribution is" characterized as the hyperspherical analog for statistical purposes. Maximum-likelihood estimates of its parameter matrices and single-sample significance tests of uniformity and rotational symmetry are presented. In terms of probability theory and stochastic processes, the hyperspherical Browniandiffusion distribution is the analog of the normal distribution in Euclidean space. An orientation distribution derived by Savelova lind. Lab. (USSR), (1984), 50, 468-474] and erroneously related to a central-limit-type theorem for rotations by Matthies, Muller & Vinel [Textures Mierostruct. (1988), 10, 77-96], with the Brownian-diffusion distribution on $4+ -H 3, recalls the fact that no simple analog of the central limit theorem in Euclidean space exists for arbitrary spaces and manifolds, e.g. for S p, SP+, SO(p), p > 2. Therefore, interpretations of preferred orientations in terms of mechanisms of texture development related to an imaginary simple analog of the central limit theorem would generally be misleading. Individual crystal-orientation measurements may initially require different nonparametric methods of analysis. As a counterpart of pole-to-orientation density inversion of diffraction data, hyperspherical kernel orientation density estimation is suggested for individual orientation measurements. With the application of non-negative kernels, the estimated orientation density is always non-negative. The critical smoothing parameter of this method is determined on information-theoretical grounds. For patterns of preferred orientation too complex to be sufficiently approximated by the parametric Bingham model distribution, a clustering of individual orientation 0021-8898/93/010112-10506.00 measurements into disjunct classes is suggested that relates similarity of orientations to multimodality of the estimated orientation density; each class of orientations may then be individually analyzed by parametric methods.