asprine 发表于 2014-5-18 13:58 
acos(x) 的值域是[0,pi]。
呵呵,可能是我没有说清楚。我有一组x值和y值,它们分属于不同的象限,而这些数值是对应于0~2*pi角度,所以需要根据它们的数值来推算对应的theta。
我找到了解决方法了:
theta1 <- acos(x/sqrt(x^2+y^2))[y>=0]
theta2 <- 2*pi - acos(x/sqrt(x^2+y^2))[y<0]
theta <- c(theta1, theta2)
theta <- sort(theta)
对应的x和y的值也应该根据重新调整后的theta来做相应的调整。