msund/gist:faf8f09c67bcd3d6ccca

## gistfile1.txt
{
 "metadata": {
  "name": "",
  "signature": "sha256:57e3953da6a8dc30a7f1354c37962b066401ca7fca37594140131e09720176cc"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "install.packages(\"devtools\")  # so we can install from github\n",
      "library(\"devtools\")\n",
      "library(reshape)\n",
      "install_github(\"ropensci/plotly\")  # plotly is part of ropensci\n",
      "install_github(\"takluyver/IRdisplay\")\n",
      "library(plotly)\n",
      " \n",
      "py <- plotly(username=\"r_user_guide\", key=\"mw5isa4yqp\")  # open plotly connection"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Installing package into \u2018/Users/matthewsundquist/Library/R/3.1/library\u2019\n",
        "(as \u2018lib\u2019 is unspecified)\n"
       ]
      },
      {
       "ename": "ERROR",
       "evalue": "Error in contrib.url(repos, \"source\"): trying to use CRAN without setting a mirror\n",
       "output_type": "pyerr",
       "traceback": [
        "Error in contrib.url(repos, \"source\"): trying to use CRAN without setting a mirror\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Downloading github repo ropensci/plotly@master\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Installing plotly\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "'/Library/Frameworks/R.framework/Resources/bin/R' --vanilla CMD INSTALL  \\\n",
        "  '/private/var/folders/vm/4bcx1bc16rd_566rnf2nwypw0000gn/T/Rtmpg0rkDG/devtools26b640b1d4cb/ropensci-plotly-c747575'  \\\n",
        "  --library='/Users/matthewsundquist/Library/R/3.1/library' --install-tests \n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Downloading github repo takluyver/IRdisplay@master\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Installing IRdisplay\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "'/Library/Frameworks/R.framework/Resources/bin/R' --vanilla CMD INSTALL  \\\n",
        "  '/private/var/folders/vm/4bcx1bc16rd_566rnf2nwypw0000gn/T/Rtmpg0rkDG/devtools26b65bbee5cf/takluyver-IRdisplay-91cc8d2'  \\\n",
        "  --library='/Users/matthewsundquist/Library/R/3.1/library' --install-tests \n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Loading required package: RCurl\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Loading required package: bitops\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Loading required package: RJSONIO\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Loading required package: ggplot2\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "dat <- data.frame(xx = c(runif(100,20,50),runif(100,40,80),runif(100,0,30)),yy = rep(letters[1:3],each = 100))\n",
      "plot <- ggplot(dat, aes(x=xx, fill=yy)) + geom_histogram(alpha=0.2, position=\"identity\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 4
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "plot"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 5,
       "text": []
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHgCAYAAAB91L6VAAAEJGlDQ1BJQ0MgUHJvZmlsZQAAOBGFVd9v21QUPolvUqQWPyBYR4eKxa9VU1u5GxqtxgZJk6XtShal6dgqJOQ6N4mpGwfb6baqT3uBNwb8AUDZAw9IPCENBmJ72fbAtElThyqqSUh76MQPISbtBVXhu3ZiJ1PEXPX6yznfOec7517bRD1fabWaGVWIlquunc8klZOnFpSeTYrSs9RLA9Sr6U4tkcvNEi7BFffO6+EdigjL7ZHu/k72I796i9zRiSJPwG4VHX0Z+AxRzNRrtksUvwf7+Gm3BtzzHPDTNgQCqwKXfZwSeNHHJz1OIT8JjtAq6xWtCLwGPLzYZi+3YV8DGMiT4VVuG7oiZpGzrZJhcs/hL49xtzH/Dy6bdfTsXYNY+5yluWO4D4neK/ZUvok/17X0HPBLsF+vuUlhfwX4j/rSfAJ4H1H0qZJ9dN7nR19frRTeBt4Fe9FwpwtN+2p1MXscGLHR9SXrmMgjONd1ZxKzpBeA71b4tNhj6JGoyFNp4GHgwUp9qplfmnFW5oTdy7NamcwCI49kv6fN5IAHgD+0rbyoBc3SOjczohbyS1drbq6pQdqumllRC/0ymTtej8gpbbuVwpQfyw66dqEZyxZKxtHpJn+tZnpnEdrYBbueF9qQn93S7HQGGHnYP7w6L+YGHNtd1FJitqPAR+hERCNOFi1i1alKO6RQnjKUxL1GNjwlMsiEhcPLYTEiT9ISbN15OY/jx4SMshe9LaJRpTvHr3C/ybFYP1PZAfwfYrPsMBtnE6SwN9ib7AhLwTrBDgUKcm06FSrTfSj187xPdVQWOk5Q8vxAfSiIUc7Z7xr6zY/+hpqwSyv0I0/QMTRb7RMgBxNodTfSPqdraz/sDjzKBrv4zu2+a2t0/HHzjd2Lbcc2sG7GtsL42K+xLfxtUgI7YHqKlqHK8HbCCXgjHT1cAdMlDetv4FnQ2lLasaOl6vmB0CMmwT/IPszSueHQqv6i/qluqF+oF9TfO2qEGTumJH0qfSv9KH0nfS/9TIp0Wboi/SRdlb6RLgU5u++9nyXYe69fYRPdil1o1WufNSdTTsp75BfllPy8/LI8G7AUuV8ek6fkvfDsCfbNDP0dvRh0CrNqTbV7LfEEGDQPJQadBtfGVMWEq3QWWdufk6ZSNsjG2PQjp3ZcnOWWing6noonSInvi0/Ex+IzAreevPhe+CawpgP1/pMTMDo64G0sTCXIM+KdOnFWRfQKdJvQzV1+Bt8OokmrdtY2yhVX2a+qrykJfMq4Ml3VR4cVzTQVz+UoNne4vcKLoyS+gyKO6EHe+75Fdt0Mbe5bRIf/wjvrVmhbqBN97RD1vxrahvBOfOYzoosH9bq94uejSOQGkVM6sN/7HelL4t10t9F4gPdVzydEOx83Gv+uNxo7XyL/FtFl8z9ZAHF4bBsrEwAAMqRJREFUeAHt3XtsHNX58PHH6/XasR07cRLIhSQQUgiEUkhTWmhp0xZKWwQRChQEqARxSYBCUZGgUi/qBcEfVfsHFVKoaJFS0QJtKZCqhZRLuL0VUAL6JSRQyK00gVxI4tiO7/brZ8I6Y2cvs+vdc87s+Q4ymZ2Znec5n3N2n53Z3dmqwaFJmBBAAAEEEEDAqEDCaDSCIYAAAggggEAgQAFmICCAAAIIIGBBgAJsAZ2QCCCAAAIIUIAZAwgggAACCFgQoABbQCckAggggAACFGDGAAIIIIAAAhYEKMAW0AmJAAIIIIAABZgxgAACCCCAgAUBCrAFdEIigAACCCBAAWYMIIAAAgggYEGAAmwBnZAIIIAAAghQgBkDCCCAAAIIWBCgAFtAJyQCCCCAAALJOBDs2bPHSJrJZFIGBgaCPyMBKyRIKpWSnp6eCmmNmWZUVVWJjrfe3l4zASskiprpD7j19/dXSIvMNEMfozrWbP/43eTJk800OCZRYlGA29vbjXBOmTJFOjs7xVS8QhpVV1cnXV1dhdzF2LazZ8+W999/38kXLq66aV6NjY2yb98+Y/0UNVAikZDq6monXxy0tLQExdfFx2htbW3wQtR2kcvUzzNmzAie12y/UKYAj+wdTkGP9OAWAggggAACRgQowEaYCYIAAggggMBIAQrwSA9uIYAAAgggYESAAmyEmSAIIIAAAgiMFKAAj/TgFgIIIIAAAkYEKMBGmAmCAAIIIIDASAHnv4akH+nXr2yYmPSrFzU1NcbiFdIm/f6jKYdC8kpvq7npd6hdm1x10+9l6td9XOxT/Y6y/unjwbVJc9LcXHRL5+aameajZukx52J+vubkfAHWgWPq+6/65X79srqpeIUMOn3CcTGvdBs0NxcLsMtu6uVin7r8PeD6+vrge8Auurn8PWA9kNHvANv+HnD6+YJ/DwlwCpqRgAACCCCAgAUBCrAFdEIigAACCCBAAWYMIIAAAgggYEGAAmwBnZAIIIAAAghQgBkDCCCAAAIIWBCgAFtAJyQCCCCAAAIUYMYAAggggAACFgQowBbQCYkAAggggAAFmDGAAAIIIICABQEKsAV0QiKAAAIIIOD8pSjpIgR8Fkju2B6p+X3TZ0Tajo0QQMAdAY6A3ekLMkEAAQQQ8EiAAuxRZ9NUBBBAAAF3BCjA7vQFmSCAAAIIeCRAAfaos2kqAggggIA7AhRgd/qCTBBAAAEEPBKgAHvU2TQVAQQQQMAdAQqwO31BJggggAACHglQgD3qbJqKAAIIIOCOAAXYnb4gEwQQQAABjwQowB51Nk1FAAEEEHBHgALsTl+QCQIIIICARwIUYI86m6YigAACCLgjQAF2py/IBAEEEEDAIwEKsEedTVMRQAABBNwRoAC70xdkggACCCDgkQAF2KPOpqkIIIAAAu4IUIDd6QsyQQABBBDwSIAC7FFn01QEEEAAAXcEKMDu9AWZIIAAAgh4JEAB9qizaSoCCCCAgDsCFGB3+oJMEEAAAQQ8EqAAe9TZNBUBBBBAwB2BpDupVGYm7+9NRW7YzJaeyNuyIQIIIIBAvAU4Ao53/5E9AggggEBMBSjAMe040kYAAQQQiLcABTje/Uf2CCCAAAIxFaAAx7TjSBsBBBBAIN4CFOB49x/ZI4AAAgjEVIACHNOOI20EEEAAgXgLUIDj3X9kjwACCCAQUwEKcEw7jrQRQAABBOItQAGOd/+RPQIIIIBATAUowDHtONJGAAEEEIi3AAU43v1H9ggggAACMRWgAMe040gbAQQQQCDeAhTgePcf2SOAAAIIxFSAAhzTjiNtBBBAAIF4C1CA491/ZI8AAgggEFMBCnBMO460EUAAAQTiLUABjnf/kT0CCCCAQEwFKMAx7TjSRgABBBCItwAFON79R/YIIIAAAjEVSJrI+4knnpA5c+bIKaecEoRbvXq1rFu3LphvamqS6667zkQaxEAAAQQQQMAZgbIW4O7ubnn44YflnXfekVmzZg03+q233pLly5dLKpWSqqqq4eXMIIAAAggg4ItAWQtwR0eHLFy4UCZOnDjsOTAwILp8/fr1Mjg4KAsWLBhepzMvv/yyvPTSS8PLZs6cKUuWLBm+Xc6ZZDIZvCgYP358ycJ0JqK/wJh29GDWuIlEQtTO1enoo492MjVX3fSFp463adOm5XQb7OnOuT69sirPftLbRfk3/aJYH5+uTWqmedXX17uWmrg61hRK3SZPnhzYOQfncUJlLcAtLS2if5s2bRombmtrk8bGxuDvgw8+kBUrVsgtt9wyvH7evHnBQEkv0Afavn370jfL+u+ECRNEj9o7OztLFufAgejE+1J9WePq2YKenp6s622u0OLb2trq5AsEV900L32hl29sJw4ciNS1AyV8jGgB1mLS398fKbbJjdQs/SLeZNwosWpqaqSvr8/JIqfF98DQWNL8bE5Tp061Gd652NGrQ4lSb25ulltvvTXY2/z582Xt2rXBk7cu12nSpEnBX3Dj4/9t3bo1fLNs8/qE09vbK11dXSWL0dOTiryvrq7cBbaUeUVOKuKGmpurR+iuuqlXvtySEV909ZVwzGrxra6uDh4LEbvf2Gb6glwfp/ncjCUUCqRH5voi2cUzB+ncXH0RH2L0atb4p6D37NkjK1euDJD11ZgOiFKe8vWq92gsAggggEBsBYwfAeupkIaGBrn//vtl9+7dcs455wSnu2IrSOIIIIAAAggUIWCkAF9wwQUjUtMPVemp3vSprhEruYEAAggggIAHAkYKcCZH/cACEwIIIIAAAr4KGH8P2Fdo2o0AAggggEBYgAIc1mAeAQQQQAABQwIUYEPQhEEAAQQQQCAsQAEOazCPAAIIIICAIQEKsCFowiCAAAIIIBAWoACHNZhHAAEEEEDAkAAF2BA0YRBAAAEEEAgLUIDDGswjgAACCCBgSIACbAiaMAgggAACCIQFKMBhDeYRQAABBBAwJEABNgRNGAQQQAABBMICFOCwBvMIIIAAAggYEqAAG4ImDAIIIIAAAmEBCnBYg3kEEEAAAQQMCVCADUETBgEEEEAAgbAABTiswTwCCCCAAAKGBCjAhqAJgwACCCCAQFiAAhzWYB4BBBBAAAFDAhRgQ9CEQQABBBBAICxAAQ5rMI8AAggggIAhAQqwIWjCIIAAAgggEBagAIc1mEcAAQQQQMCQAAXYEDRhEEAAAQQQCAtQgMMazCOAAAIIIGBIgAJsCJowCCCAAAIIhAUowGEN5hFAAAEEEDAkQAE2BE0YBBBAAAEEwgIU4LAG8wgggAACCBgSoAAbgiYMAggggAACYQEKcFiDeQQQQAABBAwJUIANQRMGAQQQQACBsAAFOKzBPAIIIIAAAoYEKMCGoAmDAAIIIIBAWIACHNZgHgEEEEAAAUMCFGBD0IRBAAEEEEAgLEABDmswjwACCCCAgCEBCrAhaMIggAACCCAQFqAAhzWYRwABBBBAwJAABdgQNGEQQAABBBAIC1CAwxrMI4AAAgggYEiAAmwImjAIIIAAAgiEBSjAYQ3mEUAAAQQQMCRAATYETRgEEEAAAQTCAhTgsAbzCCCAAAIIGBKgABuCJgwCCCCAAAJhAQpwWIN5BBBAAAEEDAlQgA1BEwYBBBBAAIGwAAU4rME8AggggAAChgQowIagCYMAAggggEBYgAIc1mAeAQQQQAABQwJJQ3GKDjM4OCjV1dVF37+QO1ZVVUkikShpPN1f1ClXOzW3XOujxijXdpqb5uja5Kpbelzk69P0dvlc8+0n3/3D6zWm/pVyn+H9j2W+HI/RseQTvm86t/AyV+bTubnYp64Y2cjD+QKsKDp4TE6ljFfIvnJtq+tyrTfpkymWq/m5nle+Ps23Pt0XUbdLb5/vX5fdNPdStzefR5T1rpqlc3c9v3SePv3rfAHWQdPX12ekT/Roe2BgoKTx+vujHwHnamcymSxpXqUG1dzVzrXJVTfNS8dbrj4PLPv7I5Hm3U+kvRzaKH3UXcp9FhA+56Y6xkr9GM0ZsICVenTZP9Rf2q+uTZqT5uZin7pmZTKf6NXBZFbEQgABBBBAoMIFKMAV3sE0DwEEEEDATQEKsJv9QlYIIIAAAhUuQAGu8A6meQgggAACbgpQgN3sF7JCAAEEEKhwAQpwhXcwzUMAAQQQcFOAAuxmv5AVAggggECFC1CAK7yDaR4CCCCAgJsCFGA3+4WsEEAAAQQqXMD5K2FVuD/NK7FAcsf2EXtMpFKS7OkZsSx9o2/6jPQs/yKAAALGBTgCNk5OQAQQQAABBEQowIwCBBBAAAEELAhQgC2gExIBBBBAAAEKMGMAAQQQQAABCwIUYAvohEQAAQQQQIACzBhAAAEEEEDAggAF2AI6IRFAAAEEEKAAMwYQQAABBBCwIEABtoBOSAQQQAABBCjAjAEEEEAAAQQsCFCALaATEgEEEEAAAQowYwABBBBAAAELAhRgC+iERAABBBBAgALMGEAAAQQQQMCCAAXYAjohEUAAAQQQoAAzBhBAAAEEELAgQAG2gE5IBBBAAAEEKMCMAQQQQAABBCwIUIAtoBMSAQQQQAABCjBjAAEEEEAAAQsCFGAL6IREAAEEEECAAswYQAABBBBAwIJA0kLMigj5/t5UyduRa5+pVFJ6eg7FnNnSU/LYpd5hcsf2SLvsmz4j0nZslFvg3W3/L/cGH6/9xOyzIm1Xjo3ePvhepN3Oq58baTs2QiDuAhwBx70HyR8BBBBAIJYCFOBYdhtJI4AAAgjEXYACHPceJH8EEEAAgVgKUIBj2W0kjQACCCAQdwEKcNx7kPwRQAABBGIpQAGOZbeRNAIIIIBA3AUowHHvQfJHAAEEEIilAAU4lt1G0ggggAACcRegAMe9B8kfAQQQQCCWAhTgWHYbSSOAAAIIxF2AAhz3HiR/BBBAAIFYClCAY9ltJI0AAgggEHcBCnDce5D8EUAAAQRiKUABjmW3kTQCCCCAQNwFKMBx70HyRwABBBCIpQAFOJbdRtIIIIAAAnEXoADHvQfJHwEEEEAglgIU4Fh2G0kjgAACCMRdgAIc9x4kfwQQQACBWAokTWT9xBNPyJw5c+SUU04Jwm3cuFHWrFkjfX19cskll8jUqVNNpEEMBBBAAAEEnBEo6xFwd3e3rFy5Ul555ZWg2Gqru7q6ZNWqVXLVVVfJRRddJA899JAzGCSCAAIIIICAKYGyHgF3dHTIwoULZeLEicPt2blzp8ycOVPq6+uDv87OzqA4J5OHUjlw4IDoX3qqra2VVCqVvlnWfxOJhGgeUeLV1NSUNZfRO9e8BgcHg8WGOEankPO2mg0MDAxvUx3RJ1HixoyOG3YbTu7jmVLHHr3/XLd1/FRVVeUda6Pbk22fycFoD+UoY1vzqq6uDvLLFq+Y5TW90R4zuXLUx6hOubYpJrdS3CfXWCvF/seyD+1TzY/JLYGy9khLS4vo36ZNm4ZbvX///qDwpheMGzdOtFA3NzcHi9auXSvPPfdcerUce+yxcuWVVw7fLueMPunoE2NDQ0PeMG0DVXm3KdcGU6YcKsSl2v/mndHaMufo7HEnTZo0Ip3Bgx0jbme7sWXgqGyrRizPFTu8YdS4ep+qKVPCdzU6r0+IWkym5MkhanuaOvdGyj9fPN2J5qZT+gVfcCPL/9a1bsyy5sjFE1OHX4gfufbwkinN2fslXYD1ecO1qRA307nrc1v4QMh0fOJlFihrAc4Usq6uTvTUdHrq7e0dUfAWLVok+heetm7dGr5Ztnl9ctIj8vb29rwxdu01c1SeTkRf8ff09AQ3a/sO/ZteN9Z/o7YlW9zZs2fLBx98MOIIOLlrV6S0dtXlf7GjO8oWe3SQ0XHDbqO37UtGOyIbfb9S3NbHgT4hqluuaXR7sm27tzdaAd6+fXu2XQwv1yKnT9j62Mw37ToYrZ/z7Se8fnt7U/jmiHl9Qd/f3y+tra0jlrtwQ8/W6WM0ygsX0/nOmDFDdu/ePfwcYjp+Op4eUDEdFijre8CHwxyemzZtmuiTgA5SLXb6L6dGDvswhwACCCDgh4DxI+CmpqbgfeEVK1ZIW1ubLF682A9pWokAAggggEBIwEgBvuCCC0IhRc4++2w588wzg/fA0u/pjNiAGwgggAACCFS4gJECnMmQ086ZVFiGAAIIIOCLgPH3gH2BpZ0IIIAAAqURuPPOO2Xu3LnBzj788EPRDzHec889oh9827JlS7D81FNPlZ/97GelCWhoLxRgQ9CEQQABBBAoTuCKK66QzZs3Bxd1evTRR+Woo46Sm266SfTT3Y888ohs2LBB1q1bJ5dffnlxASzdy9opaEvtJSwCCCCAQMwEjjvuuODrqVps33jjjeDaEPpVuaVLlwYFWL/aqp8rSh8lx6V5HAHHpafIEwEEEPBY4Oqrr5YHH3xQXnjhBfn2t78dSOgljd9880257777hpfFiYgCHKfeIlcEEEDAU4ElS5bIwYMHZcGCBTJv3rxAQS8C9OUvf1k++ugjufTSS2Mnwyno2HUZCSOAAAJ+CegVxvQ9YL1wkx716qTXntcPZOlvB+hXXeN4qU2OgP0ax7QWAQQQiJ3AM888I6eddlpw9KunonXSo2E9AtZLbP7kJz8JlsXtfxwBx63HyBcBBBDwTOAb3/jG8E/appve2NgY6Xrl6e1d/JcC7GKvkBMCCCCAwLBAqX7gIv2LVcM7tjxDAbbcAYRHAAEEEMgt0L9lc3DKOfdWudfqz8xWHzcn90aG11KADYMTDgEEEECgcIFSHQUXHrl89+BDWOWzZc8IIIAAAghkFaAAZ6VhBQIIIIAAAuUToACXz5Y9I4AAAgggkFWgqAKsVx0ZPe3atSu4Gsno5dxGAAEEEEAAgSMFCirAXV1don+LFi0K/k3f1i9E33LLLbJq1aojI7AEAQQQQAABBwXa2tqsZlVQAT7//PNl3Lhxsn79+uBfndc//UL0yy+/LGeddZbVxhAcAQQQQACBfALPP/+8fOELXwguazl//nzZtGlTvruUZX1BX0N6+umnpb+/P/jViZUrVw4nlEgkRP+YKlvg/b2pjA3sqxXZObROr82anqrbGtKzuf+ty726nGuTO7ZH2n3f9BmRtovDRlHarI/lwZmz4tAcckSgKAG9fKX+stKsWbPktttuk2effVaOP/74ovY1ljsVVID1KiLJZFL+8Ic/jCUm90UAAQQQQMCaQEtLiyxbtkxqa2tl37598olPfMJKLkUdtr700kvBT0CddNJJcuKJJw7/Pfnkk1YaQVAEEEAAAQSiCtx6663ym9/8Rh5//HE5+eSTgzO7Ue9byu0KOgJOB166dKlcc801smjRouCIOL187ty56Vn+RQABBBBAwEmByy+/XG688UbRI+H29nbZsWOHlTwLLsB6OTD9/cXvf//74tqFra0IEhQBBBBAIFYCWr86OzulpqZmxEGk6UYUfApai+55550nv/vd70R/JJkJAQQQQACBuAnoN3j0M002p4ILsCar35269tprZdKkSTJv3rzhP94DttmVxEYAAQQQiJNAUeX/pz/9qfzgBz84op28B3wECQsQQAABBBDIKFBUAb7jjjvkww8/PGKHv/jFL+Tcc889YjkLEEAAAQQQKFYgcfRUaawq9t6H7lc19Pkl16aiCrAeAaff/9V/9cpYjzzyiCxcuNC19pEPAggggEDMBf6ve4t0d3ePqRV1dXXyqYb5Y9pHqe9cVAH+7Gc/OyKPr371q/K///1PnnnmGbn44otHrOMGAggggAACYxXo6+sb6y6cu39RH8Ia3Qr9atK2bduktbV19CpuI4AAAggggEAGgaKOgPV93p07dwa70+KrP0+o36e69957M4RgEQIIIIAAAgiMFiiqAN99990jzsfr96lOOeUUSaUyX6x/dFBuI4AAAggg4LtAUaeg9cNWZ555pkyfPl0mT54sp512GsXX95FE+xFAAIGYCNx+++3yxhtvWM+2qAL87rvvyumnnx78gsQZZ5wh48ePl/vuu896Y0gAAQQQQACBfAL6c4T6qeqxfrI6X5x864sqwNdff71ccMEFwfvAe/fuldWrV4t+NWnDhg354rEeAQQQQAAB6wJ63YpLL71U9CAy03UtTCRYVAHWQvvjH/84uBRldXW1fP7zn5err75aXnzxRRM5EwMBBBBAAIExCZx99tny2GOPif4y0j333DOmfRV756IKsH7g6s033xwR84UXXgjeDx6xkBsIIIAAAgg4KKAHjjrpZ5q2b99uJcOiCvANN9wQXHJy8eLFcvPNN8sJJ5wQnEu/8MILrTSCoAgggAACCBQisG7dumDzV155JfiFv0LuW6ptiyrAerUrverVZz7zmeDTz/rbii+99FLwXeBSJcZ+EEAAAQQQKJfA008/LUuWLJHXX389+ExTueLk2m9R3wM+cOCAPProo/Lzn/9c9D3gZcuWBd8D1jezmRBAAAEEEHBZ4IEHHgjS098ysHn9iqKOgPVT0HrZyYGBgaARX/nKV4JPk+knopkQQAABBBCIg4DN4qs+RR0B6znzzZs3S1XVod+H0o9y6yejn3jiCVm6dGkc3MkRAQRyCGzseFcq8eL3OZrMKgSMCxR1BNzY2CivvfbaiGSff/55aWpqGrGMGwgggAACCCCQWaCoI+A777xTzjnnHFmwYIHMnDlT1qxZIyeddJLwKejMyCxFAAEEECheYFJti0yoHtsBXiJZXXwCZbpnUQVYv37073//O7jwRkdHR/BhrGOPPbZMKbJbBBBAAAGfBboPTBvzWyLJ5FC5a3BLsagCrE3Q7/7qHxMCCCCAAALlFNCfve3q6hpTCH3r1LWpqPeAXWsE+SCAAAIIIBA3AQpw3HqMfBFAAAEEKkKAAlwR3UgjEEAAAQTiJkABjluPkS8CCCCAQEUIUIArohtpBAIIIIBA3AQowHHrMfJFAAEEECiJQHd3d0n2U+xOiv4aUrEBC72ffvy8rq6u0LsVtb3+sERNTU2keKmUWbrgO2wft6qurrSvm0rRFr2mavra4JpmYsgxylQ9dL8oU9Q2J0btL+wWJU6mbZJlGH/qlUgk8o610e3JlJ8uqxn6L8oU5dq3eonZZLI/yC/KPku9Ta7Huz5GNb9c25Q6n6j7S+cWdXuT26lZesyZjOtqrE2bNslll10mxxxzTODyxz/+0cp4N1tFiugNHThj/f5X1LD9/f3S29sbKV5PT7TCETV2lO30lzt06uo69G+U+0TZphRt0dzCBbh6yDHK1P9xm/JtG7XNyQz7S7vli5Ftfd8Yv3+Ybb/qlW9sZ2pPpv3puI0y9Uj+saMvDPQ60LauBZ3LpL6+XvRxmmubKA7l2Ka2tlZ0rOlBg2uT5qS5jfWx4Fq7is3njjvukLvuuiv4XfsHH3xQtm3bJscdd1yxuyv6fs4X4KJbxh0RQAABBBDIIKBHwOmfz73iiisybGFmUWnPZZrJmSgIIIAAAggULbBo0SJZu3ZtcP9f/epX8uqrrxa9r7HckQI8Fj3uiwACCCAQO4Fly5bJr3/9a/nWt74VFOKFCxdaaQOnoK2wExQBBBBAwJbAvHnz5NFHHw0+S2DzA30cAdsaAcRFAAEEELAqYLP4asMpwFa7n+AIIIAAAr4KcAra156n3QgggECMBPQrqZU2UYArrUdpDwIIIFBhAnOnafFtqLBWiVCAK65LaRACCCBQWQKVePSrPcR7wJU1TmkNAggggEBMBDgCjklHkab7Askd2yMl2Td9RqTt2Ci3wNsH38u9wcdr59XPjbRdIRtlip3qSwWXsi3mUpTlyLGQ9rCtHQGOgO24ExUBBBBAwHMBCrDnA4DmI4AAAgjYEaAA23EnKgIIIICA5wIUYM8HAM1HAAEEELAjQAG2405UBBBAAAHPBSjAng8Amo8AAgggYEeAAmzHnagIIIAAAp4LUIA9HwA0HwEEEEDAjgAF2I47URFAAAEEPBegAHs+AGg+AggggIAdAQqwHXeiIoAAAgh4LkAB9nwA0HwEEEAAATsCFGA77kRFAAEEEPBcgALs+QCg+QgggAACdgQowHbciYoAAggg4LkABdjzAUDzEUAAAQTsCFCA7bgTFQEEEEDAcwEKsOcDgOYjgAACCNgRoADbcScqAggggIDnAhRgzwcAzUcAAQQQsCNAAbbjTlQEEEAAAc8FKMCeDwCajwACCCBgRyBpJyxREShMYMc7+yPdYdb4SJuxEQJOCbx98L1I+cyrnxtpOzaKhwBHwPHoJ7JEAAEEEKgwAQpwhXUozUEAAQQQiIcABTge/USWCCCAAAIVJkABrrAOpTkIIIAAAvEQoADHo5/IEgEEEECgwgQowBXWoTQHAQQQQCAeAhTgePQTWSKAAAIIVJgABbjCOpTmIIAAAgjEQ4ACHI9+IksEEEAAgQoToABXWIfSHAQQQACBeAhQgOPRT2SJAAIIIFBhAhTgCutQmoMAAgggEA8BCnA8+oksEUAAAQQqTIACXGEdSnMQQAABBOIhQAGORz+RJQIIIIBAhQlQgCusQ2kOAggggEA8BCjA8egnskQAAQQQqDABCnCFdSjNQQABBBCIhwAFOB79RJYIIIAAAhUmkLTRntWrV8u6deuC0E1NTXLdddfZSIOYCCCAAAIIWBOwUoDfeustWb58uaRSKamqqrLWeAIjgAACCCBgS8B4AR4YGJCOjg5Zv369DA4OyoIFC0a0XYvzhg0bhpdNmTLliG2GV5Z4Rl8QJBKJ4IVBvl239pk9e19dXS11dXVBWi0tA/nSK2h9Kdoyfvz4oD+HAx88ODybc2boDEikKeL+mpqqR+wu7DZiRSE3Wlqibd3eFm27of1pXslkUlry7TviPhu6GiPFbqrL760vihsba0Qfqzamlqbs3voY0LzUrymZvy2af679Fdu+TLE1p/7+/mJ3Gel+xbYl8Bp6rNnq00iN83Aj4wW4ra1t6MHdGPx98MEHsmLFCrnllluG6Wtra0WfzNNTfX192Qd1Opa+INj0oUhvb3pJrn/NPjnpk2L6waM5RpmOnRLtySDq82zVzsyBB4b6aLCzc2QBjpLg0DaDEYNXDUbzThulw4fd0ssK/nfrloLvkusO6/dvCAqIjvWDeV5YnDIwLteuhtcNRvRZ17lp+D7ZZoKzUs3Thsdbtu3Ktfz/hnyyTWqmj9Oenp5smxyxvJCi+FbbO0fcP+oCddPc9K9cUyFtCeegOeljo9j7h/fFfOkEjBfg5uZmufXWW4MWzJ8/X9auXSutra2iy3WaO3du8Bfc+Ph/W7duDd8s27weAXd39+Z9UixbAjl2rLkV8qSju2pNRXuSam9P5Yh8eFV1lmKhL5I6hwrw6OJ3+J7Z5/rb27OvDK3JFju0STDbXt0xYlExbiN2UIYb7cnO4CxLTU2NtOdpf3t7tBdRB3sjnnGI0B49CzQ41Nd9fX0Rtja7ieam4yyfWzir1oHW8M2c8+0Ho43HTDvRsdY79Oq9nAW4kLaEc9SDHjUr9DkkvI9SzE+cOLEUu6mYfZg9jzrEtmfPHlm5cmUAqA9wHRDhI96KkaUhCCCAAAII5BAwfgQ8efJkaWhokPvvv192794t55xzTvC+a44cWYUAAggggEDFCRgvwCq4ZMmS4FSNnk7SDwcwIYAAAggg4JuAlQKsyPr+FxMCCCCAAAK+Chh/D9hXaNqNAAIIIIBAWIACHNZgHgEEEEAAAUMCFGBD0IRBAAEEEEAgLEABDmswjwACCCCAgCEBCrAhaMIggAACCCAQFqAAhzWYRwABBBBAwJAABdgQNGEQQAABBBAIC1CAwxrMI4AAAgggYEiAAmwImjAIIIAAAgiEBSjAYQ3mEUAAAQQQMCRAATYETRgEEEAAAQTCAhTgsAbzCCCAAAIIGBKgABuCJgwCCCCAAAJhAQpwWIN5BBBAAAEEDAlQgA1BEwYBBBBAAIGwAAU4rME8AggggAAChgQowIagCYMAAggggEBYIBm+wXxlCby/N+V8g6p37yppjv9taxixv5qaGuntrRmxLH1j1viO9KzRf7XNiaG8BjvapXr//pyxN+Zcy8ooAm8ffC/KZrHYppC2zKufG4s2+ZwkR8A+9z5tRwABBBCwJkABtkZPYAQQQAABnwUowD73Pm1HAAEEELAmQAG2Rk9gBBBAAAGfBSjAPvc+bUcAAQQQsCZAAbZGT2AEEEAAAZ8FKMA+9z5tRwABBBCwJkABtkZPYAQQQAABnwUowD73Pm1HAAEEELAmQAG2Rk9gBBBAAAGfBSjAPvc+bUcAAQQQsCZAAbZGT2AEEEAAAZ8FKMA+9z5tRwABBBCwJkABtkZPYAQQQAABnwUowD73Pm1HAAEEELAmQAG2Rk9gBBBAAAGfBSjAPvc+bUcAAQQQsCZAAbZGT2AEEEAAAZ8FKMA+9z5tRwABBBCwJkABtkZPYAQQQAABnwUowD73Pm1HAAEEELAmQAG2Rk9gBBBAAAGfBSjAPvc+bUcAAQQQsCZAAbZGT2AEEEAAAZ8FKMA+9z5tRwABBBCwJkABtkZPYAQQQAABnwUowD73Pm1HAAEEELAmQAG2Rk9gBBBAAAGfBSjAPvc+bUcAAQQQsCZAAbZGT2AEEEAAAZ8FKMA+9z5tRwABBBCwJkABtkZPYAQQQAABnwUowD73Pm1HAAEEELAmkLQWOWLgwcFBqauri7j12Darrq6WZHJQUqnU2HZUhnsnk/a7KlFTk7VlNUPrBgYGsq63tUL7NNtkq59rpGZonCUlkUiIurk2VVVVyeDH+bmWm/an5mer73J55Bprue5XrnXh5820mY45JncE7D+r57HQgdPV1ZVnq9Ks7u/vl76+Punp6SnNDku8F9t5Vff2Zm1R79A6FwuwJqy5ZZpseabzUa/0fKb8bC3TJ+nBoceBPhZcm7SoqJutvsvloS8KtD/1oMGFKfy8qTmpmYtuLljZyoGXQ7bkiYsAAggg4LUABdjr7qfxCCCAAAK2BCjAtuSJiwACCCDgtQAF2Ovup/EIIIAAArYEKMC25ImLAAIIIOC1AAXY6+6n8QgggAACtgQowLbkiYsAAggg4LUABdjr7qfxCCCAAAK2BCjAtuSJiwACCCDgtQAF2Ovup/EIIIAAArYEKMC25ImLAAIIIOC1AAXY6+6n8QgggAACtgQowLbkiYsAAggg4LUABdjr7qfxCCCAAAK2BCjAtuSJiwACCCDgtQAF2Ovup/EIIIAAArYEKMC25ImLAAIIIOC1AAXY6+6n8QgggAACtgQowLbkiYsAAggg4LUABdjr7qfxCCCAAAK2BCjAtuSJiwACCCDgtQAF2Ovup/EIIIAAArYEKMC25ImLAAIIIOC1AAXY6+6n8QgggAACtgQowLbkiYsAAggg4LUABdjr7qfxCCCAAAK2BCjAtuSJiwACCCDgtQAF2Ovup/EIIIAAArYEKMC25ImLAAIIIOC1AAXY6+6n8QgggAACtgQowLbkiYsAAggg4LUABdjr7qfxCCCAAAK2BCjAtuSJiwACCCDgtUDS69bTeAQQQKBCBd4++N5wy/a2HpB9Hfukt7d3eFl6Zl793PQs/xoW4AjYMDjhEEAAAQQQUAEKMOMAAQQQQAABCwIUYAvohEQAAQQQQIACzBhAAAEEEEDAggAF2AI6IRFAAAEEEKAAMwYQQAABBBCwIEABtoBOSAQQQAABBCjAjAEEEEAAAQQsCFCALaATEgEEEEAAAQowYwABBBBAAAELAhRgC+iERAABBBBAgALMGEAAAQQQQMCCAAXYAjohEUAAAQQQoAAzBhBAAAEEELAgQAG2gE5IBBBAAAEEKMCMAQQQQAABBCwIUIAtoBMSAQQQQAABCjBjAAEEEEAAAQsCFGAL6IREAAEEEEAgaYNg48aNsmbNGunr65NLLrlEpk6daiMNYiKAAAIIIGBNwPgRcFdXl6xatUquuuoqueiii+Shhx6y1ngCI4AAAgggYEvA+BHwzp07ZebMmVJfXx/8dXZ2BkfCyeShVJ566inRv/Q0Z84cufHGG9M3y/pvVVWVNDUNDv01lTVOMTvX3AYHB4u5a8nuMzAwkHVfEye2ZF1nc8UQ25Bb5gymtozPvKLMS3cd7AgiaG6TJk0uc7Tidl89ZYr18ZYp8/TjoLGxMdNqq8vSuVlNIktwzW3SpEkZ186aMCvjchaWX8B4Ad6/f39QeNNNGzdunHR0dEhzc3OwaNGiRfK5z30uvVq0MG/fvn34djlndIBOqusK8ilnnGL2XVtbK93d3cXctXT3OTrzro455hjZsWOH5CrQme9Z/qW53D6QVPkTyBBhSvNJonnpmN+1a1eGLewu0ifr6q7q4IWx3UyOjD5hwoRgnB04cODIlZaXpFIp6e3tdfKFy7Rp02TPnj1BfqOZtneYeX7VuHrwxXRYwHgBrqurG1FIdMA2NDQMZ6RPTPoXnnbv3h2+WbZ5PcLUItLf31+2GMXu2NW80u1RMxcLsKtuaSsXx1oikRAtwi7mFofHqO0zVenHZPhfl93Cefo2b/w9YH0lpke0OiD09LP+mz797Bs+7UUAAQQQ8FfA+BGwvr+6cOFCWbFihbS1tcnixYv91aflCCCAAALeChgvwCp99tlny5lnnil6qkv/mBBAAAEEEPBNwEoBVmROO/s21GgvAggggEBYgMPPsAbzCCCAAAIIGBKgABuCJgwCCCCAAAJhAQpwWIN5BBBAAAEEDAlQgA1BEwYBBBBAAIGwAAU4rME8AggggAAChgQowIagCYMAAggggEBYgAIc1mAeAQQQQAABQwIUYEPQhEEAAQQQQCAsQAEOazCPAAIIIICAIQEKsCFowiCAAAIIIBAWoACHNZhHAAEEEEDAkAAF2BA0YRBAAAEEEAgLUIDDGswjgAACCCBgSKBqcGgyFMv5MCtXrpR58+bJGWec4XyuLiV4xx13yI9+9CNpbGx0KS2nc/nPf/4j//jHP+S73/2u03m6ltxjjz0mDQ0Ncu6557qWmtP53HXXXbJ06VKZPn2603n6lhxHwKEe7+npkb6+vtASZqMIdHd3C6/jokgd3mZgYEB0vDEVJtDb28tjtDCyYGt9jOqYY3JLgALsVn+QDQIIIICAJwIUYE86mmYigAACCLglwHvAof7473//K+PHj5eJEyeGljKbT2D9+vXBe+fJZDLfpqz/WKC9vV127twpxx9/PCYFCOzYsUN0nB111FEF3ItN33nnHZk1a5aMGzcODIcEKMAOdQapIIAAAgj4I8ApaH/6mpYigAACCDgkQAF2qDNIBQEEEEDAHwHetPu4r5977jlZt26dNDU1yWWXXSZ1dXX+jIICW/r3v/9d9P3y5uZmOe+886SlpUX27Nkj+h3Njo4O+dKXviSnnXZagXv1Y3N9L+7111+Xyy+/PGjwxo0bZc2aNcFXay655BKZOnWqHxARW9nW1iYvvvhiMN5OPvlk+eIXvxhY/elPfwreQz/hhBPkm9/8ZsS9+bOZfi7jX//6l9TU1MjXvva14Pu/PEbd63+OgIf6ZMuWLbJ582a56aab5LjjjpN//vOf7vWUIxlt2LBBDhw4IMuXLw+K7JNPPhlk9sgjjwTF+Nprr5WnnnpKDh486EjG7qShJo8//rhoUdGpq6tLVq1aJVdddZVcdNFF8tBDD7mTrCOZPPzww8EH/K6//vrgMaofXnv22WeDFyo333xzUITffvttR7J1I43+/n7529/+Flx44+tf/3ow5jQzHqNu9E84CwrwkMamTZvkU5/6lFRXV8vChQtFr1LElFlAX6AsXrw4WKmvrltbW4N5LcozZ84MrlI0Z84c2bZtW+YdeLz0z3/+c3A0kibQT0GrWX19vRxzzDHS2dnJRSbSOEP/6oUj1EjH2dq1a4OzBnq1NX28LliwIHi8nn766TxeQ2Y6q89jemGc7du3B2cOamtrgy14jI6CcuAmBXioE/bv3x88CWp/6Mf09TQqU2YB9dE/NdOCoqf/9Mgu/BUkLSgYjvR79dVXZcaMGSMuBRged7o1Y2+kmZ4pSJ+C1iPfX/7yl6JHd2E3xtpIM72lxffUU08Njnj1cqef/vSneYweyeTEEgrwUDfoE59eqk0nvdSdfheYKbuAvpe0YsUKufjii2X27Nmir7DDl1XUea4LfdhPX6DoKUF9kfLaa6/JRx99JPpesH7OID3udGsde3qdY6ZDAqlUKigml156qSxatCj4zrS6hR+vjLUjR8v7778fHPnefvvt8sMf/lD++te/SlVVFY/RI6msL6EAD3WBnv7bunVr0Bn6fjAXLM8+LvWU829/+1u58sorZe7cucGGespLnyz1qFdffesHtPgw0WFD9bnwwguD4qpFN+01bdq04DShmunpZ/03fCbh8B78nNNCqx+K1PfKddKxp36jH696ZoHpsIA+FtVOp0QiEcynxxyP0cNOLsxxIY6hXtD3mvRTlXpqS98nWbZsWfDAd6GDXMtBPxTz5ptvBp+A1twmT54s+sGrd999V1avXh28h6mfVuXXajL33K5du4IjEh1jOuknfPUTq3qqVd9bP/HEEzPf0dOl7733XvBpXj0FPWHChOAbCvo4/ctf/hKcMdD3h6+55prgCM9ToozN1g/76QsWfWH3yU9+Us466yweoxml7C6kAIf89XSWvnpkKl5AT6PqkyJTdAH9BS49UtE/pswCmR6bmZZlvrefS/WxqEe+o8cVj1F3xgMF2J2+IBMEEEAAAY8EeMntUWfTVAQQQAABdwQowO70BZkggAACCHgkQAH2qLNpKgIIIICAOwIUYHf6gkw8ENCvuemnUvXTvTo98MADoteA1kuhZlquX01iQgCByhTgQ1iV2a+0ymGB733ve8FXQu67777gikV6PW29BGq25Q43hdQQQGAMAhTgMeBxVwSKEdCLIcyfPz/4rvn5558vd999d7CbbMuLicF9EEDAfQFOQbvfR2RYYQJ6uckbbrghuADHd77zneHWZVs+vAEzCCBQUQIcAVdUd9KYOAjolZz0amH6p5ej/P3vfx+knW15HNpEjgggULgAR8CFm3EPBMYkcNttt4n+TqteTvHpp5+W9G8qZ1s+pmDcGQEEnBVIOpsZiSFQgQL6Y/J6nV79Efnm5ubgJ/aWL18u9957b8blep1oflmqAgcCTUJgSIBT0AwDBBBAAAEELAhwCtoCOiERQAABBBCgADMGEEAAAQQQsCBAAbaATkgEEEAAAQQowIwBBBBAAAEELAhQgC2gExIBBBBAAAEKMGMAAQQQQAABCwIUYAvohEQAAQQQQIACzBhAAAEEEEDAggAF2AI6IRFAAAEEEPj/NyswL0xPlWAAAAAASUVORK5CYII="
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "py$ggplotly(plot, session=\"notebook\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "Loading required package: IRdisplay\n"
       ]
      },
      {
       "html": [
        "<iframe height=\"525\" id=\"igraph\" scrolling=\"no\" seamless=\"seamless\"\n",
        "\t\t\t\tsrc=\"https://plot.ly/~r_user_guide/989\" width=\"100%\" frameBorder=\"0\"></iframe>"
       ],
       "metadata": {},
       "output_type": "display_data"
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "df <- structure(c(106487, 495681, 1597442,\n",
      "     2452577, 2065141, 2271925, 4735484, 3555352,\n",
      "     8056040, 4321887, 2463194, 347566, 621147,\n",
      "     1325727, 1123492, 800368, 761550, 1359737,\n",
      "     1073726, 36, 53, 141, 41538, 64759, 124160,\n",
      "     69942, 74862, 323543, 247236, 112059, 16595,\n",
      "     37028, 153249, 427642, 1588178, 2738157,\n",
      "     2795672, 2265696, 11951, 33424, 62469,\n",
      "     74720, 166607, 404044, 426967, 38972, 361888,\n",
      "     1143671, 1516716, 160037, 354804, 996944,\n",
      "     1716374, 1982735, 3615225, 4486806, 3037122,\n",
      "     17, 54, 55, 210, 312, 358, 857, 350, 7368,\n",
      "     8443, 6286, 1750, 7367, 14092, 28954, 80779,\n",
      "     176893, 354939, 446792, 33333, 69911, 53144,\n",
      "     29169, 18005, 11704, 13363, 18028, 46547,\n",
      "     14574, 8954, 2483, 14693, 25467, 25215,\n",
      "     41254, 46237, 98263, 185986), .Dim = c(19,\n",
      "     5), .Dimnames = list(c(\"1820-30\", \"1831-40\",\n",
      "     \"1841-50\", \"1851-60\", \"1861-70\", \"1871-80\",\n",
      "     \"1881-90\", \"1891-00\", \"1901-10\", \"1911-20\",\n",
      "     \"1921-30\", \"1931-40\", \"1941-50\", \"1951-60\",\n",
      "     \"1961-70\", \"1971-80\", \"1981-90\", \"1991-00\",\n",
      "     \"2001-06\"), c(\"Europe\", \"Asia\", \"Americas\",\n",
      "     \"Africa\", \"Oceania\")))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 18
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "df.m <- melt(df)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 21
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "df.m <- rename(df.m, c(X1 = \"Period\", X2 = \"Region\"))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 22
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "a <- ggplot(df.m, aes(x = Period, y = value/1e+06,\n",
      "     fill = Region)) "
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 24
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "b <- a + geom_bar(stat = \"identity\", position = \"stack\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 25
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "py$ggplotly(b, session=\"notebook\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "html": [
        "<iframe height=\"525\" id=\"igraph\" scrolling=\"no\" seamless=\"seamless\"\n",
        "\t\t\t\tsrc=\"https://plot.ly/~r_user_guide/991\" width=\"100%\" frameBorder=\"0\"></iframe>"
       ],
       "metadata": {},
       "output_type": "display_data"
      }
     ],
     "prompt_number": 26
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "set.seed(1234)\n",
      "df <- data.frame(cond = factor( rep(c(\"A\",\"B\"), each=200) ), \n",
      "                   rating = c(rnorm(200),rnorm(200, mean=.8)))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 37
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "box <- ggplot(df, aes(x=cond, y=rating, fill=cond)) + geom_boxplot()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 38
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "py$ggplotly(box, session=\"notebook\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "html": [
        "<iframe height=\"525\" id=\"igraph\" scrolling=\"no\" seamless=\"seamless\"\n",
        "\t\t\t\tsrc=\"https://plot.ly/~r_user_guide/993\" width=\"100%\" frameBorder=\"0\"></iframe>"
       ],
       "metadata": {},
       "output_type": "display_data"
      }
     ],
     "prompt_number": 39
    }
   ],
   "metadata": {}
  }
 ]
}