Skip to content

Instantly share code, notes, and snippets.

View DSLituiev's full-sized avatar

Dmytro Lituiev DSLituiev

  • UCSF
  • San Francisco, CA, USA
View GitHub Profile
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@DSLituiev
DSLituiev / find_first_nonnegative.py
Last active September 11, 2019 03:22
find index of first non-negative value in a sorted list
def split(x):
"""return the index of the first non-negative value in a sorted list:
legend for comments:
"-" stands for negative
"+" stands for non-negative
"""
# check for corner cases:
if len(x)==0:
return 0
elif len(x)==1:
@DSLituiev
DSLituiev / install_gdb_with_conda_on_mac.sh
Created April 12, 2018 14:59
installing gdb on MacOS with conda python support
# install gettext; after installation it outputs LDFLAGS and CPPFLAGS values.
brew reinstall gettext
# export gettext and conda-specific flags
export LDFLAGS="-Wl,-rpath,/Applications/anaconda3/lib/ -L/Applications/anaconda3/lib/ -L/usr/local/opt/gettext/lib"
export CPPFLAGS="-I/usr/local/opt/gettext/include"
# extract gdb
tar xvf gdb-8.1.tar.gz
cd gdb-8.1
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Wed Oct 25 13:44:51 2017
@author: dlituiev
"""
import os
import numpy as np
from uuid import uuid1
@DSLituiev
DSLituiev / segm_metric_testcase.py
Created October 31, 2017 02:44
segmentation metric with tensorflow
"""
Created on Wed Oct 25 13:44:51 2017
@author: dlituiev
"""
import numpy as np
import keras
from keras import backend as K
from keras.layers import (InputLayer, Conv2D, Dense, Activation,
AveragePooling2D, GlobalAveragePooling2D)
@DSLituiev
DSLituiev / keras_constant_init_bug.py
Created October 9, 2017 23:43
Initialization with constant results in error when chaining layers
import sys
import numpy as np
import keras
from keras.models import load_model, Model
from keras.layers import Conv2D, MaxPool2D, GlobalAveragePooling2D, Dense, Dropout
from keras import backend as K
from keras.applications.inception_v3 import InceptionV3
def get_model(
dropout=0.5,
import pandas as pd
all_cols = ['0_AnatomicRegionSequence_CodeMeaning',
'0_AnatomicRegionSequence_CodeValue',
'0_AnatomicRegionSequence_CodingSchemeDesignator',
'0_ContributingEquipmentSequence_ContributionDateTime',
'0_ContributingEquipmentSequence_ContributionDescription',
'0_ContributingEquipmentSequence_InstitutionName',
'0_ContributingEquipmentSequence_Manufacturer',
'0_ContributingEquipmentSequence_StationName',
'0_ContributingEquipmentSequence__0_PurposeOfReferenceCodeSequence_CodeMeaning',
@DSLituiev
DSLituiev / gist:31063fc862b35a729142d3949e803c73
Last active October 5, 2016 21:26
ggplot error on colouring points
dfstr = '''{"FDR":{"0":0.0000001434,"1":0.0974543339,"2":0.5637882089,"3":0.7592318769,"4":0.82376489,"5":0.876232133,"6":0.891150062,"7":0.9234070548,"8":0.9421697501,"9":0.9421697501,"10":0.9534575004,"11":0.9686912495,"12":0.9792339965,"13":0.9800272087,"14":0.9814635178,"15":0.9814635178,"16":0.9822809925,"17":0.9832209543,"18":0.9832209543,"19":0.9832209543,"20":0.0079716299,"21":0.6863477759,"22":0.7846481266,"23":0.8596431336,"24":0.8940558933,"25":0.9071788194,"26":0.9269123036,"27":0.9323909829,"28":0.936505473,"29":0.9477474057,"30":0.9479119432,"31":0.9479119432,"32":0.9479119432,"33":0.9479119432,"34":0.9479179113,"35":0.9577011312,"36":0.9589141116,"37":0.9601000647,"38":0.9614078093,"39":0.9657393252,"40":0.0000000308,"41":0.1179302537,"42":0.6074557222,"43":0.7534315164,"44":0.8290719564,"45":0.8917851883,"46":0.9325332126,"47":0.9495263504,"48":0.9532127733,"49":0.9552555496,"50":0.9594234241,"51":0.9675706204,"52":0.9737916651,"53":0.9767958687,"54":0.9767958687,"55":0.9767958687,"56":0.97679
@DSLituiev
DSLituiev / fast_lmm_corr_error.py
Last active June 5, 2016 05:42
an error which indicates a modification of the X array by reference.
from __future__ import division, print_function
import numpy as np
import pandas as pd
from pysnptools.snpreader import SnpData
from pysnptools.kernelreader import KernelData
from fastlmm.inference.lmm import LMM
def snp_from_pandas(df, pos = None):
indiv_ids = zip(df.index.tolist(), ["."] * len(df.index))
snp_ids = [("_".join(["%s"% y for y in x]) if type(x) in (tuple, list) else "%s"% x) for x in df.columns.tolist()]
<Test-sra>
<Ngs>
<Latest>
<ncbi-vdb version=""/>
<ngs-sdk version=""/>
</Latest>
<LibManager.properties/>
<ncbi-vdb>
</ncbi-vdb>
<ngs-sdk>