Scatterplot matrix in julia

pairs.jl

# Thanks to Joe Kington
#   http://stackoverflow.com/questions/7941207/is-there-a-function-to-make-scatterplot-matrices-in-matplotlib

using PyPlot

function pairs(data)
  (nobs, nvars) = size(data)
  (fig, ax) = subplots(nvars, nvars, figsize=(8,8))
  subplots_adjust(hspace=0.05, wspace=0.05)

  # Plot data
  for i = 1:nvars
  	for j = 1:nvars
  		if i != j
  			ax[i,j][:plot](data[:,j],data[:,i],"ob",mfc="none")
  		else
  			ax[i,j][:hist](data[:,i])
  		end
  		ax[i,j][:xaxis][:set_visible](false)
  		ax[i,j][:yaxis][:set_visible](false)
  	end
  end

  # Set tick positions
  for i = 1:nvars
  	ax[i,1][:yaxis][:set_ticks_position]("left")
  	ax[i,end][:yaxis][:set_ticks_position]("right")
  	ax[1,i][:xaxis][:set_ticks_position]("top")
  	ax[end,i][:xaxis][:set_ticks_position]("bottom")
  end

  # Turn ticks on
  cc = repmat([nvars, 1],int(ceil(nvars/2)),1)
  for i = 1:nvars
  	ax[i,cc[i]][:yaxis][:set_visible](true)
  	ax[cc[i],i][:xaxis][:set_visible](true)
  end
end

# Example
d = rand(100,6)
pairs(d)