If the vector is called v, then the built-in function that returns the
angle is R▶Pθ(v[1,1],v[1,2]). To avoid using the matrix
representation of the vector you could also use R▶Pθ(mat▶list(v)
[1],mat▶list(v)[2]). Either method is cumbersome and this has always
seemed like a good place for some improvements to be made.
Wayne