planets: Planet orbits

In sphunif: Uniformity Tests on the Circle, Sphere, and Hypersphere

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.

Source

Table 2a in https://ssd.jpl.nasa.gov/planets/approx_pos.html

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")

sphunif documentation built on May 29, 2024, 4:19 a.m.