NEWS PRESENCE OF ASTEROID NAMES RELATES TO ANGLES TO SUN
NEWS PRESENCE OF ASTEROID NAMES RELATES TO ANGLES TO SUN
-- COLLECTION STARTED ON FEB 15, 2022-- EVALUATED DAILY UNTIL CIRCA FEB 15, 2027 -- HERE ARE DAILY AGGREGATE RESULTS
In[]:=
(*Today'sDateandTimeofCollection&ReportingUTC*)Style[DateString[TimeZoneConvert[DateObject[Date[]],"UTC"]],"Subtitle"]
Out[]=
Sun 4 Dec 2022 16:50:41
In[]:=
SetDirectory["/home/rko/Documents/Wolfram Desktop/Asteroids1211"];
(*useyourowndirectoryhere*)In[]:=
sundata=Import["MinorPlanetSunData.csv"][[All;;-2]];
In[]:=
sundata[[-1]]
Out[]=
In[]:=
Import["asteroidentities1211.m"]
Out[]=
Normalize and then average the daily article count across all individual minor planets
Normalize and then average the daily article count across all individual minor planets
Extract Dates from Imported Data
Extract Dates from Imported Data
In[]:=
dates=Table[DateObject[ToString[sundata[[2;;All,1]][[i]]]],{i,1,Length[sundata[[2;;All,1]]]}]//Quiet;
In[]:=
dates
Out[]=
Check whether dates are sampled regularly day by day and check number of dates
Check whether dates are sampled regularly day by day and check number of dates
In[]:=
RegularlySampledQ[dates]
Out[]=
True
In[]:=
Length[dates]
Out[]=
290
Extract Article Counts and Angles of Namesakes to Sun in Radians
Extract Article Counts and Angles of Namesakes to Sun in Radians
In[]:=
toanalyse=sundata[[2;;All,2;;All]];
In[]:=
toanalyse[[All,-1]]
Out[]=
{150.7,149.48,148.27,147.06,145.86,144.67,143.48,142.3,141.13,139.96,138.8,137.65,136.51,135.37,134.24,133.12,132.01,130.9,129.8,128.71,127.63,126.55,125.49,124.43,123.38,122.33,121.29,120.31,119.29,118.27,117.27,116.27,115.28,114.29,113.31,112.34,111.38,110.42,109.47,108.53,107.59,106.66,105.73,104.82,103.91,103.,102.1,101.21,100.33,99.45,98.57,97.7,96.84,95.99,95.13,94.29,93.45,92.617,91.79,90.96,90.15,89.33,88.52,87.72,86.92,86.13,85.34,84.55,83.77,83.,82.23,81.46,80.7,79.94,79.18,78.43,77.69,76.94,76.2,75.47,74.74,74.01,73.29,72.57,71.85,71.14,70.43,69.72,69.02,68.32,67.62,66.93,66.24,65.55,64.87,64.18,63.5,62.83,62.15,61.48,60.82,60.15,59.49,58.82,58.17,57.51,56.86,56.21,55.56,54.91,54.27,53.62,52.98,52.35,51.71,51.08,50.44,49.81,49.19,48.56,47.94,47.32,46.69,46.08,45.46,44.84,44.23,43.61,43.,42.39,41.78,41.18,40.57,39.97,39.36,38.76,38.16,37.56,36.97,36.37,35.77,35.18,34.59,34.,33.4,32.81,32.23,31.64,31.05,30.47,29.88,29.3,28.71,28.13,27.55,26.97,26.39,25.81,25.23,24.65,24.07,23.5,22.92,22.34,21.77,21.19,20.62,20.05,19.47,18.9,18.33,17.75,17.18,16.61,16.04,15.47,14.9,14.33,13.76,13.19,12.62,12.05,11.49,10.92,10.35,9.78,9.21,8.64,8.08,7.51,6.94,6.37,5.8,5.23,4.66,4.1,3.53,2.96,2.39,1.82,1.25,0.68,0.11,359.54,358.98,358.41,357.84,357.27,356.7,356.12,355.55,354.98,354.41,353.84,353.26,352.69,352.12,351.54,350.97,350.39,349.82,349.24,348.67,348.09,347.51,346.93,346.35,345.77,345.19,344.61,344.03,343.45,342.87,342.29,341.7,341.12,340.53,339.95,339.36,338.78,338.19,337.6,337.01,336.42,335.83,335.24,334.64,334.05,333.45,332.86,332.26,331.66,331.06,330.46,329.86,329.26,328.66,328.05,327.45,326.84,326.24,325.63,325.02,324.41,323.8,323.16,322.55,321.93,321.32,320.7,320.08,319.46,318.84,318.22,317.59,316.97,316.34,315.71,315.08,314.45,313.82,313.19,312.55,311.92,311.28,310.64,310.,309.35,308.71,308.06}
In[]:=
MatrixPlot[toanalyse]
Out[]=
In[]:=
articles=Table[toanalyse[[All,i]],{i,1,Length[toanalyse[[1]]],2}];
In[]:=
anglesrads=Table[UnitConvert[Quantity[toanalyse[[All,i]],"AngularDegrees"],"Radians"],{i,2,Length[toanalyse[[1]]],2}];
In[]:=
Length[articles[[1]]]
Out[]=
290
In[]:=
Length[anglesrads[[1]]]
Out[]=
290
For Each Minor Planet Studied, Construct Time Series for the Normalized Article Counts
For Each Minor Planet Studied, Construct Time Series for the Normalized Article Counts
In[]:=
articlesTS=Table[TimeSeries[articles[[i]],{dates}],{i,1,Length[articles]}];
In[]:=
normarticlesTS=Table[Normalize[articlesTS[[i]]],{i,1,Length[articlesTS]}];
Look at an Example, the Second One, for a Time Series of the Number of Articles Per Day
Look at an Example, the Second One, for a Time Series of the Number of Articles Per Day
In[]:=
normarticlesTS[[2]]
Out[]=
TimeSeries
In[]:=
ListPlot[%]
Out[]=
Similarly, Construct a Time Series for the Rounded Angular Degrees to the Sun for Each Minor Planet
Similarly, Construct a Time Series for the Rounded Angular Degrees to the Sun for Each Minor Planet
In[]:=
anglesTS=Table[TimeSeries[UnitConvert[anglesrads[[i]],"AngularDegrees"],{dates}],{i,1,Length[anglesrads]}];
In[]:=
justangles=Flatten[Table[Mod[Round[QuantityMagnitude[Values[anglesTS[[i]]]]],360],{i,1,Length[anglesTS]}]];
In[]:=
justvalues=Flatten[Table[Values[normarticlesTS[[i]]],{i,Length[normarticlesTS]}]];
In[]:=
Histogram[justvalues,Automatic,"ProbabilityDensity"]
Out[]=
Now, combine angles and article counts.
Now, combine angles and article counts.
In[]:=
just=Transpose[{justangles,justvalues}];
In[]:=
Length[just](*thisisthenumberofdatapointstodate*)
Out[]=
351190
Do some basic descriptive statistics. Some will be used for calculating sampling errors.
Do some basic descriptive statistics. Some will be used for calculating sampling errors.
Calculate the average normalized article count per rounded degree.
Calculate the average normalized article count per rounded degree.
A periodogram is useful for seeing periodicities in the data.
A periodogram is useful for seeing periodicities in the data.
Compare to a random simulation to see some idea of unusual behavior.
Compare to a random simulation to see some idea of unusual behavior.
There is a peak at 2 and 3. Let’s get more particular about pulling out one Fourier Transform peak at 30. It corresponds to a frequency of ~12, i.e. this may be same basis as for the astrology signs. Also, there is an even stronger signal at 93, corresponding to a frequency of ~4. This may be the basis for squares.
There is a peak at 2 and 3. Let’s get more particular about pulling out one Fourier Transform peak at 30. It corresponds to a frequency of ~12, i.e. this may be same basis as for the astrology signs. Also, there is an even stronger signal at 93, corresponding to a frequency of ~4. This may be the basis for squares.
A spectrogram visually confirms this frequency and period. A random data spectrogram follows for comparison.
A spectrogram visually confirms this frequency and period. A random data spectrogram follows for comparison.
We can see more particularly the spikes directly in the Fourier transform of the data.
We can see more particularly the spikes directly in the Fourier transform of the data.
Going back to the data values, let’s fit the results. They look like they fit well to a normal curve.
Going back to the data values, let’s fit the results. They look like they fit well to a normal curve.
The match to the normal curve is seen in a probability plot.
The match to the normal curve is seen in a probability plot.
Let’s look closer at the sampling errors of the average normalized article counts. At most they are less than 1 percent of their corresponding values.
Let’s look closer at the sampling errors of the average normalized article counts. At most they are less than 1 percent of their corresponding values.
A plot of the values with error bars
A plot of the values with error bars
Let’s look at the peaks and the base of the peaks.
Let’s look at the peaks and the base of the peaks.
The base as a moving average with a window of 7 degrees may show 12 peaks, consistent with the Fourier transform peak above at 30 and perhaps substantiating the idea of (unequal) houses.
The base as a moving average with a window of 7 degrees may show 12 peaks, consistent with the Fourier transform peak above at 30 and perhaps substantiating the idea of (unequal) houses.
Peaks of peaks may land in house 1, 4, 7, 10; i.e. the square houses a.k.a. kendras.
Peaks of peaks may land in house 1, 4, 7, 10; i.e. the square houses a.k.a. kendras.
Idealized Theoretical Result
Idealized Theoretical Result
Visualize Angles
Visualize Angles
Play with sliders.
Play with sliders.
More documentation will be added in the coming days ... .
More documentation will be added in the coming days ... .
Reasons why we can’t go down to one hundredths of a degree: I suspect that precision largely does not reflect the multi-day process of the news publication cycle, also sampling errors would be too high. Both are good reasons to avoid it.
Phi Golden Angles (Landscheidt’s grand crosses, etc.)
Phi Golden Angles (Landscheidt’s grand crosses, etc.)
Monte Carlo Simulation to get a sense of the p-value of such successful intersections
Monte Carlo Simulation to get a sense of the p-value of such successful intersections