Cube Antics

I’m Moving!

2010-02-09T21:18:24Z

After much thought, I’ve decided to move this blog to creatingcode.com.

The new feed location is here – Creating Code

Here are the following changes I am going to make when switching over to Creating Code

Create more programming centric articles.
Focus on quality over quantity
Keep the same light-hearted attitude towards the delivery of the content, while taking the production of the content more seriously.

There are a couple of articles I’ve written for Cube Antics that I think have some real promise, but are lacking in either depth or overall writing quality. I plan on taking a few of them and revamping them for Creating Code.

Next time you visit cubeantics.com you should get redirected to creatingcode.com

It’s been fun, thanks for reading and I’ll see you on the flip-side.

Random Posts

Starting a WordPress Theme From Scratch

2010-02-05T20:59:45Z

So much of web design today is duplicating effort. This is especially true for creating WordPress themes because of the extra blog code that needs to be hooked into each theme.

I started a new project on GitHub called MinWP that will contain all of the source files necessary to get started with a new WordPress theme. The goal is to create a fully functional 2-column theme with absolutely no additional design elements added to it.

Now, to create a new WordPress theme, you can simply download the files in MinWP, add CSS and blog specific code, and upload your new theme.

This is a project that will continue to evolve over time, so if you want to contribute, feel free!

Turning Off Post Revisions in Wordpress

Using .NET to Monitor a Directory for Changes

2010-02-02T02:35:57Z

.NET’s FileSystemWatcher class can soothe the paranoid control freak in all of us by monitoring a specified folder for different types of file system changes.

The Code

using System;
using System.Text;
using System.IO;

namespace FolderWatcher {
  class Watcher {
    static void Main(string[] args) {
      FileSystemWatcher watcher = new FileSystemWatcher(@"D:\test");
      watcher.IncludeSubdirectories = true;
      watcher.Filter = "";
      watcher.Renamed += new RenamedEventHandler(renamed);
      watcher.Deleted += new FileSystemEventHandler(changed);
      watcher.Changed += new FileSystemEventHandler(changed);
      watcher.Created += new FileSystemEventHandler(changed);
      watcher.EnableRaisingEvents = true;

      Console.ReadKey();
   }

    private static void renamed(object sender, RenamedEventArgs e) {
      Console.WriteLine(DateTime.Now + ": " +
          e.ChangeType + " " + e.FullPath);
    }

    private static void changed(object sender, FileSystemEventArgs e) {
      Console.WriteLine(DateTime.Now + ": " +
          e.ChangeType + " " + e.FullPath);
    }
  }
}

Explanation

FileSystemWatcher watcher = new FileSystemWatcher(@”D:\test”);

The core of is application is .NET’s FileSystemWatcher class which monitors directories for file system changes.

watcher.IncludeSubdirectories = true;
watcher.Filter = “”;

Tells the FileSystemWatcher object to include changes in subdirectories and raise events on every type of file.

An empty string for the filter will match every file and folder in the watched directory. File types and names can also be added to the filter to create a more fine grained watch.

watcher.Renamed += new RenamedEventHandler(renamed);
watcher.Deleted += new FileSystemEventHandler(changed);
watcher.Changed += new FileSystemEventHandler(changed);
watcher.Created += new FileSystemEventHandler(changed);

Sets the types of file system events the FileSystemWatcher object will respond to. The four possible (and self-explanatory) events are: Renamed, Deleted, Changed, and Created. Each EventHandler takes the method to be executed on that particular event as a paramater in its constructor.

watcher.EnableRaisingEvents = true;

The final piece of the puzzle is to finally tell the FileSystemWatcher to start raising events on the specified directory. Once this statement is executed, we will have a live directory watcher.

Conclusion

I wrote this code as a proof of concept for a new project I am working on. My goal is to create a simple/straightforward directory monitor that anyone can get and start using with very minimal effort.

301 Redirect with ASP .NET

2010-01-25T18:47:37Z

When redirecting URL’s it is very important not to lose all of your search engine mojo. The best way to do this, is by using a 301 redirect which notifies the requester of a permanent URL move.

<script runat="server">
private void Page_Load(object sender, System.EventArgs e) {
  Response.Status = "301 Moved Permanently";
  Response.AddHeader("Location","http://www.newdomain.com");
}
</script>

Should you have to redirect a URL to a new location, using a status 301 redirect is the safest way to do it.

Creating a Fluid Two-Column Layout With a Single Fixed Width Column

2010-01-20T22:03:54Z

Ever wonder if you can create a fluid layout with a fixed witdh column? Well, you can!

layout.html

<html>
  <head>
    <title>Fixed Left Column</title>
    <link rel="stylesheet" href="test.css" type="text/css" />
  </head>
  <body>
    <div>Left div</div>
    <div>Right div</div>
  </body>
</html>

test.css

.left {
	width: 300px;
	float: left;
	background-color: #CCCCCC;
}
.right {
	margin-left: 300px;
}

There are very few times in life that you can have things both ways, this is one of them.

Using #if and [Conditional()] to Keep Debug Code Under Control

2010-01-20T19:36:59Z

The #if directive provides a perfect way to ensure your released executable doesn’t execute any code that was only intended for debug mode.

#if (<symbol>) will return true whenever the symbol being checked has been defined by default, or manually using #define.

#if DEBUG
	//Perform debug only code
#else
	//Release-level code
#endif

Remember, this is only one application of #if, you can use it any time you need to check to see if a symbol has been defined.

You can also use the [Conditional("<symbol>")] attribute on any method that has a return type of void. This will cause the compiler to remove all calls to the given method if the specified symbol is not defined.

/* using System.Diagnostics; */
[Conditional("DEBUG")]
public void printDebug() {
	Console.WriteLine("Debug Information");
}

Now printDebug() will not be executed unless the application is in debug mode.

Either method will save you time and effort by allowing you to leave debug code in your application, as well as provide you a more standardized and safe way to include debug code.

Sending Email in ASP .NET

2010-01-17T17:15:26Z

I can’t think of too many websites that don’t send automated email. Unfortunately, for ASP .NET programmers, there are several incorrect code snippets online which try to explain how to send email.

Here’s the code I use on the contact form for my personal website which does the trick quite nicely.

/* using System.Net.Mail; */
MailMessage newMail = new MailMessage();
newMail.To.Add(toAddress);
newMail.Subject = subject;
newMail.Body = body;
newMail.From = new MailAddress(fromAddress, name);
newMail.IsBodyHtml = true;

SmtpClient SmtpSender = new SmtpClient();
SmtpSender.Port = 25; //or, whatever port your SMTP server operates on
SmtpSender.Host = "mail.yourdomain.com";
SmtpSender.Send(newMail);

Here, I’m using MailMessage to build the actual email message to be sent and SmtpClient to do the actual sending.

Add this code to your Web.config file before the final </configuration> tag:

<system.net>
  <mailSettings>
    <smtp>
      <network defaultCredentials="false"
           host="mail.yourdomain.com" port="25"
           userName="name-to-send-email@yourdomain.com"
           password="*****"/>
     </smtp>
   </mailSettings>
</system.net>

Now, you should be able to send automated email messages to your heart’s content. Enjoy.

Worrying About Small Details Can Kill Your Project

2010-01-14T17:18:17Z

It’s amazing how seemingly insignificant things can completely derail a software project. I was working on my first project of the year -a simple online Regex tool- and abruptly stopped work when I realized that the name I had chosen wasn’t good enough. I then proceeded to obsess off and on for days until I came around full-circle and decided to stick with what I originally came up with.

Not only was this a horrible waste of time, it completely killed any momentum that I had built up in the first place. Now, here I am, at the exact same point in the project I was a week ago with nothing to show for it.

So, I suppose that’s my vice. For some reason, I get this weird feeling in my brain that if some small detail of the project isn’t perfect, then the it must doomed to failure. I went through the same thing with this blog, and it actually caused me to delay getting started by a few months after I finally decided to take the plunge (not to mention, the wasted $50 or so from domain name purchases that I’ll never see again.)

I don’t even want to think about how much time I’ve wasted trying to make sure all of the dominoes are in place before a project begins, when my time would have been much better spent actually producing something useful.

Unfortunately, I’ve heard that this problem is quite common among detail oriented people. This is probably why so many programmers (those who are inherently detail oriented) so frequently abandon software projects with excellent potential.

The only good thing to come of this whole realization, is that I am able to identify one of my biggest mistakes and (finally) move on. It is important to realize that getting hung up in the little details can turn your product from good to great, but it can just as easily run it into the ground when it renders you unproductive. Finding that balance between worrying about details and moving forward is absolutely vital to any project and can be the difference in the survival of a project. After all, It’s not the name, colors, or programming language selection of a product that gets people excited, it’s what your product can do to make their lives easier or better that will get them interested.

Random Posts

PDF Download: Why’s (poignant) Guide to Ruby

2010-01-13T22:01:03Z

This is one of the best (and most entertaining) introductions to a programming language that I have ever read. If you are even remotely interested in learning Ruby, _why’s guide is a great place to start.

And, wherever you are _why, thanks for such a fantastic read.

PDF download: Why’s (poignant) Guide to Ruby

Honing Your Programming Skills with Project Euler

2010-01-10T18:24:21Z

There’s nothing quite like that feeling you get when you finally figure out the solution to a problem. Whether it is solving a complex math equation, beating your friend at Tic-Tac-Toe, or even successfully programming the DVR for the first time. As programmers, we are constantly looking for new problems to solve and interesting ways to improve our programming skills.

Project Euler is a programming challenge site that hosts (currently) 272 programming problems of varying difficulty. Users can log in, attempt to solve each problem, and discuss each solution.

Here is an example of one of the Project Euler problems.

Problem #1

If we list all the natural numbers below 10 that are multiples of 3 or 5, we get 3, 5, 6 and 9. The sum of these multiples is 23.

Find the sum of all the multiples of 3 or 5 below 1000.

The Solution

This problem instantly reminds me of the FizzBuzz problem, so I think I’ll start with that for my solution.

One of my goals for this year is to become more proficient in Ruby and C#, so I’ll solve each problem using a mixture of both languages.

The Ruby Solution

max = 1000
sum = 0

(1..(max-1)).each do |i|
	 if ((i % 3 == 0) || (i % 5 == 0))
		sum += i
	 end
end

puts sum

The C# Solution

int max = 1000;
int sum = 0;

for (int x = 1; x < max; x++) {
	if ((x % 3 == 0) || (x % 5 == 0)) {
		sum += x;
	}
}

Console.WriteLine(sum);

Pretty straightforward, right? Simply sum up all of the values that are evenly divisible by 3 or 5 and print the final sum. This is probably the solution 99% of programmers would come up with on their first attempt. But, what about other ways of solving the same problem? For instance, can this problem be solved without looping through each number?

A better way to solve the problem

One of the best parts about Project Euler is the large community of smart developers all working on the same problem sets. Once you submit a correct answer to a problem, you will gain access to a forum thread for that particular problem and will be allowed to discuss your solution with everyone else who has solved the same problem.

This is one of the more interesting solutions I saw on the problem (converted to Ruby)

x = 1000;
puts 1.5*((x-1)/3)*((x+2)/3) + 2.5*((x-1)/5)*((x+4)/5) - 7.5*((x-1)/15)*((x+14)/15);

This solutions utalizes arithmetic progression to create a more elegant solution to the original problem.

Would you think of solving this problem using arithmetic progression? Probably not (at least for the first attempt anyway), I sure didn't.

Conclusion

Project Euler is a great way to improve your programming and problem solving skills. I hope this brief introduction was interesting enough to motivate you to check Project Euler out for yourself and try to solve some other problems on your own. Practicing these types of problems will (at the very least) better prepare you for your next job interview, which is never a bad thing.

Image From: Mathematician Pictures

3 Ways to Suck Less in 2010

2010-01-01T22:11:59Z

My favorite bit of advice on becoming a better programmer is to suck less every year. Programming is a very complex field that is constantly changing and requires a great deal of dedication from developers. Since 2010 officially started today, I thought it would be useful to share 3 goals which will help improve the programming skills of anyone who follows them.

My 3 ways to suck less by 2011

Read – stuffing your brain with great programming knowledge will be invaluable to us as developers. Reading about programming will not only improve your programming skills but it will introduce you to different ways of thinking about programming.

It’s also important to read source code. Writing code is one thing, but being able to read and understand source code that another programmer wrote is a great way to understand how other people solve similar problems. Three books I highly recommend are: Code Complete, The Pragmatic Programmer, and Joel on Software.
Write – Practice, practice, practice. The only way to truly get good at programming is to write code and lots of it. Reading or talking about programming is good up to a point, but nothing will improve your skill-set like actually going through the programming process. One thing I really want to do this year, is contribute to open source projects. My professional career has benefited so much from the hard work of programmers that make their software available for free, I feel like it’s my turn to give back a little. Contributing to open source is a great way to write code, work with other people, and improve the overall quality of software in the world.
Teach – teaching is a way to help yourself by helping others. You may think you know the way regular expressions work until you have to explain them to someone. Helping others not only allows you to give back to the development community, it will also help solidify your knowledge of technical subjects.

Don’t know where to get started? Create a StackOverflow account and start answering questions from other programmers just like you. Stick with it, and you’ll be surprised with how much knowledge you have gained over time.

Already have a StackOverflow account? Not a programmer? Start a blog on your hobby/trade and write about topics that interest you. You’ll learn more by going through the process of writing about interesting things and, hopefully, other people will find your work useful and learn from it.

Conclusion

Remember, great programmers didn’t become awesome in just a day. They constantly worked on improving their skills year after year until they achieved programming mastery.

Here’s to a great 2010 and I hope yours is full of awesome adventures in programming.

Happy New Year,
-Robert

Honing Your Programming Skills with Project Euler

B&N Nook: Unboxing the World’s Most Advanced eBook Reader

2009-12-09T00:30:01Z

I got the Barnes & Noble nook as a late birthday present. I’m one of the first people in the world to get one, so I thought I’d share my experience. I will write a more detailed review once I get a chance to actually use the nook. Enjoy.

First Impressions

Sexy Packaging: B&N took a page right out of the Apple playbook.
The tactile feeling on the back of the device is both tacky (feeling) and comfortable
The eInk screen is an absolute pleasure to look at
Registration was a breeze
The screen saver shows portraits of authors
Easy connection to WiFi

Here are some pictures from the unboxing.

Update: I tried subscribing to PC Magazine ($1.95) using my free $10.00 gift card from Barnes & Noble, but couldn’t because of an error. This isn’t a huge deal (to me) by any means, but it would have been nice to have everything work out of the box.

That’s it for now, I’m going to go read.

Random Posts

Getting Started With Git and TortoiseGit on Windows

2009-12-10T22:18:07Z

Version control is essential to the success of any software project. It provides the ability for multiple developers to work on the same codebase simultaneously and allows projects to be versioned for release. However, a great deal of programmers fail to leverage the great benefits of version control for their personal projects, typically due to laziness or complacency (myself included).

Luckily, for us lazy programmers, Git and TortoiseGit are extremely easy to install and configure on Windows. Now, there’s no excuse not to version control your code!

Here are the topics we will cover today:

Install Git
Install Tortoise Git
Create a new Git repository
Add an initial codebase to the Git repository
Change a file and commit
Clone the repository

Install Git

I installed the Git package available through Cygwin because I already had Cygwin installed.

Simply select the git packages and proceed with the normal Cygwin installation. If you need help with this, visit the Cygwin site for more details.

Another commonly used option is msysgit.

NOTE: msysgit and Cygwin conflict with each other, so if you already have one installed, I wouldn’t recommend using the other.

Lastly, be sure you point your system path to wherever the Git installation resides on your machine.

Install Tortoise Git

Download, Install, and restart your machine.

Create a new Git repository

Create a folder — I’m using D:\repo — to be used for your Git repository.
Right-click the folder and select Git Create Repository Here

You should now be greeted with the following message

Add an initial codebase to the Git repository

Copy your source files to the git folder.
Right click on the folder and select Git Commit
In this case, I added two files. I’m going to check both of them
click OK.

Now we have a legitimate repository with actual files in it.

Change a file and commit

Now you can change, add, or delete files within the repository.
Once you are ready to commit back to the repository, right click and select Git Commit just like above.

Clone the repository

Cloning is one of the great features of Git, and while it isn’t necessarily important for solo projects, it is important enough to mention here.

Tortoise allows you to clone projects by right clicking on the desired clone folder and selecting Clone
In the dialog, select the parent repository.
Click OK
Now, the cloned repository will be synced up with the main repository. This is useful for multi-person teams.

Conclusion

Git and TortoiseGit make it extremely easy to get personal version control up and running on your machine. Now you can develop your projects with the piece of mind that version control provides with little hassle and setup.

Random Posts

Working With Custom Time Strings in C#

2009-12-07T19:16:14Z

Let’s say we need to process some external data where the format and contents of the files are out of our control. To make matters worse, what if those same external files store tiemstamps in a foreign format that DateTime.Parse() does not understand?

It would be ill-advised to keep the dates in string format due to the level of robust date processing capabilities present in .NET.

Luckily, .NET provides a nice way to create custom date formats so any date string can be easily converted into a DateTime object.

//locale information
CultureInfo provider = CultureInfo.InvariantCulture;

//The timestamp from the external file
string dateString = "Dec 07 01:32:25 2009"; 

//the new format
string format = "MMM dd HH':'mm':'ss yyyy";

//Convert the time string into a legitimate DateTime object
DateTime result = DateTime.ParseExact(dateString, format, provider);

Now, we have the benefit of leveraging both custom timestamp values and .NET’s built in DateTime Functionality.

Do you think this might be useful on your current project? Learn more about DateTime.ParseExact() here here.

Image From: kenyonreview.org

Climategate Code Analysis Part 2

2009-12-05T22:55:53Z

There are three common issues that have been raised in my previous post that I would like to officially address concerning the CRU’s source code.

If you only get one thing from this post, please get this. I am only making a statement about the research methods of the CRU and trying to show proof that they had the means and intent to falsify data. And, until the CRU’s research results can be verified by a 3rd party, they cannot be trusted.

Here are the four most frequent concerns dealing with the CRU’s source code:

The source code that actually printed the graph was commented out and, therefore, is not valid proof.
No proof exists that shows this code was used in publishing results.
Interpolation is a normal part of dealing with large data sets, this is no different.
You need the raw climate data to prove that foul play occurred.

If anyone can think of something I missed, please let me know.

The source code that actually printed the graph was commented out and, therefore, is not valid proof.

Had I done a better job with my source analysis, I would have found a later revision of the briffa_sep98_d.pro source file (linked to in my previous post) contained in a different working tree which shows the fudge-factor array playing a direct result in the (uncommented) plotting of the data.

Snippit from: harris-tree/briffa_sep98_e.pro (see the end of the post for the full source listing)

;
; APPLY ARTIFICIAL CORRECTION
;
yearlyadj=interpol(valadj,yrloc,x)
densall=densall+yearlyadj
  ;
  ; Now plot them
  ;
  filter_cru,20,tsin=densall,tslow=tslow,/nan
  cpl_barts,x,densall,title='Age-banded MXD from all sites',$
    xrange=[1399.5,1994.5],xtitle='Year',/xstyle,$
    zeroline=tslow,yrange=[-7,3]
  oplot,x,tslow,thick=3
  oplot,!x.crange,[0.,0.],linestyle=1
  ;

Now, we can finally put this concern to rest.

Interpolation is a normal part of dealing with large data sets, this is no different.

This is partially true, the issue doesn’t lie in the fact that the CRU researchers used interpolation. The issue is the weight of the valadj array with respect to the raw data. valadj simply introduces too large of an influence to the original data to do anything productive with it.

Here is the graph I plotted of the valadj array. When we’re talking about trying to interpret temperature data that grows on the scale of one-tenths of a degree over a period of time, “fudging” a value by 2.5 is going to have a significant impact on the data set.

No proof exists that shows this code was used in publishing results.

Correct! That’s why I am (and always have) taken the following stand: Enough proof exists that the CRU had both the means and intent to intentionally falsify data. This means that all of their research results cannot be trusted until they are verified. Period.

The fact that the “fudge-factor” source code exists in the first place is reason enough for alarm. Hopefully, they didn’t use fudged results in the CRU research results, but the truth is, we just don’t know.

You need the raw climate data to prove that foul play occurred.

This is assuming the raw data are valid, which I maintain that it probably is. Several people question the validity of the climate data gathering methods used by the different climate research institutions, but I am not enough of a climate expert to have an opinion one way or the other. Furthermore, It simply doesn’t matter if the raw climate data are correct or not to demonstrate the extreme bias the valadj array forces on the raw data.

So, the raw data could actually be temperature data or corporate sales figures, the result is the same; a severe manipulation of data.

Full Source Listing

As promised, here is the entire source listing for: harris-tree/briffa_sep98_e.pro

;
; PLOTS 'ALL' REGION MXD timeseries from age banded and from hugershoff
; standardised datasets.
; Reads Harry's regional timeseries and outputs the 1600-1992 portion
; with missing values set appropriately.  Uses mxd, and just the
; "all band" timeseries
;****** APPLIES A VERY ARTIFICIAL CORRECTION FOR DECLINE*********
;
yrloc=[1400,findgen(19)*5.+1904]
valadj=[0.,0.,0.,0.,0.,-0.1,-0.25,-0.3,0.,-0.1,0.3,0.8,1.2,1.7,2.5,2.6,2.6,$
  2.6,2.6,2.6]*0.75         ; fudge factor
if n_elements(yrloc) ne n_elements(valadj) then message,'Oooops!'
;
loadct,39
def_1color,20,color='red'
plot,[0,1]
multi_plot,nrow=4,layout='large'
if !d.name eq 'X' then begin
  window, ysize=800
  !p.font=-1
endif else begin
  !p.font=0
  device,/helvetica,/bold,font_size=18
endelse
;
; Get regional tree lists and rbar
;
restore,filename='reglists.idlsave'
harryfn=['nwcan','wnam','cecan','nweur','sweur','nsib','csib','tib',$
  'esib','allsites']
;
rawdat=fltarr(4,2000)
for i = nreg-1 , nreg-1 do begin
  fn='mxd.'+harryfn(i)+'.pa.mean.dat'
  print,fn
  openr,1,fn
  readf,1,rawdat
  close,1
  ;
  densadj=reform(rawdat(2:3,*))
  ml=where(densadj eq -99.999,nmiss)
  densadj(ml)=!values.f_nan
  ;
  x=reform(rawdat(0,*))
  kl=where((x ge 1400) and (x le 1992))
  x=x(kl)
  densall=densadj(1,kl)     ; all bands
  densadj=densadj(0,kl)     ; 2-6 bands
  ;
  ; Now normalise w.r.t. 1881-1960
  ;
  mknormal,densadj,x,refperiod=[1881,1960],refmean=refmean,refsd=refsd
  mknormal,densall,x,refperiod=[1881,1960],refmean=refmean,refsd=refsd
;
; APPLY ARTIFICIAL CORRECTION
;
yearlyadj=interpol(valadj,yrloc,x)
densall=densall+yearlyadj
  ;
  ; Now plot them
  ;
  filter_cru,20,tsin=densall,tslow=tslow,/nan
  cpl_barts,x,densall,title='Age-banded MXD from all sites',$
    xrange=[1399.5,1994.5],xtitle='Year',/xstyle,$
    zeroline=tslow,yrange=[-7,3]
  oplot,x,tslow,thick=3
  oplot,!x.crange,[0.,0.],linestyle=1
  ;
endfor
;
; Restore the Hugershoff NHD1 (see Nature paper 2)
;
xband=x
restore,filename='../tree5/densadj_MEAN.idlsave'
; gets: x,densadj,n,neff
;
; Extract the post 1600 part
;
kl=where(x ge 1400)
x=x(kl)
densadj=densadj(kl)
;
; APPLY ARTIFICIAL CORRECTION
;
yearlyadj=interpol(valadj,yrloc,x)
densadj=densadj+yearlyadj
;
; Now plot it too
;
filter_cru,20,tsin=densadj,tslow=tshug,/nan
cpl_barts,x,densadj,title='Hugershoff-standardised MXD from all sites',$
  xrange=[1399.5,1994.5],xtitle='Year',/xstyle,$
  zeroline=tshug,yrange=[-7,3],bar_color=20
oplot,x,tshug,thick=3,color=20
oplot,!x.crange,[0.,0.],linestyle=1
;
; Now overplot their bidecadal components
;
plot,xband,tslow,$
  xrange=[1399.5,1994.5],xtitle='Year',/xstyle,$
  yrange=[-6,2],thick=3,title='Low-pass (20-yr) filtered comparison'
oplot,x,tshug,thick=3,color=20
oplot,!x.crange,[0.,0.],linestyle=1
;
; Now overplot their 50-yr components
;
filter_cru,50,tsin=densadj,tslow=tshug,/nan
filter_cru,50,tsin=densall,tslow=tslow,/nan
plot,xband,tslow,$
  xrange=[1399.5,1994.5],xtitle='Year',/xstyle,$
  yrange=[-6,2],thick=3,title='Low-pass (50-yr) filtered comparison'
oplot,x,tshug,thick=3,color=20
oplot,!x.crange,[0.,0.],linestyle=1
;
; Now compute the full, high and low pass correlations between the two
; series
;
perst=1400.
peren=1992.
;
openw,1,'corr_age2hug.out'
thalf=[10.,30.,50.,100.]
ntry=n_elements(thalf)
printf,1,'Correlations between timeseries'
printf,1,'Age-banded vs. Hugershoff-standardised'
printf,1,'     Region    Full   <10   >10   >30   >50  >100'
;
kla=where((xband ge perst) and (xband le peren))
klh=where((x ge perst) and (x le peren))
ts1=densadj(klh)
ts2=densall(kla)
;
r1=correlate(ts1,ts2)
rall=fltarr(ntry)
for i = 0 , ntry-1 do begin
  filter_cru,thalf(i),tsin=ts1,tslow=tslow1,tshigh=tshi1,/nan
  filter_cru,thalf(i),tsin=ts2,tslow=tslow2,tshigh=tshi2,/nan
  if i eq 0 then r2=correlate(tshi1,tshi2)
  rall(i)=correlate(tslow1,tslow2)
endfor
;
printf,1,'ALL SITES',r1,r2,rall,$
  format='(A11,2X,6F6.2)'
;
printf,1,' '
printf,1,'Correlations carried out over the period ',perst,peren
;
close,1
;
end

What do you think?

Join the conversation! Don’t just accept what one person says, gather the evidence and make an informed decision.

Image From: britsattheirbest

Cube Antics

I’m Moving!

Random Posts

Starting a WordPress Theme From Scratch

Related Posts

Using .NET to Monitor a Directory for Changes

The Code

Explanation

Conclusion

Related Posts

301 Redirect with ASP .NET

Related Posts

Creating a Fluid Two-Column Layout With a Single Fixed Width Column

layout.html

test.css

Related Posts

Using #if and [Conditional()] to Keep Debug Code Under Control

Related Posts

Sending Email in ASP .NET

Related Posts

Worrying About Small Details Can Kill Your Project

Random Posts

PDF Download: Why’s (poignant) Guide to Ruby

Related Posts

Honing Your Programming Skills with Project Euler

Problem #1

The Solution

The Ruby Solution

The C# Solution

A better way to solve the problem

Conclusion

Related Posts

3 Ways to Suck Less in 2010

My 3 ways to suck less by 2011

Conclusion

Related Posts

B&N Nook: Unboxing the World’s Most Advanced eBook Reader

First Impressions

Random Posts

Getting Started With Git and TortoiseGit on Windows

Here are the topics we will cover today:

Install Git

Install Tortoise Git

Create a new Git repository

Add an initial codebase to the Git repository

Change a file and commit

Clone the repository

Conclusion

Random Posts

Working With Custom Time Strings in C#

Related Posts

Climategate Code Analysis Part 2

The source code that actually printed the graph was commented out and, therefore, is not valid proof.

Interpolation is a normal part of dealing with large data sets, this is no different.

No proof exists that shows this code was used in publishing results.

You need the raw climate data to prove that foul play occurred.

Full Source Listing

What do you think?

Related Posts