2b:null 2c:["$","title",null,{"children":"Working with the Washington Post ARCOS database API | Aaron Williams | Observable"}] 31:T5d0,chart = svg``32:[] 2f:{"id":"22e3089a867e9eb2","slug":"arcos-api-visualization","title":"Working with the Washington Post ARCOS database API","thumbnail":"2fd9df8b9a6c14584bc6c671239f1ece0d065b44d4d4199e584738a0c9e363da","access_level":"public","create_time":"$D2019-08-15T20:34:34.530Z","update_time":"$D2020-08-10T16:43:06.886Z","view_time":null,"first_public_time":"$D2019-08-20T13:42:54.912Z","fork":false,"stars":8,"starred":false,"owner":{"type":"team","id":"eb234b21118ba2f1","login":"aboutaaron","name":"Aaron Williams","avatar_url":"https://avatars-next.observableusercontent.com/d3a33f280ffec5c8283a9aae244b575dced553a6f821571324a77d92f0ac5839","tier":"starter_2024"},"version":410,"publish_version":410,"role":null,"authored":false,"authors":[{"type":"individual","id":"c8854dae9d7f94e1","login":"aboutaaron","name":"Aaron Williams","avatar_url":"https://avatars-next.observableusercontent.com/d3a33f280ffec5c8283a9aae244b575dced553a6f821571324a77d92f0ac5839","tier":null}],"members":[{"type":"individual","id":"c8854dae9d7f94e1","login":"aboutaaron","name":"Aaron Williams","avatar_url":"https://avatars-next.observableusercontent.com/d3a33f280ffec5c8283a9aae244b575dced553a6f821571324a77d92f0ac5839","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":3,"collection":null,"type":"notebook","body":{"title":"Working with the Washington Post ARCOS database API","cells":[{"id":0,"value":"md`# Working with the Washington Post ARCOS database API\n\nLast week, my team at the Washington Post published an [API](https://arcos-api.ext.nile.works/__swagger__/) for the DEA pain pill database we published earlier this month. We also published a [wrapper written in R](https://github.com/wpinvestigative/arcos). The API exposes the raw data as well as bunch of helpful summary datasets like pharmacy locations and annual pill distributions. \n\nThis map shows cumulative number of prescription opioid pills distributed to pharmacies in Los Angeles from 2006 to 2012. Tileset data is provided by Mapbox.\n\nSee a more robust use of this data [here](https://www.washingtonpost.com/graphics/2019/investigations/pharmacies-pain-pill-map/).`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":81,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":161,"value":"path = d3.geoPath(projection)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":118,"value":"get_tiles = t => {\n const access_token = 'pk.eyJ1IjoiYWJvdXRhYXJvbiIsImEiOiJsRTRpMGJnIn0.aprlJ6wE1JQqBx4EH1lkMQ'\n \n return `https://a.tiles.mapbox.com/v4/mapbox.mapbox-streets-v8/${t[2]}/${t[0]}/${t[1]}.mvt?access_token=${access_token}`\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":116,"value":"tiles = Promise.all(tile().map(async t => {\n let data = new vt.VectorTile(new Protobuf(await d3.buffer(get_tiles(t)))).layers\n \n let keys = Object.keys(data)\n t.data = {}\n keys.map(key => {\n let collection = {\n type: \"FeatureCollection\",\n features: []\n }\n\n for (let i = 0; i < data[key].length; i++) {\n let feature = data[key].feature(i).toGeoJSON(t[0], t[1], t[2])\n if (key === 'water') {\n feature.properties['class'] = 'water'\n }\n collection.features.push(feature)\n }\n \n t.data[key] = collection\n })\n \n return t;\n}))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":250,"value":"color = scaleCluster()\n .domain(data.map(d => d.properties.cumulative_dosage_unit))\n .range(d3.schemeBuPu[9])","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":295,"value":"size = scaleCluster()\n .domain(data.map(d => d.properties.cumulative_dosage_unit))\n .range([1, 2, 3, 4, 5, 6, 7, 8, 9, 10])","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":383,"value":"data = {\n const locations = await d3.json(routes.loc);\n const information = await d3.json(routes.data);\n const grouped = d3.groups(information, key => key.BUYER_DEA_NO);\n\n return grouped\n .filter(([key, _]) => locations.find(d => d.BUYER_DEA_NO === key))\n .map(arr => {\n const [BUYER_DEA_NO, values] = arr;\n const obj = locations.find(d => d.BUYER_DEA_NO === BUYER_DEA_NO);\n const cumulative_dosage_unit = values\n .map(d => d.DOSAGE_UNIT)\n .reduce((a, b) => a + b, 0);\n\n return turf.feature(toPoint(obj), {\n BUYER_DEA_NO,\n cumulative_dosage_unit\n });\n });\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":41,"value":"routes = ({\n data:\n 'https://arcos-api.ext.nile.works/v1/combined_buyer_annual?state=CA&county=Los%20Angeles&key=WaPo',\n loc:\n 'https://arcos-api.ext.nile.works/v1/pharmacy_latlon?state=CA&county=Los%20Angeles&key=WaPo'\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":95,"value":"tile = d3.tile()\n .size([width, height])\n .scale(projection.scale() * 2 * Math.PI)\n .translate(projection([0,0]))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":206,"value":"toPoint = d => ({\n \"type\": \"Point\",\n \"coordinates\": [d.lon, d.lat]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":105,"value":"projection = d3.geoMercator()\n .center([-118.2436, 34.0522])\n .scale(400000 / (2 * Math.PI))\n .translate([width / 2, height / 2])","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":97,"value":"height = 600","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":89,"value":"Protobuf = require('https://bundle.run/pbf@3.1.0')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":92,"value":"vt = require('https://bundle.run/@mapbox/vector-tile@1.3.1')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":35,"value":"d3 = require(\n \"d3-fetch\", \n \"d3-array\", \n \"d3-geo\", \n \"d3-geo-projection\", \n \"d3-selection\", \n \"d3-tile@1\",\n \"d3-scale-chromatic\"\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":268,"value":"scaleCluster = require('https://bundle.run/d3-scale-cluster@2.0.1')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":201,"value":"turf = require('@turf/turf')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":166,"value":"import {filter, isLayer} from '@aboutaaron/mapbox-map-maker'","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"}],"files":"$W32","stdlib":"1","titling":"auto"}} 30:"$2f:body"