sphunif (version 1.4.3)

planets: Planet orbits

Description

Planet orbits data from the JPL Keplerian Elements for Approximate Positions of the Major Planets. The normal vector of a planet orbit represents is a vector on \(S^2\).

Usage

planets

Format

A data frame with 9 rows and 3 variables:

planet

names of the planets and Pluto.

i

inclination; the orbit's plane angle with respect to the ecliptic plane, in radians in \([0, \pi]\).

om

longitude of the ascending node; the counterclockwise angle from the vector pointing to the First Point of Aries and that pointing to the ascending node (the intersection between orbit and ecliptic plane), in radians in \([0, 2\pi)\). (Both vectors are heliocentric and within the ecliptic plane.)

Details

The normal vector to the ecliptic plane of the planet with inclination \(i\) and longitude of the ascending node \(\omega\) is $$(\sin(i) \sin(\omega), -\sin(i) \cos(\omega), \cos(i))'.$$

The script performing the data preprocessing is available at planets.R. The data was retrieved on 2020-05-16.

Examples

# Load data
data("planets")

# Add normal vectors
planets$normal <- cbind(sin(planets$i) * sin(planets$om),
                       -sin(planets$i) * cos(planets$om),
                       cos(planets$i))

# Tests to be performed
type_tests <- c("PCvM", "PAD", "PRt")

# Tests with Pluto
unif_test(data = planets$normal, type = type_tests, p_value = "MC")

# Tests without Pluto
unif_test(data = planets$normal[-9, ], type = type_tests, p_value = "MC")