SUMMARYWe study bootstrap con®dence intervals for correlation functions in nearest-neighbour Markov point processes, where the neighbours are characterized by an interaction of bounded radius r. In forestry statistics, the points are tree locations belonging to a region (forest) A, and the marks are qualitative or quantitative tree variables, such as tree species, the stem diameter, crown length or tree height. Estimating and analysing correlation functions between locations and marks, cross-correlations between different species and their marks, is typically a key step in statistical interpretation of mapped data sets from a forest stand. In order to de®ne the original sample, we propose coding schemes, which are ®xed and divide the observed region A of the point process into regular, conditionally independent subregions fB g located at Euclidean distance d ! 2r. Bootstrap con®dence intervals are then obtained directly, by considering kernel density estimates from all subregions fB g as conditionally independent replicates.