3D equilibria with an anisotropic pressure component in the large helical device are analyzed with respect to their magnetic axis locations. The anisotropic extension of the 3D equilibrium solver variational moments equilibrium code, anisotropic neumann inverse moments equilibrium code, is used to compute fixed-boundary plasma equilibria based on a bi-Maxwellian distribution function describing the anisotropic particles. Different heating scenarios are assessed by means of parallel and perpendicular pressure anisotropies with different radial anisotropic pressure profiles imposed. A theoretical predicted scaling of the magnetic axis location with the auxiliary parameter βeq as predicted for classical stellarators and heliotrons by Hitchon [Nucl. Fusion 23, 383 (1983)] is found to be applicable to the large helical device in the case of a flat hot-particle profile for parallel or weak perpendicular dominated anisotropies with β ⊥ / β ∥ ≤ 2. For strong perpendicular or non-flat hot-particle profiles, a deviation from the predicted scaling of the magnetic axis location is found. Whereas center-peaked profiles show a stronger shift of the magnetic axis, edge-peaked profiles show no significant change of its radial location. High critical magnetic fields are identified as a necessary condition for strong perpendicular anisotropies. The observed deviations are ascribed to the magnetic field structure and negative pressure gradients. The invalidity of the theoretical predictions in the case of certain configurations is found to be caused by higher-order terms in the pressure components, which are not accounted for by the ordering on which the theory is based.