|
3 | 3 | { |
4 | 4 | "cell_type": "code", |
5 | 5 | "execution_count": null, |
| 6 | + "id": "cfdd679e", |
6 | 7 | "metadata": {}, |
7 | 8 | "outputs": [], |
8 | 9 | "source": [ |
|
13 | 14 | { |
14 | 15 | "cell_type": "code", |
15 | 16 | "execution_count": null, |
| 17 | + "id": "37383869", |
16 | 18 | "metadata": {}, |
17 | 19 | "outputs": [], |
18 | 20 | "source": [ |
|
21 | 23 | }, |
22 | 24 | { |
23 | 25 | "cell_type": "markdown", |
| 26 | + "id": "f5a73851", |
24 | 27 | "metadata": {}, |
25 | 28 | "source": [ |
26 | 29 | "## Synced heatmaps" |
|
29 | 32 | { |
30 | 33 | "cell_type": "code", |
31 | 34 | "execution_count": null, |
| 35 | + "id": "ba0c9d3e", |
32 | 36 | "metadata": {}, |
33 | 37 | "outputs": [], |
34 | 38 | "source": [ |
|
65 | 69 | { |
66 | 70 | "cell_type": "code", |
67 | 71 | "execution_count": null, |
| 72 | + "id": "4b6c2629", |
68 | 73 | "metadata": {}, |
69 | 74 | "outputs": [], |
70 | 75 | "source": [ |
|
75 | 80 | { |
76 | 81 | "cell_type": "code", |
77 | 82 | "execution_count": null, |
| 83 | + "id": "1453f49f", |
78 | 84 | "metadata": {}, |
79 | 85 | "outputs": [], |
80 | 86 | "source": [ |
|
85 | 91 | { |
86 | 92 | "cell_type": "code", |
87 | 93 | "execution_count": null, |
| 94 | + "id": "4432d5d1", |
88 | 95 | "metadata": {}, |
89 | 96 | "outputs": [], |
90 | 97 | "source": [ |
|
103 | 110 | }, |
104 | 111 | { |
105 | 112 | "cell_type": "markdown", |
| 113 | + "id": "637bcad2", |
106 | 114 | "metadata": {}, |
107 | 115 | "source": [ |
108 | 116 | "## Value scale syncing" |
|
111 | 119 | { |
112 | 120 | "cell_type": "code", |
113 | 121 | "execution_count": null, |
| 122 | + "id": "e89a03d7", |
114 | 123 | "metadata": {}, |
115 | 124 | "outputs": [], |
116 | 125 | "source": [ |
|
122 | 131 | "# Creates a hg.RemoteTileset object\n", |
123 | 132 | "remote_tileset = hg.remote(\n", |
124 | 133 | " uid=\"CQMd6V_cRw6iCI_-Unl3PQ\",\n", |
125 | | - " server=\"https://higlass.io/api/v1/\", # optional, \"http://higlass.io/api/v1/\" default\n", |
| 134 | + " server=\"https://higlass.io/api/v1/\", # optional, \"http://higlass.io/api/v1/\"\n", |
126 | 135 | " name=\"Rao et al. (2014) GM12878 MboI (allreps) 1kb\",\n", |
127 | 136 | ")\n", |
128 | 137 | "\n", |
|
152 | 161 | }, |
153 | 162 | { |
154 | 163 | "cell_type": "markdown", |
| 164 | + "id": "ef880b6c", |
155 | 165 | "metadata": {}, |
156 | 166 | "source": [ |
157 | 167 | "## Remote heatmaps" |
|
160 | 170 | { |
161 | 171 | "cell_type": "code", |
162 | 172 | "execution_count": null, |
| 173 | + "id": "4579c459", |
163 | 174 | "metadata": {}, |
164 | 175 | "outputs": [], |
165 | 176 | "source": [ |
|
196 | 207 | }, |
197 | 208 | { |
198 | 209 | "cell_type": "markdown", |
| 210 | + "id": "74b22d9c", |
199 | 211 | "metadata": {}, |
200 | 212 | "source": [ |
201 | 213 | "## Extract track from another view config" |
|
204 | 216 | { |
205 | 217 | "cell_type": "code", |
206 | 218 | "execution_count": null, |
| 219 | + "id": "875260a4", |
207 | 220 | "metadata": {}, |
208 | 221 | "outputs": [], |
209 | 222 | "source": [ |
|
218 | 231 | { |
219 | 232 | "cell_type": "code", |
220 | 233 | "execution_count": null, |
| 234 | + "id": "f22d06ae", |
221 | 235 | "metadata": {}, |
222 | 236 | "outputs": [], |
223 | 237 | "source": [ |
|
230 | 244 | { |
231 | 245 | "cell_type": "code", |
232 | 246 | "execution_count": null, |
| 247 | + "id": "73c925a5", |
233 | 248 | "metadata": {}, |
234 | 249 | "outputs": [], |
235 | 250 | "source": [ |
|
242 | 257 | }, |
243 | 258 | { |
244 | 259 | "cell_type": "markdown", |
| 260 | + "id": "0a9ad2fc", |
245 | 261 | "metadata": {}, |
246 | 262 | "source": [ |
247 | 263 | "## Remote bigWig tiles\n", |
|
251 | 267 | { |
252 | 268 | "cell_type": "code", |
253 | 269 | "execution_count": null, |
| 270 | + "id": "3d664bb6", |
254 | 271 | "metadata": {}, |
255 | 272 | "outputs": [], |
256 | 273 | "source": [ |
|
260 | 277 | { |
261 | 278 | "cell_type": "code", |
262 | 279 | "execution_count": null, |
| 280 | + "id": "0e22f81d", |
263 | 281 | "metadata": {}, |
264 | 282 | "outputs": [], |
265 | 283 | "source": [ |
|
278 | 296 | { |
279 | 297 | "cell_type": "code", |
280 | 298 | "execution_count": null, |
| 299 | + "id": "4133c1c3", |
281 | 300 | "metadata": {}, |
282 | 301 | "outputs": [], |
283 | 302 | "source": [ |
284 | 303 | "import hg\n", |
285 | 304 | "\n", |
286 | 305 | "# get mounted path\n", |
287 | 306 | "bwpath = hg.fuse.path(\n", |
288 | | - " f\"http://hgdownload.cse.ucsc.edu/goldenpath/hg19/encodeDCC/\"\n", |
289 | | - " f\"wgEncodeSydhTfbs/wgEncodeSydhTfbsGm12878InputStdSig.bigWig\"\n", |
| 307 | + " \"http://hgdownload.cse.ucsc.edu/goldenpath/hg19/encodeDCC/\"\n", |
| 308 | + " \"wgEncodeSydhTfbs/wgEncodeSydhTfbsGm12878InputStdSig.bigWig\"\n", |
290 | 309 | ")\n", |
291 | 310 | "\n", |
292 | 311 | "track = hg.bigwig(bwpath).track(\"horizontal-bar\")\n", |
|
295 | 314 | }, |
296 | 315 | { |
297 | 316 | "cell_type": "markdown", |
298 | | - "metadata": { |
299 | | - "tags": [] |
300 | | - }, |
| 317 | + "id": "5844a303", |
| 318 | + "metadata": {}, |
301 | 319 | "source": [ |
302 | 320 | "## Local cooler files\n", |
303 | 321 | "\n", |
|
307 | 325 | { |
308 | 326 | "cell_type": "code", |
309 | 327 | "execution_count": null, |
| 328 | + "id": "0439622f", |
310 | 329 | "metadata": {}, |
311 | 330 | "outputs": [], |
312 | 331 | "source": [ |
313 | 332 | "tileset = hg.cooler(\"../test/data/Dixon2012-J1-NcoI-R1-filtered.100kb.multires.cool\")\n", |
314 | 333 | "hg.view(tileset.track(\"heatmap\"))" |
315 | 334 | ] |
316 | | - }, |
317 | | - { |
318 | | - "cell_type": "markdown", |
319 | | - "metadata": {}, |
320 | | - "source": [ |
321 | | - "## Local bigWig files (with chromsizes)\n", |
322 | | - "\n", |
323 | | - "### TODO: inline `chromsizes` tilesets are not yet implemented" |
324 | | - ] |
325 | | - }, |
326 | | - { |
327 | | - "cell_type": "code", |
328 | | - "execution_count": null, |
329 | | - "metadata": {}, |
330 | | - "outputs": [], |
331 | | - "source": [ |
332 | | - "chromsizes = [\n", |
333 | | - " (\"chr1\", 249250621),\n", |
334 | | - " (\"chr2\", 243199373),\n", |
335 | | - " (\"chr3\", 198022430),\n", |
336 | | - " (\"chr4\", 191154276),\n", |
337 | | - " (\"chr5\", 180915260),\n", |
338 | | - " (\"chr6\", 171115067),\n", |
339 | | - " (\"chr7\", 159138663),\n", |
340 | | - " (\"chr8\", 146364022),\n", |
341 | | - " (\"chr9\", 141213431),\n", |
342 | | - " (\"chr10\", 135534747),\n", |
343 | | - " (\"chr11\", 135006516),\n", |
344 | | - " (\"chr12\", 133851895),\n", |
345 | | - " (\"chr13\", 115169878),\n", |
346 | | - " (\"chr14\", 107349540),\n", |
347 | | - " (\"chr15\", 102531392),\n", |
348 | | - " (\"chr16\", 90354753),\n", |
349 | | - " (\"chr17\", 81195210),\n", |
350 | | - " (\"chr18\", 78077248),\n", |
351 | | - " (\"chr20\", 63025520),\n", |
352 | | - " (\"chr19\", 59128983),\n", |
353 | | - " (\"chr21\", 48129895),\n", |
354 | | - " (\"chr22\", 51304566),\n", |
355 | | - " (\"chrX\", 155270560),\n", |
356 | | - " (\"chrY\", 59373566),\n", |
357 | | - " (\"chrM\", 16571),\n", |
358 | | - "]\n", |
359 | | - "\n", |
360 | | - "bigwig_fp = (\n", |
361 | | - " \"../test/data/wgEncodeCaltechRnaSeqHuvecR1x75dTh1014IlnaPlusSignalRep2.bigWig\"\n", |
362 | | - ")\n", |
363 | | - "ts = hg.bigwig(bigwig_fp, chromsizes=chromsizes)\n", |
364 | | - "cs = hg.ChromSizes(chromsizes)\n", |
365 | | - "\n", |
366 | | - "view1 = View(\n", |
367 | | - " [\n", |
368 | | - " Track(\"top-axis\"),\n", |
369 | | - " Track(\"horizontal-bar\", tileset=ts),\n", |
370 | | - " Track(\"horizontal-chromosome-labels\", position=\"top\", tileset=cs),\n", |
371 | | - " ]\n", |
372 | | - ")\n", |
373 | | - "display, server, viewconf = higlass.display([view1])\n", |
374 | | - "\n", |
375 | | - "display" |
376 | | - ] |
377 | | - }, |
378 | | - { |
379 | | - "cell_type": "markdown", |
380 | | - "metadata": {}, |
381 | | - "source": [ |
382 | | - "## Local bedlike data\n", |
383 | | - "\n", |
384 | | - "### TODO: inline `bedlike` tilesets are not yet implemented" |
385 | | - ] |
386 | | - }, |
387 | | - { |
388 | | - "cell_type": "code", |
389 | | - "execution_count": null, |
390 | | - "metadata": {}, |
391 | | - "outputs": [], |
392 | | - "source": [ |
393 | | - "import hg\n", |
394 | | - "\n", |
395 | | - "bed = [\n", |
396 | | - " [\"chr1\", 1000, 2000, \"item #1\", \".\", \"+\"],\n", |
397 | | - " [\"chr2\", 3000, 3500, \"item #1\", \".\", \"-\"],\n", |
398 | | - "]\n", |
399 | | - "\n", |
400 | | - "chroms = [\n", |
401 | | - " [\"chr1\", 2100],\n", |
402 | | - " [\"chr2\", 4000],\n", |
403 | | - "]\n", |
404 | | - "\n", |
405 | | - "\n", |
406 | | - "data = bedtiles(bed, chroms)\n", |
407 | | - "track = Track(\n", |
408 | | - " track_type=\"bedlike\",\n", |
409 | | - " position=\"top\",\n", |
410 | | - " height=50,\n", |
411 | | - " data=data,\n", |
412 | | - " options={\"minusStrandColor\": \"red\"},\n", |
413 | | - ")\n", |
414 | | - "\n", |
415 | | - "\n", |
416 | | - "d, s, v = higlass.display([[track]])\n", |
417 | | - "d" |
418 | | - ] |
419 | | - }, |
420 | | - { |
421 | | - "cell_type": "markdown", |
422 | | - "metadata": {}, |
423 | | - "source": [ |
424 | | - "## Custom data" |
425 | | - ] |
426 | | - }, |
427 | | - { |
428 | | - "cell_type": "code", |
429 | | - "execution_count": null, |
430 | | - "metadata": {}, |
431 | | - "outputs": [], |
432 | | - "source": [ |
433 | | - "import numpy as np\n", |
434 | | - "\n", |
435 | | - "dim = 2000\n", |
436 | | - "I, J = np.indices((dim, dim))\n", |
437 | | - "data = -(J + 47) * np.sin(np.sqrt(np.abs(I / 2 + (J + 47)))) - I * np.sin(\n", |
438 | | - " np.sqrt(np.abs(I - (J + 47)))\n", |
439 | | - ")" |
440 | | - ] |
441 | | - }, |
442 | | - { |
443 | | - "cell_type": "code", |
444 | | - "execution_count": null, |
445 | | - "metadata": {}, |
446 | | - "outputs": [], |
447 | | - "source": [ |
448 | | - "import clodius.tiles.npmatrix\n", |
449 | | - "from hg.tilesets import LocalTileset\n", |
450 | | - "\n", |
451 | | - "ts = hg.server.add(\n", |
452 | | - " LocalTileset(\n", |
453 | | - " uid=\"my-custom-tileset\",\n", |
454 | | - " info=lambda: clodius.tiles.npmatrix.tileset_info(data),\n", |
455 | | - " tiles=lambda tids: clodius.tiles.npmatrix.tiles_wrapper(data, tids),\n", |
456 | | - " )\n", |
457 | | - ")\n", |
458 | | - "\n", |
459 | | - "hg.view(\n", |
460 | | - " (hg.track(\"top-axis\"), \"top\"),\n", |
461 | | - " (hg.track(\"left-axis\"), \"left\"),\n", |
462 | | - " (ts.track(\"heatmap\", height=250).opts(valueScaleMax=0.5), \"center\"),\n", |
463 | | - ")" |
464 | | - ] |
465 | | - }, |
466 | | - { |
467 | | - "cell_type": "code", |
468 | | - "execution_count": null, |
469 | | - "metadata": {}, |
470 | | - "outputs": [], |
471 | | - "source": [] |
472 | | - }, |
473 | | - { |
474 | | - "cell_type": "code", |
475 | | - "execution_count": null, |
476 | | - "metadata": {}, |
477 | | - "outputs": [], |
478 | | - "source": [] |
479 | 335 | } |
480 | 336 | ], |
481 | 337 | "metadata": { |
|
0 commit comments