Divertor tiles from the COMPASS tokamak have been examined with ion beam analysis for the determination of deuterium retention and co-deposits. A novel approach was used to draw 2D deuterium retention maps at different depths and correlated to the concentrations of other elements: boron, carbon, oxygen and metals. This approach, which employs all three spatial dimensions, revealed that the deuterium retention pattern is not toroidally symmetric at different depths while the overall deuterium pattern is toroidally symmetric. Analysis of visible-light camera records disclosed that parts of the divertor were not wetted by plasma. The radial profile of deuterium in both wetted and non-wetted ( "shadowed ") divertor regions was compared to strike point positions, angle of incidence, electron temperature T-e, parallel heat flux Q(par) and ion saturation current Isat in order to find the influence of these parameters on deuterium retention. There are certain series of consecutive shots for which correlations between deuterium retention and separate plasma parameters (T-e, I-sat, Q(par)) are high. Combining the three parameters in a linear model or a power law model enabled reconstruction of the radial deuterium retention pattern. After fitting the model to one third of the data points, the deuterium retention for shadowed tiles was reproduced at the other two thirds of the data points within 20% deviation. This study shows the need for detailed tile analysis in three dimensions, the great benefits of broad statistical analysis in solving challenges in fusion, and motivates the investigation of COMPASS tiles after dismantling.