This Vulcan web site presents a way to compute Vulcan’s orbital parameters by using an Blavetsky's theosophy and an IRAS (Infrared Astronomical Survey) object cataloged as 1732+239. This results in a theoretical value for Vulcan period to be 4969.0 +/- 5.75 (one sigma) years.
However, there are three other methods that lead directly to four independent computations of Vulcan’s orbital period. The first was preformed by Professor Forbes in 1880 and used the aphelion’s of long period comets. The second was the reconstruction of the path of a giant comet, CR105 over a period of that estimated to be two Vulcan’s period and measuring the time between three successive positions on the comet’s orbit. This comet has been shown to be one of the ones in a 3:2 resonance with Vulcan’s orbital period and the total of three consecutive periods must exactly match that of two of Vulcan’s periods. The third way involves measuring the time interval between weather changes believed to be cause by comet impacts (or possibly by their simply exploding in the atmosphere). The specific clusters of comets within two comet swarms were used in this calculation, Swarm A: Cl-1 and Swarm B:Cl-2. Where only the average comet period and its one sigma error was obtained, these values were multiplied by 3/2 to translated them back into the values reflecting the period of Vulcan. These are listed below.
The combined one sigma error is obtained by adding the inverse sum of the squares and taking the square root. Thus:
1/(one combined sigma)2 = 1/2002 + 1/29.252 + 1/102 + 1/162 = 0.01510007
Thus, one combined sigma = 8.14 years
The Combined Period is obtained by adding the sums of the separate periods, each divided by its separate one sigma value, Then dividing this quantity by 1/(one combined sigma)2. Thus:
(5000/2002 + 4979/29.252 + 4969.5/102 + 4959.5/162)/0.01510007 = 75.012607/0.01510007 = 4967.7
Thus, the estimated Combined Period = 4967.7 years +/- 8.14 (one sigma) years
A statistical estimate can be made to determine how well the combined period data compares to the theoretical value for Vulcan period to the Combined Period value. This method measures how far they are apart from a Gaussian perspective and that can be translated into a probability of overlap.
Gaussian Separation = [Absolute Value of the (Theoretical Period - Combined Period)]/[( one sigma)2 +(one combined sigma)2]1/2
Gaussian Separation = Absolute Value(4969 - 4967.7)/[(5.752 + 8.142)]1/2
Gaussian Separation = 1.3/9.964 = 0.13
This translates into a 90% chance that the combined measured values (4967.7 years) is related to the theoretical value (4969.0 years).
A similar calculation can be made for the case where Vulcan is not constrained to be at the IRAS point. Here the theoretical estimate is 4972 Years +/- 27.35 (one sigma) years. In this case, the Gaussian Separation = 0.105 which corresponds to a 91% chance that the theoretical and combined values are associated. Thus this statistical analysis cannot resolve whether or not Vulcan is exactly at the specified IRAS point, or is slightly removed from it (as if a typo in recording the data was made).