Wednesday, 17 April 2013


17 April 2013

Local Regression v the Gaussian Mask



Introduction



On my previous post local-regression-loess-and-sun-spot I went to great lengths to explain how local regression works in practice. To my great delight I discovered that the UK Met office is using their own version of local regression. 

But, instead of calling it local regression they call it a Gaussian mask.

See if you can spot a difference or a similarity?

Met office gaussian mask

So, instead of referring to this process as a local regression, and acknowledging the work done by the Americans, a new term to describe the technique was created



Thursday, 28 March 2013

28 March 2013

Plotting of Sun Spot Cycles: Where are we now?

In my first blog relating to the sun spot cycle count I discussed the application of local regression (LOESS) on the sunspot cycle recordings.  I showed the underlying logic of the LOESS regression and posted the first results using this method, so now I can just continue to post results of the exercises I've done. 

Sun Spot Cycle 23


Raw sun spot count daily data



MA 390 day moving average



















LOESS 390 v MA 390


















Sun Spot Cycle 22



Raw sun spot count daily data


















MA 390 day moving average

















LOESS 390 v MA 390

















Sun Spot Cycle 21



Raw sun spot count daily data




















MA 390 day moving average



















LOESS 390 v MA 390



















Sun Spot Cycle 19


Raw sun spot count daily data

















MA 390 day moving average


















LOESS 390 v MA 390


















Sun Spot Cycle 18



Raw sun spot count daily data



MA 390 day moving average


LOESS 390 v MA 390


Sun Spot Cycle 17


Raw sun spot count daily data










MA 390 day moving average

LOESS 390 v MA 390

Sun Spot Cycle 16



Raw sun spot count daily data
















MA 390 day moving average
















LOESS 390 v MA 390
















Sun Spot Cycle 15


Raw sun spot count daily data





MA 390 day moving average

LOESS 390 v MA 390

Sun Spot Cycle 14



Raw sun spot count daily data





MA 390 day moving average


LOESS 390 v MA 390

Sun Spot Cycle 13



Raw sun spot count daily data





MA 390 day moving average







LOESS 390 v MA 390














Sun Spot Cycle 12


Raw sun spot count daily data

















MA 390 day moving average
















LOESS 390 v MA 390













Sun Spot Cycle 11


Raw sun spot count daily data















MA 390 day moving average
















LOESS 390 v MA 390













Sun Spot Cycle 10


Raw sun spot count daily data
















MA 390 day moving average
















LOESS 390 v MA 390














Sun Spot Cycle 9



Raw sun spot count daily data

















MA 390 day moving average

















LOESS 390 v MA 390















Sun Spot Cycle 8


Raw sun spot count daily data


















MA 390 day moving average

















LOESS 390 v MA 390















Sun Spot Cycle 7


Raw sun spot count daily data
















MA 390 day moving average
















LOESS 390 v MA 390













Sun Spot Cycle 6


Raw sun spot count daily data
















MA 390 day moving average
















LOESS 390 v MA 390
















Updated 29 March 2013

FIN