{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "### Reanalysis of Chen & Wyble, 2015\n", "\n", "Chen and Wyble published an interesting paper (2015) where they demonstrate that participants cannot report attributes of attended stimuli unless the participants are previously informed that this attribute is important. For instance, you wouldn't remember the color of the apple if you had had just told someone the shape. I would have expected the opposite, so ... cool!\n", "\n", "After reading the paper (you can check it out at http://wyblelab.com/publications), I became curious whether participants might unconsciously retain some information about these forgotten attributes. Chen and Wyble posted their data to databrary.com (https://nyu.databrary.org/volume/79), so I downloaded the data and did some quick analyses that you see here! I want to commend Chen and Wyble for sharing their data. This is something everyone should start doing (including me). \n", "\n", "Below, I will start by showing I can replicate Chen and Wyble's analyses, then I will investigate whether there's a trace of unconscious memory for the \"forgotten\" features. \n", "\n", "\n", "##### References\n", "Chen, H., & Wyble, B. (2015). Amnesia for Object Attributes Failure to Report Attended Information That Had Just Reached Conscious Awareness. Psychological science, 26(2),203-210.\n", "\n", "Wyble, B. (2014). Amnesia for object attributes: Failure to report attended information that had just reached conscious awareness. Databrary. Retrieved November 22, 2015 from http://doi.org/10.17910/B7G010" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Load relevant libraries and write analysis functions" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "I'll start by loading the python libraries that I'll use throughout analyses. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy as np, sys, scipy.stats, pandas as pd, os, os.path, csv #loading useful libraries\n", "\n", "import matplotlib as mpl\n", "import matplotlib.pyplot as plt\n", "import pylab as pl\n", "%matplotlib inline\n", "pd.options.display.mpl_style = 'default' #load matplotlib for plotting\n", "plt.style.use('ggplot') #im addicted to ggplot. so pretty.\n", "mpl.rcParams['font.family'] = ['Bitstream Vera Sans']" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here are some quick functions I wrote for running different statistical tests and plotting the data. I won't explain this code, but encourage you to look through it later if you're wondering how I did any of the analyses. " ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def print_tests(series1, series2): #this function just presents normality and t-tests. \n", " import scipy, numpy as np #the function does these tests with arcsin(sqrt(acc)) to help with normality\n", " normTest2 = scipy.stats.shapiro(np.arcsin(np.sqrt(series1))-np.arcsin(np.sqrt(series2)))\n", " Test2 = scipy.stats.ttest_rel(np.arcsin(np.sqrt(series1)), np.arcsin(np.sqrt(series2)))\n", " Test3 = scipy.stats.wilcoxon(np.arcsin(np.sqrt(series1)), y=np.arcsin(np.sqrt(series2)))\n", " print '\\t normality test adj. Test value: %s P-value: %s' % (str(np.round(normTest2[0],2)), \n", " str(np.round(normTest2[1],4)))\n", " if normTest2[1] > 0.1: print '\\t T-test adj. Test value: %s P-value: %s' % (str(np.round(Test2[0],2)), \n", " str(np.round(Test2[1],4)))\n", " if normTest2[1] <= 0.1: print '\\t Wilcoxon. Test value: %s P-value: %s' % (str(np.round(Test3[0],2)),\n", " str(np.round(Test3[1],2)))\n", " \n", "def print_tests_ChiSq(series): #this function just presents normality and t-tests. \n", " import scipy, numpy as np #the function does these tests with arcsin(sqrt(acc)) to help with normality\n", " Test1 = scipy.stats.chisquare(series[1], f_exp = 0.25)\n", " Test2 = scipy.stats.chisquare(series[2], f_exp = 0.25)\n", " print '\\t Surprise Test. Comparison to Chance: %s P-value: %s' % (str(np.round(Test1[0],4)),\n", " str(np.round(Test1[1],4)))\n", " print '\\t After Surprise Test. Comparison to Chance: %s P-value: %s' % (str(np.round(Test2[0],4)),\n", " str(np.round(Test2[1],4)))\n", " x = scipy.stats.chi2_contingency([[sum(series[1]==1),sum(series[2]==1)], [sum(series[1]==0),sum(series[2]==0)]],\n", " correction=False) \n", " print '\\t Chi-Square Comparison: %s P-value: %s' % (str(np.round(x[0],4)),str(np.round(x[1],4)))\n", " \n", "def Analysis_and_Plot(tableRT2, CIs): #this function plots the data and writes the results, including stats tests\n", " PlotFrame = pd.DataFrame([tableRT2.mean()])\n", " PlotFrame.columns = ['Presurprise', 'Surprise', 'Post surprise']\n", " PlotFrame2 = pd.DataFrame([CIs])\n", " PlotFrame2.columns = ['Presurprise', 'Surprise', 'Post surprise']\n", " PlotFrame.plot(ylim = [0, 1], yerr=PlotFrame2, kind='bar')\n", " plt.xticks(range(1), ['Trial Type'], rotation=0);\n", "\n", " print '---------------------------------'\n", " print 'Mean Presurprise: %s' % (str(round(np.mean(tableRT2[0]),2)))\n", " print 'Mean Surprise: %s' % (str(round(np.mean(tableRT2[1]),2)))\n", " print 'Mean Postsurprise: %s' % (str(round(np.mean(tableRT2[2]),2)))\n", " print 'Presurprise - Surprise: %s' % (str(round(np.mean(tableRT2[0])-np.mean(tableRT2[1]),2)))\n", " print 'Postsurprise - Surprise: %s' % (str(round(np.mean(tableRT2[2])-np.mean(tableRT2[1]),2)))\n", " print 'Postsurprise - Presurprise: %s' % (str(round(np.mean(tableRT2[0])-np.mean(tableRT2[1]),2)))\n", " print '---------------------------------'\n", " print 'Presurprise vs Surprise'\n", " print_tests(tableRT2[1],tableRT2[0])\n", " print 'Postsuprise vs Surprise'\n", " print_tests(tableRT2[2],tableRT2[1])\n", " print 'Presurprise vs Postsurprise'\n", " print_tests(tableRT2[0],tableRT2[2])\n", " \n", "def Analysis_and_Plot_2(tableRT2, CIs): #this function also plots the data and prints results.\n", " PlotFrame = pd.DataFrame([tableRT2.mean()]) #I could probably consolidate these functions, but whatever. This works.\n", " PlotFrame.columns = ['Surprise', 'Postsurprise']\n", " PlotFrame.plot(ylim = [0, 1], kind='bar')#yerr=PlotFrame2, kind='bar')\n", " plt.xticks(range(1), ['Trial Type'], rotation=0);\n", "\n", " print '---------------------------------'\n", " print 'Mean Surprise: %s' % (str(round(np.mean(tableRT2[1]),2)))\n", " print 'Mean Postsurprise: %s' % (str(round(np.mean(tableRT2[2]),2)))\n", " print 'Postsurprise - Surprise: %s' % (str(round(np.mean(tableRT2[2])-np.mean(tableRT2[1]),2)))\n", " print '---------------------------------'\n", " print 'Postsurprise vs Surprise'\n", " print_tests_ChiSq(tableRT2)\n", " \n", "def Analysis_and_Plot_3(tableRT2, CIs): #another plotting function...i should really consolidate these.\n", " PlotFrame = pd.DataFrame([tableRT2.mean()])\n", " PlotFrame.columns = ['Match', 'Mismatch']\n", " PlotFrame2 = pd.DataFrame()\n", " PlotFrame2['Match'] = pd.DataFrame([CIs])\n", " PlotFrame2['Mismatch'] = pd.DataFrame([CIs])\n", " PlotFrame.plot(ylim = [0, 1], yerr=PlotFrame2, kind='bar')\n", " plt.xticks(range(1), ['Trial Type'], rotation=0);\n", "\n", " #disp_tab = np.round(tableRT2,2)\n", " #disp_tab['Match-Mismatch'] = disp_tab[1] - disp_tab[2]\n", " #print disp_tab\n", "\n", " print '---------------------------------'\n", " print 'Mean match: %s' % (str(round(np.mean(tableRT2[1]),4)))\n", " print 'Mean mismatch: %s' % (str(round(np.mean(tableRT2[2]),4)))\n", " print 'Match - Mismatch: %s' % (str(round(np.mean(tableRT2[1])-np.mean(tableRT2[2]),4)))\n", " print '---------------------------------'\n", " print 'Match vs Mismatch'\n", " print_tests(tableRT2[1],tableRT2[2])\n", " \n", "def OneWayConfInterval(table): #Calculates confidence intervals for a one way anova, this is Cousineau, Morey ect\n", " import scipy.stats, numpy\n", " ParticipantsMeans = []\n", " STEs = []\n", " CIs = []\n", " for participant in table.index:\n", " mean = []\n", " for condition in xrange(numpy.shape(table)[1]): #there's definitely a better way to do this, but this works...\n", " mean.append(np.array(table[table.index==participant][condition]))\n", " ParticipantsMeans.append(sum(mean)/len(mean))\n", " ConfMeans = numpy.zeros(shape=numpy.shape(table))\n", " for counter, participant in enumerate(table.index):\n", " for condition in xrange(numpy.shape(table)[1]):\n", " ConfMeans[counter][condition] = table[table.index==participant][condition]-\\\n", " ParticipantsMeans[counter]+numpy.array(ParticipantsMeans).mean()\n", " for counter, column in enumerate(ConfMeans.T):\n", " STEs.append(numpy.std(column, ddof=1)/numpy.sqrt(len(column)))\n", " CIs.append(STEs[counter]*scipy.stats.t.isf(0.025, len(ConfMeans)-1))\n", " return CIs\n", "\n", "def SimpleComparisonCI(table): #confidence interval for pairwise comparisons - masson & loftus, Baguley (2012)\n", " import scipy.stats, math\n", " ttest = scipy.stats.ttest_rel(table[1], table[2])\n", " MeanDiff_byT = abs((table[1].mean()-table[2].mean())/ttest[0])\n", " CI = MeanDiff_byT*scipy.stats.t.isf(0.025, len(table)-1)*math.pow(2,0.05)/2\n", " return CI" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Experiment 1" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Next, load Experiment 1 data" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "filename = 'Exp1.csv' #looking at Exp1 data first. \n", "\n", "if sys.platform == 'linux2': #is this my linux laptop\n", " path = '/home/dan-laptop/Documents/Databrary/Wyble_PsychSci'\n", "elif sys.platform == 'darwin': #is this my mac work comp \n", " path = '/Users/danvatterott/Dropbox Encore/Dropbox/Databrary/Wyble_PsychSci'\n", "\n", "os.chdir(path)\n", "\n", "df = pd.read_csv(filename)\n", "\n", "df.columns = ['Sub#', 'Block', 'Trial#', 'TarCol', 'Tar_Iden','Tar_Loc', 'Col_Resp', 'Iden_Resp', 'Loc_Resp', \n", " 'Col_Acc', 'Iden_Acc', 'Loc_Acc'] #naming the columns of the data file." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The data is loaded, lets just take a quick look at the data after loading it in. " ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
Sub#BlockTrial#TarColTar_IdenTar_LocCol_RespIden_RespLoc_RespCol_AccIden_AccLoc_Acc
0612124002000
1613343003001
2614131001001
3615314004001
4616221001001
\n", "
" ], "text/plain": [ " Sub# Block Trial# TarCol Tar_Iden Tar_Loc Col_Resp Iden_Resp \\\n", "0 6 1 2 1 2 4 0 0 \n", "1 6 1 3 3 4 3 0 0 \n", "2 6 1 4 1 3 1 0 0 \n", "3 6 1 5 3 1 4 0 0 \n", "4 6 1 6 2 2 1 0 0 \n", "\n", " Loc_Resp Col_Acc Iden_Acc Loc_Acc \n", "0 2 0 0 0 \n", "1 3 0 0 1 \n", "2 1 0 0 1 \n", "3 4 0 0 1 \n", "4 1 0 0 1 " ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df[0:5]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "I want to create a new variable. \n", "\n", "Before explaining the new variable, I should explain a little about Chen and Wyble's experiment. Half the participants were instructed to find the letter among numbers and the other half were instructed to find the number among letters. 4 items were briefly flashed on the screen (150 ms) then participants reported where the target item had been. Each of the 4 items was a different color. \n", "\n", "Participants reported target location for 155 trials. On the 156th trial, the participants reported the target location then (to their surprise) reported what specific letter/number the target was and what color it was. Even though participants knew where the target was, they had no idea what either the target's letter/number or color were. They had \"forgotten\" what the target was (even though they must have known what the target was when they initially located it). \n", "\n", "The new variable will code whether the trial is a \"pre-surprise\" trial (trials 1-155), a surprise trial (trial 156) or the trial after the surprise (trial 157).\n", "\n", "I'll call this variable \"TrialType\"\n", "\n", "TrialType: 0=Presurprise; 1=Surprise Trial; 2=Postsurprise trials." ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": true }, "outputs": [], "source": [ "df['TrialType'] = np.where(df['Trial#']<156, 0, np.where(df['Trial#']==156, 1, 2))\n", "df['TrialType2'] = np.where(df['Trial#']<156, 0, np.where(df['Trial#']==156, 1, np.where(df['Trial#']==157, 2, -1)))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Experiment 1: Replicating Chen & Wyble's (2015) analyses" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Lets just take a quick look at overall accuracy. Make sure everyone is doing ok on the task. Below I plot the mean accuracy of each participant...looks like participant 23 struggled a little. Chen and Wyble (2015) notes that no participants were excluded or replaced. \n", "\n", "I might have replaced participant 23 since his/her accuracy is easily 2.5 standard deviations below the mean accuracy (I print this value below)...seems like participant 23 was doing something different in this task." ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "mean accuracy\n", "0.89\n", "standard deviation of accuracies\n", "0.07\n", "2.5 standard deviations below mean accuracy\n", "0.71\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAEXCAYAAABf36TeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHRVJREFUeJzt3XtwVOX9x/HP7oYkAglkI5dcRMVgRRRQwkUiYjCx4oBm\ntFWrVrRcBNEhOvpDCChYWkVABASKGC6ibQdaGzvWG0VsJSpGIYoVC4GiNiGQsBBWlkSyOb8/HHYI\nCTkHsoHl6fs148zunuf55nvOho9nz9lz4rIsyxIA4KzmPtMNAACajzAHAAMQ5gBgAMIcAAxAmAOA\nAQhzADBAlN2AxYsXa9OmTWrXrp1mz57d6Jhly5apuLhYMTExGj9+vC644IJw9wkAaILtnnlmZqby\n8vJOuHzz5s3as2eP5s+frzFjxmjp0qVhbRAAYM82zC+55BK1adPmhMuLioo0ePBgSVK3bt0UCAR0\n4MCB8HUIALDV7GPmPp9PiYmJoeder1c+n6+5ZQEAJ4EToABggGaHudfr1b59+0LP9+3bJ6/X29yy\nAICTYPttFkmyLEsnuh9Xenq63nnnHQ0cOFDbtm1TmzZt1L59e9uaZWVlTS6Pi4uT3+930l6LzDep\nRiT0EI4akdBDpNSIhB4ipUYk9HC6aiQnJ59wmW2Yz5s3T1999ZX8fr/GjRun2267TbW1tXK5XMrK\nytKVV16pzZs366GHHlJsbKzGjRt3amsBADhltmE+YcIE2yIjR44MSzMAgFPDCVAAMABhDgAGIMwB\nwACEOQAYwNFXEwEgUpSWxqiszBN67vFIwWDremOSk4NKSak5o32cjh6ORZgDOKuUlXmUk9P0tSwF\nBQeUknJm+zgdPRyLMEejIm2v40xiW+BsQJijUZG213EmsS1wNiDMw4y9OJiM3+/IRZiHGXtxMBm/\n35GLMIfR2JOMLMe/HxLvSbgQ5jAae5KRJVK+iWIiLhoCAAMQ5gBgAA6zAKcBx+7R0ghz4DTg2D1a\nGodZAMAAEbFnzteVwouP9GgMvxdmi4gw5+tK4cVHejSG3wuzcZgFAAxAmAOAAQhzADBARBwzDwdO\n7gD4X2ZMmHNyB8D/Mg6zAIABCHMAMABhDgAGMOaYOQCcTcJ95TthDgBnQLivfOcwCwAYgDAHAANw\nmAUtprkXctnNd1IDaIyJFxkS5mgxzb2Qi7tpoqWYeJEhh1kAwADsmR/DxI9eAP43EObHMPGjF4D/\nDY7CvLi4WCtWrJBlWcrMzFROTk695YFAQAsWLFBlZaXq6uo0fPhwXXvttS3RLwCgEbZhXldXp/z8\nfD3xxBNKSEjQpEmT1LdvX6Ucs4v6zjvv6LzzztPEiRN18OBB5ebmatCgQfJ4PE1UBgCEi22Yl5SU\nKCkpSR06dJAkZWRkqKioqF6Yu1wuHT58WJJUXV2tuLg4gvwU8cetAZwK2zD3+XxKTEwMPfd6vSop\nKak35oYbbtDMmTN1//33q7q6Wrm5ueHv9H8EX8cDcCrCcgK0uLhYF154oZ588kmVl5drxowZmj17\ntmJjY084Jy4uLvTYyU68x+OpN6fh8ubNj5QabAvn88NR42zZFseLjo4+qfHh6uFseE/DUeNs2RbH\nsg1zr9erysrK0HOfzyev11tvzPvvvx86Kdq5c2d17NhRpaWluuiii05Y1+/3hx4ff1VfY4LBoPz+\nQBPLm65hNz9SarAtnM8PR42zZVscLy4urt6/ISciYT1Ox3sajhqRui2aCnbbi4bS0tJUXl6uiooK\n1dbWqrCwUOnp6fXGnHvuudqyZYsk6cCBA9q9e7c6depk2ygAIDxs98zdbrdGjhypGTNmyLIsDRky\nRKmpqVq7dq1cLpeysrJ06623atGiRXr00UclSXfddZfatm3b4s0DAH7k6Jh57969NW/evHqvZWdn\nhx4nJCQoLy8vvJ0BABzj3iwAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAA\nwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDM\nAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwQNSZbgCAvdLSGJWVeULPPR4pGGxdb0xyclAp\nKTWnuzVECMIcOAuUlXmUk9O+yTEFBQeUknKaGkLE4TALABiAMAcAAxDmAGAAwhwADECYA4ABHH2b\npbi4WCtWrJBlWcrMzFROTk6DMf/617+0cuVKBYNBxcfH68knnwx7swCAxtmGeV1dnfLz8/XEE08o\nISFBkyZNUt++fZVyzHegAoGA8vPzNWXKFHm9Xh08eLBFmwYA1Gd7mKWkpERJSUnq0KGDoqKilJGR\noaKionpjNmzYoP79+8vr9UqS4uPjW6ZbAECjbPfMfT6fEhMTQ8+9Xq9KSkrqjSkrK1MwGNT06dNV\nXV2toUOH6pprrgl/twCARoXlCtC6ujr95z//0RNPPKGamhpNmTJFF198sTp37hyO8gAAG7Zh7vV6\nVVlZGXru8/lCh1OOHRMXF6fo6GhFR0ere/fu2rVrV5NhHhcXF3rs8Zxw2DFjPPXmNFzevPmRUoNt\n4Xx+OGqwLZzPD0eN07Ee4ahxtmyLY9mGeVpamsrLy1VRUaGEhAQVFhZqwoQJ9cb07dtXy5YtU11d\nnY4cOaLt27dr2LBhTdb1+/2hx8ffMKgxwWBQfn+gieVN17CbHyk12BbO54ejBtvC+fxw1Dgd6xGO\nGpG6LZoKdtswd7vdGjlypGbMmCHLsjRkyBClpqZq7dq1crlcysrKUkpKinr16qVHH31UbrdbWVlZ\nSk1NtW0UABAejo6Z9+7dW/Pmzav3WnZ2dr3nN910k2666abwdQYAcIwrQAHAAIQ5ABiAMAcAAxDm\nAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4A\nBiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAA\nwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAEdhXlxcrNzcXE2YMEEFBQUnHFdS\nUqJf/OIX2rhxY9gaBADYsw3zuro65efnKy8vT3PmzFFhYaFKS0sbHff73/9evXr1apFGAQAnZhvm\nJSUlSkpKUocOHRQVFaWMjAwVFRU1GPf2229rwIABio+Pb5FGAQAnZhvmPp9PiYmJoeder1c+n6/B\nmKKiIl1//fXh7xAAYCssJ0BXrFihu+66K/TcsqxwlAUAOBRlN8Dr9aqysjL03Ofzyev11huzc+dO\nPf/887IsS36/X5s3b1ZUVJTS09NPWDcuLi702OOxb9Tj8dSb03B58+ZHSg22hfP54ajBtnA+Pxw1\nTsd6hKPG2bItjmUb5mlpaSovL1dFRYUSEhJUWFioCRMm1BvzwgsvhB4vWrRIffr0aTLIJcnv94ce\nB4OtbRsNBoPy+wNNLG+6ht38SKnBtnA+Pxw12BbO54ejxulYj3DUiNRt0VSw24a52+3WyJEjNWPG\nDFmWpSFDhig1NVVr166Vy+VSVlaWbUMAgJZlG+aS1Lt3b82bN6/ea9nZ2Y2OfeCBB5rfFQDgpHAF\nKAAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgD\ngAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAY\ngDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMEOVkUHFx\nsVasWCHLspSZmamcnJx6yzds2KDXX39dkhQbG6vRo0erS5cu4e8WANAo2z3zuro65efnKy8vT3Pm\nzFFhYaFKS0vrjenYsaOmT5+uWbNm6dZbb9WSJUtarGEAQEO2YV5SUqKkpCR16NBBUVFRysjIUFFR\nUb0xF198sVq3bi1J6tatm3w+X8t0CwBolG2Y+3w+JSYmhp57vd4mw3rdunXq3bt3eLoDADji6Ji5\nU19++aXef/99PfXUU7Zj4+LiQo89HvvaHo+n3pyGy5s3P1JqsC2czw9HDbaF8/nhqHE61iMcNc6W\nbXEs2zD3er2qrKwMPff5fPJ6vQ3GffPNN3rxxRc1efJktW3b1vYH+/3+0ONgsLXt+GAwKL8/0MTy\npmvYzY+UGmwL5/PDUYNt4Xx+OGqcjvUIR41I3RZNBbvtYZa0tDSVl5eroqJCtbW1KiwsVHp6er0x\nlZWVmjNnjh588EF17tzZtkEAQHjZ7pm73W6NHDlSM2bMkGVZGjJkiFJTU7V27Vq5XC5lZWXpT3/6\nk77//nvl5+fLsix5PB49/fTTp6N/AIAcHjPv3bu35s2bV++17Ozs0OOxY8dq7Nix4e0MAOAYV4AC\ngAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAY\ngDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEI\ncwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGiHIyqLi4WCtWrJBl\nWcrMzFROTk6DMcuWLVNxcbFiYmI0fvx4XXDBBeHuFQBwArZ75nV1dcrPz1deXp7mzJmjwsJClZaW\n1huzefNm7dmzR/Pnz9eYMWO0dOnSFmsYANCQbZiXlJQoKSlJHTp0UFRUlDIyMlRUVFRvTFFRkQYP\nHixJ6tatmwKBgA4cONAyHQMAGrANc5/Pp8TExNBzr9crn8930mMAAC2HE6AAYACXZVlWUwO2bdum\nNWvWKC8vT5JUUFAgSfVOgr744ou67LLLNHDgQElSbm6upk2bpvbt27dU3wCAY9jumaelpam8vFwV\nFRWqra1VYWGh0tPT641JT0/XP/7xD0k/hn+bNm0IcgA4jWz3zKUfv5q4fPlyWZalIUOGKCcnR2vX\nrpXL5VJWVpYkKT8/X8XFxYqNjdW4cePUtWvXFm8eAPAjR2EOAIhsnAAFAAMQ5gBgAMIcAAxAmAOA\nARzdaKslbd++XSkpKWrdurV++OEHFRQUaOfOnUpNTdUtt9yi1q1b29Z488031a9fP5177rmn3MfR\nr10mJCSoZ8+e2rBhg/79738rJSVFWVlZioqy31R79uzRxo0btW/fPrndbiUlJenqq692tA7A6VZV\nVaV27dqd0R78fr/i4uLOaA+mOOPfZnnkkUc0a9YseTweLVmyRDExMRowYIC2bNmib775Ro8++qht\njREjRig2NladOnVSRkaGrrrqKsXHx59UH/Pnz1cwGFRNTY3atGmj6upq9e/fX1u2bJFlWXrwwQeb\nnP/mm29q06ZN6t69uzZv3qwLLrhAbdq00SeffKJRo0apR48eJ9UPGhcJAXTUb3/7W02ePLnJMYFA\nQAUFBdq3b5+uuOIKXX311aFlL730kkaNGmX7cw4cOKA1a9bI5XLp9ttv11tvvaWNGzcqJSVF9913\nnxISEmxrfP/99/WeW5alxx9/XDNnzpQktW3b1rZGcXGxevfuHVqvlStXaseOHTrvvPM0YsQI22tL\nXn31VQ0fPlzx8fHasWOH5s6dK5fLpWAwqAcffFCXXnqpbQ8TJ05Uv379lJGRoc6dO9uOb8yOHTv0\nyiuvKCEhQXfeeacWL16skpISJScna8yYMbrwwgtta1RXV+v1118P7bxFRUWpc+fOys7O1rXXXms7\nPxAI6C9/+YuKiopUVVUll8uldu3aKT09XTk5OWrTps1Jr9cZ3zO3LEsej0eStHPnztAv1yWXXKLH\nHnvMUY1OnTrpmWee0ZYtW/Thhx9q9erV6tq1qzIyMtS/f3+dc845tjW+/fZbzZ49W8FgUGPHjtWS\nJUvkdrs1aNAgR32sW7dOs2bNktvt1rBhw/T0009r2rRpys7O1rPPPqtnn33W0bq0xJt8lJPwOdpD\npAbQ5MmTHQdQc8NH+vF38kR27dplO3/RokVKSkpS//79tX79en388ceaMGGCWrVqpe3bt9vOl6SF\nCxfqyiuvVE1NjaZPn66rr75akyZNUlFRkZYuXar/+7//s60xcuTIBp9cfT6fJk6cKJfLpRdeeMG2\nxh/+8IfQ9nz55ZeVkJCgiRMnauPGjXrxxRdt+9i0aZPuuusuSdIrr7yi3NxcpaWlqaysTPPnz9cz\nzzxj28P333+vQ4cOafr06Wrfvr0yMjI0cOBAeb1e27lHvfTSS7rtttt06NAhTZ06VSNGjNDUqVO1\nZcsWvfTSS/rNb35jW2P+/Pnq16+f8vLy9NFHH6m6uloZGRn685//rLKyMt15551Nzp87d6569OhR\n70r5AwcO6P3339fcuXM1ZcoUx+tz1BkP8/POO0/r169XZmamzj//fO3YsUMXXXSRysrKHB3akCSX\nyyW3261evXqpV69eqq2tVXFxsTZs2KBVq1YpPz/ftoZlWaqtrVV1dbVqamoUCATUtm1bHTlyRMFg\n0FEfwWBQbrdbR44cUXV1tSTp3HPPdTxfav6b3NzwkcwJoOaGjyRNmjTphHuMhw4dsp2/Z8+e0KfL\nfv366bXXXtNTTz3l6GcfVVVVpaFDh0qS3nnnndCtNIYOHar33nvPUY27775bX3zxhX75y1+qS5cu\nkqTx48dr4cKFjvs41o4dOzRr1ixJ0rBhw0JXgDelrq5OwWBQHo9HP/zwg9LS0iRJycnJOnLkiKOf\n27ZtW91zzz265557tHXrVhUWFmrixIlKTU1VRkZG6CLGpgSDQV1xxRWSfvy0MGDAAEnS5ZdfrlWr\nVjnqo6KiIrQHPmzYME2aNEk/+9nP9MADD+iRRx6xDfO9e/eGbpFyVPv27ZWTk6P169c76uF4ZzzM\nx44dq+XLl+u1115TXFycpkyZosTERCUmJur+++93VOP4I0VRUVFKT09Xenq6ampqHNXIzMxUbm6u\n6urqdMcdd+i5555Tx44dtX379tA9Z5py3XXXadKkSUpLS9PXX3+tm2++WZJ08OBBRx9hj2rum9zc\n8JHMDKBTCR9JSk1N1ZgxY5SUlNRg2bhx42zn19bWqq6uTm73j981uOWWW+T1evXkk0+G/odv59jf\n76O3mj6qrq7OUY3hw4dr4MCBWrlypRITE3XbbbfJ5XI5mntUVVWV3njjDVmWpUAgIMuyQjWcHK29\n/vrr9fTTTysnJ0e9evXS8uXL1b9/f3355Zen9Mdsunfvru7du+tXv/qVvvjiC3344YeOwrxVq1b6\n/PPPFQgE5HK59Mknn6hfv3766quvQu+TnZiYGH399de65JJL9Omnn4b+jbvdbkfbokOHDnr99dc1\nePDgBjttp3ru74yHeevWrTV+/HgFAgHt3btXdXV18nq9J3Vvl9zc3BMui4mJcVRj2LBhodD2er0a\nPHiwtmzZoqysrNAeRFNuvPFGXX755SotLdXw4cOVkpIiSYqPj9f06dMd9SA1/01ubvhI5gRQc8NH\nkn7+85+fcOx9991nO79Pnz768ssv1bNnz9Br1157rdq3b69ly5Y56iE9PV3V1dWKjY3VHXfcEXq9\nvLxcycnJjmpIUmJioh555BF9+umnmjFjhuMdnaOuu+46HT58OLQOfr9f8fHxOnDggKMwHjp0qLp0\n6aJ3331Xu3fvVjAYVHl5ufr27atbbrnFUQ+N/V673W717t079CnMzujRo/Xqq6/K5XIpLy9P7777\nrhYtWiSv1+t4B3LUqFFasmSJysvLlZqaGvq3dfDgQf30pz+1nZ+bm6uCggJNmzZNVVVVkn7caevT\np48efvhhRz00YCGi+P1+a9WqVdaECROse++917r33nut3Nxca9WqVZbf77ed/9FHH1mlpaWNLtu4\ncaOjHlatWmV9/vnnDV7fvHmz9dBDDzmq8cc//tE6fPhwg9d3795tzZ4921GNYxUVFVmTJ0+2Ro0a\n5XjO6tWr6/1XVVVlWZZl7d+/31qwYIHjOv/973+tL774osH6bN68uVnzN23adNp6OL5GTU2N9c03\n3zSrxqn0Ee71ONUa3333XVhqfP7556dcY/v27db27dsty7Ksb7/91vrrX/9qffbZZ45//vE806ZN\nm3Zq/xtAS4iOjlbPnj11ww03KCcnRzk5ObrhhhvUs2dPFRYW2p5pT01NPeFXvbZt2+boTH3Pnj3V\nqVOnBq937txZsbGxjmpcdtlljZ7zaNu2rWpqahzVOFZycrIGDRqkPn36qF27dlq/fr1tjR49etT7\n7+intNjYWB06dMhRD2+++aZefvll7dmzR6tXr1bHjh1Dn7rmzp2r7OzsJue/9dZbWrlyZaPzn3/+\nedv5djWc9NDYeiQlJYUOxzmt0dw+mrstw9HD0T5WrVp1RmusWbNGb7/9tj777DPt3btX7733nrxe\nr/75z39q//796t69u20Pxzvjh1ng3OrVq5WZmXnG5p/pGtHR0aFj56drW6xbt04zZ85UbGys9u7d\nq+eee04VFRW68cYbHR2q+fvf/96s+eGq0dz1CEcfkdDD0T6eeeaZM1rj448/1qxZs3TkyBGNGTNG\nixcvVuvWrXXTTTdp8uTJjg87HYswjzAn+l69ZVmhY2stOd+kGuHowbIsxcbGSpI6duyoadOmac6c\nOaqoqHD0j7a5802qEQk9REoNj8cjt9utmJgYderUKXRhYXR09EmfmD6KMI8wVVVVysvLa/B9csuy\nNHXq1Bafb1KNcPTQrl077dq1K3SCLzY2Vo8//rgWL16sb7/9tsXnm1QjEnqIlBpRUVGqqalRTExM\nve/XBwIBx9+oaeCUj7ajRSxatMjaunVro8uef/75Fp9vUo1w9FBZWWnt37+/0WUnqh3O+SbViIQe\nIqXGDz/80OjrVVVVoRPTJ+uMX84PAGg+7poIAAYgzAHAAIQ5ABiAMAeOMX36dMf3jgEiCV9NhJG+\n/vprvfrqq/ruu+/k8XiUkpKie++9V127dg1L/erqaj388MNavHix1q9fr++++0733HNPWGoDp4Iw\nh3EOHz6smTNnavTo0brqqqtUW1urrVu3Or6lshO7du0K3Q5g586d+slPfhK22sCp4DALjLN7925J\n0sCBA+VyudSqVSv17NlTXbp00Zo1a7RgwYLQ2IqKCt1+++317uRYXl6uyZMna8SIEZo1a1ajtw7e\nsWNHaC9/586dp3QLVyCc2DOHcZKSkuR2u7Vw4UJlZGSoW7du9a4Ctbtc+oMPPlBeXp46duyoBQsW\naNmyZXrooYckSb/73e+0ceNG1dTUqFWrVvrb3/6m6upqTZ06VV6vV3PmzGnRdQNOhD1zGOecc87R\nr3/9a7lcLi1ZskSjRo3Ss88+6/h+LIMGDVJqaqqio6N1++236+OPPw7db2Ps2LF64YUX1LFjR+Xn\n52vkyJHKzs7W8uXLCXKcUeyZw0jJycl64IEHJEllZWVasGCBVqxY4eiPORz7R0A6dOig2tpa+f1+\nbdu2TQsXLgz98Y7Ro0erurpa0dHR+uCDDzR16tSwnWAFThZhDuMlJydr8ODBWrdunbp27VrvL+zs\n37+/wfjKysrQ44qKCkVFRSkuLk7p6elavny5li5dqh49emjgwIEaM2aMFi1aFNaTq8Cp4DALjFNW\nVqY33nhDPp9P0o/hXFhYqG7duun888/X1q1bVVlZqUAgoIKCggbzP/jgA5WWlqqmpkarV6/WgAED\n6h1n37lzpy688ELt3btXCQkJBDkiAr+FME5sbKy2b9+uN954Q4FAQG3atFGfPn109913KzY2VgMH\nDtRjjz2m+Ph43Xzzzfrss8/qzb/mmmu0cOFClZWV6dJLL9WYMWNCy4LBoCorK5WUlKSNGzee9F9M\nAloKd00EAANwmAUADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAzw/2/Y6HVu\nT4m7AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "tableAcc = df.pivot_table(values='Loc_Acc', index='Sub#', aggfunc=np.mean)\n", "#print tableAcc\n", "print 'mean accuracy'\n", "print np.round(np.mean(tableAcc),2)\n", "print 'standard deviation of accuracies'\n", "print np.round(np.std(tableAcc),2)\n", "\n", "print '2.5 standard deviations below mean accuracy'\n", "print np.round(np.mean(tableAcc)-2.5*np.std(tableAcc),2)\n", "\n", "tableAcc.plot(ylim = [0, 1], kind='bar');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Lets look at participants' performance when asked to identify the target's location. I will plot performance as mean accuracy in the presurprise,surprise, and postsurprose trials. \n", "\n", "I will also run some quick statistical tests. For these tests, I take the arcsine of the square root of the accuracies (Rao, 1960) to increase the accuracies' normality (I use adj. to indiciate that the tested data is transformed). I test whether this worked with a Shapiro-Wilk test of normality. If the p-value of the Shapiro-Wilk test is greater than 0.1, I run a t test to see if the accuracy in the two conditions is significantly different. If the p-value of the Shapiro-Wilk test is less than or equal to 0.1, then I run a Wilcoxon signed rank test since this test does not care about normality. " ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "---------------------------------\n", "Mean Presurprise: 0.89\n", "Mean Surprise: 0.8\n", "Mean Postsurprise: 0.79\n", "Presurprise - Surprise: 0.09\n", "Postsurprise - Surprise: -0.01\n", "Postsurprise - Presurprise: 0.09\n", "---------------------------------\n", "Presurprise vs Surprise\n", "\t normality test adj. Test value: 0.64 P-value: 0.0\n", "\t Wilcoxon. Test value: 74.0 P-value: 0.25\n", "Postsuprise vs Surprise\n", "\t normality test adj. Test value: 0.8 P-value: 0.001\n", "\t Wilcoxon. Test value: 33.0 P-value: 0.63\n", "Presurprise vs Postsurprise\n", "\t normality test adj. Test value: 0.94 P-value: 0.2857\n", "\t T-test adj. Test value: 0.92 P-value: 0.3695\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAECCAYAAAAMxDf2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHvhJREFUeJzt3XtU1HX+P/DnXAQEQRxEYUCPF1CBLFMslZRL6GYJTYta\n1nbZZdUkXa3UjothuiRpYtH1ZAdyXXNdU8PVzrEoqj1SuvNNxlVYL6OhMgiCE4GNF5j5/P5wnR8j\nDDPADIPvno9zOofPZ96fz/s12Dznw/tzecskSZJARES3NbmnCyAioq5jmBMRCYBhTkQkAIY5EZEA\nGOZERAJgmBMRCUDpqMH777+Pw4cPo2/fvtiwYUObbQoKCqDT6eDt7Y3nnnsOQ4YMcXWdRETUDodH\n5omJicjMzLT7emlpKWpqavDWW29h3rx5+PDDD11aIBEROeYwzEeNGgU/Pz+7r2u1WsTHxwMAIiMj\nYTKZUF9f77oKiYjIoS6PmRuNRgQFBVmXVSoVjEZjV3dLREQdwBOgREQC6HKYq1QqXLp0ybp86dIl\nqFSqru6WiIg6wOHVLAAgSRLsPY8rNjYWn3/+OSZNmoSTJ0/Cz88PgYGBDvdZVVXVsUqJuom/vz8a\nGxs9XQZRK2q12u5rMkdPTczLy0N5eTkaGxvRt29fzJ49G83NzZDJZEhOTgYA5OfnQ6fTwcfHBwsW\nLMCwYcMcFsUwp56KYU49VZfC3F0Y5tRTMcypp2ovzHkClIhIAAxzIiIBOHUClOjUqVPQ6/UAAL1e\nj4iICABAREQEIiMjPVkaEYFhTk6KjIy0hvaSJUuwaNEiD1dERC1xmIWISAAMcyIiATDMiYgEwDAn\nIhIAw5yISAAMcyIiATDMiYgEwDAnIhIAw5yISAAMcyIiAfB2fiLYPnvmxx9/xNChQwHw2TN0+2CY\nE6H1s2cyMjI8XBFRx3CYhYhIAAxzIiIBMMyJiATAMXMBGAzeqKpSdGufWq1vt/WlVpsRFnat2/oj\nuh0xzAVQVaWARhPYbf3NmoVu7a+wsB5hYd3WHdFticMsREQCYJgTEQmAYU5EJACGORGRAHgClEhw\nLR9VoNfrERERAYCPKhANw5xIcLc+qmDRokUerojcgcMsREQCYJgTEQmAYU5EJACGORGRABjmREQC\nYJgTEQmAYU5EJACGORGRABjmREQCcOoOUJ1Oh82bN0OSJCQmJkKj0di8bjKZ8Pbbb6Ourg4WiwUp\nKSlISEhwR73kIcHBZQgOLgMAXLwYhejoHQCA2toY1NbGeLI0IoITYW6xWJCfn4+srCz069cPK1as\nwPjx4xHWYraAzz//HIMGDcJLL72EhoYGLFmyBJMnT4ZC0b2z35D7MLSJejaHwyx6vR6hoaEIDg6G\nUqlEXFwctFqtTRuZTIYrV64AAK5evQp/f38GORFRN3IY5kajEUFBQdZllUoFo9Fo0+aBBx5AZWUl\n5s+fj2XLluGZZ55xeaFERGSfS56aqNPpMHToUKxatQrV1dXIzs7Ghg0b4OPjY3cbf39/V3RNAET/\nI8hrYAUOGyu7tc/DxsPd1le4fziGq4Z3W3/87InJYZirVCrU1dVZl41GI1QqlU2bb775xnpSNCQk\nBAMGDIDBYMDw4fb/B21sbOxszXQLs9nX0yW4VV3TeTy1K7Xb+puFWUjZldJt/RVqCjGg14Bu64+f\nvdtXe1/EDodZIiIiUF1djdraWjQ3N6OkpASxsbE2bfr374+jR48CAOrr63HhwgUMHDiwi2UTEZGz\nHB6Zy+VypKenIzs7G5IkISkpCeHh4SgqKoJMJkNycjLS0tLw3nvvYenSpQCAJ554An369HF78UQi\nGPlzL/hWaB03dBFfbff1ZVarca3FlW/kPk6NmY8ZMwZ5eXk266ZOnWr9uV+/fsjMzHRtZUS/En41\nlxA4+6nu6WzWLATecp+IO9UXFgIM827BO0CJiATAMCciEgDDnIhIAAxzIiIBMMyJiATAMCciEgDD\nnIhIAAxzIiIBMMyJiATAMCciEgDDnIhIAAxzIiIBMMyJiATAMCciEgDDnIhIAAxzIiIBuGRCZyLq\nucqCg1EWHAwAiLp4ETuiowEAMbW1iKmt9WRp5EIMcyLBMbR/HTjMQkQkAIY5EZEAGOZERAJgmBMR\nCYAnQIkABP8SjOBfblzxcdH3IqIv3rjio9avFrV+PHlIPR/DnAgMbbr9cZiFiEgADHMiIgEwzImI\nBMAwJyISAMOciEgADHMiIgEwzImIBMAwJyISAMOciEgADHMiIgEwzImIBODUs1l0Oh02b94MSZKQ\nmJgIjUbTqk1ZWRn++te/wmw2IyAgAKtWrXJ5sURE1DaHYW6xWJCfn4+srCz069cPK1aswPjx4xEW\nFmZtYzKZkJ+fj5UrV0KlUqGhocGtRRMRkS2Hwyx6vR6hoaEIDg6GUqlEXFwctFqtTZsDBw7g3nvv\nhUqlAgAEBAS4p1oiImqTwyNzo9GIoKAg67JKpYJer7dpU1VVBbPZjNWrV+Pq1auYPn06pkyZ4vpq\niYioTS55nrnFYsGPP/6IrKwsXLt2DStXrsSIESMQEhLiit0TEZEDDsNcpVKhrq7Oumw0Gq3DKS3b\n+Pv7w8vLC15eXoiKikJFRUW7Ye7v79+FsqklhcLTFbiXTCbzdAluJfL7UygU/Kx3E4dhHhERgerq\natTW1qJfv34oKSnB4sWLbdqMHz8eBQUFsFgsaGpqwqlTpzBjxox299vY2Ni1ysnKbPb1dAluJUmS\np0twK5Hfn9lshomfdZdp74vRYZjL5XKkp6cjOzsbkiQhKSkJ4eHhKCoqgkwmQ3JyMsLCwnDXXXdh\n6dKlkMvlSE5ORnh4uEvfBBER2efUmPmYMWOQl5dns27q1Kk2y6mpqUhNTXVdZURE5DTeAUpEJACG\nORGRABjmREQCYJgTEQmAYU5EJACGORGRABjmREQCYJgTEQmAYU5EJACGORGRABjmREQCYJgTEQmA\nYU5EJACGORGRABjmREQCYJgTEQmAYU5EJACGORGRABjmREQCYJgTEQmAYU5EJACGORGRABjmREQC\nYJgTEQmAYU5EJACGORGRABjmREQCYJgTEQmAYU5EJACGORGRABjmREQCYJgTEQmAYU5EJACGORGR\nABjmREQCYJgTEQnAqTDX6XRYsmQJFi9ejMLCQrvt9Ho95syZg0OHDrmsQCIicsxhmFssFuTn5yMz\nMxO5ubkoKSmBwWBos922bdtw1113uaVQIiKyz2GY6/V6hIaGIjg4GEqlEnFxcdBqta3a7d+/HxMm\nTEBAQIBbCiUiIvschrnRaERQUJB1WaVSwWg0tmqj1Woxbdo011dIREQOueQE6ObNm/HEE09YlyVJ\ncsVuiYjISUpHDVQqFerq6qzLRqMRKpXKps2ZM2fw5ptvQpIkNDY2orS0FEqlErGxsXb36+/v34Wy\nqSWFwtMVuJdMJvN0CW4l8vtTKBT8rHcTh2EeERGB6upq1NbWol+/figpKcHixYtt2rzzzjvWn997\n7z2MGzeu3SAHgMbGxk6WTLcym309XYJbif6Xnsjvz2w2w8TPusu098XoMMzlcjnS09ORnZ0NSZKQ\nlJSE8PBwFBUVQSaTITk52aXFEhFRxzkMcwAYM2YM8vLybNZNnTq1zbYZGRldr4qIiDqEd4ASEQmA\nYU5EJACGORGRABjmREQCYJgTEQmAYU5EJACGORGRABjmREQCYJgTEQmAYU5EJACGORGRABjmREQC\nYJgTEQmAYU5EJACGORGRABjmREQCYJgTEQmAYU5EJACnpo3rTgqFAr6+Yk9Q7GqxsXLs33/N02W0\nS5KAn3+WY/lyJc6dE3c2eiJP6VFhrlAo4Ofnh4aGBk+XcluRyYDRoz1dhXM+/VSFRx7xZqATuViP\nGmbx9fVlkAsuJMSI119v9nQZRMLpUWFOvw4BARZPl0AkHIY5dTsZR1iIXI5hTkQkAIb5r9CTTz6J\nnTt3eroMInKhHnU1S1sMBm9UVSnctn+12oywMOcu67v33ntRV1cHpVIJX19fJCQkYO3atejdu7fb\n6nOHv/3tb54ugYhcrMeHeVWVAhpNoNv2X1hYj7Aw59rKZDJs2bIFcXFxqKmpweOPP44333wTK1as\nsGknSRJk3TAwbDaboVB07Iuuu2ojou7FYZYOkiQJADBw4EAkJibi+PHjmDlzJtatWweNRoOIiAic\nO3cOjY2NePHFFzF27FjExsZi/fr11m0rKiowc+ZMREVF4c4770RGRgYAoLKyEuHh4bBY/v/VHjNn\nzsT27dsBADt27IBGo8Err7yCO+64Axs3brSuW7lyJaKiopCQkIADBw7YbH9rbS33aa8WANDr9Zgz\nZw5iYmIQHx+PvXv3uveXS0Sd1uOPzHsqg8GA4uJiPPjggzh48CB2796Njz/+GMOGDYPFYsH8+fMx\nYMAAfPfddzCZTHjqqacQFhaGJ554Aq+//jri4+Oxc+dOXL9+HUeOHLHu19FRc2lpKTQaDf7zn/+g\nqakJe/bsQWlpKVJSUnDs2DF89tlnmDt3Lg4ePIi+ffsCQKvaWrJXy5UrVzBnzhwsX74c27ZtQ3l5\nOebMmYOoqChERES4+LdJRF3FI/MOSk9PR0xMDNLS0jBp0iQsWrQIADB79mxERERALpejvr4eX3/9\nNV555RX4+PhApVJh7ty52LNnDwBAqVSisrISFy5cgJeXF8aPH+90/yEhIXjmmWcgl8vh7e0NAOjf\nvz/S09OhUCiQmpqK4cOH46uvvrJu07I2pdL2+9teLUVFRRg8eDBmzZoFmUyGmJgYTJ8+nUfnRD0U\nj8w7qKCgAHFxca3Wq9Vq68+VlZVoamrC2LFjAdwYmpEkCWH/G5x/+eWXsX79esyYMQOBgYGYN28e\nHn30Uaf6b9nPTaGhoTbLYWFhqKmpaXebm+zVYjAYcPjwYcTExFjfg9lsRlpamlN1ElH3Yph30M1x\n71u1HB5Rq9Xw9vbGsWPH2hw26d+/P9avXw8A0Gq1eOyxxzBhwgT4+/sDuDHE4efnBwCora21289N\nFy5csFk2GAz4zW9+0+42jmpRq9WYOHEitm3bZndbIuo5OMziBgMGDEB8fDxWrVqFy5cvQ5IknD17\nFgcPHgQA7Nu3zxrAAQEBkMlkkMvlUKlUCAkJwa5du2CxWLB9+3acPXvWYX+XLl1CQUEBmpubsXfv\nXpw+fRr333+/U7XaqyU5ORlnzpzBrl270NzcjKamJhw5cgR6vb6TvxUicqcef2SuVptRWFjv1v07\ny94Rblvr8/Ly8OqrryIhIQEmkwmDBw+2Xily5MgRa9D3798ff/nLXzBo0CAAN05IrlixAuvWrcNj\njz2G2NhYh3Xdfffd+PHHHzF69GgEBwdj06ZN1pOfbdXWcl17tWzbtg2vvPIKVq9eDUmSEB0djVWr\nVjmsh4i6n0yyN27gZlVVVa3W+fv7o7Gx0QPV3L527NiB7du3Y/fu3Z4uxWlHjwbhgQe8nW6/5V97\n8VRxqhsr8qxz0VswaPZTni7DLeoLC2HqwAl+al975784zEJEJACnhll0Oh02b94MSZKQmJgIjUZj\n8/qBAwesl935+Phg7ty5GDx4sOurJSKiNjk8MrdYLMjPz0dmZiZyc3NRUlICg8Fg02bAgAFYvXo1\nXn/9daSlpeGDDz5wW8Fka/bs2bfVEAsRuYfDMNfr9QgNDUVwcDCUSiXi4uKg1Wpt2owYMcI6b2dk\nZCSMRqN7qiUiojY5DHOj0YigoCDrskqlajesv/rqK4wZM8Y11RERkVNcemnisWPH8M0332DNmjUO\n2968Qaaljj4BkH4dRH/Ko8jvT6FQtPlZJ9dzGOYqlQp1dXXWZaPRCJVK1ard2bNnsWnTJvz5z39G\nnz59HHbc1iWI/Eentnjo6tluI/L7M5vNMPFyY5dpLyMdDrNERESguroatbW1aG5uRklJSasbWerq\n6pCbm4uFCxciJCSk6xUTEVGHODwyl8vlSE9PR3Z2NiRJQlJSEsLDw1FUVASZTIbk5GTs3LkTly9f\nRn5+PiRJgkKhQE5OTnfU/6thMBiQlJSE48ePC/1nORF1jlNj5mPGjEFeXp7NuqlTp1p/fvbZZ/Hs\ns8+6trL/MZgMqPql9d2irqL2UyPM17mphv7973/j1VdfxcmTJ6FQKBAZGYnVq1fjzjvvdFt9N4WF\nheHEiRNu74eIbk89/tksVb9UQVOocdywkwo1hU6F+eXLl/HMM8/gtddeQ0pKCq5fv45Dhw7By8ur\nw312dOq2zkwPR0S/Lryd30lnzpyBTCZDamoqZDIZvL29MWXKFIwaNQobN260TlIBtJ7+zd7UbTk5\nOZgxYwZGjRqF9PR0/Pzzzzbbb9++Hffccw8effTRVvv8xz/+gUmTJmHkyJGYNGkSCgsLrf1v374d\nCQkJiImJwe9+97tWN3kRkXgY5k4aNmwY5HI5lixZgq+//toavDfdeqR96/Lu3buxYcMGnDx50jpJ\nxa5du/DGG2+gtLQUcrkcK1eutNnm4MGD+Pbbb/Hxxx/b7PPKlStYtWoVPv74Y5w4cQJ79uyxTiLx\n+eef45133kF+fj6OHj2Ke+65x2ZeTyISE8PcSX369MGnn34KuVyO5cuX484778Qf/vAHm8s229PW\n1G1paWmIjIxE7969sXz5cuzbt896mZpMJsPSpUvRu3dv6/RwLSkUChw/fhxXr15FcHAwIiMjAQBb\nt27FokWLMHz4cMjlcixcuBBlZWU8OicSHMO8AyIiIrBx40ZotVoUFxejpqbG6ed7t/XoypbrwsPD\n0dTUZHN37a3Twd3Uu3dvvP/++9iyZQvGjh2Lp59+GqdPnwZwY4gmKysLMTExiImJwR133AGZTIbq\n6uqOvFUius0wzDtp+PDhmDVrFk6cOAFfX19cuXLF+lrL+TdvauuEZ8tnuldWVsLLy8vmhqz2TpJO\nmTIFf//731FaWorhw4dj+fLlAG58Qaxbtw5lZWUoKytDeXk5Tp06hXHjxnXqfRLR7YFh7iS9Xo8P\nPvjAOsWawWBAYWEhxo0bh+joaBw6dAgGgwENDQ149913ndrn7t27odfrceXKFWzYsAEPPfSQNcDb\nuivw5rq6ujp88cUXuHLlCnr16gU/Pz/rdk8++STefvttnDx5EgDQ0NCAffv2dfn9E1HP1uMvTVT7\nqVGoKXTcsAv7d0afPn1QWlqKTZs2obGxEQEBAZg6dSpWrlwJPz8/pKamYurUqVCpVHjuuedQVFRk\n3dbeEXZaWhqWLFmC06dPY+LEiXjttdfa3ebmOovFgk2bNmHJkiWQyWSIjo62bvvAAw/AZDIhIyMD\nBoMB/v7+mDJlCmbMmOH074SIbj+cNs5DZs6ciZkzZ+Kxxx7zdCndjtPG2eK0ceQsThtHRCQ4hrmH\n8PkqRORKPX7MXFSffPKJp0sgIoHwyJyISAAMcyIiATDMiYgEwDAnIhIAw5yISAAMc2pXUlISDh48\n6OkyiMiBHn9porfBAEUbd4u6ilmtxrUw56aNu/fee1FXVwelUglfX18kJCRg7dq16N27d6f6rqys\nxIQJE3Du3DnI5T3ze7W4uNjTJRCRE3p8mCuqqhCocd+0cfWFhYCTYS6TybBlyxbExcWhpqYGjz/+\nON58802sWLGiU33fnD6uO56o0NGp5zhVHdHtpWceDvZgN4N34MCBSExMtE6yXFNTg9///veIiYnB\nfffdh23btlm30el0ePDBBzFq1CjcfffdWLNmDYAbD9oCgKioKIwcORKHDx9u1Z+9bb///nvExsba\ntJ0wYQIOHDgAANi4cSPmzZuHRYsWISoqCp988ol13YIFCzBy5EhMnz4d5eXlNtu/9957SE5OxogR\nI2A2m232aa8WAPjhhx/w8MMPIzo6GtOmTcP333/ftV80EXUIw7yTDAYDiouLMXr0aADAggULEBYW\nBp1Ohw8++ACvvfYavvvuOwBAVlYW/vjHP+L48eP47rvvkJKSAuDGI3AB4MSJEzhx4gTGjh3bqh97\n2wKOHwlQVFSElJQU/Pe//8UjjzxiXZeamory8nI8/PDDSE9Ph9lstm6zZ88ebN26FeXl5a2OzO3V\nUl1djaeffhrPP/88ysvL8fLLL2Pu3Lk2E20QkXsxzDsoPT0dMTExSEtLw6RJk7Bw4UJUVVXhhx9+\nQGZmJnr16oWYmBjMmTMHO3fuBAD06tULFRUVMBqN6N27N+6++26bfbY3zOLl5dXutu0ZN24cpk2b\nBgDWqedGjx6N6dOnQ6FQYP78+bh27ZrNXwTp6ekICQlpc6o6e7Xs3r0b999/PxISEgAAkydPxl13\n3cXxdqJuxDDvoIKCApSVleHgwYPIzs6Gt7c3ampqEBgYaHMiNDw83DpVW25uLk6fPo34+HjMmDED\nX375pdP9bdiwodPbtjXtXMtHaMpkMoSGhtpMKWdvqrr2aqmsrMS+ffusU9VFR0dDq9Xi4sWLTtdK\nRF3T40+A9jRtHUUPHDgQ9fX1MJlM8PX1BXBjGCYkJAQAMGTIEOvsQ5999hnmz5+PY8eOOfXkRHvb\n3jpVndlsxqVLl2y2dTRVnSRJuHDhgrVOe9s4qkWtViMtLQ3r1693+H6IyD14ZO4CarUasbGxyMnJ\nwbVr11BeXo7t27dbT3Du3r3bOn7s7+8PAJDL5VCpVJDL5aioqLC7b3vbDhs2DNeuXUNxcTGam5uR\nl5eHpqYmh7UePXoU+/fvh9lsxqZNm+Dt7d3mWH1Havntb3+LL7/8Et9++y0sFguuXr2K77//npNI\nE3WjHn9kblarb1w+6Mb9O6u9o9Z3330XL730EsaOHYvAwEAsW7YMcXFxAICvv/4aq1evxtWrVxEe\nHo7333/fOib9pz/9CRqNBmazGVu3bm01Jm5vW29vb6xduxZLly6FxWLBggUL2h0iuWnatGn45z//\nicWLF2Po0KH48MMPrSc625uqrr1a1Go1CgoKkJ2djYyMDCiVSowZMwY5OTmOf6lE5BKcNu5XZOPG\njaioqMBbb73l0To4bZwtThtHzuK0cUREgmOYExEJoMePmZPrvPDCC54ugYjchEfm1O08c5aGSGwM\nc+p2DQ38347I1XrUp8pkMiEgIMDTZZAbVVersGwZR/eIXK1HfarMZjN++eUX6w0p5JyGBjmOHetR\n/5StSNKNOpctU+LcOcd3vhJRxziVADqdDps3b4YkSUhMTISmjeeLFxQUQKfTwdvbG8899xyGDBnS\nqYLMZjOvNe+g//s/X2g0fp4ug4g8yOEwi8ViQX5+PjIzM5Gbm4uSkhIYDAabNqWlpaipqcFbb72F\nefPm4cMPP3RbwURE1JrDMNfr9QgNDUVwcDCUSiXi4uKg1Wpt2mi1WsTHxwMAIiMjYTKZUF9f756K\niYioFYdhbjQaERQUZF1WqVStJh1wpg0REblPj7qahYiIOsfhCVCVSoW6ujrrstFohEqlatWm5bO0\nL1261KrNrdp7YAx1zMMPi34jTgqenCz0GwSkJz1dgVsE/u8/cj+HR+YRERGorq5GbW0tmpubUVJS\n0moi4djYWHz77bcAgJMnT8LPzw+BgfwnJCLqLk49Alen0+Gjjz6CJElISkqCRqNBUVERZDIZkpOT\nAQD5+fnQ6XTw8fHBggULMGzYMLcXT0REN3jseeZEROQ6PAFKRCQAhjkRkQB69gM9iG5x+fJlrFmz\nBjKZDD/99BPkcjkCAgIgk8mwdu1a63ymALB27Vq88MIL8PHxsbu/BQsWIDc3F76+vtZ1K1asgMVi\nQWNjI65fvw6VSgWZTIaXXnrJ4VVaRJ7CMXO6be3cuRM+Pj6YMWNGq9ckSWp3Au6b2grzm4qLi3H+\n/Hk8/fTTLqmXyJ14ZE63rZbHIdXV1Vi/fj2GDBmCs2fPIjMz0/o8IV9fX6xbtw719fW4fv06Hnro\nISQlJXW4vy+++AJ1dXV4/PHHAQD79+/HTz/9hPj4eLzxxhtQq9U4d+4chg4dioyMDCiVSpw6dQpb\nt27F9evX0bdvX2RkZPAxz+QWHDMnYVRVVSElJQW5ubmthkMWLlyInJwc5OTkYN++fTCZTB3e/333\n3YdDhw5Zv0S++eYbJCYmAgDOnz8PjUaDN954AwDw1VdfoampCVu2bMHy5cuRk5ODiRMnYseOHV18\nl0Rt45E5CWPgwIEYOnRom6/t3bsXP/zwA4AbdzFXV1d3+F4IX19fREZG4siRI+jXrx+8vb0REhKC\nqqoqhISEWPuePHky/vWvfyEyMhKVlZVYs2YNJEmCJEkIDg7u2psksoNhTsKwd6Lz6NGjOHHiBHJy\ncqBUKpGVlYWmpqZO9ZGUlIQvvvgCgYGB1qPytshkMkiShKFDhyIrK6tTfRF1BIdZSBj2zuWbTCb4\n+flBqVTi/PnzOH36dKf7iI6OhsFggFarxcSJE63rL1y4gIqKCgBASUkJRo0ahUGDBqG2thZnzpwB\nADQ3N6OysrLTfRO1h0fmJAx7V6+MHTsWX375JV588UWEhoYiMjKyS/3cc889qK2thbe3t3Xd4MGD\nUVhYiLNnz2LIkCFISkqCUqnE888/j48++ghXr16FJElITU1FeHh4l/onagsvTSTqoOzsbMyePRsj\nRowAcOPE69tvv42cnBwPV0a/ZhxmIXJSQ0MDFi9ejMDAQGuQE/UUPDInIhIAj8yJiATAMCciEgDD\nnIhIAAxzIiIBMMyJiATAMCciEsD/AyLFK7HnHMauAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "Loc_Acc = df.pivot_table(values='Loc_Acc', index='Sub#', columns='TrialType', aggfunc=np.mean)\n", "CIs = np.array(OneWayConfInterval(Loc_Acc))\n", "Analysis_and_Plot(Loc_Acc, CIs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Replicating Chen and Wyble, participants perform no worse in the surprise and post surprise trials, indicating that they succesfully found the target. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now lets look at participants' ability to report the target's color in the surprise trial and the trial immediately following the surprise test.\n", "\n", "Below I plot the percent of participants that correctly identified the target's color in the surprise and post-surprise trials" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "---------------------------------\n", "Mean Surprise: 0.3\n", "Mean Postsurprise: 0.7\n", "Postsurprise - Surprise: 0.4\n", "---------------------------------\n", "Postsurprise vs Surprise\n", "\t Surprise Test. Comparison to Chance: 17.0 P-value: 0.5899\n", "\t After Surprise Test. Comparison to Chance: 33.0 P-value: 0.024\n", "\t Chi-Square Comparison: 6.4 P-value: 0.0114\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAECCAYAAAAMxDf2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGZZJREFUeJzt3XtwVOX9x/HP2SwEEhPCQiQ3HS6JXGK1SEAhCCTG1kug\nOwUp2DoypepwsdBWcSg0CKVGlKj1OuIEqQpFCzQOOKNGvPRnZqRpSSiXEQwImg3JJCyUYMIl2fP7\ngx/7M+ayCWTZ5PH9+ouz+5znfHdJPvvkPGfPY9m2bQsA0K05Ql0AAODSEeYAYADCHAAMQJgDgAEI\ncwAwAGEOAAZwBmrw0ksvaefOnerTp49Wr17dYpu1a9eqtLRU4eHhmjdvngYOHNjZdQIA2hBwZJ6R\nkaElS5a0+nxJSYmqqqr07LPP6v7779crr7zSqQUCAAILGObDhg1TZGRkq88XFxdr4sSJkqSUlBTV\n1dXpxIkTnVchACCgSz5n7vV61a9fP/+2y+WS1+u91G4BAB3ABCgAGOCSw9zlcunYsWP+7WPHjsnl\ncl1qtwCADgh4NYsk2bat1u7HlZaWpvfee0/jxo3TgQMHFBkZqZiYmIB9VlRUdKxS4DKJiopSbW1t\nqMsAmklISGj1OSvQXRP//Oc/a9++faqtrVWfPn00ffp0NTQ0yLIsZWVlSZLy8/NVWlqqXr16ac6c\nORo8eHDAoghzdFWEObqqSwrzYCHM0VUR5uiq2gpzJkABwACEOQAYoF0ToAC6tsjISDkcjM1M4PP5\n9M0333R4P8Ic6OYiIyN15swZNTQ0hLoUdAKn06nIyMgOBzof5UA353A4CHKDNDQ0XNRfWYQ5ABiA\nMAcAAxDmAIzl8Xg0dOjQVr/BbhImQAHDeDzhqqgIC1r/CQmNSkw80+72//znP/WnP/1JBw4cUFhY\nmFJSUrR8+XJdd911QavxgsTERO3fvz/ox+kKCHPAMBUVYXK7A98f6WIVFJxQYmL72p46dUqzZs3S\n448/rsmTJ+vs2bPasWOHevbs2eHj2rYty7La3b6xsVFhYcH7UOtqOM0CIGgOHToky7I0ZcoUWZal\n8PBwTZgwQcOGDdNTTz2lBx980N+2vLxcSUlJ8vl8kqRp06Zp1apVcrvdSk5O1ldffaVp06YpNzdX\n2dnZGjZsmGbPnq3//ve/TfbfuHGjxowZo5/97GfN+nzzzTc1btw4DR06VOPGjVNBQYH/+Bs3btSk\nSZOUmpqqX/ziF/J4PJfxnbp0hDmAoBk8eLAcDocWLlyojz76yB+8F3x3pP3d7S1btmj16tU6cOCA\nEv/vz4HNmzfr6aefVklJiRwOh5YuXdpkn88++0yffPKJ1q9f36TP+vp6LVu2TOvXr9f+/fv19ttv\nKzU1VZL03nvv6fnnn1d+fr52796tMWPGaO7cuZ33RlwGhDmAoLniiiv097//XQ6HQ4sWLdJ1112n\nX/7yl6qpqWnX/tOnT1dycrIcDoeczvNnhadOnaqUlBT17t1bixYt0rZt2/wTnJZl6aGHHlLv3r0V\nHh7erL+wsDB9/vnnOn36tGJjY5WSkiJJeuONN/Tggw9qyJAhcjgcmj9/vvbu3dutRueEOYCgSk5O\n1lNPPaXi4mJ9+OGHqqqq0rJly9q1b0t3Cfz2Y0lJSTp37lyTpSrj4+Nb7Kt379566aWX9Nprr+mG\nG27Qvffeq4MHD0o6f4omJydHqampSk1N1bXXXivLslRZWdmRlxpShDmAy2bIkCG66667tH//fkVE\nRKi+vt7/XFVVVbP2LU14fvv22eXl5erZs2eT1c3amiSdMGGC/vrXv6qkpERDhgzRokWLJJ3/gFi1\napX27t2rvXv3at++ffriiy80atSoi3qdoUCYAwiasrIyvfzyyzp69Kik89d9FxQUaNSoURoxYoR2\n7Nghj8ejkydP6oUXXmhXn1u2bFFZWZnq6+u1evVq3Xnnnf4Ab+l68guP1dTU6P3331d9fb169Oih\nyMhI/3733HOPnnvuOR04cECSdPLkSW3btu2SX//lxKWJgGESEhpVUHAiqP231xVXXKGSkhKtWbNG\ntbW1io6O1q233qqlS5cqMjJSU6ZM0a233iqXy6V58+apsLDQv29rI+ypU6dq4cKFOnjwoMaOHavH\nH3+8zX0uPObz+bRmzRotXLhQlmVpxIgR/n1vu+021dXVae7cufJ4PIqKitKECROUnZ3d7tcaaqw0\nBHxHd1tpqLvVeymmTZumadOmacaMGaEuJaha+z9lpSEAMBxhDqDb6Mg3QL9vOGcOoNv429/+FuoS\nuixG5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAeAi3XPPPdq0aVOoy5DEpYmAcTx1HlV8E7xvWCdE\nJigxon1LDd14442qqamR0+lURESEJk2apMcee0y9e/fu8HHLy8t100036auvvpLD0TXGoa+//nqo\nS/AjzAHDVHxTIXeBO2j9F7gL2h3mlmXptddeU3p6uqqqqnT33XfrmWee0eLFizt83AvLxgXjDiQX\ns8RcR5exC7au8fEGwFgXwnfAgAHKyMjQ/v37VVVVpVmzZik1NVXjx4/Xhg0b/O1LS0t1xx13aNiw\nYRo5cqRWrFgh6fwNtiRp+PDhGjp0qHbu3KnDhw9r2rRpGj58uK677jr/6kDfXS5OOn9fl40bN0qS\n3nrrLbndbj366KO69tpr9dRTT/kfW7p0qYYPH65Jkybp008/bbJ/S8vYXeiztVqk83ePnDlzplJT\nUzVx4kRt3bq1099nRuYALguPx6MPP/xQd9xxh+bOnathw4bplVde0YEDBzRz5kwNHDhQ48aNU05O\njn71q1/ppz/9qerr6/X5559LOn/r27Fjx2r//v3+EfG8efM0ceJEbdq0SWfPntWuXbv8xws0ai4p\nKZHb7dZ//vMfnTt3Tm+//bZKSko0efJk7dmzR++8847uu+8+ffbZZ+rTp4+/hvXr12vw4MFNPigk\n6cknn2yxlvr6es2cOVOLFi3Shg0btG/fPs2cOVPDhw9XcnJyp72/jMwBBNXs2bOVmpqqqVOnaty4\ncbr77rv1r3/9S0uXLlWPHj2UmpqqmTNn+icSe/ToocOHD8vr9ap3794aOXJkk/6+fZrF6XSqvLxc\nR48eVc+ePTV69Oh21xUXF6dZs2bJ4XD4l5jr37+/Zs+erbCwME2ZMkVDhgzR9u3b/fu0tIxdoFoK\nCwt19dVX66677pJlWUpNTdXtt9/e6aNzwhxAUK1du1Z79+7VZ599ppUrV6qqqkoxMTFNJkGTkpL8\nS7Tl5eXp4MGDmjhxorKzs/XBBx+02vcf/vAH2bat7Oxs3XLLLXrzzTfbXVdLt5P97pJziYmJTVZA\nausWtK3V4vF4tHPnTv+SdCNGjFBBQYGqq6vbXWt7cJoFQFB9d8JywIABOnHihOrq6hQRESHpfODF\nxcVJkgYOHOhfdeidd97RAw88oD179rR42qR///564oknJEnFxcWaMWOGbrrpJkVFRUk6f4ojMjJS\nkpqFZ0v9XVgR6QKPx6Mf//jHbe4TqJaEhASNHTu2ybxAMDAyB3BZJSQkKC0tTbm5uTpz5oz27dun\njRs3+ic4t2zZ4l+g+UIoOxwOuVwuORwOHT582N/Xtm3b/AEcHR0ty7L8bePi4rR582b5fD5t3LhR\nR44cCVjbsWPHtHbtWjU0NGjr1q06ePCgbrnllna9rtZqycrK0qFDh7R582Y1NDTo3Llz2rVrl8rK\nytr9nrUHI3PAMAmRCSpwFwS1//ZqbST7wgsv6JFHHtENN9ygmJgYPfzww0pPT5ckffTRR1q+fLlO\nnz6tpKQkvfTSS/5z2r/+9a/ldrvV2NioN954Q7t27dKyZct06tQp9e/fX3/84x911VVXSTo/Ibl4\n8WKtWrVKM2bMUFpaWsB6R44cqS+//FI/+MEPFBsbqzVr1vgnP9takk5Sm7Vs2LBBjz76qJYvXy7b\ntjVixAgtW7as3e9je7Rr2bjS0lKtW7dOtm0rIyNDbnfTa1jr6ur03HPPqaamRj6fT5MnT9akSZPa\n7JNl49BVdbdl2LpbvV3VW2+9pY0bN2rLli2hLuWilo0LODL3+XzKz89XTk6O+vbtq8WLF2v06NFK\nTPz/Lw289957uuqqq/TII4/o5MmTWrhwoW6++eYOX4QPALg4Ac+Zl5WVKT4+XrGxsXI6nUpPT1dx\ncXGTNpZlqb6+XpJ0+vRpRUVFEeQAcBkFDHOv16t+/fr5t10ul39y4oLbbrtN5eXleuCBB/Twww9r\n1qxZnV4oAATT9OnTu8QplovVKROgpaWlGjRokJYtW6bKykqtXLlSq1evVq9evVrd58IsNdDV9OzZ\ns1v9fPJXsHnCwsI6/DMYMMxdLpdqamr8216vVy6Xq0mbjz/+2D8pGhcXpyuvvFIej0dDhgxptV8m\nbNBVdbcJxe70wYP2aWxsbPFnsK3/64CnWZKTk1VZWanq6mo1NDSoqKio2SU+/fv31+7duyVJJ06c\n0NGjRzVgwICO1g/gIvh8vmZfLUf35XQ6m933pT3afWniq6++Ktu2lZmZKbfbrcLCQlmWpaysLB0/\nflwvvviijh8/Lklyu90aP358m31yaSK6qu42MpekyMjILnOPb1wan8+nb775psXn2ro0sV1hHgyE\nObqq7hjm+H5oK8z5KAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBg\nAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYABWgQW+46D3oI6cOBLqMoyQEJmg\nxIjEUJfxvUCYA99RXlsud4E71GUYocBdQJhfJpxmAQADEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYg\nzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwADtup95aWmp1q1bJ9u2lZGR\nIbe7+b2e9+7dq7/85S9qbGxUdHS0li1b1unFAgBaFjDMfT6f8vPzlZOTo759+2rx4sUaPXq0EhP/\n/4bzdXV1ys/P19KlS+VyuXTy5MmgFg0AaCrgaZaysjLFx8crNjZWTqdT6enpKi4ubtLm008/1Y03\n3iiXyyVJio6ODk61AIAWBRyZe71e9evXz7/tcrlUVlbWpE1FRYUaGxu1fPlynT59WrfffrsmTJjQ\n+dUCAFrUKWuA+nw+ffnll8rJydGZM2e0dOlSXXPNNYqLi+uM7gEAAQQMc5fLpZqaGv+21+v1n075\ndpuoqCj17NlTPXv21PDhw3X48OE2wzwqKuoSygaCxzpuhboEY4SFhfG7fpkEDPPk5GRVVlaqurpa\nffv2VVFRkRYsWNCkzejRo7V27Vr5fD6dO3dOX3zxhbKzs9vst7a29tIqB4LEtu1Ql2CMxsZGftc7\nUVsfjAHD3OFwaPbs2Vq5cqVs21ZmZqaSkpJUWFgoy7KUlZWlxMREXX/99XrooYfkcDiUlZWlpKSk\nTn0RAIDWWXaIhiEVFRWhOCwQ0E7vTk3ePDnUZRihwF2g0bGjQ12GMRISElp9jm+AAoABCHMAMABh\nDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYA\nYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAG\nIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAA7Qrz0tJSLVy4UAsWLFBB\nQUGr7crKyjRz5kzt2LGj0woEAAQWMMx9Pp/y8/O1ZMkS5eXlqaioSB6Pp8V2GzZs0PXXXx+UQgEA\nrQsY5mVlZYqPj1dsbKycTqfS09NVXFzcrN27776rm266SdHR0UEpFADQuoBh7vV61a9fP/+2y+WS\n1+tt1qa4uFg/+tGPOr9CAEBAnTIBum7dOv385z/3b9u23RndAgDayRmogcvlUk1NjX/b6/XK5XI1\naXPo0CE988wzsm1btbW1KikpkdPpVFpaWqv9RkVFXULZQPBYx61Ql2CMsLAwftcvk4BhnpycrMrK\nSlVXV6tv374qKirSggULmrR5/vnn/f9+8cUXNWrUqDaDXJJqa2svsmQguPjLsvM0Njbyu96J2vpg\nDBjmDodDs2fP1sqVK2XbtjIzM5WUlKTCwkJZlqWsrKxOLRYA0HGWHaJhSEVFRSgOCwS007tTkzdP\nDnUZRihwF2h07OhQl2GMhISEVp/jG6AAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcA\nAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAwQcA1QdH0eT7gq\nKsJCXYYxzra+MhfQZRHmBqioCJPbHRPqMozx+v+EugKg4zjNAgAGIMwBwACEOQAYgDAHAAMQ5gBg\nAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwQLsWpygt\nLdW6detk27YyMjLkdrubPP/pp5/q7bffliT16tVL9913n66++urOrxYA0KKAI3Ofz6f8/HwtWbJE\neXl5KioqksfjadLmyiuv1PLly/Xkk09q6tSpevnll4NWMACguYBhXlZWpvj4eMXGxsrpdCo9PV3F\nxcVN2lxzzTWKiIiQJKWkpMjr9QanWgBAiwKGudfrVb9+/fzbLperzbDevn27fvjDH3ZOdQCAdunU\nBZ337Nmjjz/+WCtWrAjYNioqqjMP/b0WFhbqCoCWhYWF8bt+mQQMc5fLpZqaGv+21+uVy+Vq1u7I\nkSNas2aNfv/73+uKK64IeODa2toOlorWNDZGhLoEoEWNjY38rneitj4YA55mSU5OVmVlpaqrq9XQ\n0KCioiKlpaU1aVNTU6O8vDzNnz9fcXFxl14xAKBDAo7MHQ6HZs+erZUrV8q2bWVmZiopKUmFhYWy\nLEtZWVnatGmTTp06pfz8fNm2rbCwMOXm5l6O+gEAkizbtu1QHLiioiIUhzVScXGE3O6YUJdhjNf/\nZ5vu2T451GUYocBdoNGxo0NdhjESEhJafY5vgAKAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYA\nYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAG\nIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADC\nHAAMQJgDgAEIcwAwAGEOAAZwtqdRaWmp1q1bJ9u2lZGRIbfb3azN2rVrVVpaqvDwcM2bN08DBw7s\n7FoBAK0IODL3+XzKz8/XkiVLlJeXp6KiInk8niZtSkpKVFVVpWeffVb333+/XnnllaAVDABoLmCY\nl5WVKT4+XrGxsXI6nUpPT1dxcXGTNsXFxZo4caIkKSUlRXV1dTpx4kRwKgYANBMwzL1er/r16+ff\ndrlc8nq9HW4DAAgeJkABwAABJ0BdLpdqamr8216vVy6Xq1mbY8eO+bePHTvWrM13JSQkdLRWtOIn\nP5FsO9RVmCRbvxjPG4ruJeDIPDk5WZWVlaqurlZDQ4OKioqUlpbWpE1aWpo++eQTSdKBAwcUGRmp\nmJiY4FQMAGjGsu3AY7rS0lK9+uqrsm1bmZmZcrvdKiwslGVZysrKkiTl5+ertLRUvXr10pw5czR4\n8OCgFw8AOK9dYQ4A6NqYAAUAAxDmAGCAdn2dH+gqTp06pRUrVsiyLB0/flwOh0PR0dGyLEuPPfaY\nwsLC/G0fe+wx/fa3v1WvXr1a7W/OnDnKy8tTRESE/7HFixfL5/OptrZWZ8+elcvlkmVZeuSRRwJe\npQWECufM0W1t2rRJvXr1UnZ2drPnbNuWZVkB+2gpzC/48MMP9fXXX+vee+/tlHqBYGJkjm7r2+OQ\nyspKPfHEExo4cKCOHDmiJUuW+O8nFBERoVWrVunEiRM6e/as7rzzTmVmZnb4eO+//75qamp09913\nS5LeffddHT9+XBMnTtTTTz+thIQEffXVVxo0aJDmzp0rp9OpL774Qm+88YbOnj2rPn36aO7cuYqO\nju609wC4gHPmMEZFRYUmT56svLy8ZqdD5s+fr9zcXOXm5mrbtm2qq6vrcP/jx4/Xjh07/B8iH3/8\nsTIyMiRJX3/9tdxut55++mlJ0vbt23Xu3Dm99tprWrRokXJzczV27Fi99dZbl/gqgZYxMocxBgwY\noEGDBrX43NatW/Xvf/9b0vlvMVdWVnb4uxARERFKSUnRrl271LdvX4WHhysuLk4VFRWKi4vzH/vm\nm2/WP/7xD6WkpKi8vFwrVqyQbduybVuxsbGX9iKBVhDmMEZrE527d+/W/v37lZubK6fTqZycHJ07\nd+6ijpGZman3339fMTEx/lF5SyzLkm3bGjRokHJyci7qWEBHcJoFxmhtLr+urk6RkZFyOp36+uuv\ndfDgwYs+xogRI+TxeFRcXKyxY8f6Hz969KgOHz4sSSoqKtKwYcN01VVXqbq6WocOHZIkNTQ0qLy8\n/KKPDbSFkTmM0drVKzfccIM++OAD/e53v1N8fLxSUlIu6ThjxoxRdXW1wsPD/Y9dffXVKigo0JEj\nRzRw4EBlZmbK6XTqN7/5jV599VWdPn1atm1rypQpSkpKuqTjAy3h0kSgg1auXKnp06frmmuukXR+\n4vW5555Tbm5uiCvD9xmnWYB2OnnypBYsWKCYmBh/kANdBSNzADAAI3MAMABhDgAGIMwBwACEOQAY\ngDAHAAMQ5gBggP8FPFXvS7npkwQAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "Trial_Trimmer = df['TrialType2'] > 0\n", "Col_Acc = df[Trial_Trimmer].pivot_table(values='Col_Acc', index='Sub#', columns='TrialType2', aggfunc=np.mean)\n", "CIs = SimpleComparisonCI(Col_Acc)\n", "Analysis_and_Plot_2(Col_Acc, CIs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We perfectly replicate Chen and Wyble; participants respond more accurarely in the post-surprise trial than in the surprise trial. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The next cell examines participants' ability to report the target's identity on the surprise trial and the trial immediately following the surprise trial. Remember, the participants locate the target based on whether its a letter or number, so they know the broad category of the target. Nonetheless, they cannot report the target's identity on the surprise trial" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "---------------------------------\n", "Mean Surprise: 0.25\n", "Mean Postsurprise: 0.75\n", "Postsurprise - Surprise: 0.5\n", "---------------------------------\n", "Postsurprise vs Surprise\n", "\t Surprise Test. Comparison to Chance: 15.0 P-value: 0.7226\n", "\t After Surprise Test. Comparison to Chance: 35.0 P-value: 0.014\n", "\t Chi-Square Comparison: 10.0 P-value: 0.0016\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAECCAYAAAAMxDf2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGZFJREFUeJzt3X9wVNX9//HX3awEEhPCQiS/dPiRIBCLRQMKQSAxtv4I\nfHYKUrB1ZErVEbTQVnEoNAilRJSo9eeIE6SoFC3QMOCMGEHsh8xA05JQhREMCJINySQslGDCj2Tv\n5w++7NeYH5uQLElOn4+/uLtnz33vkrz25J5777Fs27YFAOjWHJ1dAACg/QhzADAAYQ4ABiDMAcAA\nhDkAGIAwBwADOAM1eOONN7R371717t1bK1eubLLN6tWrVVxcrNDQUM2ZM0cDBgzo6DoBAC0IODJP\nS0vTwoULm32+qKhIFRUVevnll/XII4/orbfe6tACAQCBBQzzoUOHKjw8vNnnCwsLNWHCBElSUlKS\nampqdPr06Y6rEAAQULuPmXu9XvXt29e/7XK55PV629stAKANmAAFAAO0O8xdLpdOnjzp3z558qRc\nLld7uwUAtEHAs1kkybZtNXc/rpSUFG3btk1jx47VoUOHFB4erqioqIB9lpWVta1S4CqJiIhQdXV1\nZ5cBNBIXF9fsc1aguyb+6U9/0oEDB1RdXa3evXtr2rRpqqurk2VZysjIkCTl5uaquLhYPXv21GOP\nPaZBgwYFLIowR1dFmKOraleYBwthjq6KMEdX1VKYMwEKAAYgzAHAAK2aAAXQtYWHh8vhYGxmAp/P\np2+//bbNryPMgW4uPDxc58+fV11dXWeXgg7gdDoVHh7e5kDnqxzo5hwOB0FukLq6uiv6K4swBwAD\nEOYAYADCHICxPB6PbrzxxmavYDcJE6CAYTyeUJWVhQSt/7i4esXHn291+3/84x/64x//qEOHDikk\nJERJSUlasmSJRowYEbQaL4uPj9fBgweDvp+ugDAHDFNWFiK3O/D9ka5UXt5pxce3ru3Zs2c1c+ZM\nPfvss5o0aZIuXLigPXv2qEePHm3er23bsiyr1e3r6+sVEhK8L7WuhsMsAILmyJEjsixLkydPlmVZ\nCg0N1fjx4zV06FC98MILeuKJJ/xtS0tLlZCQIJ/PJ0maOnWqVqxYIbfbrcTERH3zzTeaOnWqsrOz\nlZmZqaFDh2rWrFn6z3/+0+D169ev1+jRo/XTn/60UZ/vv/++xo4dqxtvvFFjx45VXl6ef//r16/X\nxIkTlZycrJ///OfyeDxX8ZNqP8IcQNAMGjRIDodD8+bN06effuoP3su+P9L+/vamTZu0cuVKHTp0\nSPH/78+BjRs36sUXX1RRUZEcDocWLVrU4DW7d+/WZ599pvfee69Bn7W1tVq8eLHee+89HTx4UJs3\nb1ZycrIkadu2bXr11VeVm5urzz//XKNHj9bs2bM77oO4CghzAEFz7bXX6m9/+5scDofmz5+vESNG\n6Be/+IWqqqpa9fpp06YpMTFRDodDTuelo8JTpkxRUlKSevXqpfnz52vr1q3+CU7LsvTkk0+qV69e\nCg0NbdRfSEiIvvzyS507d07R0dFKSkqSJL377rt64oknNHjwYDkcDj3++OPav39/txqdE+YAgiox\nMVEvvPCCCgsLtWPHDlVUVGjx4sWtem1Tdwn87mMJCQm6ePFig6UqY2Njm+yrV69eeuONN7R27Vrd\ncssteuihh3T48GFJlw7RZGVlKTk5WcnJybrppptkWZbKy8vb8lY7FWEO4KoZPHiw7r//fh08eFBh\nYWGqra31P1dRUdGofVMTnt+9fXZpaal69OjRYHWzliZJx48fr7/85S8qKirS4MGDNX/+fEmXviBW\nrFih/fv3a//+/Tpw4IC++uor3XrrrVf0PjsDYQ4gaEpKSvTmm2/qxIkTki6d952Xl6dbb71Vw4cP\n1549e+TxeHTmzBm99tprrepz06ZNKikpUW1trVauXKn77rvPH+BNnU9++bGqqip9/PHHqq2t1TXX\nXKPw8HD/6x588EG98sorOnTokCTpzJkz2rp1a7vf/9XEqYmAYeLi6pWXdzqo/bfWtddeq6KiIq1a\ntUrV1dWKjIzUXXfdpUWLFik8PFyTJ0/WXXfdJZfLpTlz5ig/P9//2uZG2FOmTNG8efN0+PBhjRkz\nRs8++2yLr7n8mM/n06pVqzRv3jxZlqXhw4f7X3v33XerpqZGs2fPlsfjUUREhMaPH6/MzMxWv9fO\nxkpDwPd0t5WGulu97TF16lRNnTpV06dP7+xSgqq5/1NWGgIAwxHmALqNtlwB+t+GY+YAuo2//vWv\nnV1Cl8XIHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5AFyhBx98UBs2bOjsMiRxaiJgHE+NR2XfBu8K\n67jwOMWHtW6podtuu01VVVVyOp0KCwvTxIkTtXz5cvXq1avN+y0tLdXtt9+ub775Rg5H1xiHvvPO\nO51dgh9hDhim7NsyufPcQes/z53X6jC3LEtr165VamqqKioq9MADD+ill17SggUL2rzfy8vGBeMO\nJFeyxFxbl7ELtq7x9QbAWJfDt3///kpLS9PBgwdVUVGhmTNnKjk5WePGjdO6dev87YuLi3Xvvfdq\n6NChGjlypJYuXSrp0g22JGnYsGG68cYbtXfvXh09elRTp07VsGHDNGLECP/qQN9fLk66dF+X9evX\nS5I++OADud1uPfPMM7rpppv0wgsv+B9btGiRhg0bpokTJ2rXrl0NXt/UMnaX+2yuFunS3SNnzJih\n5ORkTZgwQVu2bOnwz5mROYCrwuPxaMeOHbr33ns1e/ZsDR06VG+99ZYOHTqkGTNmaMCAARo7dqyy\nsrL0y1/+Uj/5yU9UW1urL7/8UtKlW9+OGTNGBw8e9I+I58yZowkTJmjDhg26cOGC9u3b599foFFz\nUVGR3G63/v3vf+vixYvavHmzioqKNGnSJH3xxRf68MMP9fDDD2v37t3q3bu3v4b33ntPgwYNavBF\nIUnPP/98k7XU1tZqxowZmj9/vtatW6cDBw5oxowZGjZsmBITEzvs82VkDiCoZs2apeTkZE2ZMkVj\nx47VAw88oH/+859atGiRrrnmGiUnJ2vGjBn+icRrrrlGR48eldfrVa9evTRy5MgG/X33MIvT6VRp\naalOnDihHj16aNSoUa2uKyYmRjNnzpTD4fAvMdevXz/NmjVLISEhmjx5sgYPHqzt27f7X9PUMnaB\nasnPz9cNN9yg+++/X5ZlKTk5Wffcc0+Hj84JcwBBtXr1au3fv1+7d+/WsmXLVFFRoaioqAaToAkJ\nCf4l2nJycnT48GFNmDBBmZmZ+uSTT5rt+/e//71s21ZmZqbuvPNOvf/++62uq6nbyX5/ybn4+PgG\nKyC1dAva5mrxeDzau3evf0m64cOHKy8vT5WVla2utTU4zAIgqL4/Ydm/f3+dPn1aNTU1CgsLk3Qp\n8GJiYiRJAwYM8K869OGHH+rRRx/VF1980eRhk379+um5556TJBUWFmr69Om6/fbbFRERIenSIY7w\n8HBJahSeTfV3eUWkyzwej3784x+3+JpAtcTFxWnMmDEN5gWCgZE5gKsqLi5OKSkpys7O1vnz53Xg\nwAGtX7/eP8G5adMm/wLNl0PZ4XDI5XLJ4XDo6NGj/r62bt3qD+DIyEhZluVvGxMTo40bN8rn82n9\n+vU6duxYwNpOnjyp1atXq66uTlu2bNHhw4d15513tup9NVdLRkaGjhw5oo0bN6qurk4XL17Uvn37\nVFJS0urPrDUYmQOGiQuPU547L6j9t1ZzI9nXXntNTz/9tG655RZFRUXpqaeeUmpqqiTp008/1ZIl\nS3Tu3DklJCTojTfe8B/T/tWvfiW32636+nq9++672rdvnxYvXqyzZ8+qX79++sMf/qDrr79e0qUJ\nyQULFmjFihWaPn26UlJSAtY7cuRIff311/rBD36g6OhorVq1yj/52dKSdJJarGXdunV65plntGTJ\nEtm2reHDh2vx4sWt/hxbo1XLxhUXF2vNmjWybVtpaWlyuxuew1pTU6NXXnlFVVVV8vl8mjRpkiZO\nnNhinywbh66quy3D1t3q7ao++OADrV+/Xps2bersUq5o2biAI3Ofz6fc3FxlZWWpT58+WrBggUaN\nGqX4+P9/0cC2bdt0/fXX6+mnn9aZM2c0b9483XHHHW0+CR8AcGUCHjMvKSlRbGysoqOj5XQ6lZqa\nqsLCwgZtLMtSbW2tJOncuXOKiIggyAHgKgoY5l6vV3379vVvu1wu/+TEZXfffbdKS0v16KOP6qmn\nntLMmTM7vFAACKZp06Z1iUMsV6pDJkCLi4s1cOBALV68WOXl5Vq2bJlWrlypnj17Nvuay7PUQFfT\no0ePbvXzyV/B5gkJCWnzz2DAMHe5XKqqqvJve71euVyuBm127tzpnxSNiYnRddddJ4/Ho8GDBzfb\nLxM26Kq624Rid/riQevU19c3+TPY0v91wMMsiYmJKi8vV2Vlperq6lRQUNDoFJ9+/frp888/lySd\nPn1aJ06cUP/+/dtaP4Ar4PP5Gl1aju7L6XQ2uu9La7T61MS3335btm0rPT1dbrdb+fn5sixLGRkZ\nOnXqlF5//XWdOnVKkuR2uzVu3LgW++TURHRV3W1kLknh4eFd5h7faB+fz6dvv/22yedaOjWxVWEe\nDIQ5uqruGOb479BSmPNVDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAzAlQbA9xz2Htax04EXMkBg\nceFxig+LD9wQ7UaYA99TWl0qd547cEMElOfOI8yvEg6zAIABCHMAMABhDgAGIMwBwACEOQAYgDAH\nAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAw\nAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGcLamUXFxsdasWSPbtpWW\nlia3292ozf79+/XnP/9Z9fX1ioyM1OLFizu8WABA0wKGuc/nU25urrKystSnTx8tWLBAo0aNUnx8\nvL9NTU2NcnNztWjRIrlcLp05cyaoRQMAGgp4mKWkpESxsbGKjo6W0+lUamqqCgsLG7TZtWuXbrvt\nNrlcLklSZGRkcKoFADQp4Mjc6/Wqb9++/m2Xy6WSkpIGbcrKylRfX68lS5bo3LlzuueeezR+/PiO\nrxYA0KRWHTMPxOfz6euvv1ZWVpbOnz+vRYsWaciQIYqJiemI7gEAAQQMc5fLpaqqKv+21+v1H075\nbpuIiAj16NFDPXr00LBhw3T06NEWwzwiIqIdZQPBY52yOrsEY4SEhPC7fpUEDPPExESVl5ersrJS\nffr0UUFBgebOndugzahRo7R69Wr5fD5dvHhRX331lTIzM1vst7q6un2VA0Fi23Znl2CM+vp6ftc7\nUEtfjAHD3OFwaNasWVq2bJls21Z6eroSEhKUn58vy7KUkZGh+Ph43XzzzXryySflcDiUkZGhhISE\nDn0TAIDmWXYnDUPKyso6Y7dAQHu9ezVp46TOLsMIee48jYoe1dllGCMuLq7Z57gCFAAMQJgDgAEI\ncwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAH\nAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAw\nAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGaFWYFxcXa968eZo7d67y\n8vKabVdSUqIZM2Zoz549HVYgACCwgGHu8/mUm5urhQsXKicnRwUFBfJ4PE22W7dunW6++eagFAoA\naF7AMC8pKVFsbKyio6PldDqVmpqqwsLCRu0++ugj3X777YqMjAxKoQCA5gUMc6/Xq759+/q3XS6X\nvF5vozaFhYX60Y9+1PEVAgAC6pAJ0DVr1uhnP/uZf9u27Y7oFgDQSs5ADVwul6qqqvzbXq9XLper\nQZsjR47opZdekm3bqq6uVlFRkZxOp1JSUprtNyIioh1lA8FjnbI6uwRjhISE8Lt+lQQM88TERJWX\nl6uyslJ9+vRRQUGB5s6d26DNq6++6v/366+/rltvvbXFIJek6urqKywZCC7+suw49fX1/K53oJa+\nGAOGucPh0KxZs7Rs2TLZtq309HQlJCQoPz9flmUpIyOjQ4sFALSdZXfSMKSsrKwzdgsEtNe7V5M2\nTursMoyQ587TqOhRnV2GMeLi4pp9jitAAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwAGEO\nAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBg\nAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAM4OzsAtB+Hk+oyspCOrsM\nY1yI6+wKgLYjzA1QVhYitzuqs8swxjv/29kVAG3HYRYAMABhDgAGIMwBwACEOQAYgDAHAAO06myW\n4uJirVmzRrZtKy0tTW63u8Hzu3bt0ubNmyVJPXv21MMPP6wbbrih46sFADQp4Mjc5/MpNzdXCxcu\nVE5OjgoKCuTxeBq0ue6667RkyRI9//zzmjJlit58882gFQwAaCxgmJeUlCg2NlbR0dFyOp1KTU1V\nYWFhgzZDhgxRWFiYJCkpKUlerzc41QIAmhQwzL1er/r27evfdrlcLYb19u3b9cMf/rBjqgMAtEqH\nXgH6xRdfaOfOnVq6dGnAthERER256/9qIVzJjy4qJCSE3/WrJGCYu1wuVVVV+be9Xq9cLlejdseO\nHdOqVav0u9/9Ttdee23AHVdXV7exVDSnvj6ss0sAmlRfX8/vegdq6Ysx4GGWxMRElZeXq7KyUnV1\ndSooKFBKSkqDNlVVVcrJydHjjz+umJiY9lcMAGiTgCNzh8OhWbNmadmyZbJtW+np6UpISFB+fr4s\ny1JGRoY2bNigs2fPKjc3V7ZtKyQkRNnZ2VejfgCAJMu2bbszdlxWVtYZuzVSYWEYd03sQO/871Y9\nuH1SZ5dhhDx3nkZFj+rsMowRF9f8/Zm5AhQADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ\n5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwAGEO\nAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBg\nAMIcAAxAmAOAAQhzADCAszWNiouLtWbNGtm2rbS0NLnd7kZtVq9ereLiYoWGhmrOnDkaMGBAR9cK\nAGhGwJG5z+dTbm6uFi5cqJycHBUUFMjj8TRoU1RUpIqKCr388st65JFH9NZbbwWtYABAYwHDvKSk\nRLGxsYqOjpbT6VRqaqoKCwsbtCksLNSECRMkSUlJSaqpqdHp06eDUzEAoJGAYe71etW3b1//tsvl\nktfrbXMbAEDwMAEKAAYIOAHqcrlUVVXl3/Z6vXK5XI3anDx50r998uTJRm2+Ly4urq21ohn/8z+S\nbXd2FSbJ1M/H8YGiewk4Mk9MTFR5ebkqKytVV1engoICpaSkNGiTkpKizz77TJJ06NAhhYeHKyoq\nKjgVAwAasWw78JiuuLhYb7/9tmzbVnp6utxut/Lz82VZljIyMiRJubm5Ki4uVs+ePfXYY49p0KBB\nQS8eAHBJq8IcANC1MQEKAAYgzAHAAK26nB/oKs6ePaulS5fKsiydOnVKDodDkZGRsixLy5cvV0hI\niL/t8uXL9Zvf/EY9e/Zstr/HHntMOTk5CgsL8z+2YMEC+Xw+VVdX68KFC3K5XLIsS08//XTAs7SA\nzsIxc3RbGzZsUM+ePZWZmdnoOdu2ZVlWwD6aCvPLduzYoePHj+uhhx7qkHqBYGJkjm7ru+OQ8vJy\nPffccxowYICOHTumhQsX+u8nFBYWphUrVuj06dO6cOGC7rvvPqWnp7d5fx9//LGqqqr0wAMPSJI+\n+ugjnTp1ShMmTNCLL76ouLg4ffPNNxo4cKBmz54tp9Opr776Su+++64uXLig3r17a/bs2YqMjOyw\nzwC4jGPmMEZZWZkmTZqknJycRodDHn/8cWVnZys7O1tbt25VTU1Nm/sfN26c9uzZ4/8S2blzp9LS\n0iRJx48fl9vt1osvvihJ2r59uy5evKi1a9dq/vz5ys7O1pgxY/TBBx+0810CTWNkDmP0799fAwcO\nbPK5LVu26F//+pekS1cxl5eXt/laiLCwMCUlJWnfvn3q06ePQkNDFRMTo7KyMsXExPj3fccdd+jv\nf/+7kpKSVFpaqqVLl8q2bdm2rejo6Pa9SaAZhDmM0dxE5+eff66DBw8qOztbTqdTWVlZunjx4hXt\nIz09XR9//LGioqL8o/KmWJYl27Y1cOBAZWVlXdG+gLbgMAuM0dxcfk1NjcLDw+V0OnX8+HEdPnz4\nivcxfPhweTweFRYWasyYMf7HT5w4oaNHj0qSCgoKNHToUF1//fWqrKzUkSNHJEl1dXUqLS294n0D\nLWFkDmM0d/bKLbfcok8++US//e1vFRsbq6SkpHbtZ/To0aqsrFRoaKj/sRtuuEF5eXk6duyYBgwY\noPT0dDmdTv3617/W22+/rXPnzsm2bU2ePFkJCQnt2j/QFE5NBNpo2bJlmjZtmoYMGSLp0sTrK6+8\nouzs7E6uDP/NOMwCtNKZM2c0d+5cRUVF+YMc6CoYmQOAARiZA4ABCHMAMABhDgAGIMwBwACEOQAY\ngDAHAAP8HyIf70sQsw2rAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "Trial_Trimmer = df['TrialType2'] > 0\n", "Iden_Acc = df[Trial_Trimmer].pivot_table(values='Iden_Acc', index='Sub#', columns='TrialType2', aggfunc=np.mean)\n", "CIs = SimpleComparisonCI(Iden_Acc)\n", "Analysis_and_Plot_2(Iden_Acc, CIs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Experiment 1 - Intertrial analyses" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "So far, I've perfectly replicated Chen & Wyble (which is good since this is their data). \n", "\n", "Now I want to see if the target's color or identity on the previous trial influences the current trial's performance in the location task. I am only examining presurprise trials, so this should be trials when the participants don't \"remember\" the target's color or identity.\n", "\n", "First I want to make some variables representing whether the target's color and identity repeat across trials." ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": true }, "outputs": [], "source": [ "df['Prev_TarCol'] = df['TarCol'].shift(periods=1)\n", "df['Prev_TarCol_match'] = np.where(df['Prev_TarCol']==df['TarCol'], 1, 2)\n", "df['Prev_Iden'] = df['Tar_Iden'].shift(periods=1)\n", "df['Prev_Iden_match'] = np.where(df['Prev_Iden']==df['Tar_Iden'], 1, 2)\n", "df['Prev_Col+Iden_match'] = np.where((df['Prev_Iden_match']==1) & (df['Prev_TarCol_match']==1), 1, 2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Lets see what happens when the target's color and identity repeat. " ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "---------------------------------\n", "Mean match: 0.918\n", "Mean mismatch: 0.8925\n", "Match - Mismatch: 0.0255\n", "---------------------------------\n", "Match vs Mismatch\n", "\t normality test adj. Test value: 0.92 P-value: 0.0821\n", "\t Wilcoxon. Test value: 51.0 P-value: 0.04\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAECCAYAAAAMxDf2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF/9JREFUeJzt3X1wVNX9x/HPJguExBCyMZBHBJvwEKQiJhiMFpMJtEwx\nvyhiRX4MrVQpYAu0UqVgGGmGtNWAQKtVG0gVbSc6NhanKpGKSrR0BxILUpGoPCUGEhYkEHnI7v39\nwc8d1jxsAhuSHN+vv3L3ntz73TD72cO5595jsyzLEgCgRwvq6gIAAJeOMAcAAxDmAGAAwhwADECY\nA4ABCHMAMIDdX4Mnn3xSO3bsUEREhB577LEW26xbt06VlZXq06eP5s2bp8GDBwe6TgBAG/z2zDMz\nM7VkyZJW91dUVOjw4cNas2aN7rvvPj3zzDMBLRAA4J/fMB8+fLjCwsJa3e90OjV+/HhJUnJyshob\nG3X8+PHAVQgA8OuSx8xdLpeioqK82w6HQy6X61IPCwDoAC6AAoABLjnMHQ6Hjh496t0+evSoHA7H\npR4WANABfmezSJJlWWrteVypqal64403dOONN+rjjz9WWFiY+vfv7/eYNTU1HasUuEzCw8PV0NDQ\n1WUAzcTFxbW6z+bvqYmrV6/W7t271dDQoIiICN15551qamqSzWZTdna2JKmoqEiVlZUKCQnRnDlz\ndPXVV/stijBHd0WYo7u6pDDvLIQ5uivCHN1VW2HOBVAAMABhDgAGaNcFUABmCwsLU1AQfbvuwOPx\n6NSpUx3+PcIc+IYLCwvTmTNn1NTU1NWlQJLdbldYWFiHA52vYuAbLigoiCDvRpqami7qf0mEOQAY\ngDAHAAMQ5gDw/9LT07V169auLuOicAEUgI/q6j6qqQnutOPHxbkVH3+mXW1vuOEG1dXVafv27YqM\njPS+PnHiRO3evVvbtm1TfHx8q79/6NAhpaen68CBA8bP1iHMAfioqQlWbq7/5ytdrNLS42ojf33Y\nbDYlJibqlVde0Q9/+ENJ0kcffaTTp0/LZrP5/X3LsmSz2Vp9tpRJzP6qQpfbu3evXnvtNb322mta\nu3at9+e9e/d2dWnoIaZMmaIXX3zRu/3iiy9q6tSp3u3Nmzfru9/9roYPH66xY8dq5cqVPr8rSSNG\njNCwYcO0Y8cOSdLzzz+vW265RcOGDVNWVpZ27drl/Z1du3YpOztbKSkpmjt3rs6ePdvZbzEgCHN0\nquTkZE2aNEmTJk3SJ5984v05OTm5q0tDDzFmzBidPHlSVVVV8ng8+vvf/67bb7/d29sOCwvTmjVr\n9NFHH+nZZ5/Vc889p02bNkmSXn75ZUnSnj17tGfPHo0ZM0YbN27UqlWrtHbtWu3Zs0fr16/3GcJ5\n9dVX9Ze//EXvv/++du/erZKSksv/pi8CwywAur2veufjxo1TcnKyYmJivPvS09O9Pw8fPlw5OTl6\n//33NXHiRO/rXw23SNJf//pXzZ07V6NGjZIkXXXVVT7n+vGPf6zo6GhJ0oQJE/Thhx922vsKJMIc\nQLc3ZcoU3X777Tp48KDuuOMOn307duxQQUGB9uzZo3Pnzuns2bOaPHlyq8eqqalpFuAXuvLKK70/\n9+3bV4cPH770N3AZEOYG6OzZB4HkdIZ2dQl+XXWVpQEDuroKXCg+Pl6JiYl66623VFhYKEnenvZP\nf/pT3XPPPXrhhRfUq1cvLVu2TMeOHfNpc6G4uDjt37//8hV/mRDmBujs2QeBMnWqekSdGzc2EObd\n0MqVK/XFF1+ob9++crvd3jHzU6dOKSIiQr169VJFRYVKS0s1fvx4SeeXtQwKCtK+ffu8i+ZMmzZN\ny5cvV1pamkaNGqV9+/apV69ebU5x7AkIc6CH2Lt3r6qqqiRJVVVVSkpKkiQlJSUF9IJyXJxbpaXH\nA3a8lo7fXhf2rAcNGtTivhUrVuiRRx7R0qVLlZ6erpycHH3xxReSzg+T/OxnP1Nubq7cbrc2bNig\nyZMn6/jx45o3b54OHz6sxMRErV69WvHx8e2a7thdsdKQAZzO0B7R45069U69+GL3nxmwcWODxozp\n3isNLViwQI8//nhAjsXKSt1Pa/8mba00RM8cnSo6+kNFR5+fDXDkyAilpJwP87q6kaqrG9mVpQFG\nIczRqQht4PIgzIGvCb7yEznruv9sB2edMyDHSe2TKpt67lgxziPMga852nRIM0pzu7qMNk3VVOUG\nqMbX73pdo8JHBeRY6Drczg8ABiDMAcAADLMAPUT0qWhFnzr/zJAjoUeUciRFklQXVqe6sLquLA3d\nAGEO9BCENtrCMAuAHu2hhx7S6tWru7qMVpWUlOi2227r9PPQMwfgo7qxWjWnOu8O7biwOMWHtu85\nKO1ZNu43v/lNZ5XaqoULFyouLk6LFi1qV/vL8ZgAwhyAj5pTNQGb9tiS0tzSdof5pS4b903CMAuA\nbs3fsnELFy7Uo48+KklyuVyaOXOmUlJSNHLkSO+ycdL5RSz++Mc/Kjs7W0OHDtWiRYtUX1+vGTNm\naNiwYZo2bZpOnDjhbT979mxdd911SklJ0R133OFd6vD555/X3/72Nz355JMaNmyYfvSjH0k6/7yp\ne++9V9/+9rc1atQoPfzww95jWZalX//61xo5cqRuvPFGvfXWWwH/OxHmALq11paNa8lTTz2luLg4\n7dq1S//5z3/00EMP+ez/xz/+oZKSEr377rvatGmTZsyYocWLF2vnzp1yu90qKiryts3KytJ7772n\nDz74QNdcc43mzZsnSZo+fbpuu+02zZkzx7vsnMfj0cyZM5WYmKh///vf2r59u3JycrzHqqioUHJy\nsnbt2qWf/OQneuCBBwL+dyLMAXR7X/XO33nnHe+ycS098LVXr146cuSIDhw4oODgYKWlpfnsv+ee\ne+RwODRw4ECNHTvW2/Pu3bu3Jk2a5LNE3A9+8AP17dtXvXr10sKFC7V7926dPHmyxfoqKip05MgR\nLV26VCEhIerdu7fPuRMSEnTXXXfJZrNp6tSpOnLkiOrr6wP01zmPMAfQ7U2ZMkWlpaUqKSnxLhvX\n0pj5nDlzdNVVV+nuu+9WRkaG/vCHP/jsv3BJuJCQEO9an19tnzp1SpLk8Xi0YsUKZWRkaMSIERo3\nbpxsNptcLleL9dXU1CghIUFBQS1H6oALVjvp27evLMvynitQCHMA3d6Fy8ZNmjSp1XZhYWHKy8vT\ne++9p/Xr1+vpp59WeXl5h8/38ssvq6ysTCUlJfrvf/+rf/3rX7Isy/u/ga9/kcTFxam6uloej6fD\n5woUwhxAj7By5UqVlJSob9++ktTiMMubb76pffv2STof7Ha7XcHBHV8f99SpU+rdu7ciIiLU2Nio\ngoICnwCPjo7WgQMHvNvXXXedBgwYoBUrVujLL7/UmTNn5HQG5qmW7cXURAA+4sLiVJpb2qnHb6/2\nLBt3oc8++0xLly6Vy+VSRESEZs6cqfT09BbbtzW1cerUqXr77bd1/fXXKzIyUosWLdKGDRu8+++6\n6y7Nnj1bI0eO1Lhx4/SnP/1JxcXFevjhh5WWlqagoCDl5uY2G7Nvz7kvVruWjausrFRxcbEsy1Jm\nZqZyc33noDY2Nmrt2rWqr6+Xx+PRrbfeqltuuaXNY7JsXOD0lGXjeorn3n1VMzbf2tVlXDY8Arf7\n6ZRl4zwej4qKipSXl6fIyEgtXrxYaWlpPitZv/HGG0pMTNSDDz6oEydOaMGCBbr55psv6r83AICO\n8ztmXlVVpdjYWEVHR8tutysjI6PZWJDNZtOXX34pSTp9+rTCw8MJcgC4jPyGucvlUlRUlHfb4XA0\nm57zve99T4cOHdLs2bO1aNEi7223AIDLIyAXQCsrKzVkyBAtW7ZMtbW1ys/P12OPPaaQkJBWfyc8\nPDwQp4Yk/hMEmCU4OLjDGek3zB0Oh8+dSi6XSw6Hw6fNli1bvBdFY2JiNGDAAFVXV+tb3/pWq8dt\naXAfF8ftDu3qEtCDtWMOBC4zt9vdYka2FfB+h1mSkpJUW1ururo6NTU1qby8XKmpqT5trrzySu3c\nuVOSdPz4cX3++ecaOHBgR+sH0AV2H90tj63rbnaBL7vdflE3H/ntmQcFBWnWrFnKz8+XZVnKyspS\nQkKCysrKZLPZlJ2drSlTpuiJJ57wPjxm+vTpuuKKKzr+LgBcdg+8+YCULaVEpQR8/vM1V16jfr37\nBfSYpvN4PBd1q3+75pl3BuaZBw7zzAPrmzbPvDOV5pYqLbrlG2fQcW3NM+d2fgAwAGEOAAYgzAHA\nAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxA\nmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5\nABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwAD29jSqrKxUcXGxLMtSZmamcnNzm7X5\n8MMP9ec//1lut1v9+vXTsmXLAl4sAKBlfsPc4/GoqKhIeXl5ioyM1OLFi5WWlqb4+Hhvm8bGRhUV\nFWnp0qVyOBw6ceJEpxYNAPDld5ilqqpKsbGxio6Olt1uV0ZGhpxOp0+brVu36oYbbpDD4ZAk9evX\nr3OqBQC0yG/P3OVyKSoqyrvtcDhUVVXl06ampkZut1uPPPKITp8+rUmTJuk73/lO4KsFALSoXWPm\n/ng8Hn322WfKy8vTmTNntHTpUg0dOlQxMTGBODwAwA+/Ye5wOFRfX+/ddrlc3uGUC9uEh4erd+/e\n6t27t0aMGKF9+/a1Gebh4eGXUDYuFBzc1RUALQsODuazfpn4DfOkpCTV1taqrq5OkZGRKi8v1/z5\n833apKWlad26dfJ4PDp37pz27t2ryZMnt3nchoaGS6scXm53aFeXALTI7XbzWQ+gtr4Y/YZ5UFCQ\nZs2apfz8fFmWpaysLCUkJKisrEw2m03Z2dmKj4/XtddeqwceeEBBQUHKzs5WQkJCQN8EAKB17Roz\nHz16tFavXu3z2oQJE3y2c3JylJOTE7jKAADtxh2gAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACE\nOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgD\ngAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAY\ngDAHAAMQ5gBgAMIcAAxAmAOAAQhzADBAu8K8srJSCxYs0Pz581VaWtpqu6qqKk2bNk3btm0LWIEA\nAP/8hrnH41FRUZGWLFmiwsJClZeXq7q6usV2L7zwgq699tpOKRQA0Dq/YV5VVaXY2FhFR0fLbrcr\nIyNDTqezWbvXX39d6enp6tevX6cUCgBond8wd7lcioqK8m47HA65XK5mbZxOpyZOnBj4CgEAfgXk\nAmhxcbGmT5/u3bYsKxCHBQC0k91fA4fDofr6eu+2y+WSw+HwafPpp5/q8ccfl2VZamhoUEVFhex2\nu1JTU1s9bnh4+CWUjQsFB3d1BUDLgoOD+axfJn7DPCkpSbW1taqrq1NkZKTKy8s1f/58nza///3v\nvT8/8cQTuv7669sMcklqaGi4yJLxdW53aFeXALTI7XbzWQ+gtr4Y/YZ5UFCQZs2apfz8fFmWpays\nLCUkJKisrEw2m03Z2dkBLRYA0HF+w1ySRo8erdWrV/u8NmHChBbbzp0799KrAgB0CHeAAoABCHMA\nMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwAD\nEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABh\nDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAA9vY0qqysVHFxsSzL\nUmZmpnJzc332b926Va+88ookKSQkRPfee68GDRoU+GoBAC3y2zP3eDwqKirSkiVLVFhYqPLyclVX\nV/u0GTBggB555BE9+uijmjJlip566qlOKxgA0JzfMK+qqlJsbKyio6Nlt9uVkZEhp9Pp02bo0KEK\nDQ2VJCUnJ8vlcnVOtQCAFvkNc5fLpaioKO+2w+FoM6w3b96s0aNHB6Y6AEC7tGvMvL127dqlLVu2\naPny5X7bhoeHB/LU32jBwV1dAdCy4OBgPuuXid8wdzgcqq+v9267XC45HI5m7fbv36+nn35av/rV\nr3TFFVf4PXFDQ0MHS0Vr3O7Qri4BaJHb7eazHkBtfTH6HWZJSkpSbW2t6urq1NTUpPLycqWmpvq0\nqa+vV2Fhoe6//37FxMRcesUAgA7x2zMPCgrSrFmzlJ+fL8uylJWVpYSEBJWVlclmsyk7O1svvfSS\nTp48qaKiIlmWpeDgYBUUFFyO+gEAkmyWZVldceKampquOK2RnM5Q5eb27+oyjPHcu69qxuZbu7oM\nI5TmliotOq2ryzBGXFxcq/u4AxQADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIc\nAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHA\nAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxA\nmAOAAQhzADCAvT2NKisrVVxcLMuylJmZqdzc3GZt1q1bp8rKSvXp00fz5s3T4MGDA10rAKAVfnvm\nHo9HRUVFWrJkiQoLC1VeXq7q6mqfNhUVFTp8+LDWrFmj++67T88880ynFQwAaM5vmFdVVSk2NlbR\n0dGy2+3KyMiQ0+n0aeN0OjV+/HhJUnJyshobG3X8+PHOqRgA0IzfMHe5XIqKivJuOxwOuVyuDrcB\nAHQeLoACgAH8XgB1OByqr6/3brtcLjkcjmZtjh496t0+evRoszZfFxcX19Fa0Yr/+R/Jsrq6CpNM\n1v/exB8UPYvfnnlSUpJqa2tVV1enpqYmlZeXKzU11adNamqq3n77bUnSxx9/rLCwMPXv379zKgYA\nNGOzLP99usrKSq1fv16WZSkrK0u5ubkqKyuTzWZTdna2JKmoqEiVlZUKCQnRnDlzdPXVV3d68QCA\n89oV5gCA7o0LoABgAMIcAAzQrtv5ge7i5MmTWr58uWw2m44dO6agoCD169dPNptNK1asUHBwsLft\nihUr9POf/1whISGtHm/OnDkqLCxUaGio97XFixfL4/GooaFBZ8+elcPhkM1m04MPPuh3lhbQVRgz\nR4/10ksvKSQkRJMnT262z7Is2Ww2v8doKcy/8s9//lMHDx7UzJkzA1Iv0JnomaPHurAfUltbq9/9\n7ncaPHiw9u/fryVLlnifJxQaGqrf/va3On78uM6ePavvf//7ysrK6vD5Nm3apPr6et19992SpNdf\nf13Hjh3T+PHjtWrVKsXFxenAgQMaMmSI5s6dK7vdrr1792rDhg06e/asIiIiNHfuXPXr1y9gfwPg\nK4yZwxg1NTW69dZbVVhY2Gw45P7771dBQYEKCgr06quvqrGxscPHv+mmm7Rt2zbvl8iWLVuUmZkp\nSTp48KByc3O1atUqSdLmzZt17tw5Pfvss/rlL3+pgoICjRs3TiUlJZf4LoGW0TOHMQYOHKghQ4a0\nuG/jxo3avn27pPN3MdfW1nb4XojQ0FAlJyfrgw8+UGRkpPr06aOYmBjV1NQoJibGe+6bb75Z77zz\njpKTk3Xo0CEtX75clmXJsixFR0df2psEWkGYwxitXejcuXOn9uzZo4KCAtntduXl5encuXMXdY6s\nrCxt2rRJ/fv39/bKW2Kz2WRZloYMGaK8vLyLOhfQEQyzwBitXctvbGxUWFiY7Ha7Dh48qE8++eSi\nz5GSkqLq6mo5nU6NGzfO+/rnn3+uffv2SZLKy8s1fPhwJSYmqq6uTp9++qkkqampSYcOHbrocwNt\noWcOY7Q2e2XMmDF688039Ytf/EKxsbFKTk6+pPOMHTtWdXV16tOnj/e1QYMGqbS0VPv379fgwYOV\nlZUlu92uhQsXav369Tp9+rQsy1JOTo4SEhIu6fxAS5iaCHRQfn6+7rzzTg0dOlTS+Quva9euVUFB\nQRdXhm8yhlmAdjpx4oTmz5+v/v37e4Mc6C7omQOAAeiZA4ABCHMAMABhDgAGIMwBwACEOQAYgDAH\nAAP8H/CrGY+sA/5dAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "Trial_Trimmer = df['TrialType'] == 0\n", "ColandIden_Acc1 = df[Trial_Trimmer].pivot_table(values='Loc_Acc', index='Sub#', columns='Prev_Col+Iden_match',\n", " aggfunc=np.mean)\n", "CIs = SimpleComparisonCI(ColandIden_Acc1)\n", "Analysis_and_Plot_3(ColandIden_Acc1, CIs)" ] }, { "cell_type": "markdown", "metadata": { "collapsed": false }, "source": [ "Looks like a 2.5% increase in accuracy. This time significant. Now, this wasn't really a planned comparison, so please take this result with a grain of salt." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "As a sanity check, lets look at how repetitions in the target's location (the reported feature) effect performance. \n", "\n", "We have to quickly create a new variable coding target location repetitions" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": true }, "outputs": [], "source": [ "df['Prev_Loc'] = df['Tar_Loc'].shift(periods=1)\n", "df['Prev_Loc_match'] = np.where(df['Prev_Loc']==df['Tar_Loc'], 1, 2)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "---------------------------------\n", "Mean match: 0.9101\n", "Mean mismatch: 0.8883\n", "Match - Mismatch: 0.0218\n", "---------------------------------\n", "Match vs Mismatch\n", "\t normality test adj. Test value: 0.93 P-value: 0.1812\n", "\t T-test adj. Test value: 2.62 P-value: 0.0168\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAECCAYAAAAMxDf2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF8lJREFUeJzt3XtwVOX9x/HPJguExBCyMZArgk24BK2ICQZTi8kEKlPM\nL4pYkTK0pUoBLdBKlYJxoBnSViMFWq8NpEovEx0bq9NRIhWVaOkOJBZEkajcEgMJa0ogckn2/P7g\n5/7YZpNNyOb2+H79lbPn2XO+GyaffXjOc85jsyzLEgCgXwvq7QIAAF1HmAOAAQhzADAAYQ4ABiDM\nAcAAhDkAGMDur8ETTzyh3bt3KyIiQo8++qjPNps2bVJlZaUGDRqkxYsXa+TIkYGuEwDQDr8988zM\nTK1cubLN/RUVFTp27Jg2bNige+65R88880xACwQA+Oc3zMeOHauwsLA29zudTk2ZMkWSlJycrKam\nJjU0NASuQgCAX10eM3e5XIqKivJsOxwOuVyurh4WANAJXAAFAAN0OcwdDodOnDjh2T5x4oQcDkdX\nDwsA6AS/s1kkybIstfU8rtTUVL322mu64YYb9NFHHyksLExDhw71e8yamprOVQr0kPDwcDU2NvZ2\nGUArcXFxbe6z+Xtq4vr167Vv3z41NjYqIiJCd9xxh5qbm2Wz2ZSdnS1JKioqUmVlpUJCQrRw4UJd\neeWVfosizNFXEeboq7oU5t2FMEdfRZijr2ovzLkACgAGIMwBwAAdugAKwGxhYWEKCqJv1xe43W6d\nPn260+8jzIGvuLCwMJ09e1bNzc29XQok2e12hYWFdTrQ+SoGvuKCgoII8j6kubn5kv6XRJgDgAEI\ncwAwAGEOAP8nPT1dO3bs6O0yLgkXQAF4qa4epJqa4G47flxci+Ljz3ao7fXXX6+6ujrt2rVLkZGR\nntenTZumffv2aefOnYqPj2/z/UePHlV6eroOHz5s/GwdwhyAl5qaYOXm+n++0qUqLW1QO/nrxWaz\nKTExUS+99JK+973vSZI+/PBDnTlzRjabze/7LcuSzWZr89lSJjH7qwpAvzdz5kw9//zznu3nn39e\ns2bN8mxv27ZN3/rWtzR27FhNmjRJjz32mNd7JWncuHEaM2aMdu/eLUn64x//qJtuukljxoxRVlaW\n9u7d63nP3r17lZ2drZSUFC1atEjnzp3r7o8YEIQ5gD5t4sSJOnXqlKqqquR2u/W3v/1Nt912m6e3\nHRYWpg0bNujDDz/Us88+q+eee05bt26VJL344ouSpP3792v//v2aOHGiXn75Za1bt04bN27U/v37\ntXnzZq8hnFdeeUV//vOf9e6772rfvn0qKSnp+Q99CRhmAdDnfdk7nzx5spKTkxUTE+PZl56e7vl5\n7NixysnJ0bvvvqtp06Z5Xv9yuEWS/vKXv2jRokW6+uqrJUlXXHGF17l++MMfKjo6WpI0depUvf/+\n+932uQKJMEe3OnDggKqqqiRJVVVVSkpKkiQlJSUpOTm5N0tDPzJz5kzddtttOnLkiG6//Xavfbt3\n71ZBQYH279+v8+fP69y5c5oxY0abx6qpqWkV4Be7/PLLPT8PHjxYx44d6/oH6AEMs6BbJScna/r0\n6Zo+fbo+/vhjz88EOTojPj5eiYmJeuONNzR9+nRJ8vS077vvPt18883atWuXPvjgA333u9/1DMH4\nukgaFxenQ4cO9VzxPYQwB9AvPPbYYyopKdHgwYMlyRPYp0+fVkREhAYMGKCKigqVlpZ63uNwOBQU\nFKSDBw96Xps9e7aefPJJ7dmzR5J08OBBVVdX99wH6SYMsxigu+cFB5LTGdrbJfh1xRWWhg3r7Sp6\nT1xci0pLG7r1+B11cc96xIgRPvetXbtWq1ev1qpVq5Senq6cnBz95z//kXRhmOTHP/6xcnNz1dLS\noi1btmjGjBlqaGjQ4sWLdezYMSUmJmr9+vWKj4/v0HTHvoqVhgzgdIZ267zgQJk16w49/3zfnxnw\n8suNmjjxq7PSECsr9T1t/Zu0t9IQPXOgn+BiMtpDmKNbRUe/r+joC1O7jh8fp5SUCz3zurrxqqsb\n35ul9TvJycme0F66dKnuu+++Xq4IfQlhjm5FaAM9gzAH/kvw5R/LWdf3p64565wBOU7qoFTZ1H8v\n/OECwhz4Lyeaj2puaW5vl9GuWZql3ADV+Oqdr+rq8KsDciz0HsIc6CeiT0cr+vSF28yPhx5XyvEU\nSVJdWJ3qwup6szT0AYQ50E8Q2mgPd4ACgAEIcwD92oMPPqj169f3dhltKikp0a233trt52GYBYCX\n6qZq1Zzuvju048LiFB/asaWGOrJs3C9/+cvuKrVNy5YtU1xcnJYvX96h9j3xmADCHICXmtM1AZsp\n40tpbmmHw7yry8Z9lTDMAqBP87ds3LJly/TII49Iklwul+bNm6eUlBSNHz/es2ycdGERiyeffFLZ\n2dkaPXq0li9frvr6es2dO1djxozR7NmzdfLkSU/7BQsW6Nprr1VKSopuv/12HThwQNKFJef++te/\n6oknntCYMWP0/e9/X9KF503dfffd+vrXv66rr75aDz30kOdYlmXpF7/4hcaPH68bbrhBb7zxRsB/\nT4Q5gD6trWXjfHnqqacUFxenvXv36t///rcefPBBr/1///vfVVJSorfffltbt27V3LlztWLFCu3Z\ns0ctLS0qKirytM3KytI777yj9957T1dddZUWL14sSZozZ45uvfVWLVy40LPsnNvt1rx585SYmKh/\n/etf2rVrl3JycjzHqqioUHJysvbu3asf/ehHuv/++wP+eyLMAfR5X/bO33rrLc+ycb4e+DpgwAAd\nP35chw8fVnBwsNLS0rz2/+AHP5DD4dDw4cM1adIkT8974MCBmj59utcScd/5znc0ePBgDRgwQMuW\nLdO+fft06tQpn/VVVFTo+PHjWrVqlUJCQjRw4ECvcyckJOjOO++UzWbTrFmzdPz4cdXX1wfot3MB\nYQ6gz5s5c6ZKS0tVUlLiWTbO15j5woULdcUVV+iuu+5SRkaGfve733ntv3hJuJCQEM9an19unz59\nWpLkdru1du1aZWRkaNy4cZo8ebJsNptcLpfP+mpqapSQkKCgIN+ROuyiB+QPHjxYlmV5zhUohDmA\nPs/XsnG+hIWFKS8vT++88442b96sp59+WuXl5Z0+34svvqiysjKVlJTogw8+0D//+U9ZltXmcnRx\ncXGqrq6W2+3u9LkChTAH0C+0tWzcxV5//XXPEnFhYWGy2+0KDu78KlynT5/WwIEDFRERoaamJhUU\nFHgFeHR0tA4fPuzZvvbaazVs2DCtXbtWX3zxhc6ePSunMzAPQusopiYC8BIXFqfS3FL/Dbtw/I7q\nyLJxF/v000+1atUquVwuRUREaN68eUpPT/fZvr2pjbNmzdKbb76p6667TpGRkVq+fLm2bNni2X/n\nnXdqwYIFGj9+vCZPnqzf//73Ki4u1kMPPaS0tDQFBQUpNze31Zh9R859qTq0bFxlZaWKi4tlWZYy\nMzOVm+s9B7WpqUkbN25UfX293G63brnlFt10003tHpNl4wKnvywb11889/Yrmrvtlt4uo8fw1MS+\np1uWjXO73SoqKlJeXp4iIyO1YsUKpaWlKT7+/yf9v/baa0pMTNQDDzygkydPaunSpbrxxhsv6b83\nAIDO8ztmXlVVpdjYWEVHR8tutysjI6PVWJDNZtMXX3whSTpz5ozCw8MJcgDoQX7D3OVyKSoqyrPt\ncDhaTc+5+eabdfToUS1YsEDLly/33HYLAOgZAbkAWllZqVGjRunhhx9WbW2t8vPz9eijjyokJKTN\n94SHhwfi1JDEf4IAswQHB3c6I/2GucPh8LpTyeVyyeFweLXZvn2756JoTEyMhg0bpurqan3ta19r\n87i+BvdxaVpaQnu7BPRjHZgDgR7W0tLiMyPbC3i/wyxJSUmqra1VXV2dmpubVV5ertTUVK82l19+\nufbs2SNJamho0Geffabhw4d3tn4AvWDfiX1y23rvZhd4s9vtl3Tzkd+eeVBQkObPn6/8/HxZlqWs\nrCwlJCSorKxMNptN2dnZmjlzph5//HHPw2PmzJmjyy67rPOfAkCPu//1+6VsKSUqJeDzn6+6/CoN\nGTgkoMc0ndvtvqRb/Ts0z7w7MM88cJhnHlhftXnm3ak0t1Rp0b5vnEHntTfPnNv5AcAAhDkAGIAw\nBwADEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMA\nMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwAD\nEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAA9g70qiyslLFxcWyLEuZmZnKzc1t\n1eb999/XH/7wB7W0tGjIkCF6+OGHA14sAMA3v2HudrtVVFSkvLw8RUZGasWKFUpLS1N8fLynTVNT\nk4qKirRq1So5HA6dPHmyW4sGAHjzO8xSVVWl2NhYRUdHy263KyMjQ06n06vNjh07dP3118vhcEiS\nhgwZ0j3VAgB88tszd7lcioqK8mw7HA5VVVV5tampqVFLS4tWr16tM2fOaPr06frmN78Z+GoBAD51\naMzcH7fbrU8//VR5eXk6e/asVq1apdGjRysmJiYQhwcA+OE3zB0Oh+rr6z3bLpfLM5xycZvw8HAN\nHDhQAwcO1Lhx43Tw4MF2wzw8PLwLZeNiwcG9XQHgW3BwMH/rPcRvmCclJam2tlZ1dXWKjIxUeXm5\nlixZ4tUmLS1NmzZtktvt1vnz53XgwAHNmDGj3eM2NjZ2rXJ4tLSE9nYJgE8tLS38rQdQe1+MfsM8\nKChI8+fPV35+vizLUlZWlhISElRWViabzabs7GzFx8frmmuu0f3336+goCBlZ2crISEhoB8CANC2\nDo2ZT5gwQevXr/d6berUqV7bOTk5ysnJCVxlAIAO4w5QADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYA\nYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAG\nIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADC\nHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABigQ2FeWVmppUuXasmSJSotLW2zXVVVlWbPnq2dO3cG\nrEAAgH9+w9ztdquoqEgrV65UYWGhysvLVV1d7bPdn/70J11zzTXdUigAoG1+w7yqqkqxsbGKjo6W\n3W5XRkaGnE5nq3avvvqq0tPTNWTIkG4pFADQNr9h7nK5FBUV5dl2OBxyuVyt2jidTk2bNi3wFQIA\n/ArIBdDi4mLNmTPHs21ZViAOCwDoILu/Bg6HQ/X19Z5tl8slh8Ph1eaTTz7Rb37zG1mWpcbGRlVU\nVMhutys1NbXN44aHh3ehbFwsOLi3KwB8Cw4O5m+9h/gN86SkJNXW1qqurk6RkZEqLy/XkiVLvNr8\n9re/9fz8+OOP67rrrms3yCWpsbHxEkvGf2tpCe3tEgCfWlpa+FsPoPa+GP2GeVBQkObPn6/8/HxZ\nlqWsrCwlJCSorKxMNptN2dnZAS0WANB5fsNckiZMmKD169d7vTZ16lSfbRctWtT1qgAAncIdoABg\nAMIcAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYg\nzAHAAIQ5ABiAMAcAAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIc\nAAxAmAOAAQhzADAAYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwgL0jjSorK1Vc\nXCzLspSZmanc3Fyv/Tt27NBLL70kSQoJCdHdd9+tESNGBL5aAIBPfnvmbrdbRUVFWrlypQoLC1Ve\nXq7q6mqvNsOGDdPq1av1yCOPaObMmXrqqae6rWAAQGt+w7yqqkqxsbGKjo6W3W5XRkaGnE6nV5vR\no0crNDRUkpScnCyXy9U91QIAfPIb5i6XS1FRUZ5th8PRblhv27ZNEyZMCEx1AIAO6dCYeUft3btX\n27dv15o1a/y2DQ8PD+Spv9KCg3u7AsC34OBg/tZ7iN8wdzgcqq+v92y7XC45HI5W7Q4dOqSnn35a\nP//5z3XZZZf5PXFjY2MnS0VbWlpCe7sEwKeWlhb+1gOovS9Gv8MsSUlJqq2tVV1dnZqbm1VeXq7U\n1FSvNvX19SosLNS9996rmJiYrlcMAOgUvz3zoKAgzZ8/X/n5+bIsS1lZWUpISFBZWZlsNpuys7P1\nwgsv6NSpUyoqKpJlWQoODlZBQUFP1A8AkGSzLMvqjRPX1NT0xmmN5HSGKjd3aG+XYYzn3n5Fc7fd\n0ttlGKE0t1Rp0Wm9XYYx4uLi2tzHHaAAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcA\nAxDmAGAAwhwADECYA4ABCHMAMABhDgAGIMwBwACEOQAYgDAHAAMQ5gBgAMIcAAxAmAOAAQhzADAA\nYQ4ABiDMAcAAhDkAGIAwBwADEOYAYADCHAAMQJgDgAEIcwAwAGEOAAYgzAHAAIQ5ABiAMAcAAxDm\nAGAAwhwADECYA4AB7B1pVFlZqeLiYlmWpczMTOXm5rZqs2nTJlVWVmrQoEFavHixRo4cGehaAQBt\n8Nszd7vdKioq0sqVK1VYWKjy8nJVV1d7tamoqNCxY8e0YcMG3XPPPXrmmWe6rWAAQGt+w7yqqkqx\nsbGKjo6W3W5XRkaGnE6nVxun06kpU6ZIkpKTk9XU1KSGhobuqRgA0IrfMHe5XIqKivJsOxwOuVyu\nTrcBAHQfLoACgAH8XgB1OByqr6/3bLtcLjkcjlZtTpw44dk+ceJEqzb/LS4urrO1og3/8z+SZfV2\nFSaZoe9+g18o+he/PfOkpCTV1taqrq5Ozc3NKi8vV2pqqleb1NRUvfnmm5Kkjz76SGFhYRo6dGj3\nVAwAaMVmWf77dJWVldq8ebMsy1JWVpZyc3NVVlYmm82m7OxsSVJRUZEqKysVEhKihQsX6sorr+z2\n4gEAF3QozAEAfRsXQAHAAIQ5ABigQ7fzA33FqVOntGbNGtlsNn3++ecKCgrSkCFDZLPZtHbtWgUH\nB3varl27Vj/5yU8UEhLS5vEWLlyowsJChYaGel5bsWKF3G63Ghsbde7cOTkcDtlsNj3wwAN+Z2kB\nvYUxc/RbL7zwgkJCQjRjxoxW+yzLks1m83sMX2H+pX/84x86cuSI5s2bF5B6ge5Ezxz91sX9kNra\nWv3617/WyJEjdejQIa1cudLzPKHQ0FD96le/UkNDg86dO6dvf/vbysrK6vT5tm7dqvr6et11112S\npFdffVWff/65pkyZonXr1ikuLk6HDx/WqFGjtGjRItntdh04cEBbtmzRuXPnFBERoUWLFmnIkCEB\n+x0AX2LMHMaoqanRLbfcosLCwlbDIffee68KCgpUUFCgV155RU1NTZ0+/je+8Q3t3LnT8yWyfft2\nZWZmSpKOHDmi3NxcrVu3TpK0bds2nT9/Xs8++6x+9rOfqaCgQJMnT1ZJSUkXPyXgGz1zGGP48OEa\nNWqUz30vv/yydu3aJenCXcy1tbWdvhciNDRUycnJeu+99xQZGalBgwYpJiZGNTU1iomJ8Zz7xhtv\n1FtvvaXk5GQdPXpUa9askWVZsixL0dHRXfuQQBsIcxijrQude/bs0f79+1VQUCC73a68vDydP3/+\nks6RlZWlrVu3aujQoZ5euS82m02WZWnUqFHKy8u7pHMBncEwC4zR1rX8pqYmhYWFyW6368iRI/r4\n448v+RwpKSmqrq6W0+nU5MmTPa9/9tlnOnjwoCSpvLxcY8eOVWJiourq6vTJJ59Ikpqbm3X06NFL\nPjfQHnrmMEZbs1cmTpyo119/XT/96U8VGxur5OTkLp1n0qRJqqur06BBgzyvjRgxQqWlpTp06JBG\njhyprKws2e12LVu2TJs3b9aZM2dkWZZycnKUkJDQpfMDvjA1Eeik/Px83XHHHRo9erSkCxdeN27c\nqIKCgl6uDF9lDLMAHXTy5EktWbJEQ4cO9QQ50FfQMwcAA9AzBwADEOYAYADCHAAMQJgDgAEIcwAw\nAGEOAAb4X1aNDcJ2Mq6xAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "Trial_Trimmer = df['TrialType'] == 0\n", "Loc_Acc1 = df[Trial_Trimmer].pivot_table(values='Loc_Acc', index='Sub#', columns='Prev_Loc_match', aggfunc=np.mean)\n", "CIs = SimpleComparisonCI(Loc_Acc1)\n", "Analysis_and_Plot_3(Loc_Acc1, CIs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Target location repetitions lead to a 2% increase in performance. Again, this result is robust. \n", "\n", "It's a good sign that this effect is about the same size as repetitions in the unreported features." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Replicate Experiments 1 Intertrial Analyses with Experiment 1b" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Experiment 1 had some evidence that participants unconsciously knew the color and identity of the target, since they performed a little better when the color and identity repeated. The effect was small, so I am not 100% confident that it's robust. \n", "\n", "The best way to demonstrate that this effect is real would be to show that it also exists in another similar Experiment. Chen and Wyble provide a replication of Experiment 1. In this experiment, the only difference is the target and distractors appear for longer and are not masked (making them easier to see). \n", "\n", "If participants response more accurately when the target color and identity repeat in Experiment 1b, then we can be a little more confident that participants are unconsciously aware of the target's color and identity. " ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": false }, "outputs": [], "source": [ "filename = 'Exp1b.csv'\n", "df = pd.read_csv(filename)\n", "df.columns = ['Sub#', 'Block', 'Trial#', 'TarCol', 'Tar_Iden','Tar_Loc', 'Col_Resp', 'Iden_Resp', 'Loc_Resp', \n", " 'Col_Acc', 'Iden_Acc', 'Loc_Acc'] #naming the columns of the data file.\n", "\n", "df['TrialType'] = np.where(df['Trial#']<156, 0, np.where(df['Trial#']==156, 1, 2))\n", "df['TrialType2'] = np.where(df['Trial#']<156, 0, np.where(df['Trial#']==156, 1, np.where(df['Trial#']==157, 2, -1)))\n", "\n", "df['Prev_TarCol'] = df['TarCol'].shift(periods=1)\n", "df['Prev_TarCol_match'] = np.where(df['Prev_TarCol']==df['TarCol'], 1, 2)\n", "df['Prev_Iden'] = df['Tar_Iden'].shift(periods=1)\n", "df['Prev_Iden_match'] = np.where(df['Prev_Iden']==df['Tar_Iden'], 1, 2)\n", "df['Prev_Col+Iden_match'] = np.where((df['Prev_Iden_match']==1) & (df['Prev_TarCol_match']==1), 1, 2)\n", "\n", "Trial_Trimmer = df['TrialType'] == 0 #only interested in pre-surprise trials\n", "ColandIden_Acc = df[Trial_Trimmer].pivot_table(values='Loc_Acc', index='Sub#', \n", " columns='Prev_Col+Iden_match', aggfunc=np.mean)" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "---------------------------------\n", "Mean match: 0.9716\n", "Mean mismatch: 0.9644\n", "Match - Mismatch: 0.0072\n", "---------------------------------\n", "Match vs Mismatch\n", "\t normality test adj. Test value: 0.93 P-value: 0.1875\n", "\t T-test adj. Test value: 2.81 P-value: 0.0112\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAECCAYAAAAMxDf2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF7VJREFUeJzt3XtwVOXBx/HfJguExBCyMZArgk24BK2ICQajxWQClSmm\nq4gVKUNbqhTQAq1UKRgHmiFtJVKg9dpAqvQy0bGxOh0lUvESLd2BxIIokiq3xEDCmhKIXJI97x+M\n+7LNZROyYcPj9/NXzu6Tc57NTL45OXvOHptlWZYAAJe0kGBPAADQc8QcAAxAzAHAAMQcAAxAzAHA\nAMQcAAxg9zfgiSee0M6dOxUVFaU1a9a0O2bjxo2qqqrSgAEDtHDhQg0fPjzQ8wQAdMLvnnl2draW\nL1/e4fOVlZU6cuSI1q9fr3vvvVfPPPNMQCcIAPDPb8xHjx6tiIiIDp93uVyaNGmSJCk1NVXNzc1q\nbGwM3AwBAH71+Ji52+1WTEyMd9nhcMjtdvd0tQCAbuANUAAwQI9j7nA4dOzYMe/ysWPH5HA4erpa\nAEA3+D2bRZIsy1JHn8eVnp6u1157TTfccIM+/vhjRUREaPDgwX7XWVtb272ZAhdJZGSkmpqagj0N\noI2EhIQOn7P5+9TEdevWac+ePWpqalJUVJTuvPNOtbS0yGazKTc3V5JUXFysqqoqhYWFaf78+bry\nyiv9ToqYo68i5uirehTz3kLMvxr27dun6upqSVJ1dbVSUlIkSSkpKUpNTQ3m1DpEzNFXdRbzLh1m\nAS5UamqqN9qLFy/W/fffH+QZAWYi5sAl4lL8LwcXDzE3QE3NANXWhgZ7Gl3icoUHewp+XXGFpSFD\ngj2Ltnrzv5yIiAiFhHCmcl/g8Xh08uTJbn8fMTdAbW2onE7/ZxAF24wZuiTm+fLLTX0y5r0lIiJC\np0+fVktLS7CnAkl2u10RERHdDjoxB/5H6OX/kav+QLCn4Zer3hWQ9aQPSJetxRaQdaHnWlpaNHDg\nwG5/HzFHr4qN/UCxsR9Iko4eHaO0tFJJUn39WNXXjw3m1Dp0rOWwZpc5gz2NTs3QDDkDNMdX73pV\nV0deHZB1IXiIOXpVX472pSb2ZKxiT8ZKko6GH1Xa0TRJUn1Eveoj6oM5NfQBxBy4RBDt3peZmak1\na9boxhtvDPZUuo2YA/DR22dHJSS0KjHxdJfGXn/99aqvr9eOHTsUHR3tfXzKlCnas2ePtm/frsTE\nxA6///Dhw8rMzNTBgweNP1uHmAPw0dtnR5WVNaqT/vqw2WxKTk7WSy+9pO9973uSpI8++kinTp2S\nzeb/TVvLsmSz2Tr8bCmTmP2nCsAlb/r06Xr++ee9y88//7xmzJjhXd66dau++c1vavTo0ZowYYIe\ne+wxn++VpDFjxmjUqFHauXOnJOmPf/yjbr75Zo0aNUo5OTnavXu393t2796t3NxcpaWlacGCBTpz\n5kxvv8SAIOYA+rTx48frxIkTqq6ulsfj0d/+9jfdfvvt3r3tiIgIrV+/Xh999JGeffZZPffcc9qy\nZYsk6cUXX5Qk7d27V3v37tX48eP18ssva+3atdqwYYP27t2rTZs2+RzCeeWVV/TnP/9Z7733nvbs\n2aPS0tKL/6IvAIdZAPR5X+6dT5w4UampqYqLi/M+l5mZ6f169OjRysvL03vvvacpU6Z4H//ycIsk\n/eUvf9GCBQt09dXnTse84oorfLb1wx/+ULGx584amjx5sj744INee12BRMwB9HnTp0/X7bffrkOH\nDumOO+7weW7nzp0qLCzU3r17dfbsWZ05c0bTpk3rcF21tbVtAn6+yy+/3Pv1wIEDdeTIkZ6/gIuA\nwywA+rzExEQlJyfrjTfe0NSpUyXJu6d9//3365ZbbtGOHTv04Ycf6rvf/a73EEx7b5ImJCTowIG+\nf4VvdxFzAJeExx57TKWlpd5L3b8M9smTJxUVFaV+/fqpsrJSZWVl3u9xOBwKCQnR/v37vY/NnDlT\nTz75pHbt2iVJ2r9/v2pqai7eC+klHGYB4CMhoVVlZY29uv6uOn/PetiwYe0+t3r1aq1cuVIrVqxQ\nZmam8vLy9N///lfSucMkP/7xj+V0OtXa2qrNmzdr2rRpamxs1MKFC3XkyBElJydr3bp1SkxM7NLp\njn0VdxoygMsVfkl8GuGl4rm3X9HsrbcGexoXDZ/N0vd0dLerzu40xGEWADAAMQcAAxBzADAAMQcA\nAxBzADAAMQcAAxBzADAAMQcAAxBzAJe0hx56SOvWrQv2NDpUWlqq2267rde3w+X8AHzUNNeo9mTv\nXaGdEJGgxPCu3WqoK7eN++Uvf9lbU+3QkiVLlJCQoKVLl3Zp/MX4mABiDsBH7claOcucvbb+MmdZ\nl2Pe09vGfZVwmAVAn+bvtnFLlizRo48+Kklyu92aM2eO0tLSNHbsWO9t46RzN7F48sknlZubq5Ej\nR2rp0qVqaGjQ7NmzNWrUKM2cOVPHjx/3jp83b56uvfZapaWl6Y477tC+ffsknbvl3F//+lc98cQT\nGjVqlL7//e9LOvd5U/fcc4++/vWv6+qrr9bDDz/sXZdlWfrFL36hsWPH6oYbbtAbb7wR8J8TMQfQ\np3V027j2PPXUU0pISNDu3bv173//Ww899JDP83//+99VWlqqt99+W1u2bNHs2bO1bNky7dq1S62t\nrSouLvaOzcnJ0bvvvqv3339fV111lRYuXChJmjVrlm677TbNnz/fe9s5j8ejOXPmKDk5Wf/617+0\nY8cO5eXleddVWVmp1NRU7d69Wz/60Y/0wAMPBPznRMwB9Hlf7p2/9dZb3tvGtfeBr/369dPRo0d1\n8OBBhYaGKiMjw+f5H/zgB3I4HBo6dKgmTJjg3fPu37+/pk6d6nOLuO985zsaOHCg+vXrpyVLlmjP\nnj06ceJEu/OrrKzU0aNHtWLFCoWFhal///4+205KStJdd90lm82mGTNm6OjRo2poaAjQT+ccYg6g\nz5s+fbrKyspUWlrqvW1ce8fM58+fryuuuEJ33323srKy9Lvf/c7n+fNvCRcWFua91+eXyydPnpQk\neTwerV69WllZWRozZowmTpwom80mt9vd7vxqa2uVlJSkkJD2kzpkyBDv1wMHDpRlWd5tBQoxB9Dn\ntXfbuPZEREQoPz9f7777rjZt2qSnn35aFRUV3d7eiy++qPLycpWWlurDDz/UP//5T1mW1eHt6BIS\nElRTUyOPx9PtbQUKMQdwSejotnHne/311723iIuIiJDdbldoaGi3t3Xy5En1799fUVFRam5uVmFh\noU/AY2NjdfDgQe/ytddeqyFDhmj16tX64osvdPr0ablcrm5vtyc4NRGAj4SIBJU5y/wP7MH6u6or\nt40736effqoVK1bI7XYrKipKc+bMUWZmZrvjOzu1ccaMGXrzzTd13XXXKTo6WkuXLtXmzZu9z991\n112aN2+exo4dq4kTJ+r3v/+9SkpK9PDDDysjI0MhISFyOp1tjtl3ZdsXqku3jauqqlJJSYksy1J2\ndracTt9zUJubm7VhwwY1NDTI4/Ho1ltv1c0339zpOrltXOBw27jA4rZxCLYLuW2c3z1zj8ej4uJi\n5efnKzo6WsuWLVNGRoYSE///pP/XXntNycnJevDBB3X8+HEtXrxYN9100wX9ewMA6D6/x8yrq6sV\nHx+v2NhY2e12ZWVltTkWZLPZ9MUXX0iSTp06pcjISEIOABeR35i73W7FxMR4lx0OR5vTc2655RYd\nPnxY8+bN09KlS72X3QIALo6AvAFaVVWlESNG6JFHHlFdXZ0KCgq0Zs0ahYWFdfg9kZGRgdg0JPFP\nEGCW0NDQbjfSb8wdDofPlUput1sOh8NnzLZt27xvisbFxWnIkCGqqanR1772tQ7X297BfVyY1tbw\nYE8Bl7AunAOBi6y1tbXdRnYWeL+HWVJSUlRXV6f6+nq1tLSooqJC6enpPmMuv/xy7dq1S5LU2Nio\nzz77TEOHDu3u/AEEwZ5je+SxBe9iF/iy2+0XdPGR3z3zkJAQzZ07VwUFBbIsSzk5OUpKSlJ5ebls\nNptyc3M1ffp0Pf74494Pj5k1a5Yuu+yy7r8KABfdA68/IOVKaTFpAT//+arLr9Kg/oMCuk7TeTye\nC7rUv0vnmfcGzjMPHM4zD6yv2nnmvanMWaaM2PYvnEH3dXaeOZfzA4ABiDkAGICYA4ABiDkAGICY\nA4ABiDkAGICYA4ABiDkAGICYA4ABiDkAGICYA4ABiDkAGICYA4ABiDkAGICYA4ABiDkAGICYA4AB\niDkAGICYA4ABiDkAGICYA4ABiDkAGICYA4ABiDkAGICYA4ABiDkAGICYA4ABiDkAGICYA4ABiDkA\nGICYA4ABiDkAGICYA4ABiDkAGICYA4AB7F0ZVFVVpZKSElmWpezsbDmdzjZjPvjgA/3hD39Qa2ur\nBg0apEceeSTgkwUAtM9vzD0ej4qLi5Wfn6/o6GgtW7ZMGRkZSkxM9I5pbm5WcXGxVqxYIYfDoePH\nj/fqpAEAvvweZqmurlZ8fLxiY2Nlt9uVlZUll8vlM+add97R9ddfL4fDIUkaNGhQ78wWANAuv3vm\nbrdbMTEx3mWHw6Hq6mqfMbW1tWptbdXKlSt16tQpTZ06Vd/4xjcCP1sAQLu6dMzcH4/Ho08//VT5\n+fk6ffq0VqxYoZEjRyouLi4QqwcA+OE35g6HQw0NDd5lt9vtPZxy/pjIyEj1799f/fv315gxY7R/\n//5OYx4ZGdmDaeN8oaHBngHQvtDQUH7XLxK/MU9JSVFdXZ3q6+sVHR2tiooKLVq0yGdMRkaGNm7c\nKI/Ho7Nnz2rfvn2aNm1ap+ttamrq2czh1doaHuwpAO1qbW3ldz2AOvvD6DfmISEhmjt3rgoKCmRZ\nlnJycpSUlKTy8nLZbDbl5uYqMTFR11xzjR544AGFhIQoNzdXSUlJAX0RAICOdemY+bhx47Ru3Tqf\nxyZPnuyznJeXp7y8vMDNDADQZVwBCgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABi\nDgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAG\nIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAGIOYAYABiDgAGIOYA\nYABiDgAGIOYAYIAuxbyqqkqLFy/WokWLVFZW1uG46upqzZw5U9u3bw/YBAEA/vmNucfjUXFxsZYv\nX66ioiJVVFSopqam3XF/+tOfdM011/TKRAEAHfMb8+rqasXHxys2NlZ2u11ZWVlyuVxtxr366qvK\nzMzUoEGDemWiAICO+Y252+1WTEyMd9nhcMjtdrcZ43K5NGXKlMDPEADgV0DeAC0pKdGsWbO8y5Zl\nBWK1AIAusvsb4HA41NDQ4F12u91yOBw+Yz755BP95je/kWVZampqUmVlpex2u9LT0ztcb2RkZA+m\njfOFhgZ7BkD7QkND+V2/SPzGPCUlRXV1daqvr1d0dLQqKiq0aNEinzG//e1vvV8//vjjuu666zoN\nuSQ1NTVd4JTxv1pbw4M9BaBdra2t/K4HUGd/GP3GPCQkRHPnzlVBQYEsy1JOTo6SkpJUXl4um82m\n3NzcgE4WANB9fmMuSePGjdO6det8Hps8eXK7YxcsWNDzWQEAuoUrQAHAAMQcAAxAzAHAAMQcAAxA\nzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHA\nAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQc\nAAxAzAHAAMQcAAxAzAHAAMQcAAxAzAHAAMQcAAxg78qgqqoqlZSUyLIsZWdny+l0+jz/zjvv6KWX\nXpIkhYWF6Z577tGwYcMCP1sAQLv87pl7PB4VFxdr+fLlKioqUkVFhWpqanzGDBkyRCtXrtSjjz6q\n6dOn66mnnuq1CQMA2vIb8+rqasXHxys2NlZ2u11ZWVlyuVw+Y0aOHKnw8HBJUmpqqtxud+/MFgDQ\nLr8xd7vdiomJ8S47HI5OY71161aNGzcuMLMDAHRJl46Zd9Xu3bu1bds2rVq1yu/YyMjIQG76Ky00\nNNgzANoXGhrK7/pF4jfmDodDDQ0N3mW32y2Hw9Fm3IEDB/T000/r5z//uS677DK/G25qaurmVNGR\n1tbwYE8BaFdrayu/6wHU2R9Gv4dZUlJSVFdXp/r6erW0tKiiokLp6ek+YxoaGlRUVKT77rtPcXFx\nPZ8xAKBb/O6Zh4SEaO7cuSooKJBlWcrJyVFSUpLKy8tls9mUm5urF154QSdOnFBxcbEsy1JoaKgK\nCwsvxvwBAJJslmVZwdhwbW1tMDZrJJcrXE7n4GBPwxjPvf2KZm+9NdjTMEKZs0wZsRnBnoYxEhIS\nOnyOK0ABwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAM\nQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwB\nwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwADEHAAMQMwBwAD2rgyqqqpSSUmJ\nLMtSdna2nE5nmzEbN25UVVWVBgwYoIULF2r48OGBnisAoAN+98w9Ho+Ki4u1fPlyFRUVqaKiQjU1\nNT5jKisrdeTIEa1fv1733nuvnnnmmV6bMACgLb8xr66uVnx8vGJjY2W325WVlSWXy+UzxuVyadKk\nSZKk1NRUNTc3q7GxsXdmDABow2/M3W63YmJivMsOh0Nut7vbYwAAvYc3QAHAAH7fAHU4HGpoaPAu\nu91uORyONmOOHTvmXT527FibMf8rISGhu3NFB779bcmygj0Lk0zTd2/kB4pLi98985SUFNXV1am+\nvl4tLS2qqKhQenq6z5j09HS9+eabkqSPP/5YERERGjx4cO/MGADQhs2y/O/TVVVVadOmTbIsSzk5\nOXI6nSovL5fNZlNubq4kqbi4WFVVVQoLC9P8+fN15ZVX9vrkAQDndCnmAIC+jTdAAcAAxBwADNCl\ny/mBvuLEiRNatWqVbDabPv/8c4WEhGjQoEGy2WxavXq1QkNDvWNXr16tn/zkJwoLC+twffPnz1dR\nUZHCw8O9jy1btkwej0dNTU06c+aMHA6HbDabHnzwQb9naQHBwjFzXLJeeOEFhYWFadq0aW2esyxL\nNpvN7zrai/mX/vGPf+jQoUOaM2dOQOYL9Cb2zHHJOn8/pK6uTr/+9a81fPhwHThwQMuXL/d+nlB4\neLh+9atfqbGxUWfOnNG3vvUt5eTkdHt7W7ZsUUNDg+6++25J0quvvqrPP/9ckyZN0tq1a5WQkKCD\nBw9qxIgRWrBggex2u/bt26fNmzfrzJkzioqK0oIFCzRo0KCA/QyAL3HMHMaora3VrbfeqqKiojaH\nQ+677z4VFhaqsLBQr7zyipqbm7u9/htvvFHbt2/3/hHZtm2bsrOzJUmHDh2S0+nU2rVrJUlbt27V\n2bNn9eyzz+pnP/uZCgsLNXHiRJWWlvbwVQLtY88cxhg6dKhGjBjR7nMvv/yyduzYIencVcx1dXXd\nvhYiPDxcqampev/99xUdHa0BAwYoLi5OtbW1iouL8277pptu0ltvvaXU1FQdPnxYq1atkmVZsixL\nsbGxPXuRQAeIOYzR0Rudu3bt0t69e1VYWCi73a78/HydPXv2graRk5OjLVu2aPDgwd698vbYbDZZ\nlqURI0YoPz//grYFdAeHWWCMjt7Lb25uVkREhOx2uw4dOqT//Oc/F7yNtLQ01dTUyOVyaeLEid7H\nP/vsM+3fv1+SVFFRodGjRys5OVn19fX65JNPJEktLS06fPjwBW8b6Ax75jBGR2evjB8/Xq+//rp+\n+tOfKj4+XqmpqT3azoQJE1RfX68BAwZ4Hxs2bJjKysp04MABDR8+XDk5ObLb7VqyZIk2bdqkU6dO\nybIs5eXlKSkpqUfbB9rDqYlANxUUFOjOO+/UyJEjJZ1743XDhg0qLCwM8szwVcZhFqCLjh8/rkWL\nFmnw4MHekAN9BXvmAGAA9swBwADEHAAMQMwBwADEHAAMQMwBwADEHAAM8H8lngrVkxCiMwAAAABJ\nRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "CIs = SimpleComparisonCI(ColandIden_Acc)\n", "Analysis_and_Plot_3(ColandIden_Acc, CIs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Wow. Only a 1% change in accuracy, so again not big. Nonetheless, this result is signficant. So, Some evidence that participants perform a little better when the targets' color and identity repeat. \n", "\n", "This suggests that participants retain some information about the targets' color and identity even though they cannot explicitly report these attributes. \n", "\n", "Now, I would probably want to replicate this result again before trusting it, but I'm relatively confident that participants unconsciously retain some information about the target's color and identity.\n", "\n" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.857142857143\n", "1.0\n", "0.864864864865\n", "0.993197278912\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAECCAYAAAAYfWtSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFhtJREFUeJzt3W1sU+fBxvHr4DSkCZZd562BiKW8TAyasIq0o13VhFFp\nFWNq+qFIpaqWamhtAlvnvWUUoWnTtBUNCCmMFKmsrdZ9mJBINKoJadPsqmWVliyJFC1LKRrjKTCW\nFxPPJDipnfN84Kmfsoa82MeOw/3/fYlyfJ9zLru9feX4+Bws27ZtAQCMtGi+AwAA5g8lAAAGowQA\nwGCUAAAYjBIAAINRAgBgsJyZBrS2tqqrq0sej0f79++XJL355pv661//qpycHJWWlqqxsVH5+fmS\npLa2NgUCAblcLtXX12v9+vXpfQYAgKTNeCSwadMm7dmz56ZlVVVVOnDggH7xi1+orKxM7e3tkqSL\nFy/qvffeU3Nzs3bv3q1XX31VXIYAANlrxhJYs2aNCgoKblpWVVWlRYturLp69WoNDw9Lkjo7O/XQ\nQw/J5XKppKREZWVlOnfuXBpiAwCckPI5gUAgoPvuu0+SFAqFVFRUlHjM5/MpFAqlugsAQJqkVAIn\nT56Uy+XSww8/7FQeAEAGJV0CwWBQ3d3deuGFFxLLfD6fhoaGEr8PDw/L5/OllhAAkDYzfjtIkmzb\nvukEb09Pj373u9/pxz/+se64447E8urqar388svaunWrQqGQrly5olWrVs24/cuXLycRPbPcbrci\nkch8x5gROZ1FTucshIzSwsm5dOlSR7YzYwm0tLSor69PkUhEDQ0N2rZtm9ra2hSLxfTTn/5U0o2T\nwzt27FB5ebkefPBB+f1+5eTkaMeOHbIsy5GgAADnWdlwK2mOBJxDTmeR0zkLIaO0cHI6dSTAFcMA\nYDBKAAAMRgkAgMEoAQAwGCUAAAajBADAYJQAABiMEgAAg1ECAGAwSgAADEYJAIDBKAEAMBglAAAG\nowQAwGCUAAAYjBIAAINRAgBgMEoAAAxGCQCAwSgBADAYJQAABqMEAMBglAAAGCxnvgMAgOlcV4ek\n0ODsV/AVS0uXOrJvSgAA5ltoUBMvNc16eO4P9zm2az4OAgCDUQIAYDBKAAAMNuM5gdbWVnV1dcnj\n8Wj//v2SpGvXrunQoUMaHBxUSUmJ/H6/8vPzJUltbW0KBAJyuVyqr6/X+vXr0/sMAABJm/FIYNOm\nTdqzZ89Ny9rb21VZWamWlhatW7dObW1tkqSLFy/qvffeU3Nzs3bv3q1XX31Vtm2nJzkAIGUzlsCa\nNWtUUFBw07LOzk7V1NRIkmpra9XR0ZFY/tBDD8nlcqmkpERlZWU6d+5cGmIDAJyQ1DmBcDgsr9cr\nSfJ6vQqHw5KkUCikoqKixDifz6dQKORATABAOjhyYtiyLCc2AwDIsKQuFvN6vRoZGUn89Hg8km78\n5T80NJQYNzw8LJ/PN+P23G53MjEyKjc3l5wOIqezFkLOhZBRmp+c4665vRW75jh+OrPakm3bN53g\n3bBhg4LBoOrq6hQMBlVdXS1Jqq6u1ssvv6ytW7cqFArpypUrWrVq1Yzbj0QiScbPHLfbTU4HkdNZ\nCyHnQsgozU9OVzw2p/HxOY6fzowl0NLSor6+PkUiETU0NGjbtm2qq6tTc3OzAoGAiouL5ff7JUnl\n5eV68MEH5ff7lZOTox07dvBREQBksRlL4IUXXphy+d69e6dc/sQTT+iJJ55ILRUAICO4YhgADEYJ\nAIDBKAEAMBglAAAGowQAwGCUAAAYjBIAAINRAgBgMEoAAAxGCQCAwSgBADAYJQAABqMEAMBglAAA\nGIwSAACDUQIAYDBKAAAMRgkAgMEoAQAwGCUAAAajBADAYJQAABiMEgAAg1ECAGAwSgAADEYJAIDB\nKAEAMBglAAAGowQAwGA5qaz81ltvKRAIyLIsLV++XI2NjYpGozp06JAGBwdVUlIiv9+v/Px8p/IC\nAByU9JFAKBTS6dOntW/fPu3fv1/xeFzvvvuu2tvbVVlZqZaWFq1bt05tbW1O5gUAOCilj4MmJycV\njUYVj8c1MTEhn8+nzs5O1dTUSJJqa2vV0dHhSFAAgPOS/jjI5/Np69atamxs1OLFi1VVVaWqqiqF\nw2F5vV5JktfrVTgcdiwsAMBZSZfA6OioOjs7dfToUeXn5+vgwYN65513PjXOsqwZt+V2u5ONkTG5\nubnkdBA5nbUQci6EjNL85Bx3ze2t2DXH8dNJeku9vb0qKSnRkiVLJEkPPPCA3n//fXm9Xo2MjCR+\nejyeGbcViUSSjZExbrebnA4ip7MWQs6FkFGan5yueGxO4+NzHD+dpM8JFBUV6YMPPtDExIRs21Zv\nb6/Ky8u1YcMGBYNBSVIwGFR1dbVTWQEADkv6SGDVqlXauHGjmpqa5HK5VFFRoUcffVTRaFTNzc0K\nBAIqLi6W3+93Mi8AwEEpfbD05JNP6sknn7xp2ZIlS7R3796UQgEAMoMrhgHAYJQAABiMEgAAg1EC\nAGAwSgAADEYJAIDBKAEAMBglAAAGowQAwGCUAAAYjBIAAINRAgBgMEoAAAxGCQCAwSgBADAYJQAA\nBqMEAMBglAAAGIwSAACDUQIAYDBKAAAMRgkAgMEoAQAwGCUAAAajBADAYJQAABiMEgAAg1ECAGCw\nnFRWHhsb0yuvvKIPP/xQlmWpoaFBZWVlOnTokAYHB1VSUiK/36/8/Hyn8gIAHJRSCbz22mu67777\n9J3vfEfxeFzj4+M6efKkKisr9fjjj6u9vV1tbW16+umnncoLAHBQ0h8HjY2Nqb+/X5s2bZIkuVwu\n5efnq7OzUzU1NZKk2tpadXR0OJMUAOC4pI8EBgYG5Ha7dfToUV24cEErVqxQfX29wuGwvF6vJMnr\n9SocDjsWFgDgrKRLYHJyUufPn9fXv/51rVy5Uq+//rra29s/Nc6yrBm35Xa7k42RMbm5ueR0EDmd\ntRByLoSM0vzkHHfN7a3YNcfx00l6Sz6fT4WFhVq5cqUkaePGjWpvb5fX69XIyEjip8fjmXFbkUgk\n2RgZ43a7yekgcjprIeRcCBml+cnpisfmND4+x/HTSfqcgNfrVWFhoS5fvixJ6u3tVXl5uTZs2KBg\nMChJCgaDqq6udiQoAMB5KR1TPPvsszp8+LBisZhKS0vV2NioyclJNTc3KxAIqLi4WH6/36msAACH\npVQCFRUV+vnPf/6p5Xv37k1lswCADOGKYQAwGCUAAAajBADAYJQAABiMEgAAg1ECAGAwSgAADEYJ\nAIDBKAEAMBglAAAGowQAwGCUAAAYjBIAAINRAgBgMEoAAAxGCQCAwSgBADAYJQAABqMEAMBglAAA\nGIwSAACDUQIAYDBKAAAMRgkAgMEoAQAwGCUAAAajBADAYJQAABgsJ9UNTE5Oavfu3fL5fGpqatK1\na9d06NAhDQ4OqqSkRH6/X/n5+U5kBQA4LOUjgd///vdatmxZ4vf29nZVVlaqpaVF69atU1tbW6q7\nAACkSUolMDw8rO7ubm3evDmxrLOzUzU1NZKk2tpadXR0pJYQAJA2KZXAG2+8oWeeeUaWZSWWhcNh\neb1eSZLX61U4HE4tIQAgbZI+J9DV1SWPx6OKigr97W9/u+W4TxbErbjd7mRjZExubi45HUROZy2E\nnAshozQ/Ocddc3srds1x/HSS3lJ/f786OzvV3d2tiYkJXb9+XYcPH5bX69XIyEjip8fjmXFbkUgk\n2RgZ43a7yekgcjprIeRcCBml+cnpisfmND4+x/HTSboEtm/fru3bt0uS+vr6dOrUKX3zm9/Um2++\nqWAwqLq6OgWDQVVXVzsWFgCy3VB0UkOjs3uTLirIUVHe/H5T37ljiv9TV1en5uZmBQIBFRcXy+/3\nO70LAMhaQ6MxNZ3+x6zG7ntshYryctOcaHqOlMDatWu1du1aSdKSJUu0d+9eJzYLAEgzrhgGAINR\nAgBgMEoAAAxGCQCAwSgBADAYJQAABqMEAMBglAAAGIwSAACDUQIAYDBKAAAMRgkAgMEoAQAwGCUA\nAAajBADAYJQAABiMEgAAg1ECAGAwSgAADEYJAIDBKAEAMBglAAAGowQAwGCUAAAYjBIAAINRAgBg\nsJz5DgAA6TYUndTQaGxWY8tiY/LkzG2dooIcFeUtzL+pKQEAt72h0ZiaTv9jVmP3f2W1PF7XnNbZ\n99gKFeXlphJx3iRdAsPDwzpy5IjC4bAsy9LmzZu1ZcsWXbt2TYcOHdLg4KBKSkrk9/uVn5/vZGYA\ngEOSLgGXy6Wvfe1rqqioUDQaVVNTk9avX69AIKDKyko9/vjjam9vV1tbm55++mknMwMAHJL0h1he\nr1cVFRWSpLy8PC1btkzDw8Pq7OxUTU2NJKm2tlYdHR2OBAUAOM+RMxkDAwO6cOGCPvvZzyocDsvr\n9Uq6URThcNiJXQAA0iDlEohGozp48KDq6+uVl5f3qccty0p1FwCANEnp20HxeFwHDhzQI488ovvv\nv1/Sjb/+R0ZGEj89Hs+M23G73anEyIjc3FxyOoiczloIOeczo2tkZNZjrUWL5Ha757SOy+VKPLdk\n1ht3ze2t2DXH8dNJaUutra0qLy/Xli1bEss2bNigYDCouro6BYNBVVdXz7idSCSSSoyMcLvd5HQQ\nOZ21EHLOZ8Z4PD7rsfbkpCKRyJzWicfjieeWzHqu+OyuR/j/9eY2fjpJl0B/f7/eeecdLV++XD/4\nwQ9kWZaeeuop1dXVqbm5WYFAQMXFxfL7/Y6FBQA4K+kSWLNmjX77299O+djevXuTDgQAyJyFeZ0z\nAMARlAAAGIwSAACDUQIAYDDuIgrcpjJ5K+TZ7Ms1MqK78qzEfky5VXO2owSA21Qmb4U82319cj+m\n3Ko521ECU7DGo1p0feymZbGxiFwTH029wh13KF6Q3VdrAsBUKIEpuCIjGn/x+ZuWTUwzPnfXHqnq\n/vSGAoA0oARuxZ6cy+C0xQCAdOJMCwAYjBIAAINRAgBgMEoAAAxGCQCAwSgBADAYJQAABqMEAMBg\nlAAAGIwrhtNgZMJW38D1WY1dXZin4jtv7y6e6m6RrpGRKf9B7vm4W+R0d7P875zczRK3G0ogDcZj\ntva9/T+zGnvwKytv/xLI8rtFZns+IJ0ogSzjujokhQZnv4KvWPG7itIXKEstxNdpIWbG7Y8SyDah\nQU281DTr4bk/3CeZ+EaxEF+nhZgZt73b+3MIAMC0KAEAMBglAAAGowQAwGCUAAAYjBIAAIOl7Sui\nPT09ev3112XbtjZt2qS6urp07QoAkKS0HAlMTk7q+PHj2rNnjw4cOKAzZ87o0qVL6dgVACAFaSmB\nc+fOqaysTMXFxcrJydEXv/hFdXR0pGNXAIAUpKUEQqGQCgsLE7/7fD6FQqF07AoAkAJuGzGVO3J1\nx/bnblq0yLI0adtTDl90d7k+eT/MXJelbzxQNqtduRfznwDA/LFs+xbvbCk4e/asTpw4oT179kiS\n2tvbJYmTwwCQZdLycdCqVat05coVDQ4OKhaL6cyZM6qurk7HrgAAKUjLkYB04yuir732mmzb1pe+\n9CWOAgAgC6WtBAAA2Y8rhgHAYJQAABgsrd9PnOnWEWNjYzp8+LCGhoY0OTmpr371q6qtrU08Pjk5\nqd27d8vn86mpafb/IlMmc46NjemVV17Rhx9+KMuy1NDQoNWrV2dVxrfeekuBQECWZWn58uVqbGxU\nTk56/tPPlHN0dFStra3697//rdzcXDU0NKi8vHxW62ZDzuHhYR05ckThcFiWZWnz5s3asmVL1uX8\nWLbMoelyZmoOpZozU/OotbVVXV1d8ng82r9//5RjfvWrX6mnp0eLFy/Wzp07VVFRMavnNyU7TeLx\nuL1r1y57YGDA/uijj+zvfe979sWLF28ac/LkSfs3v/mNbdu2HQ6H7WeffdaOxWKJx0+dOmW3tLTY\nL730UrpippzzyJEj9p/+9Cfbtm07FovZo6OjWZVxeHjY3rlzp/3RRx/Ztm3bBw8etIPBoOMZZ5vz\n17/+tX3ixAnbtm370qVL9k9+8pNZr5sNOa9evWqfP3/etm3bvn79uv2tb30rK3N+LFvm0HQ5MzGH\nUs2ZyXn097//3T5//rz93e9+d8rHu7q67J/97Ge2bdv22bNn7RdffHHWz28qafs4aDa3jrAsS9ev\nX5ckRaNRud1uuVwuSdLw8LC6u7u1efPmdEVMOefY2Jj6+/u1adMmSZLL5VJ+fn5WZZRu/DUYjUYV\nj8c1Pj6uu+66y/GMs8158eJF3XvvvZKkpUuXamBgQP/5z38yequRVHJ6vd7EX115eXlatmxZ2q6G\nTyWnlF1z6FY5MzWHUs0pZW4erVmzRgUFBbd8vKOjQzU1NZKk1atXa2xsTCMjI0nPobSVwGxuHfHY\nY4/p4sWLeu655/T9739f9fX1icfeeOMNPfPMM7IsK10RU845MDAgt9uto0ePqqmpSceOHdPExERW\nZfT5fNq6dasaGxv1/PPPq6CgQFVVVY5nnG3Oz3zmM/rLX/4i6cakHBoa0vDwcEZvNZJKzk8aGBjQ\nhQsX0vbRRao5s2kO3SpnpuZQqjkzOY9mcqvnkewcmtcTwz09Pbrnnnt07Ngx7du3T8ePH1c0Gk18\nHlZRUSHbtmXP87dYb5VzcnJS58+f15e//GXt27dPixcvTlwdnS0ZR0dH1dnZqaNHj+rYsWOKRqN6\n99135yWjdOOq8WvXrqmpqUmnT5/WPffco0WLsu/7CTPljEajOnjwoOrr65WXl5d1ObNtDt0qZzbN\noelyZts8clLaTgz7fD4NDQ0lfg+FQvL5fDeNCQaDiRMXd999t0pKSnTp0iX19/ers7NT3d3dmpiY\n0PXr13XkyBHt2rUrq3IWFhaqsLBQK1eulCRt3LgxLf8Dp5JxcHBQJSUlWrJkiSTpC1/4gt5//309\n/PDD85LzzjvvVGNjY+L3nTt3qrS0VOPj4zOumw05JSkej+vAgQN65JFHdP/996clY7I5d+3apdLS\nUp05cyar5tCtXs9oNJqROZRszo9fz56enozNo5n4fL6bjko/PlKJxWJJzaG0/Qk2m1tHFBUVqbe3\nV5I0MjKif/3rXyotLdX27dvV2tqqI0eO6Nvf/rbuvffetPzPm2pOr9erwsJCXb58WZLU29t70zcz\nsiFjUVGRPvjgA01MTMi2bfX29mrZsmWOZ5xtzrGxMcViMUnSH//4R61du1Z5eXkZvdVIKjmlG9/e\nKC8vT+u3gpLN+bnPfU55eXlZN4du9Xpmag4lm/Pj1zOT80jStEdv1dXVevvttyXduE9bQUGBvF5v\n0nMorVcMT3XriD/84Q+yLEuPPvqorl69qqNHj+rq1auSbhyK/Xez9vX16dSpU2n/eluyOf/5z3/q\n2LFjisViKi0tVWNjY1pObKWS8cSJE/rzn/8sl8uliooKPf/882n9iuh0Oc+ePatf/vKXWrRokcrL\ny9XQ0JB4vTJ5q5Fkc/b39+tHP/qRli9fLsuyZFmWnnrqKX3+85/PqpyflA1zaLqcmZpDqebM1Dxq\naWlRX1+fIpGIPB6Ptm3bplgslsgoScePH1dPT4/y8vLU0NCgFStW3PL5zYTbRgCAwbLvjBwAIGMo\nAQAwGCUAAAajBADAYJQAABiMEgAAg1ECAGAwSgAADPa/pte4cbPfqJQAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "def Histograms(df1,df2):\n", " x=np.array(df1)\n", " print np.min(x)\n", " print np.max(x)\n", " pl.hist(x, 20, normed=1,histtype='bar', rwidth=0.80)\n", " x=np.array(df2)\n", " print np.min(x)\n", " print np.max(x)\n", " pl.hist(x, 20, normed=1, histtype='bar', rwidth=0.80);\n", "\n", "Histograms(ColandIden_Acc[1],ColandIden_Acc[2])" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
Prev_Col+Iden_match12
Sub#
121.0000000.992908
130.8571430.959459
141.0000000.993056
151.0000000.944444
160.9090910.965278
171.0000000.979452
180.9333330.964286
190.8750000.959184
201.0000000.993151
210.8571430.864865
531.0000000.993197
541.0000000.993197
551.0000000.979310
561.0000000.901408
571.0000000.986395
581.0000000.979592
591.0000000.944444
601.0000000.986577
611.0000000.986486
621.0000000.921986
\n", "
" ], "text/plain": [ "Prev_Col+Iden_match 1 2\n", "Sub# \n", "12 1.000000 0.992908\n", "13 0.857143 0.959459\n", "14 1.000000 0.993056\n", "15 1.000000 0.944444\n", "16 0.909091 0.965278\n", "17 1.000000 0.979452\n", "18 0.933333 0.964286\n", "19 0.875000 0.959184\n", "20 1.000000 0.993151\n", "21 0.857143 0.864865\n", "53 1.000000 0.993197\n", "54 1.000000 0.993197\n", "55 1.000000 0.979310\n", "56 1.000000 0.901408\n", "57 1.000000 0.986395\n", "58 1.000000 0.979592\n", "59 1.000000 0.944444\n", "60 1.000000 0.986577\n", "61 1.000000 0.986486\n", "62 1.000000 0.921986" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ColandIden_Acc" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "WilcoxonResult(statistic=45.0, pvalue=0.025068176000876803)" ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" } ], "source": [ "scipy.stats.wilcoxon(np.arcsin(np.sqrt(ColandIden_Acc[1])), y=np.arcsin(np.sqrt(ColandIden_Acc[2])))" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.11" } }, "nbformat": 4, "nbformat_minor": 0 }