In this paper, parameter estimation for R-dimensional (R -D) sinusoids with R >; 2 in additive white Gaussian noise is addressed. With the use of tensor algebra and principal-singular-vector utilization for modal analysis, the sinusoidal parameters at one dimension are first accurately estimated according to an iterative procedure which utilizes the linear prediction property and weighted least squares. The damping factors and frequencies in the remaining dimensions are then solved such that pairing of the R-D parameters is automatically achieved. Algorithm modification for a single R -D tone is made and it is proved that the frequency estimates are asymptotically unbiased while their variances approach Cramér-Rao lower bound at sufficiently high signal-to-noise ratio conditions. Computer simulations are also included to compare the proposed approach with conventional R -D harmonic retrieval schemes in terms of mean square error performance and computational complexity.