{ "metadata": { "name": "", "signature": "sha256:1bf5ab6e8f96a5440ab804e9eea168fa81df084c4794dca5cd0ccdcb8930b2c2" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "#R" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Lets start with a typical R command. This will throw an exception." ] }, { "cell_type": "code", "collapsed": false, "input": [ "X=c(1,4,5,7);" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'c' is not defined", "output_type": "pyerr", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mX\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mc\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m4\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m5\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m7\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m;\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mNameError\u001b[0m: name 'c' is not defined" ] } ], "prompt_number": 1 }, { "cell_type": "markdown", "metadata": {}, "source": [ "To enable the magics below, execute %load_ext rpy2.ipython." ] }, { "cell_type": "code", "collapsed": false, "input": [ "%load_ext rpy2.ipython" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "%R X=c(1,4,5,7); sd(X); mean(X)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 3, "text": [ "\n", "[4.250000]" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "X" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'X' is not defined", "output_type": "pyerr", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mX\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mNameError\u001b[0m: name 'X' is not defined" ] } ], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Again X is only defined in R." ] }, { "cell_type": "code", "collapsed": false, "input": [ "%%R\n", "Y = c(2,4,3,9)\n", "summary(lm(Y~X))" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "text": [ "\n", "Call:\n", "lm(formula = Y ~ X)\n", "\n", "Residuals:\n", " 1 2 3 4 \n", " 0.88 -0.24 -2.28 1.64 \n", "\n", "Coefficients:\n", " Estimate Std. Error t value Pr(>|t|)\n", "(Intercept) 0.0800 2.3000 0.035 0.975\n", "X 1.0400 0.4822 2.157 0.164\n", "\n", "Residual standard error: 2.088 on 2 degrees of freedom\n", "Multiple R-squared: 0.6993,\tAdjusted R-squared: 0.549 \n", "F-statistic: 4.651 on 1 and 2 DF, p-value: 0.1638\n", "\n" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "%R plot(X, Y)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHgCAYAAAB91L6VAAAEJGlDQ1BJQ0MgUHJvZmlsZQAAOBGF\nVd9v21QUPolvUqQWPyBYR4eKxa9VU1u5GxqtxgZJk6XtShal6dgqJOQ6N4mpGwfb6baqT3uBNwb8\nAUDZAw9IPCENBmJ72fbAtElThyqqSUh76MQPISbtBVXhu3ZiJ1PEXPX6yznfOec7517bRD1fabWa\nGVWIlquunc8klZOnFpSeTYrSs9RLA9Sr6U4tkcvNEi7BFffO6+EdigjL7ZHu/k72I796i9zRiSJP\nwG4VHX0Z+AxRzNRrtksUvwf7+Gm3BtzzHPDTNgQCqwKXfZwSeNHHJz1OIT8JjtAq6xWtCLwGPLzY\nZi+3YV8DGMiT4VVuG7oiZpGzrZJhcs/hL49xtzH/Dy6bdfTsXYNY+5yluWO4D4neK/ZUvok/17X0\nHPBLsF+vuUlhfwX4j/rSfAJ4H1H0qZJ9dN7nR19frRTeBt4Fe9FwpwtN+2p1MXscGLHR9SXrmMgj\nONd1ZxKzpBeA71b4tNhj6JGoyFNp4GHgwUp9qplfmnFW5oTdy7NamcwCI49kv6fN5IAHgD+0rbyo\nBc3SOjczohbyS1drbq6pQdqumllRC/0ymTtej8gpbbuVwpQfyw66dqEZyxZKxtHpJn+tZnpnEdrY\nBbueF9qQn93S7HQGGHnYP7w6L+YGHNtd1FJitqPAR+hERCNOFi1i1alKO6RQnjKUxL1GNjwlMsiE\nhcPLYTEiT9ISbN15OY/jx4SMshe9LaJRpTvHr3C/ybFYP1PZAfwfYrPsMBtnE6SwN9ib7AhLwTrB\nDgUKcm06FSrTfSj187xPdVQWOk5Q8vxAfSiIUc7Z7xr6zY/+hpqwSyv0I0/QMTRb7RMgBxNodTfS\nPqdraz/sDjzKBrv4zu2+a2t0/HHzjd2Lbcc2sG7GtsL42K+xLfxtUgI7YHqKlqHK8HbCCXgjHT1c\nAdMlDetv4FnQ2lLasaOl6vmB0CMmwT/IPszSueHQqv6i/qluqF+oF9TfO2qEGTumJH0qfSv9KH0n\nfS/9TIp0Wboi/SRdlb6RLgU5u++9nyXYe69fYRPdil1o1WufNSdTTsp75BfllPy8/LI8G7AUuV8e\nk6fkvfDsCfbNDP0dvRh0CrNqTbV7LfEEGDQPJQadBtfGVMWEq3QWWdufk6ZSNsjG2PQjp3ZcnOWW\ning6noonSInvi0/Ex+IzAreevPhe+CawpgP1/pMTMDo64G0sTCXIM+KdOnFWRfQKdJvQzV1+Bt8O\nokmrdtY2yhVX2a+qrykJfMq4Ml3VR4cVzTQVz+UoNne4vcKLoyS+gyKO6EHe+75Fdt0Mbe5bRIf/\nwjvrVmhbqBN97RD1vxrahvBOfOYzoosH9bq94uejSOQGkVM6sN/7HelL4t10t9F4gPdVzydEOx83\nGv+uNxo7XyL/FtFl8z9ZAHF4bBsrEwAAK15JREFUeAHt3Qu0VXWdB/AfCAryFA2UsDJNUVJGGzUf\n+ZhKnYjGDF+jooMZOPkksnRaVOqSamg10+RUFtrQgxp8ML7WIJGOz6YYbYIIGHPNYCr5FkEBwab/\nWYu7vHk57JP7nLPP3p+91oFz9v7f//7/Pv997/fufR631+//sISFAAECBAgQaKlA75buzc4IECBA\ngACBmoAAdiAQIECAAIE2CAjgNqDbJQECBAgQEMCOAQIECBAg0AYBAdwGdLskQIAAAQIC2DFAgAAB\nAgTaICCA24BulwQIECBAQAA7BggQIECAQBsEBHAb0O2SAAECBAgIYMcAAQIECBBog4AAbgO6XRIg\nQIAAAQHsGCBAgAABAm0QEMBtQLdLAgQIECAggB0DBAgQIECgDQICuA3odkmAAAECBASwY4AAAQIE\nCLRBQAC3Ad0uCRAgQICAAHYMECBAgACBNggI4Dag2yUBAgQIEBDAjgECBAgQINAGAQHcBnS7JECA\nAAECAtgxQIAAAQIE2iAggNuAbpcECBAgQEAAOwYIECBAgEAbBARwG9DtkgABAgQICGDHAAECBAgQ\naIOAAG4Dul0SIECAAAEB7BggQIAAAQJtEBDAbUC3SwIECBAgIIAdAwQIECBAoA0CArgN6HZJgAAB\nAgQEsGOAAAECBAi0QUAAtwHdLgkQIECAgAB2DBAgQIAAgTYICOA2oNslAQIECBAQwI4BAgQIECDQ\nBgEB3AZ0uyRAgAABAgLYMUCAAAECBNogIIDbgG6XBAgQIEBAADsGCBAgQIBAGwQEcBvQ7ZIAAQIE\nCAhgxwABAgQIEGiDgABuA7pdEiBAgAABAewYIECAAAECbRAQwG1At0sCBAgQICCAHQMECBAgQKAN\nAgK4Deh2SYAAAQIEBLBjgAABAgQItEFAALcB3S4JECBAgIAAdgwQIECAAIE2CAjgNqDbJQECBAgQ\nEMCOAQIECBAg0AYBAdwGdLskQIAAAQIC2DFAgAABAgTaICCA24BulwQIECBAQAA7BggQIECAQBsE\nBHAb0O2SAAECBAj0qRLB9ddfHxs3bqxSyWolQIAAgToCw4cPj7/4i7+o06J5m3r9/g9L87ovTs83\n3HBDfPnLX44zzzyzOIMyEgIECBBoq8BXv/rV+P73vx9/9md/1vJxVOYMOJ35Tpw4MSZPntxyZDsk\nQIAAgWIKrFixIl599dW2DM5zwG1ht1MCBAgQqLqAAK76EaB+AgQIEGiLgABuC7udEiBAgEDVBQRw\n1Y8A9RMgQIBAWwQEcFvY7ZQAAQIEqi4ggKt+BKifAAECBNoiIIDbwm6nBAgQIFB1gcq8D7jqE61+\nAgQIlFVg6dKl8cILL8Ree+0Vw4YN65gyC3cGvHbt2njiiSeiIh/Q1TEHioESIECgaAIpJ6ZPnx5X\nXHFFpE87HDVqVDzwwANFG+YWx1OIAN60aVNccsklNbyBAwfGbrvtFjvssEPst99+ccstt2xx8DYQ\nIECAQHUFrrzyykgnbXPmzImZM2fGgw8+GBdeeGH83//9X0egFOIS9KWXXhpPPfVULFiwIJ577rn4\n3Oc+V/uNpnfv3jXMdevWxYknnrhV0Lvvvjt++tOf9tjunnvuifSh21OmTOlxu5UECBAg0FkCixYt\nim9/+9tdgx49enRMmjQp7rvvvnjrW9/atb6odwoRwLfeemvtssGQIUNqThdccEFcffXVMXv27Jg1\na1Z85jOfyRTAO++8c+y77749Wt9///21kO9xo5UECBAg0HEC6YppOkF77fL888/Hjjvu+NpVhb1f\niAAeM2ZM3HnnnXH88cfXoNJZ7K677lq7ny4lbA7mrSnuueeekW49LXfccUesWrWqp03WESBAgEAH\nCnz4wx+uPX157bXXRv/+/WPevHmRrqiuXr26I6opRABfdNFFMX78+PjHf/zHWL9+fTz88MO1a/m3\n3XZb7ZJx+ju+FgIECBAg8FqBCRMm1K5sHnDAAbWTr3RGnE7aBg0a9Npmhb1fiAA+7LDD4vHHH4+b\nb745tt9++zj22GOjb9++tVv6U1HpNxsLAQIECBD4Y4Fzzz030q0Tl0IEcILr169fnHTSSd0MR4wY\n0e2xBwQIECBAoCwChXgbUlkw1UGAAAECBLIKCOCsUtoRIECAAIEcBQRwjpi6IkCAAAECWQUEcFYp\n7QgQIECAQI4CAjhHTF0RIECAAIGsAgI4q5R2BAgQIEAgRwEBnCOmrggQIECAQFYBAZxVSjsCBAgQ\nIJCjgADOEVNXBAgQIEAgq4AAziqlHQECBAgQyFFAAOeIqSsCBAgQIJBVQABnldKOAAECBAjkKCCA\nc8TUFQECBAgQyCoggLNKaUeAAAECBHIUEMA5YuqKAAECBAhkFRDAWaW0I0CAAAECOQoI4BwxdUWA\nAAECBLIKCOCsUtoRIECAAIEcBQRwjpi6IkCAAAECWQUEcFYp7QgQIECAQI4CAjhHTF0RIECAAIGs\nAgI4q5R2BAgQIEAgRwEBnCOmrggQIECAQFYBAZxVSjsCBAgQIJCjgADOEVNXBAgQIEAgq4AAziql\nHQECBAgQyFFAAOeIqSsCBAgQIJBVQABnldKOAAECBAjkKCCAc8TUFQECBAgQyCoggLNKaUeAAAEC\nBHIUEMA5YuqKAAECBAhkFRDAWaW0I0CAAAECOQoI4BwxdUWAAAECBLIKCOCsUtoRIECAAIEcBQRw\njpi6IkCAAAECWQUEcFYp7QgQIECAQI4CAjhHTF0RIECAAIGsAgI4q5R2BAgQIEAgR4HCBvCrr74a\na9asybFUXREgQIAAgeIIFCKAV69eHTNnzozx48fHnXfeGTfffHOMGDEiRo4cGR/72McEcXGOFyMh\nQIAAgZwE+uTUzxvq5gtf+EI8/PDD8YEPfCAuvPDC2LhxYy2E99prr5g6dWrMnTs3/uZv/mar+1i4\ncGHcd999Pbb76U9/GoMGDepxm5UECBAgQKDVAoUI4H/7t3+Ln/3sZzFgwID43e9+F08//XQccsgh\nNYtLL720FsJZAni33XaL3r17Pqn/5S9/Ga+88kqrfe2PAAECBAj0KFCIAN57771jwYIFcfTRR8fd\nd98dL7/8ctdgU3AecMABXY/r3Xn7298e6dbTki5rr1q1qqdN1hEgQIAAgZYLFCKAP/GJT8SkSZPi\nkUceiQsuuCBefPHFSKE8duzYuPfee+Ouu+5qOYwdEiBAgACBZgoUIoDT5ealS5fGs88+GzvuuGOs\nX78+5s+fH88//3xcd9110b9//2Ya6JsAAQIECLRcoBABnKru1atXLXzT/e222y4+9KEPpbsWAgQI\nECBQSoGeX7FUylIVRYAAAQIEiiMggIszF0ZCgAABAhUSEMAVmmylEiBAgEBxBARwcebCSAgQIECg\nQgICuEKTrVQCBAgQKI6AAC7OXBgJAQIECFRIQABXaLKVSoAAAQLFERDAxZkLIyFAgACBCgkI4ApN\ntlIJECBAoDgCArg4c2EkBAgQIFAhAQFcoclWKgECBAgUR0AAF2cujIQAAQIEKiQggCs02UolQIAA\ngeIICODizIWRECBAgECFBARwhSZbqQQIECBQHAEBXJy5MBICBAgQqJCAAK7QZCuVAAECBIojIICL\nMxdGQoAAAQIVEhDAFZpspRIgQIBAcQQEcHHmwkgIECBAoEICArhCk61UAgQIECiOgAAuzlwYCQEC\nBAhUSEAAV2iylUqAAAECxREQwMWZCyMhQIAAgQoJCOAKTbZSCRAgQKA4AgK4OHNhJAQIECBQIQEB\nXKHJVioBAgQIFEdAABdnLoyEAAECBCokIIArNNlKJUCAAIHiCAjg4syFkRAgQIBAhQQEcIUmW6kE\nCBAgUBwBAVycuTASAgQIEKiQgACu0GQrlQABAgSKIyCAizMXRkKAAAECFRIQwBWabKUSIECAQHEE\nBHBx5sJICBAgQKBCAgK4QpOtVAIECBAojoAALs5cGAkBAgQIVEigEAH8+OOPx6ZNmyrErlQCBAgQ\nqLpAIQJ43LhxccQRR8QjjzxS9flQPwECBAhURKBPUercd99949BDD41LLrkkzj777BgyZEjDQ/v3\nf//3+I//+I8ev+7ee++NoUOH9rjNSgIECBAg0GqBwgTwOeecE9OmTYtzzz03PvvZz8Zpp50Wp59+\neuy3334xePDgTC5jxozZYsims+v169dn6kcjAgQIECDQbIHCBHAqdI899ogFCxbEsmXL4tvf/nac\nccYZsXLlypg0aVJ861vf2qrFrrvuGunW0zJy5MhYtWpVT5usI0CAAAECLRcoVABvrn706NExc+bM\n2m3t2rXxzDPPbN7kfwIECBAgUAqBQgRwet53S2euAwYMiHSzECBAgACBMgkUIoBPPfXUMpmqhQAB\nAgQIbFWgEG9D2uooNSBAgAABAiUTEMAlm1DlECBAgEBnCAjgzpgnoyRAgACBkgkI4JJNqHIIECBA\noDMEBHBnzJNREiBAgEDJBARwySZUOQQIECDQGQICuDPmySgJECBAoGQCArhkE6ocAgQIEOgMAQHc\nGfNklAQIECBQMgEBXLIJVQ4BAgQIdIaAAO6MeTJKAgQIECiZgAAu2YQqhwABAgQ6Q0AAd8Y8GSUB\nAgQIlExAAJdsQpVDgAABAp0hIIA7Y56MkgABAgRKJiCASzahyiFAgACBzhAQwJ0xT0ZJgAABAiUT\nEMAlm1DlECBAgEBnCAjgzpgnoyRAgACBkgkI4JJNqHIIECBAoDMEBHBnzJNREiBAgEDJBARwySZU\nOQQIECDQGQICuDPmySgJECBAoGQCArhkE6ocAgQIEOgMAQHcGfNklAQIECBQMgEBXLIJVQ4BAgQI\ndIaAAO6MeTJKAgQIECiZgAAu2YQqhwABAgQ6Q0AAd8Y8GSUBAgQIlExAAJdsQpVDgAABAp0hIIA7\nY56MkgABAgRKJiCASzahyiFAgACBzhAQwJ0xT0ZJgAABAiUTEMAlm1DlECBAgEBnCAjgzpgnoyRA\ngACBkgkI4JJNqHIIECBAoDMEBHBnzJNREiBAgEDJBARwySZUOQQIECDQGQKFDeB169bFK6+80hmK\nRkmAAAECBBoUKEQAr1y5MiZOnBiLFi2Kp556Ks4+++zYeeedY+jQoTFp0qTYsGFDg2VpToAAAQIE\nii3QpwjDmz59erzlLW+JMWPGxIwZM2Ljxo2xZMmSWL9+fXz605+OK664onbb2lhvvfXWWLhwYY/N\n7rrrrthxxx173GYlAQIECBBotUAhAvjuu++OZcuWxbbbbhs33XRTzJs3L0aNGlWzSOE7ZcqUTC7v\nete7ur7uj7/gySefjJdffvmPV3tMgAABAgTaIlCIAN5zzz1j9uzZ8dGPfjSOOuqouP322+P888+v\ngaSz2ne84x2ZcHbZZZdIt56W4cOHx6pVq3raZB0BAgQIEGi5QCEC+Oqrr44PfvCDMWvWrNhjjz1i\n2rRpce2110bv3r1j9erVkc6QLQQIECBAoEwChQjg3XffPZYuXRoLFiyI5cuX154P3mGHHWpnvuPG\njYs+fQoxzDLNu1oIECBAoM0ChUm2Xr16xTHHHFO7tdnE7gkQIECAQNMFCvE2pKZXaQcECBAgQKBg\nAgK4YBNiOAQIECBQDQEBXI15ViUBAgQIFExAABdsQgyHAAECBKohIICrMc+qJECAAIGCCQjggk2I\n4RAgQIBANQQEcDXmWZUECBAgUDABAVywCTEcAgQIEKiGgACuxjyrkgABAgQKJiCACzYhhkOAAAEC\n1RAQwNWYZ1USIECAQMEEBHDBJsRwCBAgQKAaAgK4GvOsSgIECBAomIAALtiEGA4BAgQIVEPgdQH8\nyU9+Ml588cVqVK9KAgQIECDQJoHXBfDKlStjv/32i3vuuadNQ7JbAgQIECBQfoHXBfCPfvSjuOqq\nq2LChAlxySWXxIYNG8qvoEICBAgQINBigT497e/UU0+N973vfTFt2rQ48MAD46STTupqtvfee8cJ\nJ5zQ9dgdAgQIECBAoHGBHgM4ddOrV6/o27dvrFq1KpYsWdLV88CBA7vuu0OAAAECBAj8aQI9BvCc\nOXPiggsuiCOPPDIWL14cw4cP/9N691UECBAgQIBAjwKvC+CTTz45Fi5cGF/72tfilFNO6fGLrCRA\ngAABAgTemMDrAnjIkCHxq1/9KkaMGPHGevbVBAgQIECAwBYFXhfA11xzzRYb20CAAAECBAjkI/C6\ntyHl061eCBAgQIAAgXoCAriejm0ECBAgQKBJAgK4SbC6JUCAAAEC9QQEcD0d2wgQIECAQJMEBHCT\nYHVLgAABAgTqCQjgejq2ESBAgACBJgkI4CbB6pYAAQIECNQTEMD1dGwjQIAAAQJNEhDATYLVLQEC\nBAgQqCcggOvp2EaAAAECBJokIICbBKtbAgQIECBQT0AA19OxjQABAgQINElAADcJVrcECBAgQKCe\ngACup2MbAQIECBBokoAAbhKsbgkQIECAQD2BwgbwU089FRs3bqw3dtsIECBAgEDHChQigCdOnBjL\nli2rIS5fvjzGjRsXu+66a+y8885x3nnnxSuvvNKxwAZOgAABAgR6EujT08pWr1uyZEmsXbu2ttsZ\nM2bE6NGj47vf/W48/fTTMXXq1Ejrpk+fvtVhzZs3L+bPn99ju/vuuy/e9KY39bjNSgIECBAg0GqB\nQgTwa4tOAbpixYoYNGhQDBs2LK688spaCGcJ4MMOO6wW3q/tb/P9l156KdasWbP5of8JECBAgEBb\nBQoTwPfff3+MHDky3v3ud8czzzxTC+Aks3jx4th///0zIaUz3C2d5aYw37BhQ6Z+NCJAgAABAs0W\nKEQAn3baaXHLLbfEFVdcES+88EL069cv5syZE5/73Ofi6quvjoULFzbbQf8ECBAgQKClAoUI4E98\n4hORbml57LHHYvXq1bX7xx13XEybNi0GDhxYe+wfAgQIECBQFoFCBPBrMd/85jdHuqUlXY62ECBA\ngACBMgoU4m1IZYRVEwECBAgQqCcggOvp2EaAAAECBJokIICbBKtbAgQIECBQT0AA19OxjQABAgQI\nNElAADcJVrcECBAgQKCegACup2MbAQIECBBokoAAbhKsbgkQIECAQD0BAVxPxzYCBAgQINAkAQHc\nJFjdEiBAgACBegICuJ6ObQQIECBAoEkCArhJsLolQIAAAQL1BARwPR3bCBAgQIBAkwQEcJNgdUuA\nAAECBOoJCOB6OrYRIECAAIEmCQjgJsHqlgABAgQI1BMQwPV0bCNAgAABAk0SEMBNgtUtAQIECBCo\nJyCA6+nYRoAAAQIEmiQggJsEq1sCBAgQIFBPQADX07GNAAECBAg0SUAANwlWtwQIECBAoJ6AAK6n\nYxsBAgQIEGiSgABuEqxuCRAgQIBAPQEBXE/HNgIECBAg0CQBAdwkWN0SIECAAIF6AgK4no5tBAgQ\nIECgSQICuEmwuiVAgAABAvUEBHA9HdsIECBAgECTBARwk2B1S4AAAQIE6gkI4Ho6thEgQIAAgSYJ\nCOAmweqWAAECBAjUExDA9XRsI0CAAAECTRIQwE2C1S0BAgQIEKgnIIDr6dhGgAABAgSaJCCAmwSr\nWwIECBAgUE9AANfTsY0AAQIECDRJQAA3CVa3BAgQIECgnkCfehttI0CgfAILFiyIZ599Nnbffff4\n8z//8/IVqCICHSJQ2DPgdevWxerVqzuE0TAJdIbAlClTYu7cubUAnjBhQsyYMaMzBm6UBEooUNgA\nvuGGG2Lq1KklJFcSgfYIfOMb34ihQ4fGNddcE+eee26sWLEi7r777rjzzjvbMyB7JVBxgUJcgn7H\nO94RTz/9dLep2LBhQ2zcuDFSEB9//PFx3XXXddve04Prr78+brvttp42xc9+9rMYMWJEj9usJFAF\ngfQ98KlPfaqr1G233TbSGXFaf/TRR3etd4cAgdYIFCKAU7hOmjQpTj/99DjzzDNrlc+bNy8eeOCB\n+OIXvxgDBgzIpPHe9753i89pXX755S5pZ1LUqKwCAwcOjMceeyz22muvrhKXLVuW+fur64vcIUAg\nF4FCBPDhhx8eixYtivPOO6922fmb3/xm7LTTTpF+YLz1rW/NXOgOO+wQ6dbTMmTIkHj55Zd72mQd\ngUoITJ48uXbpediwYTF69Oi48cYb49Of/nSsWbOmEvUrkkDRBAoRwAll8ODBMXv27PjXf/3XOOKI\nI+Lggw+ObbbZpmhexkOgYwXGjBkT6Xngiy++OHr37h3pqZ8nnnjCGXDHzqiBd7pAYQJ4M+RJJ50U\nhx56aO039bFjx25e7X8CBHIQ2GeffWL+/Pk59KQLAgTeqEDhAjgVNGrUqLjlllveaG2+ngABAgQI\nFFagsG9DKqyYgREgQIAAgRwEBHAOiLogQIAAAQKNCgjgRsW0J0CAAAECOQgI4BwQdUGAAAECBBoV\nEMCNimlPgAABAgRyEBDAOSDqggABAgQINCoggBsV054AAQIECOQgIIBzQNQFAQIECBBoVEAANyqm\nPQECBAgQyEFAAOeAqAsCBAgQINCogABuVEx7AgQIECCQg4AAzgFRFwQIECBAoFEBAdyomPYECBAg\nQCAHAQGcA6IuCBAgQIBAowICuFEx7QkQIECAQA4CAjgHRF0QIECAAIFGBQRwo2LaEyBAgACBHAQE\ncA6IuiBAgAABAo0KCOBGxbQnQIAAAQI5CAjgHBB1QYAAAQIEGhUQwI2KaU+AAAECBHIQEMA5IOqC\nAAECBAg0KiCAGxXTngABAgQI5CAggHNA1AUBAgQIEGhUQAA3KqY9AQIECBDIQUAA54CoCwIECBAg\n0KiAAG5UTHsCBAgQIJCDgADOAVEXBAgQIECgUQEB3KiY9gQIECBAIAcBAZwDoi4IECBAgECjAgK4\nUTHtCRAgQIBADgICOAdEXRAgQIAAgUYFBHCjYtoTIECAAIEcBARwDoi6IECAAAECjQoI4EbFtCdA\ngAABAjkICOAcEHVBgAABAgQaFRDAjYppT4AAAQIEchDok0MfuiBAgEBhBR566KFYvHhxDBw4MP7q\nr/4qttlmm8KO1cCqJVDYM+BNmzbF6tWrqzUbqiVAIFeBa665Ji677LJYv3593HbbbfHOd74zXnrp\npVz3oTMCf6pAIQL4lVdeiRkzZsSkSZPiwQcfjB/+8IcxYsSIGDp0aJxwwgm1b54/tUBfR4BANQXS\nz5JPfvKTMXfu3DjnnHNi1qxZceKJJ8bMmTOrCaLqwgkU4hJ0+ib59a9/HQcccECcfPLJ0adPn7jh\nhhti1KhRcfHFF8e8efNq67emN2fOnLj55pt7bJa+GXfZZZcet1lJgED5BNKl53/4h3+oXXreXN1F\nF10UZ5111uaH/ifQVoFCBPDtt98eixYtisGDB0f//v3jySefjCOPPLIGc+WVV8ZnPvOZTAE8fvz4\nOProo3sEnT59ejz33HM9brOSAIHyCaTnfJcsWdKtsEcffdRzwN1EPGinQCEC+O1vf3ssW7YsDjro\noPjoRz8av/3tb7tM0osn9thjj67H9e6kb7h062kZMGBAvPjiiz1tso4AgRIKHH/88bWns6699tr4\nyEc+Ev/7v/8b6Qw4nRVbCBRBoBABPHXq1NqrE7/xjW/U/h85cmTNJr14In3z/PjHPy6ClTEQINBB\nAtttt12kp6UuvPDCuP7662PIkCFx1VVXxdixYzuoCkMts0AhAviYY46J5cuXx9q1a7tZp0vK6fLz\n9ttv3229BwQIEMgi0K9fv/jmN7+Zpak2BFouUIgATlWn53/T7bXLIYcc8tqH7hMgQIAAgdIIFOJt\nSKXRVAgBAgQIEMgoIIAzQmlGgAABAgTyFBDAeWrqiwABAgQIZBQQwBmhNCNAgAABAnkKCOA8NfVF\ngAABAgQyCgjgjFCaESBAgACBPAUEcJ6a+iJAgAABAhkFBHBGKM0IECBAgECeAgI4T019ESBAgACB\njAICOCOUZgQIECBAIE8BAZynpr4IECBAgEBGAQGcEUozAgQIECCQp4AAzlNTXwQIECBAIKOAAM4I\npRkBAgQIEMhTQADnqakvAgQIECCQUUAAZ4TSjAABAgQI5CkggPPU1BcBAgQIEMgoIIAzQmlGgAAB\nAgTyFBDAeWrqiwABAgQIZBQQwBmhNCNAgAABAnkKCOA8NfVFgAABAgQyCgjgjFCaESBAgACBPAUE\ncJ6a+iJAgAABAhkFBHBGKM0IECBAgECeAgI4T019ESBAgACBjAICOCOUZgQIECBAIE8BAZynpr4I\nECBAgEBGAQGcEUozAgQIECCQp4AAzlNTXwQIECBAIKOAAM4IpRkBAgQIEMhTQADnqakvAgQIECCQ\nUUAAZ4TSjAABAgQI5CkggPPU1BcBAgQIEMgoIIAzQmlGgAABAgTyFBDAeWrqiwABAgQIZBQQwBmh\nNCNAgAABAnkKCOA8NfVFgAABAgQyCvTJ2K6yzebPnx+PP/547LLLLnHcccdV1kHhBAgQIJCvQOHP\ngDdt2hTr16/Pt+qMvZ155plx/fXXR9++feNLX/pSnHjiifHqq69m/GrNCBAgQIDAlgUKEcCPPvpo\nTJw4MQYOHBjvf//74+GHH+4a8dy5c+OMM87oetyqO3PmzImlS5fGt771rTj99NPjJz/5SQwePDi+\n973vtWoI9kOAAAECJRYoxCXor3zlK7VLvIsWLYof/OAHccQRR8Rdd90Ve+65Z0P03/3ud+PGG2/s\n8Wt++ctfxqhRo3rc1tPKX/ziF5HG9dpl8uTJ8aMf/ei1q9wnQIAAAQJ/kkAhAvj222+Phx56KPr3\n7x+XX3557LPPPnHsscfGvffe21BR6RLx+PHje/yadCa9du3aHrf1tDKd7S5fvjwOP/zwrs1pjGm9\nhQABAgQIvFGBQgRwCtx09vue97ynVs8pp5xSe+HTX/7lX0Y668y69OvXL9KtpyUFZ3o+Oetyzjnn\nxIQJE2LnnXeOgw46KO65556YMmVKPP/881m70I4AAQIECGxRoBDPAadgS2evX/ziF7sGOnXq1PjI\nRz4SF198cde6Vt4ZPnx4zJs3r3ZJ+7zzzos777wzVq5cGUOGDGnlMOyLAAECBEoqUIgz4GOOOSZ+\n85vfxCOPPNKN+bOf/WwceeSRtW3dNrTowbBhw2LWrFkt2pvdECBAgECVBAoRwAl8wIABse+++77O\n/qijjop0sxAgQIAAgTIJFOISdJlA1UKAAAECBLIICOAsStoQIECAAIGcBQRwzqC6I0CAAAECWQQE\ncBYlbQgQIECAQM4CAjhnUN0RIECAAIEsAgI4i5I2BAgQIEAgZwEBnDOo7ggQIECAQBaBXr//w5Kl\nYae3SX9cYdy4cbH//vs3XEr6FKwtfcRlw511wBe8/PLLse2228Y222zTAaN940PcuHFj7WNKt9tu\nuzfeWYf0kD4XPb33virLSy+9VPus+V69elWi5A0bNtTqTN/HVVg2x9ihhx7acLnpA6AWLFgQb37z\nmxv+2jf6BZUJ4DcClT4IJP11pqosF154YZx99tmx3377VaLk9M3385//PC677LJK1JuKrNoxnT7W\nNv1p0fTpdlVYrrvuutov0OnPvFZhefLJJ+P888/vuL9W5xJ0FY5ONRIgQIBA4QQEcOGmxIAIECBA\noAoCArgKs6xGAgQIECicgAAu3JQYEAECBAhUQUAAV2GW1UiAAAEChRMQwIWbEgMiQIAAgSoIeBtS\nhll+4oknYpdddsnQshxNnnnmmRg0aFDtvcDlqKh+Fel9z+vXr4+hQ4fWb1iirVU7pn/3u9/Fm970\npujduxrnHGvWrKkdrQMHDizRUbvlUl599dV4+umnY/jw4VtuVMAtAriAk2JIBAgQIFB+gWr8Olj+\neVQhAQIECHSYgADusAkzXAIECBAoh4AALsc8qoIAAQIEOkxAAHfYhBkuAQIECJRDQACXYx5VQYAA\nAQIdJiCAO2zCDJcAAQIEyiEggMsxj6ogQIAAgQ4TEMAZJiz9sedNmzZlaKkJAQIECBDIJiCAt+KU\nPmHlpJNOir//+7/fSsvO37x06dI49dRTY+zYsfHe97634/64daMzkH6puvjii2P//fePffbZJ77y\nla802kXHtr/jjjsq8cfpFy1aFG95y1u63R577LGOnbetDfy5556LU045pXY8H3TQQfG9731va1/S\n0dtPOOGE2HXXXbvdxo8f3zE19emYkbZhoP/1X/8VF110UfzqV7+Kd73rXW0YQWt3mWo9/fTTY86c\nOfH444/Xgumoo46KESNGtHYgLdrb1VdfHStXrow0z+njKPfcc8845JBD4t3vfneLRtCe3aQf0tOm\nTYt0ZafsSwrg97///fFP//RPXaX279+/637Z7lx22WWx++67xw9/+MN48skn44gjjohjjz229jGc\nZas11fODH/wg0klSWtatWxcHH3xwnHbaabXHnfCPM+A6s/Qv//IvccEFF9TOCus0K8WmdBD/7d/+\nbVetI0eOrH0e9IMPPliK+noqYvLkyfGd73yn9vnA6TOCN2zYUImnGs4///zamX+vXr16YinVul/8\n4he1H8opjJ566qnYfvvto6x1pys61113XXzqU5+q/UKZPtt82bJlpQ3fdKD269evNqdpXq+66qo4\n7LDDalcAOuUgFsB1ZuqrX/1qnHjiiXValGdTCqDjjz8++vbtWytq4cKFkc6U0hlhWZftttuu9ktG\nOlt4z3veE+PGjSt1vWke586dW/uhlZ5iqMKSAnjmzJlxzDHHxNve9rZaOJW17lWrVtWO5y996Uu1\n0B0yZEhcc801ZS23W10rVqyIr3/96x33NJIA7jaNHiSBdDCfccYZ8bWvfa0SfyEonRGl54DT5cpf\n//rXpT0I0g/oyy+/vBZIpS3yjwpLTx3NmjWrdkynqznpUnQ6Ey7jkup69tln49FHH43f/va3tV+2\n0tlw+ktfZV/Sc90f+tCHYocdduioUgVwR01X8webLlml532nT5/edTm6+Xtt7x5OPvnkmD17du3y\nVfqlo6zLxz/+8Tj88MPj3nvvjR//+Me1S+633nprqX9Ap+f509WNtKQX26VLlDfeeGMppzhdck5P\nJaXngdP9D37wg7HXXntFesFd2Zfvf//7cdZZZ3VcmQK446aseQN+5JFH4n3ve1/83d/9XUyZMqV5\nOypIz9dee23853/+Z9doxowZE//zP//T9bhsd7bddtv47//+79pzZSmY0otW0vNma9euLVuptXpS\nfZ///OdrdW4u8KWXXirtc6Lpb5Zvs802tcvQm+tNc55qLvOyePHieOGFF2ovtuu0OgVwp81YE8eb\nLjunVxCmtyKlS1npll6YVNYl1ZbO9Ddu3BjPP/987Sw4nf2XdUmvbr///vtrt5tuuikGDx5cuz9s\n2LBSlpxeoPOTn/ykdgk6FZh+2XrooYdqzweXseD0mob0tpz0y1Va0tsK09Mqhx56aBnL7aopPW2U\nfnlOr2PptKXzRtxpwh0y3p///Oe1H8bpBRw77rhj1y390C7rMnHixNhpp51i7733jgMPPLD23uf0\n9hxLeQRmzJhRe6tKeovZcccdVwvjgQMHlqfAP6rky1/+cu3phVRvegtSemFSep9smZfly5fH6NGj\nO7LEXn94L2D53wzYkVNj0K0SSJdg09lDnz7eFt8q81bvJ13NSc+LduJZ0p9i9cwzz9TqTZekLcUV\nEMDFnRsjI0CAAIESC7gEXeLJVRoBAgQIFFdAABd3boyMAAECBEosIIBLPLlKI0CAAIHiCgjg4s6N\nkREgQIBAiQUEcIknV2kECBAgUFwBAVzcuTEyAgQIECixgAAu8eQqjQABAgSKKyCAizs3RkaAAAEC\nJRYQwCWeXKURIECAQHEFBHBx58bICBAgQKDEAgK4xJOrNAIECBAoroAALu7cGBkBAgQIlFhAAJd4\ncpVGgAABAsUVEMDFnRsjI0CAAIESCwjgEk+u0ggQIECguAICuLhzY2QECBAgUGIBAVziyVUaAQIE\nCBRXQAAXd26MjAABAgRKLCCASzy5SiOwWWDdunXxzne+My677LLNq2r/n3XWWfHXf/3X3dZ5QIBA\nawT6tGY39kKAQDsF+vXrF3PmzImDDz44DjzwwPjwhz8cX/jCF+KBBx6IRYsWtXNo9k2gsgICuLJT\nr/CqCey777610D3nnHOib9++ccUVV8T9998fgwYNqhqFegkUQqDX7/+wFGIkBkGAQEsEPvCBD8Qd\nd9wR//zP/xwf+9jHWrJPOyFA4PUCngN+vYk1BEotsMcee8SmTZtip512KnWdiiNQdAEBXPQZMj4C\nOQrcddddMXv27Lj88svj4x//eDz33HM59q4rAgQaEXAJuhEtbQl0sMCaNWsiPQ986aWXRnoe+Kij\njorddtstvvOd73RwVYZOoHMFBHDnzp2RE2hIYPLkyfGb3/wmFixYEL169YoVK1bE2LFj46abborj\njjuuob40JkDgjQsI4DduqAcChReYP39+TJgwIRYvXhxve9vbusab3or09a9/PZYsWeLV0F0q7hBo\njYAAbo2zvRAgQIAAgW4CXoTVjcMDAgQIECDQGgEB3BpneyFAgAABAt0EBHA3Dg8IECBAgEBrBARw\na5zthQABAgQIdBMQwN04PCBAgAABAq0REMCtcbYXAgQIECDQTUAAd+PwgAABAgQItEZAALfG2V4I\nECBAgEA3AQHcjcMDAgQIECDQGgEB3BpneyFAgAABAt0EBHA3Dg8IECBAgEBrBARwa5zthQABAgQI\ndBMQwN04PCBAgAABAq0REMCtcbYXAgQIECDQTUAAd+PwgAABAgQItEZAALfG2V4IECBAgEA3AQHc\njcMDAgQIECDQGoH/B7vCmYiHRnyRAAAAAElFTkSuQmCC\n" } ], "prompt_number": 7 }, { "cell_type": "markdown", "metadata": {}, "source": [ "The plot above is a R plot." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "From ipython to R." ] }, { "cell_type": "code", "collapsed": false, "input": [ "Z = [3, 1, 8, 2]\n", "%R -i Z mean(Z)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 24, "text": [ "\n", "[3.500000]" ] } ], "prompt_number": 24 }, { "cell_type": "markdown", "metadata": {}, "source": [ "from R to ipython" ] }, { "cell_type": "code", "collapsed": false, "input": [ "%R -o Zm Zm=mean(Z)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 28, "text": [ "\n", "[3.500000]" ] } ], "prompt_number": 28 }, { "cell_type": "code", "collapsed": false, "input": [ "Zm[0]" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 30, "text": [ "3.5" ] } ], "prompt_number": 30 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }