From 046285dc01238bc370cd4b0f60e62cc84477513d Mon Sep 17 00:00:00 2001 From: Aaron Sheldon Date: Wed, 8 Feb 2017 20:29:52 -0700 Subject: [PATCH] one time contingency table generation --- contingency.csv | 1025 ++++++++++++++++++++++++++++++++++ longitudinal.csv | 1025 ++++++++++++++++++++++++++++++++++ simple-data-generation.ipynb | 645 +++++++++++++-------- 3 files changed, 2447 insertions(+), 248 deletions(-) create mode 100644 contingency.csv create mode 100644 longitudinal.csv diff --git a/contingency.csv b/contingency.csv new file mode 100644 index 0000000..3a0dafe --- /dev/null +++ b/contingency.csv @@ -0,0 +1,1025 @@ +"sample","dimensionone","dimensiontwo","dimensionthree","dimensionfour" +1,"B","0","C","0" +2,"B","0","C","0" +3,"B","1","C","-1" +4,"B","-1","D","-1" +5,"B","1","F","0" +6,"B","0","F","-2" +7,"B","1","C","-1" +8,"B","-1","C","1" +9,"B","1","C","2" +10,"B","0","E","2" +11,"B","0","F","-2" +12,"A","0","F","-2" +13,"B","1","C","-1" +14,"A","-1","E","-2" +15,"B","1","E","1" +16,"A","-1","F","0" +17,"A","1","F","2" +18,"B","1","C","2" +19,"A","-1","F","-1" +20,"B","-1","C","-2" +21,"B","0","D","0" +22,"B","1","E","1" +23,"B","1","F","-2" +24,"B","1","E","-2" +25,"A","-1","C","-2" +26,"B","1","E","-1" +27,"B","0","C","1" +28,"B","1","F","-2" +29,"B","-1","E","2" +30,"A","-1","D","0" +31,"A","-1","E","-2" +32,"B","1","F","-2" +33,"B","1","E","2" +34,"B","1","F","-2" +35,"B","-1","F","-2" +36,"B","0","D","2" +37,"A","0","E","-2" +38,"A","-1","E","-1" +39,"B","0","D","0" +40,"B","1","C","1" +41,"B","-1","C","1" +42,"B","1","F","1" +43,"B","0","E","1" +44,"A","-1","C","0" +45,"A","-1","C","2" +46,"B","1","C","2" +47,"A","-1","F","0" +48,"B","0","E","2" +49,"A","-1","C","1" +50,"A","-1","E","-1" +51,"B","1","D","0" +52,"B","1","E","0" +53,"A","0","F","0" +54,"A","-1","E","0" +55,"B","0","F","2" +56,"B","0","C","2" +57,"B","-1","E","-1" +58,"A","0","F","-1" +59,"B","1","F","-2" +60,"A","-1","D","1" +61,"A","0","D","1" +62,"B","0","D","2" +63,"A","-1","F","-1" +64,"A","-1","C","2" +65,"A","-1","F","2" +66,"B","1","C","0" +67,"A","-1","C","-1" +68,"B","1","C","-2" +69,"A","-1","F","-1" +70,"B","0","F","-2" +71,"B","0","F","1" +72,"B","0","F","0" +73,"B","0","C","-1" +74,"A","-1","D","-1" +75,"A","0","F","2" +76,"B","0","D","0" +77,"A","-1","F","-2" +78,"A","-1","F","-1" +79,"B","1","E","-1" +80,"B","1","D","0" +81,"B","1","C","1" +82,"A","0","D","1" +83,"A","-1","F","-1" +84,"A","0","C","2" +85,"B","1","D","1" +86,"B","1","E","2" +87,"B","-1","C","0" +88,"B","1","C","0" +89,"B","1","C","-2" +90,"B","-1","C","0" +91,"A","-1","C","1" +92,"B","1","C","-2" +93,"B","1","F","-2" +94,"B","-1","C","2" +95,"A","0","C","1" +96,"A","0","C","0" +97,"A","1","C","1" +98,"B","-1","C","2" +99,"B","1","F","0" +100,"B","-1","C","1" +101,"B","1","D","1" +102,"B","0","F","-1" +103,"A","0","D","0" +104,"B","0","C","1" +105,"A","0","D","-2" +106,"A","1","C","2" +107,"A","-1","D","-1" +108,"A","-1","D","1" +109,"A","-1","F","0" +110,"B","0","C","-1" +111,"B","1","C","0" +112,"B","0","F","2" +113,"B","0","C","-1" +114,"A","-1","C","2" +115,"A","-1","F","-1" +116,"B","0","F","0" +117,"A","0","E","2" +118,"A","-1","E","-1" +119,"A","-1","E","2" +120,"B","1","E","-1" +121,"A","0","E","1" +122,"A","1","C","0" +123,"B","1","F","-2" +124,"A","-1","C","0" +125,"B","1","D","0" +126,"A","-1","E","2" +127,"A","1","F","-1" +128,"B","1","E","-2" +129,"A","-1","F","-2" +130,"B","1","E","0" +131,"B","-1","C","1" +132,"B","1","C","1" +133,"B","1","F","-1" +134,"A","0","F","-1" +135,"A","1","C","1" +136,"A","0","C","2" +137,"B","-1","D","-2" +138,"B","1","E","0" +139,"B","1","F","1" +140,"A","-1","C","0" +141,"A","-1","F","-2" +142,"A","-1","E","2" +143,"A","-1","E","-2" +144,"A","-1","F","-2" +145,"A","-1","F","-2" +146,"B","1","F","0" +147,"A","-1","F","-2" +148,"A","-1","C","0" +149,"A","-1","F","-1" +150,"B","-1","E","2" +151,"A","-1","E","-1" +152,"B","-1","F","-2" +153,"B","0","E","1" +154,"A","1","F","-1" +155,"B","0","C","2" +156,"A","-1","F","1" +157,"A","-1","F","-1" +158,"B","1","E","2" +159,"A","0","C","0" +160,"B","0","C","2" +161,"B","1","C","2" +162,"B","1","C","1" +163,"A","0","F","-2" +164,"A","-1","E","-1" +165,"A","0","D","-1" +166,"B","1","E","-2" +167,"A","0","E","-1" +168,"A","-1","C","2" +169,"A","1","C","2" +170,"A","-1","F","-2" +171,"B","0","F","0" +172,"A","-1","C","0" +173,"B","1","F","0" +174,"A","0","D","2" +175,"A","-1","C","0" +176,"A","-1","C","2" +177,"B","1","C","-2" +178,"A","-1","F","-2" +179,"B","-1","D","-2" +180,"A","-1","E","-2" +181,"B","-1","F","-1" +182,"A","-1","D","0" +183,"B","1","E","-2" +184,"B","1","E","1" +185,"A","-1","D","-1" +186,"B","1","D","1" +187,"A","1","F","-2" +188,"A","0","E","-1" +189,"A","1","F","-2" +190,"A","1","C","1" +191,"B","1","F","-1" +192,"B","0","C","1" +193,"B","0","D","0" +194,"A","-1","F","-2" +195,"B","1","F","-2" +196,"A","-1","C","-2" +197,"A","-1","C","2" +198,"B","-1","F","-1" +199,"A","0","F","1" +200,"B","1","F","1" +201,"A","-1","F","-1" +202,"B","0","F","-2" +203,"B","0","C","2" +204,"A","-1","C","0" +205,"A","0","F","-1" +206,"B","-1","F","0" +207,"B","0","F","-2" +208,"B","1","C","1" +209,"A","-1","F","-2" +210,"A","1","F","-2" +211,"B","0","C","-1" +212,"B","1","E","0" +213,"A","1","F","0" +214,"A","-1","F","0" +215,"A","0","F","-1" +216,"A","-1","C","2" +217,"B","0","D","1" +218,"A","-1","D","0" +219,"A","0","D","-1" +220,"A","1","E","2" +221,"A","-1","E","-2" +222,"A","-1","E","-1" +223,"A","0","E","0" +224,"B","1","D","-2" +225,"A","-1","C","0" +226,"B","1","F","2" +227,"B","1","F","-1" +228,"A","-1","C","1" +229,"B","0","C","1" +230,"B","0","F","0" +231,"A","0","F","0" +232,"B","-1","D","0" +233,"B","-1","C","1" +234,"A","0","C","2" +235,"A","-1","D","-1" +236,"A","1","C","2" +237,"B","-1","E","-2" +238,"B","0","E","-2" +239,"B","1","C","1" +240,"A","0","C","0" +241,"A","0","C","2" +242,"A","-1","C","0" +243,"B","-1","F","-2" +244,"B","1","C","1" +245,"A","0","C","1" +246,"A","-1","F","0" +247,"B","0","E","-2" +248,"A","-1","D","1" +249,"A","1","E","-1" +250,"B","1","C","1" +251,"B","0","E","2" +252,"B","0","F","-1" +253,"B","0","C","2" +254,"A","-1","C","-1" +255,"A","-1","F","-2" +256,"B","0","C","1" +257,"B","0","D","0" +258,"B","1","C","1" +259,"A","-1","D","1" +260,"B","0","D","-1" +261,"B","0","E","-2" +262,"A","1","E","-2" +263,"B","0","D","-1" +264,"B","1","F","-2" +265,"A","-1","F","-1" +266,"B","1","F","-2" +267,"A","0","C","-1" +268,"A","-1","F","-1" +269,"A","0","F","-2" +270,"B","1","D","0" +271,"A","0","C","1" +272,"B","1","F","0" +273,"B","1","C","1" +274,"B","1","C","1" +275,"A","-1","C","2" +276,"B","0","C","1" +277,"B","1","E","0" +278,"B","0","E","-2" +279,"B","0","F","-1" +280,"A","-1","F","1" +281,"B","0","F","-1" +282,"A","1","E","2" +283,"A","-1","F","1" +284,"B","1","C","0" +285,"B","-1","C","-1" +286,"B","1","E","-2" +287,"A","0","C","1" +288,"A","-1","C","0" +289,"B","0","E","-1" +290,"A","0","E","2" +291,"A","-1","C","2" +292,"A","-1","F","1" +293,"B","-1","F","0" +294,"B","0","D","2" +295,"B","1","E","1" +296,"A","-1","F","0" +297,"A","-1","F","-1" +298,"B","1","E","2" +299,"A","-1","F","-1" +300,"A","-1","E","-2" +301,"B","0","F","-1" +302,"B","-1","F","-1" +303,"B","0","C","0" +304,"B","1","F","1" +305,"A","-1","E","2" +306,"B","1","E","2" +307,"A","0","D","0" +308,"B","1","F","-2" +309,"A","0","E","-2" +310,"A","0","F","2" +311,"B","0","E","-2" +312,"A","0","D","0" +313,"B","1","C","2" +314,"B","0","F","1" +315,"B","1","E","2" +316,"B","0","F","0" +317,"A","-1","E","0" +318,"B","1","F","1" +319,"B","-1","C","0" +320,"A","0","F","0" +321,"B","1","F","-2" +322,"A","-1","D","-1" +323,"B","-1","C","1" +324,"A","-1","F","-1" +325,"B","1","E","-2" +326,"B","1","F","1" +327,"A","-1","C","2" +328,"B","0","E","1" +329,"A","0","C","-1" +330,"B","1","C","2" +331,"A","-1","F","-1" +332,"A","0","C","0" +333,"B","1","D","0" +334,"B","-1","F","1" +335,"B","1","E","-1" +336,"A","-1","C","0" +337,"A","-1","E","-2" +338,"B","0","E","-1" +339,"B","1","C","-1" +340,"A","0","C","1" +341,"A","-1","C","-2" +342,"B","0","C","2" +343,"B","0","F","0" +344,"B","1","D","2" +345,"B","-1","C","-1" +346,"B","1","C","0" +347,"A","-1","E","-2" +348,"B","0","C","1" +349,"A","-1","D","0" +350,"A","1","E","2" +351,"A","-1","C","2" +352,"B","1","C","2" +353,"B","1","F","-2" +354,"A","-1","C","-2" +355,"A","-1","E","2" +356,"B","1","C","2" +357,"B","0","F","-2" +358,"A","0","E","2" +359,"B","0","E","0" +360,"A","1","F","-2" +361,"A","-1","F","-2" +362,"B","-1","F","-1" +363,"A","-1","E","0" +364,"A","-1","E","0" +365,"B","1","D","0" +366,"A","0","F","-1" +367,"A","-1","C","1" +368,"A","-1","D","0" +369,"B","1","C","-1" +370,"B","1","F","-2" +371,"A","1","D","1" +372,"B","-1","D","0" +373,"A","-1","F","-2" +374,"A","-1","F","0" +375,"B","0","F","-1" +376,"B","0","F","-2" +377,"A","-1","E","1" +378,"B","1","F","-2" +379,"B","1","C","1" +380,"A","-1","F","0" +381,"A","-1","E","-2" +382,"A","-1","D","0" +383,"B","1","D","-1" +384,"A","-1","E","-1" +385,"A","-1","C","0" +386,"A","1","D","0" +387,"A","-1","D","0" +388,"B","1","F","-1" +389,"B","1","E","-2" +390,"B","1","C","1" +391,"B","1","F","-1" +392,"B","0","F","-2" +393,"A","0","E","2" +394,"B","1","F","1" +395,"A","0","E","-1" +396,"A","-1","E","0" +397,"B","1","E","2" +398,"B","-1","C","1" +399,"B","0","E","-1" +400,"A","-1","F","-2" +401,"B","-1","F","-2" +402,"A","0","C","1" +403,"A","-1","D","0" +404,"B","1","D","0" +405,"B","0","F","-2" +406,"B","1","C","1" +407,"B","0","F","0" +408,"A","0","C","2" +409,"A","-1","D","0" +410,"A","0","D","0" +411,"A","-1","E","0" +412,"A","1","E","-1" +413,"A","-1","F","-2" +414,"B","1","E","1" +415,"B","0","D","1" +416,"A","-1","D","-1" +417,"A","1","E","-2" +418,"B","1","C","1" +419,"A","-1","F","0" +420,"A","-1","C","2" +421,"B","0","C","2" +422,"B","0","C","1" +423,"A","1","F","-2" +424,"A","1","C","0" +425,"A","1","C","2" +426,"A","0","E","-1" +427,"A","-1","D","-1" +428,"A","1","F","0" +429,"B","0","E","-2" +430,"A","0","C","2" +431,"B","1","C","-1" +432,"A","0","E","1" +433,"A","1","D","1" +434,"A","1","C","-2" +435,"A","-1","F","2" +436,"A","-1","D","1" +437,"A","-1","D","2" +438,"A","0","C","0" +439,"A","1","C","2" +440,"A","-1","D","0" +441,"A","-1","F","-2" +442,"A","-1","F","0" +443,"A","-1","D","0" +444,"B","1","E","2" +445,"A","-1","F","-2" +446,"B","0","E","1" +447,"B","-1","C","1" +448,"A","-1","F","-2" +449,"A","1","C","1" +450,"B","0","F","0" +451,"B","0","E","-2" +452,"B","-1","F","0" +453,"A","-1","C","0" +454,"B","1","E","-2" +455,"A","-1","D","0" +456,"A","-1","F","0" +457,"B","1","F","1" +458,"A","0","F","1" +459,"A","-1","E","2" +460,"B","0","C","0" +461,"B","1","C","1" +462,"A","0","D","0" +463,"A","0","F","-2" +464,"B","-1","F","-2" +465,"A","-1","C","1" +466,"A","-1","C","1" +467,"B","1","F","-2" +468,"A","1","E","1" +469,"B","-1","C","2" +470,"A","-1","D","0" +471,"B","1","F","-2" +472,"A","-1","F","-1" +473,"A","0","F","-1" +474,"B","1","C","2" +475,"B","1","D","-1" +476,"B","0","D","0" +477,"B","-1","D","0" +478,"B","1","F","-1" +479,"B","1","F","2" +480,"B","1","F","-2" +481,"A","1","C","0" +482,"B","1","C","-2" +483,"A","0","C","2" +484,"A","1","E","-2" +485,"A","1","E","-2" +486,"A","-1","F","-1" +487,"B","-1","C","2" +488,"A","0","D","0" +489,"B","1","C","-1" +490,"B","1","F","-2" +491,"A","-1","E","2" +492,"B","0","D","1" +493,"B","1","C","2" +494,"A","0","C","-1" +495,"B","0","F","-1" +496,"A","-1","D","2" +497,"A","-1","F","-1" +498,"B","1","C","-1" +499,"A","0","E","-2" +500,"A","-1","E","-2" +501,"B","0","E","0" +502,"B","0","D","0" +503,"A","-1","F","0" +504,"B","1","F","-2" +505,"B","-1","D","2" +506,"B","0","F","-1" +507,"A","1","E","1" +508,"A","-1","C","1" +509,"A","-1","F","1" +510,"B","0","C","-2" +511,"B","0","F","-1" +512,"B","1","E","-2" +513,"A","0","C","2" +514,"A","-1","C","2" +515,"B","-1","C","-1" +516,"A","0","C","-1" +517,"B","1","D","0" +518,"A","-1","D","-1" +519,"B","1","C","-2" +520,"A","-1","F","0" +521,"A","1","D","2" +522,"B","1","C","-1" +523,"A","1","F","0" +524,"B","1","E","-2" +525,"A","0","C","2" +526,"A","0","E","0" +527,"A","0","F","-2" +528,"A","0","C","0" +529,"B","0","F","-1" +530,"A","-1","F","1" +531,"A","0","C","0" +532,"B","-1","F","-1" +533,"A","1","F","0" +534,"B","0","F","-2" +535,"A","0","F","1" +536,"B","0","E","-2" +537,"B","-1","F","-2" +538,"B","0","D","-1" +539,"B","1","C","0" +540,"B","-1","E","-2" +541,"A","-1","E","-2" +542,"B","1","E","-2" +543,"A","-1","E","2" +544,"B","0","F","-1" +545,"B","0","C","2" +546,"A","-1","C","2" +547,"A","0","E","2" +548,"B","0","C","2" +549,"B","0","D","1" +550,"A","1","C","-2" +551,"B","0","F","-2" +552,"B","1","F","-1" +553,"A","0","E","-1" +554,"A","1","F","0" +555,"B","1","F","1" +556,"B","1","E","2" +557,"B","1","D","-1" +558,"A","1","C","2" +559,"B","-1","E","-1" +560,"B","0","E","-2" +561,"B","1","C","-1" +562,"A","-1","E","-2" +563,"B","1","D","-1" +564,"B","1","C","1" +565,"A","0","E","2" +566,"B","1","F","-1" +567,"B","1","F","-1" +568,"B","0","F","-1" +569,"B","0","F","-1" +570,"B","0","C","1" +571,"A","-1","E","-1" +572,"B","-1","C","2" +573,"B","-1","E","-1" +574,"A","-1","C","0" +575,"B","-1","C","2" +576,"B","0","E","2" +577,"B","-1","C","1" +578,"B","1","F","-2" +579,"A","0","C","2" +580,"B","-1","E","2" +581,"A","-1","F","2" +582,"B","1","E","2" +583,"A","0","F","-2" +584,"B","-1","E","2" +585,"A","0","E","0" +586,"A","-1","F","-1" +587,"A","-1","C","0" +588,"A","0","F","2" +589,"A","0","D","-2" +590,"B","0","D","-1" +591,"A","-1","E","2" +592,"B","1","C","-2" +593,"B","1","E","2" +594,"A","-1","E","-1" +595,"B","1","D","-1" +596,"A","0","F","-2" +597,"B","-1","C","0" +598,"A","-1","C","2" +599,"A","0","F","2" +600,"B","0","F","2" +601,"A","1","C","2" +602,"A","1","E","2" +603,"A","1","F","0" +604,"A","1","C","-1" +605,"A","0","E","-2" +606,"A","-1","E","-1" +607,"B","-1","F","0" +608,"A","0","D","0" +609,"A","-1","C","2" +610,"B","-1","F","1" +611,"B","-1","F","-2" +612,"B","1","E","1" +613,"B","0","E","-2" +614,"B","1","F","-2" +615,"B","1","F","-2" +616,"B","1","C","0" +617,"A","-1","F","-2" +618,"B","0","C","-1" +619,"A","-1","D","0" +620,"B","0","C","2" +621,"B","1","D","0" +622,"B","0","D","2" +623,"A","0","F","-2" +624,"A","-1","E","0" +625,"A","-1","D","-2" +626,"A","-1","E","-2" +627,"B","-1","C","0" +628,"A","-1","C","2" +629,"A","0","C","-1" +630,"A","-1","D","2" +631,"B","0","E","1" +632,"A","-1","F","-1" +633,"A","1","E","-1" +634,"A","-1","C","0" +635,"A","-1","E","2" +636,"A","0","C","-2" +637,"A","-1","D","-1" +638,"A","-1","F","-1" +639,"B","-1","E","-2" +640,"B","1","D","-1" +641,"B","0","F","-2" +642,"A","-1","E","1" +643,"A","-1","D","1" +644,"A","-1","D","0" +645,"B","0","C","1" +646,"B","-1","E","-1" +647,"A","1","F","1" +648,"A","-1","E","2" +649,"B","-1","F","2" +650,"B","1","C","2" +651,"A","-1","F","-2" +652,"B","0","C","1" +653,"B","1","F","0" +654,"A","1","C","-1" +655,"A","0","F","-1" +656,"B","0","D","2" +657,"B","-1","F","-2" +658,"A","-1","E","-1" +659,"B","1","D","-1" +660,"A","0","D","0" +661,"A","0","D","1" +662,"B","-1","F","-2" +663,"A","-1","E","1" +664,"A","-1","C","-1" +665,"B","1","C","-1" +666,"B","-1","D","-2" +667,"A","-1","E","-1" +668,"B","1","E","1" +669,"B","-1","E","-1" +670,"A","-1","F","-2" +671,"A","0","F","-2" +672,"B","-1","F","-2" +673,"A","0","C","2" +674,"B","-1","D","0" +675,"B","-1","C","-2" +676,"B","-1","C","2" +677,"A","1","D","1" +678,"B","0","D","1" +679,"A","-1","E","1" +680,"B","1","F","-2" +681,"B","1","C","-1" +682,"B","1","C","0" +683,"A","-1","C","2" +684,"B","0","D","-1" +685,"A","-1","C","1" +686,"A","-1","D","1" +687,"B","0","C","1" +688,"B","1","E","2" +689,"A","-1","C","2" +690,"B","-1","D","0" +691,"A","-1","C","0" +692,"B","0","C","2" +693,"A","0","F","1" +694,"B","1","F","-1" +695,"B","-1","F","2" +696,"A","0","C","1" +697,"A","-1","D","1" +698,"A","-1","C","-1" +699,"B","1","F","-2" +700,"A","-1","E","-1" +701,"B","-1","E","-2" +702,"B","-1","C","2" +703,"A","0","E","0" +704,"B","1","C","1" +705,"A","1","F","-1" +706,"B","1","F","-1" +707,"B","0","F","1" +708,"A","-1","F","-1" +709,"A","0","F","-2" +710,"A","0","E","-2" +711,"B","1","D","1" +712,"B","1","F","0" +713,"A","1","D","0" +714,"A","-1","D","1" +715,"B","0","F","-1" +716,"B","1","F","-2" +717,"A","-1","C","1" +718,"A","1","C","1" +719,"A","-1","F","-2" +720,"B","0","C","1" +721,"A","1","F","-2" +722,"A","1","C","2" +723,"B","1","C","0" +724,"B","0","E","0" +725,"A","-1","E","2" +726,"B","1","D","-1" +727,"A","-1","C","1" +728,"B","0","E","0" +729,"B","0","E","-2" +730,"A","-1","F","-2" +731,"B","-1","C","2" +732,"A","1","C","2" +733,"B","1","C","2" +734,"A","-1","D","1" +735,"A","0","C","1" +736,"B","1","E","1" +737,"A","1","E","2" +738,"A","-1","F","1" +739,"A","-1","C","2" +740,"A","-1","C","1" +741,"B","1","C","2" +742,"A","0","F","1" +743,"A","-1","C","1" +744,"A","0","E","2" +745,"B","1","E","1" +746,"B","-1","E","1" +747,"A","-1","E","2" +748,"B","0","F","-2" +749,"A","0","F","-2" +750,"B","1","C","2" +751,"A","0","F","2" +752,"A","-1","C","2" +753,"A","-1","C","-1" +754,"A","1","D","1" +755,"A","-1","D","1" +756,"B","0","F","-2" +757,"B","0","C","2" +758,"B","0","C","0" +759,"A","0","C","-1" +760,"B","1","C","0" +761,"B","1","E","-1" +762,"B","-1","F","-1" +763,"B","0","C","2" +764,"A","0","F","0" +765,"A","0","E","0" +766,"B","1","C","2" +767,"A","-1","D","1" +768,"A","-1","F","-2" +769,"A","0","C","1" +770,"B","0","F","-1" +771,"B","1","F","-2" +772,"B","-1","E","2" +773,"B","1","C","1" +774,"A","-1","E","1" +775,"B","0","C","-1" +776,"B","-1","F","0" +777,"B","1","C","1" +778,"B","1","C","1" +779,"A","1","D","1" +780,"B","0","C","-2" +781,"A","0","D","2" +782,"B","1","C","1" +783,"A","0","D","0" +784,"A","-1","D","2" +785,"B","-1","D","1" +786,"A","-1","E","-2" +787,"B","1","C","0" +788,"A","-1","E","-2" +789,"B","1","C","-2" +790,"A","-1","C","-1" +791,"B","0","F","-1" +792,"A","0","C","0" +793,"A","-1","D","0" +794,"B","1","E","-1" +795,"B","1","C","-1" +796,"A","0","C","-2" +797,"A","-1","F","1" +798,"A","0","F","0" +799,"A","0","F","1" +800,"B","0","F","-1" +801,"B","0","E","2" +802,"A","-1","C","2" +803,"A","1","D","0" +804,"A","1","C","1" +805,"A","-1","C","1" +806,"B","1","D","-1" +807,"A","0","F","-2" +808,"B","-1","F","0" +809,"B","0","F","-1" +810,"A","-1","F","1" +811,"A","-1","C","2" +812,"B","0","D","1" +813,"A","-1","E","2" +814,"A","1","E","2" +815,"B","1","F","-2" +816,"B","1","F","2" +817,"A","-1","F","-2" +818,"B","1","C","1" +819,"A","-1","D","-2" +820,"B","1","C","1" +821,"B","1","F","1" +822,"B","0","D","0" +823,"A","1","C","1" +824,"A","0","E","0" +825,"B","1","E","1" +826,"A","0","E","2" +827,"A","0","E","-2" +828,"A","0","C","0" +829,"B","-1","F","-2" +830,"A","1","D","0" +831,"B","0","C","1" +832,"A","-1","D","1" +833,"B","1","C","0" +834,"A","-1","E","1" +835,"A","1","E","-2" +836,"B","0","F","-2" +837,"A","-1","F","-2" +838,"A","-1","C","2" +839,"B","-1","F","1" +840,"A","-1","E","2" +841,"A","-1","C","0" +842,"A","1","C","0" +843,"A","0","F","0" +844,"A","0","F","-2" +845,"A","0","E","-2" +846,"B","-1","E","2" +847,"B","1","E","-1" +848,"A","-1","C","1" +849,"A","-1","F","0" +850,"B","1","E","1" +851,"B","1","D","0" +852,"B","-1","E","1" +853,"B","1","F","-2" +854,"B","1","E","-1" +855,"B","-1","D","-2" +856,"A","-1","C","1" +857,"B","0","C","1" +858,"B","1","D","-1" +859,"B","0","F","1" +860,"A","-1","F","-2" +861,"A","-1","F","0" +862,"B","1","D","-1" +863,"B","1","F","-2" +864,"A","-1","C","2" +865,"A","-1","D","1" +866,"A","0","E","1" +867,"A","0","D","-1" +868,"A","0","C","-1" +869,"A","-1","F","-2" +870,"A","0","D","-1" +871,"A","0","C","2" +872,"B","0","F","1" +873,"A","-1","E","-1" +874,"B","0","C","1" +875,"A","0","F","1" +876,"B","1","E","1" +877,"A","-1","F","-2" +878,"A","-1","D","1" +879,"A","0","E","1" +880,"B","0","F","1" +881,"A","-1","E","2" +882,"B","0","F","-2" +883,"A","-1","C","2" +884,"B","0","D","2" +885,"A","-1","C","0" +886,"A","-1","C","-2" +887,"A","-1","F","-2" +888,"B","0","C","1" +889,"B","1","C","-1" +890,"A","0","F","0" +891,"A","0","F","1" +892,"A","-1","D","-2" +893,"B","0","F","2" +894,"B","1","E","-2" +895,"A","0","C","-2" +896,"A","-1","C","0" +897,"A","0","F","-2" +898,"A","-1","E","1" +899,"A","1","F","-2" +900,"A","-1","D","1" +901,"A","0","E","1" +902,"B","1","F","-1" +903,"A","-1","E","-2" +904,"A","-1","F","0" +905,"A","-1","F","-1" +906,"A","-1","F","-2" +907,"B","1","C","2" +908,"A","-1","D","2" +909,"B","-1","E","-2" +910,"A","-1","D","-1" +911,"B","0","D","1" +912,"B","0","C","2" +913,"B","1","C","2" +914,"B","1","C","-2" +915,"A","-1","E","-2" +916,"B","1","F","-2" +917,"A","-1","F","-1" +918,"B","0","D","2" +919,"B","-1","E","-2" +920,"A","0","D","-2" +921,"A","-1","F","0" +922,"A","1","F","1" +923,"B","1","C","-1" +924,"B","-1","C","-2" +925,"B","0","F","0" +926,"B","1","F","-1" +927,"B","1","F","0" +928,"A","-1","F","-2" +929,"A","-1","C","-2" +930,"A","1","C","-2" +931,"A","-1","C","1" +932,"A","-1","C","0" +933,"A","-1","D","1" +934,"A","0","C","-1" +935,"A","-1","E","2" +936,"B","-1","D","0" +937,"B","0","C","1" +938,"B","1","C","2" +939,"A","0","D","0" +940,"B","0","E","2" +941,"B","1","E","-2" +942,"A","1","F","-2" +943,"A","1","F","-1" +944,"A","0","F","-2" +945,"A","-1","E","1" +946,"A","0","C","-2" +947,"A","-1","C","2" +948,"A","-1","C","2" +949,"A","-1","F","-2" +950,"B","0","C","-1" +951,"B","-1","F","-1" +952,"A","-1","E","-2" +953,"A","0","E","-2" +954,"B","1","F","0" +955,"B","1","D","1" +956,"A","-1","C","2" +957,"A","-1","E","-1" +958,"A","0","C","2" +959,"A","0","D","-2" +960,"A","-1","F","-1" +961,"B","0","C","2" +962,"A","-1","F","1" +963,"A","-1","F","1" +964,"B","-1","D","0" +965,"A","-1","F","0" +966,"B","0","F","-1" +967,"B","-1","F","-1" +968,"A","0","C","2" +969,"A","0","C","0" +970,"B","1","F","-2" +971,"B","-1","E","-1" +972,"B","0","C","1" +973,"A","-1","C","0" +974,"B","1","C","2" +975,"B","1","F","0" +976,"A","1","E","-1" +977,"B","0","C","2" +978,"A","0","C","1" +979,"A","0","F","-1" +980,"A","0","E","-2" +981,"A","0","F","0" +982,"A","0","D","1" +983,"B","-1","C","-1" +984,"A","1","C","-2" +985,"B","-1","F","-2" +986,"A","0","F","0" +987,"A","0","F","-2" +988,"B","-1","F","0" +989,"A","0","E","-2" +990,"B","0","C","2" +991,"A","-1","C","0" +992,"A","-1","C","1" +993,"B","0","D","-1" +994,"B","0","D","1" +995,"B","-1","C","-1" +996,"B","1","F","-2" +997,"B","1","C","2" +998,"B","0","D","2" +999,"B","1","C","-1" +1000,"B","0","F","2" +1001,"B","1","D","1" +1002,"A","0","C","-2" +1003,"B","1","F","-2" +1004,"A","-1","C","1" +1005,"B","0","F","-1" +1006,"A","1","E","-2" +1007,"B","1","E","1" +1008,"B","0","E","2" +1009,"B","0","C","1" +1010,"A","1","C","2" +1011,"A","-1","F","-1" +1012,"A","0","E","-1" +1013,"A","-1","C","1" +1014,"A","0","D","-2" +1015,"A","-1","C","-1" +1016,"A","-1","C","1" +1017,"B","1","F","0" +1018,"A","-1","E","-2" +1019,"A","1","D","-2" +1020,"A","-1","C","1" +1021,"B","-1","C","2" +1022,"B","1","E","2" +1023,"A","-1","F","-1" +1024,"B","0","E","2" diff --git a/longitudinal.csv b/longitudinal.csv new file mode 100644 index 0000000..8f01b3e --- /dev/null +++ b/longitudinal.csv @@ -0,0 +1,1025 @@ +"subject","time","event" +1,5,"C" +1,7,"B" +1,15,"D" +1,20,"B" +1,24,"C" +1,32,"C" +1,40,"C" +1,45,"C" +1,51,"C" +1,57,"C" +1,65,"D" +1,70,"D" +1,74,"C" +1,75,"D" +1,83,"A" +1,86,"C" +1,89,"B" +1,91,"B" +1,92,"C" +1,94,"A" +1,99,"C" +1,106,"D" +2,7,"D" +2,10,"D" +2,15,"C" +2,22,"A" +2,27,"E" +2,28,"E" +2,36,"D" +2,41,"B" +2,44,"B" +2,48,"A" +2,51,"C" +2,56,"C" +2,57,"C" +2,59,"C" +2,61,"E" +2,69,"D" +2,77,"B" +2,85,"C" +2,90,"C" +3,1,"C" +3,3,"A" +3,7,"C" +3,11,"C" +3,14,"C" +3,22,"D" +3,28,"C" +3,29,"E" +3,30,"C" +3,32,"C" +3,39,"C" +3,45,"E" +4,1,"C" +4,7,"C" +4,15,"C" +4,22,"C" +4,30,"B" +4,32,"D" +4,35,"C" +4,40,"E" +4,46,"C" +4,53,"A" +4,55,"C" +4,63,"D" +4,69,"C" +5,7,"B" +5,14,"A" +5,19,"B" +5,24,"B" +5,26,"E" +5,31,"D" +5,37,"B" +5,43,"C" +5,48,"D" +5,52,"B" +5,53,"C" +5,54,"B" +5,61,"C" +5,68,"D" +5,75,"C" +5,83,"C" +6,5,"C" +6,9,"C" +6,15,"C" +6,16,"E" +6,21,"C" +6,26,"C" +6,34,"C" +6,42,"C" +6,43,"B" +6,45,"C" +6,50,"C" +6,58,"B" +6,59,"D" +6,67,"C" +6,74,"B" +6,82,"C" +6,90,"C" +6,98,"C" +6,105,"B" +7,2,"C" +7,4,"E" +7,11,"E" +7,16,"B" +7,21,"C" +7,25,"B" +7,27,"C" +8,3,"D" +8,4,"D" +8,9,"C" +8,14,"C" +8,20,"E" +8,25,"D" +8,30,"D" +8,32,"A" +8,34,"E" +8,40,"A" +8,48,"E" +8,56,"C" +8,62,"C" +8,70,"C" +8,78,"C" +8,80,"C" +9,6,"D" +9,7,"A" +9,10,"D" +9,16,"C" +9,18,"C" +9,21,"B" +9,27,"C" +9,31,"D" +9,36,"B" +9,39,"B" +9,44,"D" +9,52,"C" +9,58,"B" +9,61,"D" +10,7,"C" +10,13,"B" +10,17,"A" +10,23,"D" +10,28,"B" +10,35,"A" +10,36,"D" +10,42,"D" +10,45,"D" +10,53,"B" +10,57,"A" +10,62,"C" +10,69,"A" +10,77,"B" +10,80,"C" +10,83,"D" +10,88,"C" +11,4,"E" +11,11,"A" +11,14,"E" +11,15,"D" +11,22,"C" +11,30,"C" +11,32,"C" +11,39,"E" +11,42,"C" +11,45,"D" +11,46,"C" +11,47,"C" +11,48,"C" +11,55,"B" +11,63,"B" +11,64,"B" +11,71,"C" +11,72,"D" +12,5,"D" +12,9,"C" +12,10,"C" +12,14,"C" +12,15,"C" +12,19,"B" +12,27,"E" +12,35,"D" +12,37,"C" +12,42,"A" +12,45,"C" +12,47,"C" +12,54,"C" +12,56,"C" +12,57,"C" +12,63,"D" +12,66,"C" +12,70,"D" +12,72,"C" +12,73,"C" +12,75,"C" +12,80,"C" +13,5,"C" +13,9,"C" +13,14,"A" +13,18,"B" +13,23,"C" +13,27,"C" +13,30,"E" +13,37,"B" +13,40,"E" +13,43,"C" +13,46,"C" +13,49,"C" +13,52,"C" +13,59,"B" +13,63,"B" +13,64,"A" +13,68,"A" +13,76,"C" +14,7,"D" +14,10,"C" +14,15,"C" +14,20,"E" +14,23,"C" +14,24,"C" +14,27,"C" +14,30,"B" +14,31,"B" +14,34,"D" +14,36,"C" +14,38,"E" +14,46,"A" +14,48,"C" +14,54,"B" +14,56,"C" +14,63,"B" +14,64,"C" +14,65,"C" +14,69,"B" +14,72,"B" +15,7,"C" +15,15,"C" +15,16,"C" +15,18,"C" +15,22,"E" +15,30,"D" +15,33,"C" +15,35,"C" +15,43,"C" +15,44,"C" +15,48,"C" +15,49,"C" +15,53,"C" +15,57,"D" +15,63,"D" +15,68,"D" +15,70,"D" +15,78,"D" +15,83,"B" +15,84,"C" +16,1,"C" +16,2,"A" +16,5,"C" +16,7,"D" +16,9,"A" +16,12,"C" +16,13,"C" +16,16,"D" +16,24,"D" +16,27,"C" +16,32,"C" +16,40,"B" +16,42,"C" +16,50,"C" +16,58,"B" +16,62,"C" +16,69,"C" +16,77,"C" +16,82,"B" +17,3,"D" +17,10,"A" +17,13,"B" +17,21,"C" +17,27,"D" +17,28,"C" +17,33,"D" +17,36,"B" +17,39,"A" +17,40,"C" +18,6,"B" +18,9,"C" +18,14,"C" +18,22,"C" +18,25,"C" +18,33,"B" +18,39,"B" +18,40,"D" +18,44,"D" +18,52,"E" +18,56,"D" +19,7,"C" +19,9,"E" +19,13,"E" +19,15,"C" +19,20,"A" +19,23,"C" +19,28,"A" +19,33,"C" +19,34,"E" +19,39,"A" +19,43,"A" +19,46,"D" +19,51,"C" +19,53,"A" +19,56,"B" +19,64,"B" +20,6,"B" +20,10,"D" +20,11,"B" +20,12,"D" +20,19,"D" +20,22,"C" +20,25,"C" +20,32,"C" +20,38,"B" +20,44,"B" +20,45,"B" +20,52,"E" +20,55,"D" +20,58,"A" +20,65,"B" +20,66,"C" +21,6,"C" +21,10,"A" +21,17,"D" +21,23,"E" +21,24,"D" +21,27,"C" +21,35,"C" +21,38,"B" +21,40,"C" +21,42,"B" +21,47,"B" +21,51,"D" +21,55,"C" +21,62,"E" +21,67,"C" +22,5,"A" +22,7,"C" +22,9,"C" +22,17,"D" +22,20,"D" +22,21,"C" +22,27,"B" +22,35,"C" +22,41,"A" +22,45,"D" +22,47,"D" +22,52,"D" +22,57,"C" +23,5,"D" +23,10,"C" +23,17,"C" +23,23,"C" +23,31,"B" +23,37,"B" +23,42,"D" +23,43,"C" +23,48,"C" +23,51,"E" +23,58,"C" +23,61,"B" +23,65,"C" +24,4,"D" +24,5,"E" +24,10,"C" +24,12,"E" +24,18,"D" +24,23,"B" +24,28,"B" +24,29,"C" +24,31,"A" +24,33,"B" +24,35,"D" +24,36,"C" +24,42,"C" +25,4,"A" +25,7,"B" +25,13,"A" +25,18,"B" +25,25,"B" +25,29,"C" +25,34,"C" +25,37,"D" +25,40,"C" +25,43,"D" +25,45,"C" +25,48,"D" +26,2,"B" +26,7,"D" +26,15,"E" +26,17,"D" +26,21,"C" +26,26,"E" +26,32,"C" +26,34,"C" +26,36,"C" +26,37,"C" +27,6,"C" +27,14,"B" +27,16,"D" +27,23,"D" +27,28,"E" +27,30,"C" +27,36,"B" +27,40,"E" +27,42,"D" +27,43,"C" +27,48,"E" +27,56,"C" +27,58,"E" +27,62,"D" +27,65,"C" +27,72,"C" +28,2,"B" +28,7,"D" +28,10,"C" +28,12,"D" +28,18,"C" +28,22,"C" +28,24,"E" +28,25,"C" +28,30,"C" +28,38,"B" +28,44,"C" +28,49,"D" +28,55,"E" +28,58,"B" +28,64,"B" +29,1,"C" +29,9,"E" +29,12,"C" +29,16,"A" +29,24,"D" +29,29,"C" +29,30,"D" +29,38,"A" +29,46,"C" +29,49,"A" +29,52,"C" +29,57,"C" +29,59,"E" +29,63,"C" +29,64,"D" +30,3,"E" +30,11,"D" +30,18,"A" +30,23,"D" +30,29,"B" +30,37,"B" +30,43,"D" +30,51,"C" +30,55,"E" +30,59,"C" +30,66,"C" +30,70,"C" +30,78,"C" +30,84,"C" +30,85,"B" +30,93,"E" +30,98,"C" +30,101,"E" +30,109,"C" +30,110,"C" +30,117,"B" +30,118,"B" +30,120,"A" +30,127,"C" +31,2,"E" +31,9,"E" +31,13,"B" +31,14,"C" +31,22,"C" +31,28,"E" +31,33,"B" +31,39,"B" +31,45,"C" +31,46,"C" +31,53,"C" +31,58,"A" +31,64,"D" +31,67,"B" +31,70,"A" +31,74,"D" +31,78,"D" +32,3,"B" +32,4,"C" +32,7,"E" +32,12,"E" +32,18,"A" +32,24,"B" +32,31,"C" +32,37,"B" +32,45,"C" +32,47,"D" +32,49,"C" +32,50,"D" +32,51,"A" +32,55,"D" +32,62,"B" +32,68,"D" +33,7,"B" +33,14,"E" +33,22,"D" +33,29,"D" +33,33,"A" +33,40,"C" +33,41,"C" +33,49,"E" +33,56,"B" +33,59,"C" +33,63,"A" +34,4,"C" +34,6,"D" +34,12,"B" +34,13,"B" +34,21,"B" +34,29,"B" +34,35,"D" +34,41,"C" +34,48,"B" +34,53,"A" +34,54,"B" +34,55,"B" +34,57,"C" +34,62,"D" +34,66,"E" +34,74,"D" +35,7,"D" +35,15,"C" +35,23,"B" +35,29,"B" +35,30,"B" +35,38,"D" +35,45,"A" +35,52,"C" +35,59,"C" +35,65,"A" +36,8,"A" +36,12,"E" +36,13,"B" +36,20,"D" +36,25,"C" +36,29,"D" +36,37,"C" +36,45,"C" +36,51,"D" +36,53,"C" +36,58,"C" +36,59,"E" +36,65,"B" +37,6,"B" +37,8,"E" +37,12,"D" +37,13,"B" +37,14,"D" +37,19,"C" +37,24,"C" +37,28,"B" +37,33,"C" +37,36,"A" +37,37,"C" +38,4,"C" +38,6,"D" +38,14,"B" +38,18,"C" +38,19,"C" +38,26,"B" +38,29,"E" +38,30,"D" +38,34,"C" +38,41,"C" +38,49,"B" +38,53,"C" +38,60,"C" +38,62,"B" +38,65,"E" +38,73,"C" +38,81,"D" +39,1,"D" +39,4,"C" +39,6,"D" +39,9,"B" +39,14,"B" +39,20,"D" +39,26,"C" +39,29,"C" +39,32,"C" +39,33,"B" +39,39,"D" +39,41,"C" +39,42,"D" +39,50,"A" +39,58,"C" +39,62,"B" +39,69,"C" +40,3,"B" +40,10,"D" +40,16,"C" +40,23,"D" +40,30,"B" +40,37,"D" +40,44,"C" +40,52,"B" +40,59,"B" +40,60,"C" +40,61,"E" +40,64,"B" +41,3,"C" +41,11,"E" +41,14,"A" +41,22,"C" +41,24,"A" +41,32,"D" +41,34,"B" +41,35,"D" +41,39,"C" +41,44,"C" +41,51,"C" +41,52,"C" +41,59,"A" +41,61,"C" +41,64,"C" +41,65,"D" +41,68,"C" +42,3,"C" +42,8,"C" +42,12,"C" +42,15,"A" +42,18,"B" +42,26,"D" +42,29,"D" +42,33,"C" +42,39,"C" +42,47,"A" +42,48,"C" +42,52,"C" +42,55,"C" +42,58,"D" +42,61,"B" +42,62,"D" +42,68,"D" +42,70,"C" +42,77,"B" +42,84,"B" +43,6,"B" +43,10,"C" +43,18,"A" +43,19,"E" +43,23,"C" +43,30,"E" +43,31,"E" +43,34,"C" +43,40,"B" +43,46,"C" +43,53,"C" +43,60,"C" +43,67,"C" +43,74,"C" +43,82,"B" +43,86,"A" +44,3,"D" +44,11,"C" +44,19,"C" +44,25,"C" +44,33,"C" +44,40,"B" +44,43,"B" +44,47,"A" +44,54,"C" +44,62,"A" +44,68,"C" +44,76,"D" +44,83,"B" +44,88,"C" +44,89,"D" +44,90,"C" +45,8,"C" +45,16,"E" +45,21,"D" +45,23,"C" +45,27,"E" +45,34,"C" +45,40,"E" +45,42,"B" +45,44,"E" +45,50,"C" +45,52,"D" +45,59,"B" +45,62,"E" +46,4,"D" +46,7,"E" +46,10,"A" +46,14,"A" +46,15,"E" +46,16,"B" +46,23,"D" +46,27,"B" +46,29,"C" +46,32,"C" +46,37,"C" +46,38,"D" +46,39,"B" +46,45,"A" +46,47,"A" +46,53,"C" +46,54,"E" +46,61,"D" +46,67,"C" +46,69,"C" +47,5,"A" +47,10,"C" +47,15,"C" +47,21,"C" +47,26,"D" +47,32,"B" +47,37,"B" +47,42,"C" +47,48,"C" +47,55,"E" +47,58,"C" +47,66,"A" +47,73,"B" +47,80,"C" +47,88,"D" +47,96,"C" +47,103,"D" +47,107,"B" +47,113,"C" +48,3,"C" +48,9,"C" +48,16,"E" +48,21,"C" +48,27,"A" +48,28,"C" +48,34,"D" +48,39,"A" +48,42,"C" +48,46,"C" +48,48,"B" +48,53,"B" +48,55,"E" +48,61,"D" +48,68,"B" +48,74,"D" +48,75,"E" +48,79,"C" +48,84,"C" +48,89,"B" +48,95,"B" +48,97,"A" +48,104,"A" +48,107,"B" +48,112,"B" +49,7,"A" +49,14,"B" +49,21,"A" +49,22,"A" +49,30,"C" +49,35,"C" +49,42,"E" +49,44,"B" +49,52,"A" +49,59,"E" +49,63,"C" +49,68,"B" +49,71,"B" +49,73,"B" +49,76,"D" +49,78,"E" +49,86,"B" +49,91,"C" +50,7,"C" +50,9,"C" +50,13,"B" +50,21,"C" +50,23,"B" +50,24,"C" +50,30,"C" +50,38,"B" +50,46,"C" +50,50,"B" +50,51,"C" +50,57,"B" +50,61,"D" +50,62,"C" +50,63,"B" +50,69,"D" +51,4,"C" +51,5,"D" +51,10,"C" +51,18,"D" +51,24,"B" +51,31,"C" +51,39,"B" +51,47,"C" +51,51,"A" +51,54,"C" +51,61,"C" +51,68,"E" +51,72,"B" +51,79,"D" +51,86,"B" +51,94,"C" +51,95,"B" +51,96,"A" +52,5,"E" +52,9,"D" +52,12,"C" +52,13,"C" +52,18,"C" +52,26,"D" +52,29,"B" +52,36,"C" +52,40,"C" +52,45,"A" +52,51,"D" +52,56,"A" +52,64,"B" +52,68,"C" +52,75,"A" +52,78,"B" +52,79,"D" +52,85,"A" +53,7,"C" +53,11,"C" +53,17,"D" +53,21,"C" +53,24,"B" +53,31,"B" +53,38,"D" +53,39,"C" +53,40,"C" +53,47,"B" +53,48,"B" +53,50,"D" +53,52,"D" +53,56,"D" +53,58,"D" +54,2,"A" +54,3,"A" +54,6,"C" +54,7,"C" +54,14,"B" +54,22,"B" +54,29,"C" +54,33,"C" +54,40,"B" +54,41,"B" +54,45,"B" +54,47,"C" +54,49,"A" +54,50,"C" +54,56,"B" +54,64,"E" +54,68,"B" +54,73,"B" +54,79,"C" +54,80,"C" +54,82,"B" +54,90,"D" +55,7,"D" +55,15,"D" +55,20,"B" +55,26,"D" +55,28,"C" +55,35,"D" +55,40,"D" +55,44,"C" +55,46,"B" +55,49,"D" +55,56,"D" +55,64,"A" +55,69,"C" +55,70,"B" +55,72,"D" +55,75,"E" +56,8,"C" +56,9,"C" +56,17,"A" +56,25,"D" +56,28,"D" +56,32,"C" +56,38,"D" +56,42,"C" +56,48,"B" +56,54,"C" +56,60,"C" +56,63,"C" +56,68,"C" +56,69,"E" +56,72,"C" +56,80,"A" +56,84,"D" +56,88,"D" +56,89,"C" +56,95,"B" +56,103,"C" +56,111,"C" +56,112,"C" +56,113,"C" +57,3,"D" +57,11,"C" +57,14,"E" +57,22,"D" +57,29,"D" +57,30,"B" +57,33,"C" +57,39,"D" +57,43,"A" +57,48,"D" +57,55,"B" +57,61,"C" +57,62,"C" +57,68,"B" +57,69,"B" +57,73,"C" +57,77,"C" +57,80,"E" +58,7,"E" +58,8,"B" +58,16,"C" +58,18,"C" +58,26,"C" +58,28,"E" +58,30,"C" +58,33,"B" +58,34,"C" +58,36,"B" +58,41,"D" +58,49,"C" +59,1,"D" +59,9,"D" +59,12,"D" +59,18,"C" +59,22,"B" +59,24,"C" +59,25,"D" +59,27,"C" +59,34,"C" +59,39,"D" +59,41,"C" +59,42,"E" +59,49,"E" +59,54,"E" +59,61,"A" +59,64,"D" +59,66,"D" +59,71,"B" +59,75,"B" +59,76,"D" +59,78,"A" +60,7,"C" +60,13,"C" +60,16,"D" +60,18,"B" +60,26,"B" +60,31,"C" +60,33,"C" +60,35,"D" +60,40,"C" +61,4,"C" +61,8,"B" +61,14,"A" +61,17,"B" +61,18,"C" +61,24,"D" +61,32,"C" +61,34,"C" +61,36,"D" +61,43,"D" +61,50,"E" +61,56,"C" +61,57,"A" +62,6,"A" +62,12,"C" +62,15,"D" +62,21,"C" +62,27,"C" +62,35,"D" +62,38,"B" +62,39,"E" +62,47,"A" +62,49,"C" +62,50,"D" +62,57,"D" +62,58,"A" +63,4,"A" +63,10,"C" +63,14,"D" +63,18,"E" +63,26,"C" +63,32,"C" +63,37,"C" +63,39,"C" +63,47,"B" +63,49,"C" +63,57,"E" +63,64,"D" +63,70,"B" +63,77,"C" +63,82,"B" +64,2,"D" +64,6,"C" +64,11,"C" +64,15,"C" +64,21,"C" +64,23,"C" +64,27,"C" +64,30,"C" +64,36,"C" +64,41,"E" +64,43,"B" +64,45,"D" +64,50,"C" +64,57,"C" +64,62,"C" +64,70,"E" +64,77,"C" +64,85,"D" diff --git a/simple-data-generation.ipynb b/simple-data-generation.ipynb index 7cb5e68..4d208b4 100644 --- a/simple-data-generation.ipynb +++ b/simple-data-generation.ipynb @@ -2,23 +2,31 @@ "cells": [ { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ "# Simple Data Generation for Tableau" ] }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ - "We will store the data in a data frame to aid with exporting to a \\*.csv file" + "We will store the data in data frames to aid with exporting to \\*.csv files" ] }, { "cell_type": "code", - "execution_count": 1, + "execution_count": 3, "metadata": { - "collapsed": true + "collapsed": false, + "deletable": true, + "editable": true }, "outputs": [], "source": [ @@ -27,8 +35,13 @@ }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ + "## Time to First Event Data\n", + "\n", "We initialize the data frame with 3 columns and 1024 samples\n", "\n", "1. The first column contains the subject `subject`, randomly selected from 64 identifiers.\n", @@ -38,373 +51,464 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": 14, "metadata": { - "collapsed": false + "collapsed": false, + "deletable": true, + "editable": true }, "outputs": [ { "data": { "text/html": [ - "
subjecttimeevent
116A
2114A
3121D
4124A
5126B
6131C
7132B
8137C
9141D
10148C
11153B
12158C
13162B
14264C
15272C
16274B
17276D
18284A
19287D
20289B
21291D
22397E
233104C
243105D
253111D
263115A
273119E
283122B
293129C
303131C
" + "
subjecttimeevent
115C
217B
3115D
4120B
5124C
6132C
7140C
8145C
" ], "text/plain": [ - "1024×3 DataFrames.DataFrame\n", - "│ Row │ subject │ time │ event │\n", - "├──────┼─────────┼──────┼───────┤\n", - "│ 1 │ 1 │ 6 │ 'A' │\n", - "│ 2 │ 1 │ 14 │ 'A' │\n", - "│ 3 │ 1 │ 21 │ 'D' │\n", - "│ 4 │ 1 │ 24 │ 'A' │\n", - "│ 5 │ 1 │ 26 │ 'B' │\n", - "│ 6 │ 1 │ 31 │ 'C' │\n", - "│ 7 │ 1 │ 32 │ 'B' │\n", - "│ 8 │ 1 │ 37 │ 'C' │\n", - "│ 9 │ 1 │ 41 │ 'D' │\n", - "│ 10 │ 1 │ 48 │ 'C' │\n", - "│ 11 │ 1 │ 53 │ 'B' │\n", - "⋮\n", - "│ 1013 │ 64 │ 4561 │ 'C' │\n", - "│ 1014 │ 64 │ 4569 │ 'C' │\n", - "│ 1015 │ 64 │ 4570 │ 'E' │\n", - "│ 1016 │ 64 │ 4574 │ 'C' │\n", - "│ 1017 │ 64 │ 4579 │ 'D' │\n", - "│ 1018 │ 64 │ 4584 │ 'C' │\n", - "│ 1019 │ 64 │ 4588 │ 'A' │\n", - "│ 1020 │ 64 │ 4590 │ 'C' │\n", - "│ 1021 │ 64 │ 4592 │ 'D' │\n", - "│ 1022 │ 64 │ 4599 │ 'B' │\n", - "│ 1023 │ 64 │ 4601 │ 'B' │\n", - "│ 1024 │ 64 │ 4603 │ 'C' │" + "8×3 DataFrames.DataFrame\n", + "│ Row │ subject │ time │ event │\n", + "├─────┼─────────┼──────┼───────┤\n", + "│ 1 │ 1 │ 5 │ 'C' │\n", + "│ 2 │ 1 │ 7 │ 'B' │\n", + "│ 3 │ 1 │ 15 │ 'D' │\n", + "│ 4 │ 1 │ 20 │ 'B' │\n", + "│ 5 │ 1 │ 24 │ 'C' │\n", + "│ 6 │ 1 │ 32 │ 'C' │\n", + "│ 7 │ 1 │ 40 │ 'C' │\n", + "│ 8 │ 1 │ 45 │ 'C' │" ] }, - "execution_count": 25, + "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "exportdata = DataFrame(\n", + "longitudinal = DataFrame(\n", " subject = sort(rand(1:64, 1024)),\n", " time = cumsum(rand(1:8, 1024)),\n", " event = rand(['A', 'B', 'B', 'C', 'C', 'C', 'C', 'D', 'D', 'E'], 1024)\n", - ")" + ");\n", + "longitudinal[1:8, :]" ] }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ "Notice that the time is increasing across all subjects. Instead we want to be the cumulative sum strictly within each subject. To accomplish this we use logical indexing to pull the last time from each subject and subtract it from the times of the next subject." ] }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ "First find the logic index of the row before the subject changed." ] }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 15, "metadata": { - "collapsed": false + "collapsed": false, + "deletable": true, + "editable": true }, "outputs": [ { "data": { "text/html": [ - "
subjecttimeeventbefore
116Afalse
2114Afalse
3121Dfalse
4124Afalse
5126Bfalse
6131Cfalse
7132Bfalse
8137Cfalse
9141Dfalse
10148Cfalse
11153Bfalse
12158Cfalse
13162Btrue
14264Cfalse
15272Cfalse
16274Bfalse
17276Dfalse
18284Afalse
19287Dfalse
20289Bfalse
21291Dtrue
22397Efalse
233104Cfalse
243105Dfalse
253111Dfalse
263115Afalse
273119Efalse
283122Bfalse
293129Cfalse
303131Cfalse
" + "
subjecttimeeventbefore
115Cfalse
217Bfalse
3115Dfalse
4120Bfalse
5124Cfalse
6132Cfalse
7140Cfalse
8145Cfalse
" ], "text/plain": [ - "1024×4 DataFrames.DataFrame\n", - "│ Row │ subject │ time │ event │ before │\n", - "├──────┼─────────┼──────┼───────┼────────┤\n", - "│ 1 │ 1 │ 6 │ 'A' │ false │\n", - "│ 2 │ 1 │ 14 │ 'A' │ false │\n", - "│ 3 │ 1 │ 21 │ 'D' │ false │\n", - "│ 4 │ 1 │ 24 │ 'A' │ false │\n", - "│ 5 │ 1 │ 26 │ 'B' │ false │\n", - "│ 6 │ 1 │ 31 │ 'C' │ false │\n", - "│ 7 │ 1 │ 32 │ 'B' │ false │\n", - "│ 8 │ 1 │ 37 │ 'C' │ false │\n", - "│ 9 │ 1 │ 41 │ 'D' │ false │\n", - "│ 10 │ 1 │ 48 │ 'C' │ false │\n", - "│ 11 │ 1 │ 53 │ 'B' │ false │\n", - "⋮\n", - "│ 1013 │ 64 │ 4561 │ 'C' │ false │\n", - "│ 1014 │ 64 │ 4569 │ 'C' │ false │\n", - "│ 1015 │ 64 │ 4570 │ 'E' │ false │\n", - "│ 1016 │ 64 │ 4574 │ 'C' │ false │\n", - "│ 1017 │ 64 │ 4579 │ 'D' │ false │\n", - "│ 1018 │ 64 │ 4584 │ 'C' │ false │\n", - "│ 1019 │ 64 │ 4588 │ 'A' │ false │\n", - "│ 1020 │ 64 │ 4590 │ 'C' │ false │\n", - "│ 1021 │ 64 │ 4592 │ 'D' │ false │\n", - "│ 1022 │ 64 │ 4599 │ 'B' │ false │\n", - "│ 1023 │ 64 │ 4601 │ 'B' │ false │\n", - "│ 1024 │ 64 │ 4603 │ 'C' │ false │" + "8×4 DataFrames.DataFrame\n", + "│ Row │ subject │ time │ event │ before │\n", + "├─────┼─────────┼──────┼───────┼────────┤\n", + "│ 1 │ 1 │ 5 │ 'C' │ false │\n", + "│ 2 │ 1 │ 7 │ 'B' │ false │\n", + "│ 3 │ 1 │ 15 │ 'D' │ false │\n", + "│ 4 │ 1 │ 20 │ 'B' │ false │\n", + "│ 5 │ 1 │ 24 │ 'C' │ false │\n", + "│ 6 │ 1 │ 32 │ 'C' │ false │\n", + "│ 7 │ 1 │ 40 │ 'C' │ false │\n", + "│ 8 │ 1 │ 45 │ 'C' │ false │" ] }, - "execution_count": 26, + "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "exportdata[:before] = [exportdata[:subject][1:end-1] .!= exportdata[:subject][2:end];false];\n", - "exportdata" + "longitudinal[:before] = [longitudinal[:subject][1:end-1] .!= longitudinal[:subject][2:end];false];\n", + "longitudinal[1:8, :]" ] }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ "Similarly find the logical index of the row after the subject changed." ] }, { "cell_type": "code", - "execution_count": 27, + "execution_count": 16, "metadata": { - "collapsed": false + "collapsed": false, + "deletable": true, + "editable": true }, "outputs": [ { "data": { "text/html": [ - "
subjecttimeeventbeforeafter
116Afalsefalse
2114Afalsefalse
3121Dfalsefalse
4124Afalsefalse
5126Bfalsefalse
6131Cfalsefalse
7132Bfalsefalse
8137Cfalsefalse
9141Dfalsefalse
10148Cfalsefalse
11153Bfalsefalse
12158Cfalsefalse
13162Btruefalse
14264Cfalsetrue
15272Cfalsefalse
16274Bfalsefalse
17276Dfalsefalse
18284Afalsefalse
19287Dfalsefalse
20289Bfalsefalse
21291Dtruefalse
22397Efalsetrue
233104Cfalsefalse
243105Dfalsefalse
253111Dfalsefalse
263115Afalsefalse
273119Efalsefalse
283122Bfalsefalse
293129Cfalsefalse
303131Cfalsefalse
" + "
subjecttimeeventbeforeafter
115Cfalsefalse
217Bfalsefalse
3115Dfalsefalse
4120Bfalsefalse
5124Cfalsefalse
6132Cfalsefalse
7140Cfalsefalse
8145Cfalsefalse
" ], "text/plain": [ - "1024×5 DataFrames.DataFrame\n", - "│ Row │ subject │ time │ event │ before │ after │\n", - "├──────┼─────────┼──────┼───────┼────────┼───────┤\n", - "│ 1 │ 1 │ 6 │ 'A' │ false │ false │\n", - "│ 2 │ 1 │ 14 │ 'A' │ false │ false │\n", - "│ 3 │ 1 │ 21 │ 'D' │ false │ false │\n", - "│ 4 │ 1 │ 24 │ 'A' │ false │ false │\n", - "│ 5 │ 1 │ 26 │ 'B' │ false │ false │\n", - "│ 6 │ 1 │ 31 │ 'C' │ false │ false │\n", - "│ 7 │ 1 │ 32 │ 'B' │ false │ false │\n", - "│ 8 │ 1 │ 37 │ 'C' │ false │ false │\n", - "│ 9 │ 1 │ 41 │ 'D' │ false │ false │\n", - "│ 10 │ 1 │ 48 │ 'C' │ false │ false │\n", - "│ 11 │ 1 │ 53 │ 'B' │ false │ false │\n", - "⋮\n", - "│ 1013 │ 64 │ 4561 │ 'C' │ false │ false │\n", - "│ 1014 │ 64 │ 4569 │ 'C' │ false │ false │\n", - "│ 1015 │ 64 │ 4570 │ 'E' │ false │ false │\n", - "│ 1016 │ 64 │ 4574 │ 'C' │ false │ false │\n", - "│ 1017 │ 64 │ 4579 │ 'D' │ false │ false │\n", - "│ 1018 │ 64 │ 4584 │ 'C' │ false │ false │\n", - "│ 1019 │ 64 │ 4588 │ 'A' │ false │ false │\n", - "│ 1020 │ 64 │ 4590 │ 'C' │ false │ false │\n", - "│ 1021 │ 64 │ 4592 │ 'D' │ false │ false │\n", - "│ 1022 │ 64 │ 4599 │ 'B' │ false │ false │\n", - "│ 1023 │ 64 │ 4601 │ 'B' │ false │ false │\n", - "│ 1024 │ 64 │ 4603 │ 'C' │ false │ false │" + "8×5 DataFrames.DataFrame\n", + "│ Row │ subject │ time │ event │ before │ after │\n", + "├─────┼─────────┼──────┼───────┼────────┼───────┤\n", + "│ 1 │ 1 │ 5 │ 'C' │ false │ false │\n", + "│ 2 │ 1 │ 7 │ 'B' │ false │ false │\n", + "│ 3 │ 1 │ 15 │ 'D' │ false │ false │\n", + "│ 4 │ 1 │ 20 │ 'B' │ false │ false │\n", + "│ 5 │ 1 │ 24 │ 'C' │ false │ false │\n", + "│ 6 │ 1 │ 32 │ 'C' │ false │ false │\n", + "│ 7 │ 1 │ 40 │ 'C' │ false │ false │\n", + "│ 8 │ 1 │ 45 │ 'C' │ false │ false │" ] }, - "execution_count": 27, + "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "exportdata[:after] = [false; exportdata[:before][1:end-1]];\n", - "exportdata" + "longitudinal[:after] = [false; longitudinal[:before][1:end-1]];\n", + "longitudinal[1:8, :]" ] }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ "Next initialize a column to contain the shift of baseline time for each subject." ] }, { "cell_type": "code", - "execution_count": 28, + "execution_count": 17, "metadata": { - "collapsed": false + "collapsed": false, + "deletable": true, + "editable": true }, "outputs": [ { "data": { "text/html": [ - "
subjecttimeeventbeforeaftershift
116Afalsefalse0
2114Afalsefalse0
3121Dfalsefalse0
4124Afalsefalse0
5126Bfalsefalse0
6131Cfalsefalse0
7132Bfalsefalse0
8137Cfalsefalse0
9141Dfalsefalse0
10148Cfalsefalse0
11153Bfalsefalse0
12158Cfalsefalse0
13162Btruefalse0
14264Cfalsetrue0
15272Cfalsefalse0
16274Bfalsefalse0
17276Dfalsefalse0
18284Afalsefalse0
19287Dfalsefalse0
20289Bfalsefalse0
21291Dtruefalse0
22397Efalsetrue0
233104Cfalsefalse0
243105Dfalsefalse0
253111Dfalsefalse0
263115Afalsefalse0
273119Efalsefalse0
283122Bfalsefalse0
293129Cfalsefalse0
303131Cfalsefalse0
" + "
subjecttimeeventbeforeaftershift
115Cfalsefalse0
217Bfalsefalse0
3115Dfalsefalse0
4120Bfalsefalse0
5124Cfalsefalse0
6132Cfalsefalse0
7140Cfalsefalse0
8145Cfalsefalse0
" ], "text/plain": [ - "1024×6 DataFrames.DataFrame\n", - "│ Row │ subject │ time │ event │ before │ after │ shift │\n", - "├──────┼─────────┼──────┼───────┼────────┼───────┼───────┤\n", - "│ 1 │ 1 │ 6 │ 'A' │ false │ false │ 0 │\n", - "│ 2 │ 1 │ 14 │ 'A' │ false │ false │ 0 │\n", - "│ 3 │ 1 │ 21 │ 'D' │ false │ false │ 0 │\n", - "│ 4 │ 1 │ 24 │ 'A' │ false │ false │ 0 │\n", - "│ 5 │ 1 │ 26 │ 'B' │ false │ false │ 0 │\n", - "│ 6 │ 1 │ 31 │ 'C' │ false │ false │ 0 │\n", - "│ 7 │ 1 │ 32 │ 'B' │ false │ false │ 0 │\n", - "│ 8 │ 1 │ 37 │ 'C' │ false │ false │ 0 │\n", - "│ 9 │ 1 │ 41 │ 'D' │ false │ false │ 0 │\n", - "│ 10 │ 1 │ 48 │ 'C' │ false │ false │ 0 │\n", - "│ 11 │ 1 │ 53 │ 'B' │ false │ false │ 0 │\n", - "⋮\n", - "│ 1013 │ 64 │ 4561 │ 'C' │ false │ false │ 0 │\n", - "│ 1014 │ 64 │ 4569 │ 'C' │ false │ false │ 0 │\n", - "│ 1015 │ 64 │ 4570 │ 'E' │ false │ false │ 0 │\n", - "│ 1016 │ 64 │ 4574 │ 'C' │ false │ false │ 0 │\n", - "│ 1017 │ 64 │ 4579 │ 'D' │ false │ false │ 0 │\n", - "│ 1018 │ 64 │ 4584 │ 'C' │ false │ false │ 0 │\n", - "│ 1019 │ 64 │ 4588 │ 'A' │ false │ false │ 0 │\n", - "│ 1020 │ 64 │ 4590 │ 'C' │ false │ false │ 0 │\n", - "│ 1021 │ 64 │ 4592 │ 'D' │ false │ false │ 0 │\n", - "│ 1022 │ 64 │ 4599 │ 'B' │ false │ false │ 0 │\n", - "│ 1023 │ 64 │ 4601 │ 'B' │ false │ false │ 0 │\n", - "│ 1024 │ 64 │ 4603 │ 'C' │ false │ false │ 0 │" + "8×6 DataFrames.DataFrame\n", + "│ Row │ subject │ time │ event │ before │ after │ shift │\n", + "├─────┼─────────┼──────┼───────┼────────┼───────┼───────┤\n", + "│ 1 │ 1 │ 5 │ 'C' │ false │ false │ 0 │\n", + "│ 2 │ 1 │ 7 │ 'B' │ false │ false │ 0 │\n", + "│ 3 │ 1 │ 15 │ 'D' │ false │ false │ 0 │\n", + "│ 4 │ 1 │ 20 │ 'B' │ false │ false │ 0 │\n", + "│ 5 │ 1 │ 24 │ 'C' │ false │ false │ 0 │\n", + "│ 6 │ 1 │ 32 │ 'C' │ false │ false │ 0 │\n", + "│ 7 │ 1 │ 40 │ 'C' │ false │ false │ 0 │\n", + "│ 8 │ 1 │ 45 │ 'C' │ false │ false │ 0 │" ] }, - "execution_count": 28, + "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "exportdata[:shift] = zeros(Int, 1024);\n", - "exportdata" + "longitudinal[:shift] = zeros(Int, 1024);\n", + "longitudinal[1:8, :]" ] }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ "Into the baseline shift of each row following the change in subject place the difference between the successive start times of each subject." ] }, { "cell_type": "code", - "execution_count": 29, + "execution_count": 18, "metadata": { - "collapsed": false + "collapsed": false, + "deletable": true, + "editable": true }, "outputs": [ { "data": { "text/html": [ - "
subjecttimeeventbeforeaftershift
116Afalsefalse0
2114Afalsefalse0
3121Dfalsefalse0
4124Afalsefalse0
5126Bfalsefalse0
6131Cfalsefalse0
7132Bfalsefalse0
8137Cfalsefalse0
9141Dfalsefalse0
10148Cfalsefalse0
11153Bfalsefalse0
12158Cfalsefalse0
13162Btruefalse0
14264Cfalsetrue62
15272Cfalsefalse0
16274Bfalsefalse0
17276Dfalsefalse0
18284Afalsefalse0
19287Dfalsefalse0
20289Bfalsefalse0
21291Dtruefalse0
22397Efalsetrue29
233104Cfalsefalse0
243105Dfalsefalse0
253111Dfalsefalse0
263115Afalsefalse0
273119Efalsefalse0
283122Bfalsefalse0
293129Cfalsefalse0
303131Cfalsefalse0
" + "
subjecttimeeventbeforeaftershift
115Cfalsefalse0
217Bfalsefalse0
3115Dfalsefalse0
4120Bfalsefalse0
5124Cfalsefalse0
6132Cfalsefalse0
7140Cfalsefalse0
8145Cfalsefalse0
" ], "text/plain": [ - "1024×6 DataFrames.DataFrame\n", - "│ Row │ subject │ time │ event │ before │ after │ shift │\n", - "├──────┼─────────┼──────┼───────┼────────┼───────┼───────┤\n", - "│ 1 │ 1 │ 6 │ 'A' │ false │ false │ 0 │\n", - "│ 2 │ 1 │ 14 │ 'A' │ false │ false │ 0 │\n", - "│ 3 │ 1 │ 21 │ 'D' │ false │ false │ 0 │\n", - "│ 4 │ 1 │ 24 │ 'A' │ false │ false │ 0 │\n", - "│ 5 │ 1 │ 26 │ 'B' │ false │ false │ 0 │\n", - "│ 6 │ 1 │ 31 │ 'C' │ false │ false │ 0 │\n", - "│ 7 │ 1 │ 32 │ 'B' │ false │ false │ 0 │\n", - "│ 8 │ 1 │ 37 │ 'C' │ false │ false │ 0 │\n", - "│ 9 │ 1 │ 41 │ 'D' │ false │ false │ 0 │\n", - "│ 10 │ 1 │ 48 │ 'C' │ false │ false │ 0 │\n", - "│ 11 │ 1 │ 53 │ 'B' │ false │ false │ 0 │\n", - "⋮\n", - "│ 1013 │ 64 │ 4561 │ 'C' │ false │ false │ 0 │\n", - "│ 1014 │ 64 │ 4569 │ 'C' │ false │ false │ 0 │\n", - "│ 1015 │ 64 │ 4570 │ 'E' │ false │ false │ 0 │\n", - "│ 1016 │ 64 │ 4574 │ 'C' │ false │ false │ 0 │\n", - "│ 1017 │ 64 │ 4579 │ 'D' │ false │ false │ 0 │\n", - "│ 1018 │ 64 │ 4584 │ 'C' │ false │ false │ 0 │\n", - "│ 1019 │ 64 │ 4588 │ 'A' │ false │ false │ 0 │\n", - "│ 1020 │ 64 │ 4590 │ 'C' │ false │ false │ 0 │\n", - "│ 1021 │ 64 │ 4592 │ 'D' │ false │ false │ 0 │\n", - "│ 1022 │ 64 │ 4599 │ 'B' │ false │ false │ 0 │\n", - "│ 1023 │ 64 │ 4601 │ 'B' │ false │ false │ 0 │\n", - "│ 1024 │ 64 │ 4603 │ 'C' │ false │ false │ 0 │" + "8×6 DataFrames.DataFrame\n", + "│ Row │ subject │ time │ event │ before │ after │ shift │\n", + "├─────┼─────────┼──────┼───────┼────────┼───────┼───────┤\n", + "│ 1 │ 1 │ 5 │ 'C' │ false │ false │ 0 │\n", + "│ 2 │ 1 │ 7 │ 'B' │ false │ false │ 0 │\n", + "│ 3 │ 1 │ 15 │ 'D' │ false │ false │ 0 │\n", + "│ 4 │ 1 │ 20 │ 'B' │ false │ false │ 0 │\n", + "│ 5 │ 1 │ 24 │ 'C' │ false │ false │ 0 │\n", + "│ 6 │ 1 │ 32 │ 'C' │ false │ false │ 0 │\n", + "│ 7 │ 1 │ 40 │ 'C' │ false │ false │ 0 │\n", + "│ 8 │ 1 │ 45 │ 'C' │ false │ false │ 0 │" ] }, - "execution_count": 29, + "execution_count": 18, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "exportdata[:shift][exportdata[:after]] = \n", - " exportdata[:time][exportdata[:before]] - \n", - " [0; exportdata[:time][exportdata[:before]][1:end-1]];\n", - "exportdata" + "longitudinal[:shift][longitudinal[:after]] = \n", + " longitudinal[:time][longitudinal[:before]] - \n", + " [0; longitudinal[:time][longitudinal[:before]][1:end-1]];\n", + "longitudinal[1:8, :]" ] }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ "Fill in the zeroes by running a cummulative sum over the shift column." ] }, { "cell_type": "code", - "execution_count": 30, + "execution_count": 19, "metadata": { - "collapsed": false + "collapsed": false, + "deletable": true, + "editable": true }, "outputs": [ { "data": { "text/html": [ - "
subjecttimeeventbeforeaftershift
116Afalsefalse0
2114Afalsefalse0
3121Dfalsefalse0
4124Afalsefalse0
5126Bfalsefalse0
6131Cfalsefalse0
7132Bfalsefalse0
8137Cfalsefalse0
9141Dfalsefalse0
10148Cfalsefalse0
11153Bfalsefalse0
12158Cfalsefalse0
13162Btruefalse0
14264Cfalsetrue62
15272Cfalsefalse62
16274Bfalsefalse62
17276Dfalsefalse62
18284Afalsefalse62
19287Dfalsefalse62
20289Bfalsefalse62
21291Dtruefalse62
22397Efalsetrue91
233104Cfalsefalse91
243105Dfalsefalse91
253111Dfalsefalse91
263115Afalsefalse91
273119Efalsefalse91
283122Bfalsefalse91
293129Cfalsefalse91
303131Cfalsefalse91
" + "
subjecttimeeventbeforeaftershift
115Cfalsefalse0
217Bfalsefalse0
3115Dfalsefalse0
4120Bfalsefalse0
5124Cfalsefalse0
6132Cfalsefalse0
7140Cfalsefalse0
8145Cfalsefalse0
" ], "text/plain": [ - "1024×6 DataFrames.DataFrame\n", - "│ Row │ subject │ time │ event │ before │ after │ shift │\n", - "├──────┼─────────┼──────┼───────┼────────┼───────┼───────┤\n", - "│ 1 │ 1 │ 6 │ 'A' │ false │ false │ 0 │\n", - "│ 2 │ 1 │ 14 │ 'A' │ false │ false │ 0 │\n", - "│ 3 │ 1 │ 21 │ 'D' │ false │ false │ 0 │\n", - "│ 4 │ 1 │ 24 │ 'A' │ false │ false │ 0 │\n", - "│ 5 │ 1 │ 26 │ 'B' │ false │ false │ 0 │\n", - "│ 6 │ 1 │ 31 │ 'C' │ false │ false │ 0 │\n", - "│ 7 │ 1 │ 32 │ 'B' │ false │ false │ 0 │\n", - "│ 8 │ 1 │ 37 │ 'C' │ false │ false │ 0 │\n", - "│ 9 │ 1 │ 41 │ 'D' │ false │ false │ 0 │\n", - "│ 10 │ 1 │ 48 │ 'C' │ false │ false │ 0 │\n", - "│ 11 │ 1 │ 53 │ 'B' │ false │ false │ 0 │\n", - "⋮\n", - "│ 1013 │ 64 │ 4561 │ 'C' │ false │ false │ 4531 │\n", - "│ 1014 │ 64 │ 4569 │ 'C' │ false │ false │ 4531 │\n", - "│ 1015 │ 64 │ 4570 │ 'E' │ false │ false │ 4531 │\n", - "│ 1016 │ 64 │ 4574 │ 'C' │ false │ false │ 4531 │\n", - "│ 1017 │ 64 │ 4579 │ 'D' │ false │ false │ 4531 │\n", - "│ 1018 │ 64 │ 4584 │ 'C' │ false │ false │ 4531 │\n", - "│ 1019 │ 64 │ 4588 │ 'A' │ false │ false │ 4531 │\n", - "│ 1020 │ 64 │ 4590 │ 'C' │ false │ false │ 4531 │\n", - "│ 1021 │ 64 │ 4592 │ 'D' │ false │ false │ 4531 │\n", - "│ 1022 │ 64 │ 4599 │ 'B' │ false │ false │ 4531 │\n", - "│ 1023 │ 64 │ 4601 │ 'B' │ false │ false │ 4531 │\n", - "│ 1024 │ 64 │ 4603 │ 'C' │ false │ false │ 4531 │" + "8×6 DataFrames.DataFrame\n", + "│ Row │ subject │ time │ event │ before │ after │ shift │\n", + "├─────┼─────────┼──────┼───────┼────────┼───────┼───────┤\n", + "│ 1 │ 1 │ 5 │ 'C' │ false │ false │ 0 │\n", + "│ 2 │ 1 │ 7 │ 'B' │ false │ false │ 0 │\n", + "│ 3 │ 1 │ 15 │ 'D' │ false │ false │ 0 │\n", + "│ 4 │ 1 │ 20 │ 'B' │ false │ false │ 0 │\n", + "│ 5 │ 1 │ 24 │ 'C' │ false │ false │ 0 │\n", + "│ 6 │ 1 │ 32 │ 'C' │ false │ false │ 0 │\n", + "│ 7 │ 1 │ 40 │ 'C' │ false │ false │ 0 │\n", + "│ 8 │ 1 │ 45 │ 'C' │ false │ false │ 0 │" ] }, - "execution_count": 30, + "execution_count": 19, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "exportdata[:shift] = cumsum(exportdata[:shift]);\n", - "exportdata" + "longitudinal[:shift] = cumsum(longitudinal[:shift]);\n", + "longitudinal[1:8, :]" ] }, { "cell_type": "markdown", - "metadata": {}, + "metadata": { + "deletable": true, + "editable": true + }, "source": [ "Rectify the time for each subject by subtracting the baseline shift." ] }, { "cell_type": "code", - "execution_count": 31, + "execution_count": 20, + "metadata": { + "collapsed": false, + "deletable": true, + "editable": true + }, + "outputs": [ + { + "data": { + "text/html": [ + "
subjecttimeeventbeforeaftershift
115Cfalsefalse0
217Bfalsefalse0
3115Dfalsefalse0
4120Bfalsefalse0
5124Cfalsefalse0
6132Cfalsefalse0
7140Cfalsefalse0
8145Cfalsefalse0
" + ], + "text/plain": [ + "8×6 DataFrames.DataFrame\n", + "│ Row │ subject │ time │ event │ before │ after │ shift │\n", + "├─────┼─────────┼──────┼───────┼────────┼───────┼───────┤\n", + "│ 1 │ 1 │ 5 │ 'C' │ false │ false │ 0 │\n", + "│ 2 │ 1 │ 7 │ 'B' │ false │ false │ 0 │\n", + "│ 3 │ 1 │ 15 │ 'D' │ false │ false │ 0 │\n", + "│ 4 │ 1 │ 20 │ 'B' │ false │ false │ 0 │\n", + "│ 5 │ 1 │ 24 │ 'C' │ false │ false │ 0 │\n", + "│ 6 │ 1 │ 32 │ 'C' │ false │ false │ 0 │\n", + "│ 7 │ 1 │ 40 │ 'C' │ false │ false │ 0 │\n", + "│ 8 │ 1 │ 45 │ 'C' │ false │ false │ 0 │" + ] + }, + "execution_count": 20, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "longitudinal[:time] = longitudinal[:time] - longitudinal[:shift];\n", + "longitudinal[1:8, :]" + ] + }, + { + "cell_type": "markdown", + "metadata": { + "deletable": true, + "editable": true + }, + "source": [ + "We can now export the data, choosing the first 3 columns." + ] + }, + { + "cell_type": "code", + "execution_count": 21, + "metadata": { + "collapsed": false, + "deletable": true, + "editable": true + }, + "outputs": [], + "source": [ + "writetable(\"longitudinal.csv\", longitudinal[[:subject, :time, :event]])" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Contingency Data\n", + "\n", + "The goal is to generate four columns of dimensional classifiers, as two pairs of dependent dimensions. We start by initializing the two distributions to sample:" + ] + }, + { + "cell_type": "code", + "execution_count": 52, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "distributionone = [\n", + " [\"A\" -1],\n", + " [\"A\" -1],\n", + " [\"A\" -1],\n", + " [\"A\" 0],\n", + " [\"A\" 0],\n", + " [\"A\" 1],\n", + " [\"B\" -1],\n", + " [\"B\" 0],\n", + " [\"B\" 0],\n", + " [\"B\" 1],\n", + " [\"B\" 1],\n", + " [\"B\" 1]\n", + "];\n", + "distributiontwo = [\n", + " [\"C\" -2],\n", + " [\"C\" -1],\n", + " [\"C\" -1],\n", + " [\"C\" 0],\n", + " [\"C\" 0],\n", + " [\"C\" 0],\n", + " [\"C\" 1],\n", + " [\"C\" 1],\n", + " [\"C\" 1],\n", + " [\"C\" 1],\n", + " [\"C\" 2],\n", + " [\"C\" 2],\n", + " [\"C\" 2],\n", + " [\"C\" 2],\n", + " [\"C\" 2],\n", + " [\"D\" -2],\n", + " [\"D\" -1],\n", + " [\"D\" -1],\n", + " [\"D\" 0],\n", + " [\"D\" 0],\n", + " [\"D\" 0],\n", + " [\"D\" 1],\n", + " [\"D\" 1],\n", + " [\"D\" 2],\n", + " [\"E\" -2],\n", + " [\"E\" -2],\n", + " [\"E\" -2],\n", + " [\"E\" -1],\n", + " [\"E\" -1],\n", + " [\"E\" 0],\n", + " [\"E\" 1],\n", + " [\"E\" 1],\n", + " [\"E\" 2],\n", + " [\"E\" 2],\n", + " [\"E\" 2],\n", + " [\"F\" 2],\n", + " [\"F\" 1],\n", + " [\"F\" 1],\n", + " [\"F\" 0],\n", + " [\"F\" 0],\n", + " [\"F\" 0],\n", + " [\"F\" -1],\n", + " [\"F\" -1],\n", + " [\"F\" -1],\n", + " [\"F\" -1],\n", + " [\"F\" -2],\n", + " [\"F\" -2],\n", + " [\"F\" -2],\n", + " [\"F\" -2],\n", + " [\"F\" -2]\n", + "];" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "We then generate another 1024 samples to populate a data frame." + ] + }, + { + "cell_type": "code", + "execution_count": 53, "metadata": { "collapsed": false }, @@ -412,64 +516,109 @@ { "data": { "text/html": [ - "
subjecttimeeventbeforeaftershift
116Afalsefalse0
2114Afalsefalse0
3121Dfalsefalse0
4124Afalsefalse0
5126Bfalsefalse0
6131Cfalsefalse0
7132Bfalsefalse0
8137Cfalsefalse0
9141Dfalsefalse0
10148Cfalsefalse0
11153Bfalsefalse0
12158Cfalsefalse0
13162Btruefalse0
1422Cfalsetrue62
15210Cfalsefalse62
16212Bfalsefalse62
17214Dfalsefalse62
18222Afalsefalse62
19225Dfalsefalse62
20227Bfalsefalse62
21229Dtruefalse62
2236Efalsetrue91
23313Cfalsefalse91
24314Dfalsefalse91
25320Dfalsefalse91
26324Afalsefalse91
27328Efalsefalse91
28331Bfalsefalse91
29338Cfalsefalse91
30340Cfalsefalse91
" + "
samplepreonetwoprethreefour
11Any[\"B\" 0]Any[\"C\" 0]
22Any[\"B\" 0]Any[\"C\" 0]
33Any[\"B\" 1]Any[\"C\" -1]
44Any[\"B\" -1]Any[\"D\" -1]
55Any[\"B\" 1]Any[\"F\" 0]
66Any[\"B\" 0]Any[\"F\" -2]
77Any[\"B\" 1]Any[\"C\" -1]
88Any[\"B\" -1]Any[\"C\" 1]
" ], "text/plain": [ - "1024×6 DataFrames.DataFrame\n", - "│ Row │ subject │ time │ event │ before │ after │ shift │\n", - "├──────┼─────────┼──────┼───────┼────────┼───────┼───────┤\n", - "│ 1 │ 1 │ 6 │ 'A' │ false │ false │ 0 │\n", - "│ 2 │ 1 │ 14 │ 'A' │ false │ false │ 0 │\n", - "│ 3 │ 1 │ 21 │ 'D' │ false │ false │ 0 │\n", - "│ 4 │ 1 │ 24 │ 'A' │ false │ false │ 0 │\n", - "│ 5 │ 1 │ 26 │ 'B' │ false │ false │ 0 │\n", - "│ 6 │ 1 │ 31 │ 'C' │ false │ false │ 0 │\n", - "│ 7 │ 1 │ 32 │ 'B' │ false │ false │ 0 │\n", - "│ 8 │ 1 │ 37 │ 'C' │ false │ false │ 0 │\n", - "│ 9 │ 1 │ 41 │ 'D' │ false │ false │ 0 │\n", - "│ 10 │ 1 │ 48 │ 'C' │ false │ false │ 0 │\n", - "│ 11 │ 1 │ 53 │ 'B' │ false │ false │ 0 │\n", - "⋮\n", - "│ 1013 │ 64 │ 30 │ 'C' │ false │ false │ 4531 │\n", - "│ 1014 │ 64 │ 38 │ 'C' │ false │ false │ 4531 │\n", - "│ 1015 │ 64 │ 39 │ 'E' │ false │ false │ 4531 │\n", - "│ 1016 │ 64 │ 43 │ 'C' │ false │ false │ 4531 │\n", - "│ 1017 │ 64 │ 48 │ 'D' │ false │ false │ 4531 │\n", - "│ 1018 │ 64 │ 53 │ 'C' │ false │ false │ 4531 │\n", - "│ 1019 │ 64 │ 57 │ 'A' │ false │ false │ 4531 │\n", - "│ 1020 │ 64 │ 59 │ 'C' │ false │ false │ 4531 │\n", - "│ 1021 │ 64 │ 61 │ 'D' │ false │ false │ 4531 │\n", - "│ 1022 │ 64 │ 68 │ 'B' │ false │ false │ 4531 │\n", - "│ 1023 │ 64 │ 70 │ 'B' │ false │ false │ 4531 │\n", - "│ 1024 │ 64 │ 72 │ 'C' │ false │ false │ 4531 │" + "8×3 DataFrames.DataFrame\n", + "│ Row │ sample │ preonetwo │ prethreefour │\n", + "├─────┼────────┼─────────────┼──────────────┤\n", + "│ 1 │ 1 │ Any[\"B\" 0] │ Any[\"C\" 0] │\n", + "│ 2 │ 2 │ Any[\"B\" 0] │ Any[\"C\" 0] │\n", + "│ 3 │ 3 │ Any[\"B\" 1] │ Any[\"C\" -1] │\n", + "│ 4 │ 4 │ Any[\"B\" -1] │ Any[\"D\" -1] │\n", + "│ 5 │ 5 │ Any[\"B\" 1] │ Any[\"F\" 0] │\n", + "│ 6 │ 6 │ Any[\"B\" 0] │ Any[\"F\" -2] │\n", + "│ 7 │ 7 │ Any[\"B\" 1] │ Any[\"C\" -1] │\n", + "│ 8 │ 8 │ Any[\"B\" -1] │ Any[\"C\" 1] │" ] }, - "execution_count": 31, + "execution_count": 53, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "exportdata[:time] = exportdata[:time] - exportdata[:shift];\n", - "exportdata" + "contingency = DataFrame(\n", + " sample = 1:1024,\n", + " preonetwo = rand(distributionone, 1024),\n", + " prethreefour = rand(distributiontwo, 1024)\n", + ");\n", + "contingency[1:8, :]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ - "We can now export the data, choosing the first 3 columns." + "Next we extract the individual column elements into dimension columns." ] }, { "cell_type": "code", - "execution_count": 34, + "execution_count": 54, "metadata": { "collapsed": false }, + "outputs": [ + { + "data": { + "text/html": [ + "
samplepreonetwoprethreefourdimensiononedimensiontwodimensionthreedimensionfour
11Any[\"B\" 0]Any[\"C\" 0]B0C0
22Any[\"B\" 0]Any[\"C\" 0]B0C0
33Any[\"B\" 1]Any[\"C\" -1]B1C-1
44Any[\"B\" -1]Any[\"D\" -1]B-1D-1
55Any[\"B\" 1]Any[\"F\" 0]B1F0
66Any[\"B\" 0]Any[\"F\" -2]B0F-2
77Any[\"B\" 1]Any[\"C\" -1]B1C-1
88Any[\"B\" -1]Any[\"C\" 1]B-1C1
" + ], + "text/plain": [ + "8×7 DataFrames.DataFrame\n", + "│ Row │ sample │ preonetwo │ prethreefour │ dimensionone │ dimensiontwo │\n", + "├─────┼────────┼─────────────┼──────────────┼──────────────┼──────────────┤\n", + "│ 1 │ 1 │ Any[\"B\" 0] │ Any[\"C\" 0] │ \"B\" │ 0 │\n", + "│ 2 │ 2 │ Any[\"B\" 0] │ Any[\"C\" 0] │ \"B\" │ 0 │\n", + "│ 3 │ 3 │ Any[\"B\" 1] │ Any[\"C\" -1] │ \"B\" │ 1 │\n", + "│ 4 │ 4 │ Any[\"B\" -1] │ Any[\"D\" -1] │ \"B\" │ -1 │\n", + "│ 5 │ 5 │ Any[\"B\" 1] │ Any[\"F\" 0] │ \"B\" │ 1 │\n", + "│ 6 │ 6 │ Any[\"B\" 0] │ Any[\"F\" -2] │ \"B\" │ 0 │\n", + "│ 7 │ 7 │ Any[\"B\" 1] │ Any[\"C\" -1] │ \"B\" │ 1 │\n", + "│ 8 │ 8 │ Any[\"B\" -1] │ Any[\"C\" 1] │ \"B\" │ -1 │\n", + "\n", + "│ Row │ dimensionthree │ dimensionfour │\n", + "├─────┼────────────────┼───────────────┤\n", + "│ 1 │ \"C\" │ 0 │\n", + "│ 2 │ \"C\" │ 0 │\n", + "│ 3 │ \"C\" │ -1 │\n", + "│ 4 │ \"D\" │ -1 │\n", + "│ 5 │ \"F\" │ 0 │\n", + "│ 6 │ \"F\" │ -2 │\n", + "│ 7 │ \"C\" │ -1 │\n", + "│ 8 │ \"C\" │ 1 │" + ] + }, + "execution_count": 54, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "contingency[:dimensionone] = vcat(contingency[:, :preonetwo]...)[:, 1];\n", + "contingency[:dimensiontwo] = vcat(contingency[:, :preonetwo]...)[:, 2];\n", + "contingency[:dimensionthree] = vcat(contingency[:, :prethreefour]...)[:, 1];\n", + "contingency[:dimensionfour] = vcat(contingency[:, :prethreefour]...)[:, 2];\n", + "contingency[1:8, :]" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Finally we export our contingency table data." + ] + }, + { + "cell_type": "code", + "execution_count": 55, + "metadata": { + "collapsed": true + }, "outputs": [], "source": [ - "writetable(\"export.csv\", exportdata[[:subject, :time, :event]])" + "writetable(\"contingency.csv\", contingency[[:sample, :dimensionone, :dimensiontwo, :dimensionthree, :dimensionfour]])" ] } ],