Talk About Network



Register and Login
Nick
Password
Register create new account Sign up is FREE and you can post replies, new topics, bookmark posts and more!
Recover lost password


Programming > Fractals > Re: Analysis of...
Latest [ Topics | Posts ] Archive Post A New Topic Post a Reply
<< Topic < Post Post 5 of 7 Topic 186 of 206
Post > Topic >>

Re: Analysis of Sunspot number cycles in Mathematica

by Roger Bagula <rlbagula@[EMAIL PROTECTED] > Dec 7, 2007 at 05:53 AM

I compared the:
1600 to 1980 temperatures as a data file here:
http://www.ncdc.noaa.gov/paleo/ei/ei_data/latbands.dat

with Group Sunspot Numbers (Doug Hoyt re-evaluation) 1610-1995
ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUNSPOT_NUMBERS/GROUP_SUNSPOT_NUMBERS/yearrg.dat

What I get is an unrelated as far as I can analyze noise blot.

As far as I'm concerned that pretty much puts to rest the oil company 
argument that the warming is due to solar activity.
There is no real detectable correlation between the two data sets;
even their correlation dimensions are very different.
The only thing one can say is that they both appear to be increasing.
That there is a strange line effect as related to the golden mean...
The scatter is really awful.
You can get a picture of the scatter plot from me
or see it at:
http://profile.imeem.com/GUmj0c/photo/W86TQ5sS0k/
Mathematica:
Directory[]
Apply[Plus, {0.06`, -0.26`, -0.01`}]/3
raw = Flatten[ReadList["latbands2.dat", Number, RecordLists -> True]]
Dimensions[raw]
l = Length[raw]
b = Table[((raw[[n]] + raw[[n + 1]] + raw[[
      n + 2]])/3 - (-0.11900269541778977))/0.2866666666666667, {n, 2 + 
4*10, \
Length[raw], 4}]
Length[b]
Max[b]
Apply[Plus, b]/Length[b]
raw2 = Flatten[ReadList["yearrg.dat", Number, RecordLists -> True]]
Dimensions[raw2]
Length[raw2]
b2 = Table[(raw2[[n]] - 34.19611398963731)/175.1, {n, 2, Length[raw2], 2}]
Length[b2]
Max[b2]
Apply[Plus, b2]/Length[b2]
d = Sort[Table[{Log[Abs[b[[n]]]], Log[Abs[b2[[n]]]]/(1 +
Sqrt[5])/2 + 1/4}, {n, 1, Length[b]}]]
ga = ListPlot[d]
ListPlot[d, PlotJoined -> True, PlotRange -> All]
StandardDeviation[b]
StandardDeviation[b2]
Table[Length[Delete[
        Union[Table[If[Abs[b[[n]]] > m*StandardDeviation[b], {n + 1609, \
b[[n]]}, {}], {n, 1, Length[b]}]], 1]], {m, 1, 5}]
Dimensions[%]
Table[Length[Delete[Union[Table[If[Abs[b2[[n]]] > m*StandardDeviation[b2],
{
        n + 1609, b2[[n]]}, {}], {n, 1, Length[b2]}]], 1]], {m, 1, 5}]
Dimensions[%]
<< Statistics`LinearRegression`
Regress[d, {1, x}, x]
f[x_] = Fit[d, {1, x}, x]
gb = Plot[f[x], {x, -6, 1}]
Show[{ga, gb}]
Respectfully, Roger L. Bagula
 11759Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814 
:http://www.geocities.com/rlbagulatftn/Index.html
alternative email: rlbagula@[EMAIL PROTECTED]





 7 Posts in Topic:
Analysis of Sunspot number cycles in Mathematica
Roger Bagula <rlbagula  2007-12-03 05:18:16 
Re: Analysis of Sunspot number cycles in Mathematica
"Simon" <new  2007-12-04 21:11:09 
Re: Analysis of Sunspot number cycles in Mathematica
Roger Bagula <rlbagula  2007-12-04 03:43:00 
Re: Analysis of Sunspot number cycles in Mathematica
"Simon" <new  2007-12-06 21:35:48 
Re: Analysis of Sunspot number cycles in Mathematica
Roger Bagula <rlbagula  2007-12-07 05:53:05 
Re: Analysis of Sunspot number cycles in Mathematica
Roger Bagula <rlbagula  2007-12-09 12:55:33 
IPCC WG1 AR4 Report->Palaeoclimate
Roger Bagula <rlbagula  2007-12-09 03:45:40 

Post A Reply:
  Go here to Signup

AddThis Feed Button


About - Advertising - Contact - Frequently Asked Questions - Privacy Policy - Terms of Use - Signup

Contact
tan12V112 Sat May 17 1:06:27 CDT 2008.