Epidermal nerve fibre (ENF) density and summed length of ENFs per epidermal surface area are reduced, and ENFs may appear more clustered within the epidermis in subjects suffering from diabetic neuropathy compared to healthy subjects. Therefore, it is important to understand the spatial behaviour of ENFs in healthy and neuropathy subjects. By using confocal microscopy data , we study the spatial structure of epidermal nerves by regarding the nerve tree locations as realizations of marked point processes . The termination points of the fibres of a nerve tree are used to define a reactive territory which is taken as a mark for the nerve tree location. We study the differences in the spatial pattern of ENFs between healthy subjects and subjects suffering from mild diabetic neuropathy by using Ripley's K function and the mark correlation function. In addition, we propose a marked sequential point process model for the nerve tree locations. Data are replicated point patterns, where we have several patterns from each subject and from each group.

Diabetic neuropathy is a disorder characterized by impaired nerve function and reduction of the number of epidermal nerve fibers per epidermal surface. Additionally, as neuropathy related nerve fiber loss and regrowth progresses over time, the two‐dimensional spatial arrangement of the nerves becomes more clustered. These observations suggest that with development of neuropathy, the spatial pattern of diminished skin innervation is defined by a thinning process which remains incompletely characterized. We regard samples obtained from healthy controls and subjects suffering from diabetic neuropathy as realisations of planar point processes consisting of nerve entry points and nerve endings, and propose point process models based on spatial thinning to describe the change as neuropathy advances. Initially, the hypothesis that the nerve removal occurs completely at random is tested using independent random thinning of healthy patterns. Then, a dependent parametric thinning model that favors the removal of isolated nerve trees is proposed. Approximate Bayesian computation is used to infer the distribution of the model parameters, and the goodness‐of‐fit of the models is evaluated using both non‐spatial and spatial summary statistics. Our findings suggest that the nerve mortality process changes as neuropathy advances.