Replicated linear functional relationship model is often used to describe
relationships between two circular variables where both variables have error terms and
replicate observations are available. We derive the estimate of the rotation parameter
of the model using the maximum likelihood method. The performance of the proposed
method is studied through simulation, and it is found that the biasness of the estimates
is small, thus implying the suitability of the method. Practical application of the
method is illustrated by using a real data set.