Pathway analysis


In [35]:
import getpass
from transmart_api import TransmartApi

api = TransmartApi(
    host = 'http://localhost:8080',
    user = raw_input('Username:'),
    password = getpass.getpass('Password:'))

api.access()


Username:admin
Password:········
Out[35]:
'SUCCESS'

In [36]:
import pandas

pandas.set_option('max_colwidth', 1000)
pandas.set_option("display.max_rows",100)

In [37]:
from pandas.io.json import json_normalize

obs = api.get_observations(study = 'GSE8581')

obs_df1 = json_normalize(obs)
obs_df2 = obs_df1.pivot(index = 'subject.inTrialId', columns = 'label', values = 'value')
obs_df3 = obs_df2.convert_objects()
obs_df4 = obs_df3.rename(
    columns = lambda c: c.replace('\Public Studies\GSE8581\\', '')[:-1],
    inplace = False)
obs_df4


Out[37]:
label Biomarker Data\GPL570 Endpoints\Diagnosis Endpoints\FEV1 Endpoints\Forced Expiratory Volume Ratio Subjects\Age Subjects\Height (inch) Subjects\Lung Disease Subjects\Organism Subjects\Race Subjects\Sex
subject.inTrialId
GSE8581GSM210004 None non-small cell squamous cell carcinoma 2.54 58.00 63 72 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM210005 None non-small cell adenocarcinoma 1.69 83.66 84 60 control Homo sapiens Afro American female
GSE8581GSM210006 None non-small cell adenocarcinoma 1.41 51.00 65 66 chronic obstructive pulmonary disease Homo sapiens Afro American female
GSE8581GSM210007 E non-small cell adenocarcinoma 2.51 80.96 46 66 not specified Homo sapiens Caucasian male
GSE8581GSM210008 E non-small cell adenocarcinoma 1.64 57.00 53 65 chronic obstructive pulmonary disease Homo sapiens Caucasian female
GSE8581GSM210009 E non-small cell squamous cell carcinoma 2.72 74.00 53 64 control Homo sapiens Caucasian female
GSE8581GSM210010 E non-small cell adenocarcinoma 1.45 73.00 77 63 not specified Homo sapiens Caucasian female
GSE8581GSM210011 E non-small cell squamous cell carcinoma 1.87 56.00 56 72 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM210012 E non-small cell adenocarcinoma 2.76 70.58 61 69 not specified Homo sapiens Caucasian male
GSE8581GSM210014 E non-small cell adenocarcinoma 1.98 78.00 71 63 control Homo sapiens Caucasian female
GSE8581GSM210015 E non-small cell adenocarcinoma 2.59 74.00 68 65 control Homo sapiens Caucasian male
GSE8581GSM210071 E emphysema 1.16 43.00 68 67 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM210087 E non-small cell adenocarcinoma 2.59 74.00 68 65 control Homo sapiens Caucasian female
GSE8581GSM210090 E hematoma 4.08 68.00 50 72 not specified Homo sapiens Caucasian male
GSE8581GSM210188 E metastatic non-small cell adenocarcinoma 1.20 72.00 65 63 not specified Homo sapiens Caucasian female
GSE8581GSM210192 E non-small cell adenocarcinoma 2.31 75.00 50 65 control Homo sapiens Caucasian female
GSE8581GSM210193 E Unknown 1.63 80.78 67 66 not specified Homo sapiens Caucasian female
GSE8581GSM210194 E non-small cell squamous cell carcinoma 0.52 58.00 56 60 chronic obstructive pulmonary disease Homo sapiens Caucasian female
GSE8581GSM210196 E non-small cell adenocarcinoma 2.30 74.00 57 65 control Homo sapiens Caucasian female
GSE8581GSM210978 E non-small cell adenocarcinoma 1.53 74.00 73 66 not specified Homo sapiens Caucasian female
GSE8581GSM210979 E non-small cell adenocarcinoma 2.17 76.00 55 63 control Homo sapiens Caucasian female
GSE8581GSM210992 E non-small cell squamous cell carcinoma 0.90 55.90 70 71 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM210993 E giant bullae 0.40 45.00 61 73 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM210994 E carcinoid 1.69 76.00 64 63 not specified Homo sapiens Caucasian female
GSE8581GSM211007 E non-small cell squamous cell carcinoma 1.70 59.00 61 66 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM211008 E non-small cell adenocarcinoma 3.66 72.00 71 71 control Homo sapiens Caucasian male
GSE8581GSM211009 E non-small cell adenocarcinoma 1.64 41.00 75 71 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM211010 E non-small cell adenocarcinoma 1.44 86.00 61 60 not specified Homo sapiens Caucasian female
GSE8581GSM211865 E non-small cell squamous cell carcinoma 2.13 71.00 69 67 not specified Homo sapiens Caucasian male
GSE8581GSM211872 E non-small cell squamous cell carcinoma 2.61 70.00 74 70 not specified Homo sapiens Caucasian male
GSE8581GSM212067 E inflammation 4.04 79.00 55 69 control Homo sapiens Caucasian male
GSE8581GSM212068 E metastatic renal cell carcinoma 1.99 85.00 78 63 control Homo sapiens Caucasian female
GSE8581GSM212069 E no malignancy 1.60 77.00 81 66 not specified Homo sapiens Caucasian female
GSE8581GSM212070 E carcinoid 3.70 81.86 40 67 control Homo sapiens Caucasian male
GSE8581GSM212074 E non-small cell squamous cell carcinoma 2.06 54.00 64 65 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM212075 E non-small cell squamous cell carcinoma 1.30 78.00 79 63 chronic obstructive pulmonary disease Homo sapiens Afro American female
GSE8581GSM212787 E NSC-Mixed 3.66 78.00 78 70 control Homo sapiens Caucasian male
GSE8581GSM212788 E non-small cell squamous cell carcinoma 1.29 53.00 77 67 chronic obstructive pulmonary disease Homo sapiens Caucasian female
GSE8581GSM212789 E non-small cell squamous cell carcinoma 1.94 86.00 62 62 control Homo sapiens Caucasian female
GSE8581GSM212790 E non-small cell squamous cell carcinoma 2.83 75.00 54 75 control Homo sapiens Caucasian female
GSE8581GSM212809 E non-small cell adenocarcinoma 1.28 89.00 61 65 not specified Homo sapiens Caucasian female
GSE8581GSM212810 E non-small cell squamous cell carcinoma 0.88 52.00 52 66 chronic obstructive pulmonary disease Homo sapiens Caucasian female
GSE8581GSM212811 E non-small cell squamous cell carcinoma 1.37 82.00 77 58 control Homo sapiens Caucasian female
GSE8581GSM212848 E non-small cell squamous cell carcinoma 0.55 44.00 59 58 chronic obstructive pulmonary disease Homo sapiens Caucasian female
GSE8581GSM212849 E non-small cell adenocarcinoma 2.73 80.00 39 73 not specified Homo sapiens Caucasian male
GSE8581GSM212850 E carcinoid 2.32 76.00 71 69 not specified Homo sapiens Caucasian male
GSE8581GSM212852 E Giant Cell Tumor 0.99 63.00 72 60 chronic obstructive pulmonary disease Homo sapiens Caucasian female
GSE8581GSM212853 E lymphoma 2.42 80.00 71 66 control Homo sapiens Caucasian female
GSE8581GSM212854 E non-small cell adenocarcinoma 2.26 83.00 64 66 not specified Homo sapiens Caucasian male
GSE8581GSM212855 E non-small cell adenocarcinoma 2.43 69.00 73 68 not specified Homo sapiens Caucasian male
GSE8581GSM213017 E non-small cell adenocarcinoma 2.12 73.00 79 69 not specified Homo sapiens Caucasian male
GSE8581GSM213018 E non-small cell adenocarcinoma 2.34 81.00 82 68 not specified Homo sapiens Caucasian male
GSE8581GSM213019 E non-small cell adenocarcinoma 1.70 72.00 71 71 not specified Homo sapiens Caucasian female
GSE8581GSM213020 E non-small cell adenocarcinoma 0.67 42.00 75 63 chronic obstructive pulmonary disease Homo sapiens Caucasian male
GSE8581GSM213034 E non-small cell adenocarcinoma 2.56 61.00 70 69 not specified Homo sapiens Caucasian male
GSE8581GSM213035 E non-small cell squamous cell carcinoma 3.11 73.00 67 71 control Homo sapiens Caucasian male
GSE8581GSM213036 E non-small cell squamous cell carcinoma 2.68 71.00 59 69 control Homo sapiens Caucasian female
GSE8581GSM213037 E non-small cell squamous cell carcinoma 2.25 72.00 77 66 not specified Homo sapiens Caucasian male

In [38]:
control = obs_df4['Subjects\Lung Disease'] == 'control'
treatment = obs_df4['Subjects\Lung Disease'] == 'chronic obstructive pulmonary disease'

In [39]:
obs_df4[control].describe()


Out[39]:
Endpoints\FEV1 Endpoints\Forced Expiratory Volume Ratio Subjects\Age Subjects\Height (inch)
count 19.000000 19.000000 19.000000 19.000000
mean 2.618421 77.448421 64.105263 65.842105
std 0.737219 4.540892 11.503114 4.193423
min 1.370000 71.000000 40.000000 58.000000
25% 2.080000 74.000000 55.000000 63.000000
50% 2.590000 76.000000 67.000000 65.000000
75% 2.970000 80.930000 71.000000 69.000000
max 4.040000 86.000000 84.000000 75.000000

In [40]:
obs_df4[treatment].describe()


Out[40]:
Endpoints\FEV1 Endpoints\Forced Expiratory Volume Ratio Subjects\Age Subjects\Height (inch)
count 17.000000 17.000000 17.000000 17.000000
mean 1.265882 53.523529 65.058824 66.176471
std 0.594433 9.205130 8.554754 4.530939
min 0.400000 41.000000 52.000000 58.000000
25% 0.880000 45.000000 59.000000 63.000000
50% 1.290000 54.000000 64.000000 66.000000
75% 1.640000 58.000000 72.000000 71.000000
max 2.540000 78.000000 79.000000 73.000000

In [41]:
(hdHeader, hdRows) = api.get_hd_node_data(study = 'GSE8581',
                                          node_name = 'Lung',
                                          projection='zscore')

In [42]:
rows = [row.value[0].doubleValue._values for row in hdRows if row.bioMarker != 'null']
row_names = [row.bioMarker + "\\" + row.label for row in hdRows if row.bioMarker != 'null']
col_names = [assay.patientId for assay in hdHeader.assay]

In [43]:
from pandas import DataFrame

hd_df = DataFrame(rows, columns = col_names, index = row_names)
hd_df


Out[43]:
GSE8581GSM213034 GSE8581GSM212811 GSE8581GSM213036 GSE8581GSM212075 GSE8581GSM211008 GSE8581GSM210090 GSE8581GSM212855 GSE8581GSM212070 GSE8581GSM212810 GSE8581GSM210193 ... GSE8581GSM210979 GSE8581GSM213017 GSE8581GSM211007 GSE8581GSM210192 GSE8581GSM213020 GSE8581GSM210994 GSE8581GSM210010 GSE8581GSM210011 GSE8581GSM210007 GSE8581GSM213019
KIAA1377\235956_at -2.374995e-01 7.941626e-01 0.222018 1.117909 6.363410e-01 -4.181891e-01 -1.823707e+00 1.354521 -2.500000 -2.198561e-02 ... 1.234281e+00 -1.578853e-01 0.642150 -0.824503 -1.418721e-01 7.586432e-01 -1.102181e+00 -1.168079e+00 -6.064866e-01 2.864873e-01
ZNF358\226260_x_at 4.708822e-01 -2.500000e+00 0.551998 0.298824 -7.286189e-01 4.696274e-01 -2.807687e-01 0.299635 -2.500000 -3.315997e-01 ... 3.100759e-01 1.628198e-01 0.695030 0.129990 1.036752e-01 -1.037046e+00 -3.647481e-01 3.708714e-01 -3.854062e-01 1.054395e+00
GPR135\214503_x_at 6.308876e-02 -1.971560e+00 0.429403 -0.606546 -7.913499e-17 -2.291924e-01 1.167696e+00 1.319348 -2.267266 9.553151e-01 ... 5.955446e-01 6.090966e-01 1.238395 -0.078148 9.911758e-01 -7.427790e-01 -7.761581e-02 4.237702e-01 -2.863704e-01 -7.628633e-01
GPATCH2\242224_at 1.101496e-01 -3.358226e-01 -0.656001 -0.014406 -2.024513e-01 -1.037138e-02 -1.055835e+00 0.522693 -2.500000 4.395237e-01 ... 3.377338e-01 -1.245019e+00 0.809489 0.112782 1.517590e-01 5.940620e-01 -1.477183e-01 -1.779910e-01 1.171425e-01 1.539858e-01
RPL24\214143_x_at 1.868383e-01 1.476438e-02 0.003760 -0.026180 5.286380e-02 -2.765255e-15 1.267947e-01 0.125205 -2.500000 -8.126549e-02 ... 5.964594e-02 -5.714644e-02 0.044090 -0.187937 1.330189e-01 1.035764e-01 -1.438904e-01 -2.244213e-01 -3.951937e-01 3.198973e-01
PALM2\1554640_at -1.460607e-01 -2.671897e-01 -0.625874 0.321917 -2.148522e+00 -1.406787e+00 1.648780e+00 1.031431 -2.370058 -1.142304e-01 ... 7.204935e-01 6.870224e-01 1.292842 -0.257120 5.744479e-01 3.985112e-01 4.682653e-02 2.494825e-01 1.336069e-01 9.463886e-01
LOC339316\242084_at -8.695599e-02 -1.501149e-01 2.208113 -0.725581 1.596191e+00 5.343815e-02 1.648072e-01 0.454762 1.022814 -8.821926e-01 ... 1.952618e-01 -6.714563e-01 2.500000 -0.980878 0.000000e+00 4.029939e-01 -8.271509e-01 -6.114075e-01 2.050778e+00 2.693141e-01
FECH\203115_at 3.809356e-01 1.352984e-01 -0.442387 -0.219601 3.856423e-01 1.784519e-01 -1.081562e+00 0.720528 -2.500000 2.179938e-01 ... 1.106194e+00 -1.906898e-01 0.315945 -0.107265 -4.421974e-01 2.162255e-01 -9.024499e-02 -4.498444e-01 2.732963e-01 6.115749e-01
KCNJ9\207527_at 2.446787e-01 -7.501915e-01 0.221697 0.497712 -2.342015e-02 8.869240e-01 -5.943556e-01 0.204983 -2.500000 1.110048e+00 ... 7.171176e-01 9.340802e-02 0.381381 -0.367554 -1.282485e-15 5.520371e-01 3.760669e-01 -1.293106e+00 -1.034804e+00 6.951946e-01
C21orf96\220918_at 6.877719e-01 4.009709e-01 0.859799 1.189436 -1.291472e+00 -4.096633e-01 -4.467567e-01 0.340609 -2.500000 6.170242e-01 ... 1.322267e-01 -2.444516e+00 0.736984 -0.937333 1.411358e-01 -8.504788e-01 1.726166e-02 -4.542636e-01 -8.424365e-01 1.258965e+00
SAT2\225272_at 5.220952e-01 3.603069e-01 0.197867 0.709173 -1.437147e-01 6.631793e-01 -5.187783e-01 0.918461 -2.500000 -2.303941e-01 ... -4.467455e-01 -7.023133e-01 -0.510265 0.211429 -1.295257e-02 -7.275942e-01 1.857135e-01 -2.238315e-01 -3.345544e-01 5.608342e-01
EHD1\209038_s_at -1.151206e-01 7.838067e-01 0.384543 0.581821 1.054911e-01 2.080865e-01 7.126174e-01 1.657426 -2.500000 -5.027341e-02 ... 2.102385e-01 -9.875349e-01 -1.849258 -0.533168 4.666680e-01 -2.500000e+00 -6.914566e-01 -4.362436e-01 -1.043725e-01 3.922175e-01
CBLB\208348_s_at 7.738033e-01 -1.933735e+00 0.082362 0.469809 1.340361e+00 -1.081758e+00 -6.607074e-01 0.263457 -0.694266 1.307891e+00 ... 1.983327e+00 8.274782e-01 -0.140937 0.570401 -7.240150e-01 1.757159e-01 -9.647979e-01 5.032972e-01 1.024542e+00 -9.315211e-01
C20orf85\229542_at -2.908508e-03 5.987500e-01 0.254282 0.753694 -6.652577e-01 3.341133e-02 -1.511821e+00 0.860719 -2.485342 -2.410139e-01 ... -8.473255e-16 8.755296e-02 0.088931 -0.375043 2.624213e-01 -6.056148e-02 -1.464985e+00 6.121084e-02 1.220506e-01 -4.128167e-01
UGGT1\218257_s_at 1.782258e-01 -1.137344e-01 -0.617899 0.513928 3.984019e-01 6.856900e-01 -1.595056e-01 1.127709 -2.500000 4.515102e-01 ... 4.807792e-01 -1.718055e+00 0.321934 -0.424538 7.535896e-02 1.116358e-01 1.440283e-01 2.094979e-03 -8.030938e-02 6.544336e-02
LOC283585\1561442_at -6.699713e-01 -6.624370e-01 0.700004 -2.500000 1.482911e+00 8.532391e-01 4.876661e-02 -0.356320 -2.500000 4.831032e-01 ... 3.031844e-01 -1.996948e+00 0.973158 -0.977061 1.243125e+00 9.759445e-01 7.292454e-02 6.582480e-01 3.560563e-01 1.075439e+00
EMP2\204975_at -1.648862e-02 -4.951437e-01 -0.333029 0.105847 2.403381e-02 2.640703e-01 -2.462476e-01 0.435395 -2.500000 -1.813071e-01 ... 3.241176e-01 4.882364e-02 0.215725 -0.207513 4.261064e-01 1.583930e-01 -2.786914e-01 -4.835402e-01 -1.116450e-01 2.027803e-01
LHFP\218656_s_at -1.672177e-01 -8.173977e-02 0.082007 -0.039112 2.256442e-01 -1.412568e-15 -1.941996e-01 0.587759 -2.500000 6.075659e-01 ... 6.274778e-01 -2.734732e-01 0.124004 0.031477 2.239015e-01 4.204985e-01 2.847553e-02 -6.270991e-01 -2.628937e-01 3.536542e-01
SNCAIP\237834_at 4.284552e-01 -1.163575e+00 0.243658 -0.489694 -8.015225e-02 -4.781464e-02 -1.275460e-01 -0.762873 -0.175935 -2.769998e-01 ... 4.296933e-01 -2.500000e+00 1.247165 0.317115 4.454265e-01 -1.746180e+00 -6.298814e-01 -7.808986e-01 3.364997e-01 -7.405259e-01
RDH8\220683_at -3.343337e-04 2.719586e-01 -2.500000 0.331118 0.000000e+00 4.078399e-01 -1.036219e+00 0.762181 -2.500000 4.214579e-01 ... -6.652780e-01 -2.310751e+00 1.231245 0.406053 2.393105e-01 1.022069e+00 -9.697710e-01 -1.266995e+00 1.571626e-01 7.023662e-01
SRSF4\239512_at 3.433336e-02 -1.190874e+00 -0.029209 0.606425 4.751952e-01 2.473199e-01 -9.144223e-02 0.660245 -2.500000 1.593003e-02 ... 8.720908e-01 -7.960976e-01 0.753250 0.027204 5.499334e-01 7.468201e-01 1.441688e-01 3.253595e-02 5.998565e-01 7.623570e-01
ACTB\213867_x_at 3.806367e-01 -2.437849e-02 0.235402 0.203462 -8.948707e-02 7.836163e-02 3.256786e-01 0.262792 -2.500000 -6.313343e-02 ... -3.470505e-14 1.877330e-01 0.066712 -0.231704 3.914000e-01 -1.944697e-01 -3.105753e-01 -1.443515e-01 -3.226865e-01 3.269460e-01
FAM182B\215487_x_at -7.827940e-01 -1.743531e+00 1.542259 -0.815641 7.408833e-01 8.590252e-01 -5.966756e-01 0.922570 -0.327944 -1.910653e-01 ... 2.551106e-15 -1.221011e+00 0.357840 -0.717567 -1.581761e+00 5.548705e-01 1.118574e+00 4.014430e-01 6.299936e-01 -7.558542e-01
C6orf108\213871_s_at 1.283673e+00 -8.853416e-01 -0.148429 -0.091426 1.303612e+00 9.460156e-01 -9.347134e-01 -0.421391 -1.730166 2.410012e-01 ... -1.440800e-01 1.992118e-01 1.008553 1.869795 2.182196e-01 4.145624e-16 -3.242178e-02 -1.134264e-01 -3.916806e-01 -1.043065e-01
CSRNP2\225042_s_at 4.124055e-01 -8.381178e-01 0.435375 0.529975 1.396664e-01 3.393898e-02 3.863428e-01 0.942473 -2.500000 1.824226e-01 ... -2.157490e-01 -3.603765e-01 0.423193 0.130419 2.230886e-01 -1.696119e-01 -2.235140e-15 -6.171768e-03 3.178128e-01 7.780457e-01
SLCO2B1\203473_at 8.176764e-01 -2.952853e-02 0.454380 0.590203 1.299635e-01 3.416065e-01 -3.586638e-01 0.675063 -2.500000 -1.824153e-01 ... -3.013392e-01 -1.904205e-01 0.384603 0.289718 3.438642e-01 -3.036421e-01 -2.657001e-01 1.615090e-01 5.801628e-02 6.904673e-01
MYNN\237510_at -8.085378e-01 -2.362796e+00 -0.099629 -0.447579 9.998493e-01 -7.860182e-02 2.155692e-01 0.385779 -1.658656 2.851661e-01 ... 8.639631e-01 -1.698898e-01 0.028395 0.537210 2.123092e-01 2.644444e-01 -8.311871e-01 -4.848538e-02 -1.477602e-01 8.573179e-01
NAP1L1\208752_x_at 9.602593e-02 7.153991e-02 -0.160236 0.107609 2.225467e-01 1.687107e-01 -3.239409e-02 0.271276 -2.500000 1.624003e-01 ... 1.767542e-01 -3.788714e-01 0.297318 -0.034827 3.441874e-01 2.408059e-01 -1.646946e-01 -1.708129e-03 -2.785391e-01 4.158654e-01
ADAM28\205997_at 1.772528e-01 1.117742e+00 -0.731367 1.274166 1.885729e-01 2.602659e-01 -1.066456e+00 0.863145 -2.500000 -1.823998e-01 ... -1.079184e+00 3.550066e-01 -2.500000 0.006168 8.469244e-01 -7.686321e-02 -8.059085e-01 5.649578e-02 -9.739542e-01 -1.514475e-01
LOC100507248\235236_at -1.072643e+00 4.594174e-01 -0.927991 -0.655609 -7.862437e-01 -2.286303e+00 -5.365266e-02 -2.491841 -1.202169 1.742903e-15 ... -1.297602e+00 1.865324e-01 -0.693526 -1.443704 2.288546e-02 7.171869e-01 -3.005735e-01 4.729783e-02 -1.715550e+00 5.227358e-01
PHF8\238837_at -2.947776e-01 -2.915566e-01 -0.060843 -0.705987 -2.392442e-01 -4.892896e-01 -6.090763e-01 1.377218 -1.571625 1.116291e-01 ... 9.266539e-01 -3.207794e-01 1.748153 0.145365 8.748371e-03 9.995139e-01 -8.309919e-01 -4.290229e-01 -5.363052e-01 7.551622e-01
C9orf96\1554787_at 1.149944e+00 -9.559984e-01 -0.343755 -1.705164 1.327615e+00 1.246688e+00 -1.392792e-01 -0.472804 -0.277043 5.783098e-01 ... -7.369633e-01 4.219920e-01 1.099040 -0.259896 3.860776e-01 -6.389184e-01 2.974728e-01 6.941508e-01 8.204310e-01 -1.153295e+00
TNPO2\210585_s_at 8.415549e-01 -1.558058e+00 -0.108350 1.293719 2.108132e-01 1.376872e+00 6.629526e-01 1.427286 -2.259321 -4.671707e-01 ... -6.485810e-01 -7.156987e-01 -0.157029 0.966534 1.498387e-01 -6.615023e-01 -5.128144e-01 -2.662810e-02 2.944100e-01 5.898351e-01
FLJ31485\239569_at -5.540852e-01 -1.379922e+00 0.198699 -0.409129 -1.486759e+00 6.456342e-01 -3.645415e-01 0.323915 -2.500000 3.639542e-01 ... -2.603269e-15 -1.020649e+00 1.020018 0.373374 5.062770e-01 -1.502273e-01 3.419281e-01 6.765752e-02 3.784114e-01 -1.268814e-02
PUM1\201164_s_at 3.616000e-01 -8.103959e-02 -0.073435 0.156986 3.594612e-01 5.380874e-02 -4.399876e-01 0.345643 -2.500000 -2.482161e-01 ... 5.972622e-01 -2.876168e-01 0.429244 -0.091790 6.488431e-01 5.650431e-01 -5.263984e-01 -5.719179e-01 -4.366248e-01 6.478936e-01
COX4I1\202698_x_at 4.272540e-01 1.905541e-01 0.160764 0.231704 -2.588752e-01 2.185515e-01 1.461314e-01 0.416512 -2.500000 -4.172664e-02 ... -1.665304e-01 -4.477226e-02 0.132276 -0.099144 3.437716e-01 1.535771e-02 -1.547498e-01 -1.107080e-01 -1.705950e-02 4.616539e-01
IZUMO4\215734_at 9.570135e-01 -8.113009e-01 0.330161 0.648157 6.302402e-01 5.921951e-01 5.853311e-01 0.311937 -2.500000 -7.680267e-01 ... -5.166969e-01 8.629959e-02 -1.336871 0.426588 -1.975528e-02 -2.494942e-01 -2.396816e-01 1.031394e+00 7.834997e-03 8.988842e-01
TPSAB1\216485_s_at -6.078397e-01 1.184266e+00 0.983453 2.500000 1.250009e+00 -5.077744e-01 1.048912e+00 1.714357 -0.355648 -2.467265e-02 ... -4.449256e-01 -5.165284e-01 0.478172 -0.708736 -1.101193e-01 9.643122e-01 -1.608353e-01 1.718105e-01 -8.349604e-01 7.185777e-01
LMAN1L\220420_at -3.930113e-01 1.297623e+00 -0.190354 1.930165 1.658784e+00 5.384069e-03 -3.424678e-01 -0.334402 0.136873 -6.443848e-02 ... 1.369716e+00 -1.247334e+00 2.500000 -0.578041 5.389191e-01 5.576721e-01 2.955870e-01 3.838427e-01 -4.748444e-01 8.219235e-01
LY6K\223687_s_at -1.461647e+00 -2.265928e-02 0.145340 0.641209 1.941769e-01 -1.522822e+00 7.950834e-01 0.541757 -2.044659 1.781366e-01 ... 3.818696e-01 -1.102716e+00 0.495548 -0.542092 -1.535721e+00 5.965487e-01 2.741282e-01 -5.391458e-01 -1.108703e-15 -1.547624e-01
ACAD10\238655_at 3.744543e-01 -2.191149e+00 0.452272 -2.198740 -1.868530e-01 5.587617e-01 3.074604e-15 0.222280 -2.500000 7.626799e-02 ... 1.423953e-01 -4.586533e-01 0.768854 0.201477 1.869591e-01 6.070856e-01 -5.366020e-01 -2.556770e-01 2.514190e-02 4.341299e-01
TIMM22\219184_x_at 1.860355e-01 -4.626873e-01 0.259960 0.374713 5.352339e-01 -3.145453e-02 -1.147574e+00 0.528904 -2.500000 1.228076e-01 ... 7.864175e-01 -1.967618e-01 0.741393 -0.230203 4.937277e-01 8.848425e-16 -1.197899e+00 -6.914155e-01 -4.348508e-01 2.612987e-01
HUS1\217618_x_at -3.573032e-01 -7.975662e-01 1.097615 0.183515 6.254582e-01 1.400126e+00 -8.121844e-01 0.702328 -2.500000 9.160469e-01 ... 6.075843e-01 9.858510e-02 0.306931 0.394713 4.329534e-01 9.430008e-01 2.014285e-01 6.965350e-01 6.140034e-01 2.169261e-01
EPHA7\206852_at 9.453371e-01 1.055212e+00 1.330655 -1.462244 9.950828e-01 6.651014e-01 1.403666e+00 -0.917735 -2.025684 1.031804e+00 ... 1.284141e+00 1.716224e-01 -0.295550 -0.472820 5.728688e-01 -7.848282e-01 1.158740e+00 -1.160183e+00 7.239862e-01 1.487858e+00
CSH1\208069_x_at 8.291225e-01 3.785826e-01 0.508698 -0.611269 6.796478e-01 6.744793e-01 -5.521105e-01 -0.501059 -1.350985 -7.695647e-01 ... -1.704326e-01 -1.103740e-01 1.161543 -0.046680 2.500000e+00 -7.618032e-02 1.701280e+00 -4.317074e-01 -1.565720e-01 9.939125e-01
FMNL2\242665_at -4.044236e-02 -5.865016e-01 0.694424 -0.088772 4.113713e-01 8.806274e-02 -1.591566e+00 -0.115378 -2.500000 6.792634e-01 ... 3.431706e-01 -8.605117e-01 0.268740 0.598927 -6.343997e-01 3.910055e-01 1.036978e+00 1.723116e-01 3.036245e-01 -3.656789e-01
PUS3\221277_s_at 1.923533e-01 4.075534e-01 -0.082974 0.539869 2.213191e-01 4.199566e-01 -5.393975e-01 0.568746 -2.500000 -8.174332e-02 ... 4.791909e-01 -4.053745e-01 -0.244712 -0.093528 -1.603453e-01 2.567217e-01 1.749341e-01 -8.780012e-02 -4.179677e-01 6.323481e-01
CELF4\231220_at 1.376934e-02 -6.964947e-01 1.562912 1.683338 3.190191e-03 1.506210e+00 -6.010970e-01 2.389317 -0.016451 -4.306440e-01 ... -1.978729e-01 -5.503310e-01 0.073438 1.548142 3.880176e-02 2.120659e+00 1.215759e+00 -8.241560e-02 -9.420770e-01 -4.496815e-01
RBM6\228030_at 4.779255e-01 2.041114e-01 0.923409 0.414006 8.910700e-01 -1.997272e-01 -6.127431e-01 0.487421 -2.500000 1.069007e-01 ... 5.565498e-01 3.281905e-01 0.694981 -0.560279 3.989280e-01 1.095893e+00 -8.442382e-01 -8.199109e-01 -9.166730e-02 6.334560e-01
AVIL\233333_x_at -9.783464e-01 -8.856929e-01 0.766268 0.590044 8.821904e-01 6.341817e-15 -6.550214e-01 1.043906 -2.500000 -9.072674e-02 ... 9.266262e-01 -2.048787e+00 1.667175 0.164728 4.518892e-02 -2.652245e-01 -1.004089e-01 -5.809123e-01 -7.326943e-01 9.717940e-01
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
MON1B\233557_s_at 8.810125e-01 -8.945361e-01 0.859893 1.026511 -1.432378e-01 1.829283e-15 6.114901e-01 1.870331 -2.500000 -3.662924e-01 ... -5.728325e-02 -2.661090e-02 1.288685 0.026459 9.216581e-01 -1.923043e-01 1.823643e-01 7.080469e-01 -6.358982e-02 5.178159e-01
CLLU1\1558185_at 1.430754e+00 1.600067e+00 0.205752 1.844058 1.180818e+00 -7.810344e-01 1.170296e+00 1.839945 0.358849 1.499494e+00 ... -8.249987e-01 1.817452e+00 0.213771 -1.334990 -6.085725e-01 -6.868481e-01 1.461691e+00 -3.411989e-01 1.006889e+00 -4.339588e-01
KTN1\214709_s_at 3.844828e-02 -1.179759e-01 0.128939 -0.142158 6.254535e-01 3.115187e-01 -5.017092e-01 0.439360 -2.500000 1.275202e-01 ... 6.770990e-01 -4.167230e-01 0.304909 -0.190950 3.769602e-01 2.094832e-01 -1.297410e-01 -1.386177e-01 -3.966711e-01 4.634931e-01
KRTAP4-7\216810_at 1.320582e-01 -3.909804e-01 0.371147 -1.050310 2.146137e-01 -1.459212e-01 1.285379e+00 0.313965 -0.859916 2.764068e-01 ... 7.241434e-01 -6.335071e-02 1.508098 -0.823065 1.238118e+00 -1.143154e-01 -5.926425e-02 1.169225e+00 3.633507e-01 -4.482008e-01
RBKS\219222_at 8.771964e-01 6.335409e-01 0.379687 1.048153 -3.290674e-01 -6.103011e-01 2.530613e-01 1.375255 -2.500000 -6.246420e-01 ... 1.816752e-01 4.958249e-02 0.869107 -0.619181 6.231288e-01 5.711128e-01 -4.304674e-01 4.280768e-01 -9.638054e-01 8.646066e-01
LOC100506621\227061_at -9.809685e-02 1.704357e-02 1.378401 0.101114 -3.492839e-15 -1.135276e+00 2.534411e-01 0.526410 -2.500000 4.511687e-01 ... -3.947754e-01 -5.189949e-02 0.768363 -0.989644 -4.984056e-02 2.476006e-01 -9.890471e-02 -4.462636e-01 4.291749e-02 4.474606e-01
CCDC67\1554071_at 8.725676e-01 1.336801e+00 -0.277102 1.659277 -1.640883e-16 -8.522715e-01 1.393413e+00 0.394425 -2.500000 -1.011547e+00 ... 1.246248e-01 -4.246001e-01 1.143442 -0.145376 1.621148e+00 -1.311573e-01 1.093230e+00 1.401844e+00 -2.792740e-01 -7.094043e-01
MCM9\219673_at 3.667066e-01 -3.280029e-01 0.178328 0.352335 5.635915e-01 8.023715e-01 -8.550817e-01 0.681810 -2.500000 5.930072e-01 ... 1.829210e-01 -1.001783e+00 0.065792 0.150269 3.375726e-01 -1.414786e-01 1.637869e-01 -2.240233e-01 5.632103e-02 4.680923e-01
LOC730101\213248_at 1.105918e-01 6.911357e-01 0.480355 0.364861 7.344180e-01 6.011730e-01 -1.075058e+00 0.623014 -2.500000 2.793485e-01 ... 7.322634e-01 -1.684741e-01 0.452439 -0.302283 -6.829856e-02 7.417542e-01 2.116920e-01 -4.113453e-01 -3.517702e-01 1.225708e+00
LOC100508546\244037_at 2.490395e-01 4.797762e-01 0.456635 0.498514 1.499252e+00 -2.834839e-02 -1.883011e+00 0.673851 -2.500000 4.316784e-01 ... -1.763166e-01 -6.568147e-01 -0.730289 0.285510 -2.138232e+00 -2.682613e-01 -9.140803e-01 -2.113736e+00 -3.117951e-01 6.025652e-01
RNPS1\200060_s_at 1.969155e-01 8.339107e-02 0.060937 0.563833 -1.760587e-01 2.259313e-01 -1.988471e-01 0.948589 -2.500000 -1.236302e-02 ... -2.785834e-02 -7.577448e-01 0.083098 0.104838 2.969498e-01 -9.367385e-02 -5.940128e-02 -1.382211e-15 -1.962634e-01 4.175789e-01
C8orf74\1569246_a_at -1.588470e-15 -1.263879e+00 -0.568192 0.916387 1.467641e+00 1.005052e+00 4.655200e-01 1.023241 -1.648073 -1.059640e+00 ... 8.633091e-01 7.904270e-01 1.893826 1.071859 -8.684117e-03 6.454497e-01 7.726947e-01 -1.148032e+00 1.151162e+00 -9.630953e-01
UTP14C\203614_at 2.956424e-01 -4.169494e-02 -0.051720 0.094925 2.115923e-01 1.629282e-01 -4.304645e-01 0.403530 -2.500000 8.880886e-02 ... 3.864589e-01 -1.518737e-01 0.269590 0.031779 3.964468e-01 3.702738e-01 -7.403777e-02 -4.434015e-01 -1.775024e-01 4.281801e-01
RHEB\201452_at 7.536154e-01 -1.326204e+00 0.549601 -0.296880 -2.060105e+00 -9.396856e-01 7.855606e-01 0.724879 -2.500000 -7.087826e-01 ... -2.097672e+00 2.704643e-01 -1.281542 0.897462 6.127546e-01 6.354196e-01 -7.908111e-01 4.390502e-01 -8.706542e-02 5.826829e-01
CXCL16\223454_at 6.844901e-01 -2.374285e-01 0.502841 1.002893 -1.351760e-03 4.880478e-01 7.649009e-03 0.700994 -2.500000 -2.956657e-01 ... -7.074534e-01 -4.119786e-01 0.303856 0.100922 2.757123e-01 -2.456973e-01 1.392537e-15 2.917335e-01 1.385660e-01 3.440072e-01
PSPC1\218371_s_at 1.626253e-01 5.674949e-02 0.221091 0.411380 2.193392e-01 2.966490e-01 -4.952883e-01 0.757099 -2.500000 3.823820e-01 ... 1.847272e-01 -7.089064e-01 0.112730 0.307784 4.462130e-01 -4.546926e-02 -2.685357e-01 -2.925314e-01 -4.102348e-01 4.417630e-01
SIK3\213034_at 2.113850e-01 4.319275e-02 -0.160323 0.061715 2.944546e-01 1.844371e-01 3.111475e-01 0.951732 -2.500000 3.382042e-02 ... 4.818038e-01 -4.862550e-01 0.166349 -0.084095 6.242145e-01 -8.634988e-02 -3.233842e-01 -6.604433e-02 -3.057007e-01 8.177099e-01
ACOT9\221641_s_at 7.807096e-01 6.206924e-02 0.177504 0.466942 -2.468935e-01 3.845749e-01 9.261722e-01 0.874142 -2.500000 -1.452015e-01 ... -1.259537e-01 -3.770752e-01 0.300626 -0.163270 4.699680e-01 4.756456e-02 -3.029091e-01 2.857613e-01 1.896554e-01 2.485417e-01
KIAA1661\215943_at 6.726885e-01 1.706044e-01 0.382563 -1.789371 -1.751758e+00 2.479091e-01 5.441798e-01 1.195781 -0.542833 -7.176655e-01 ... -7.684039e-01 -1.853766e+00 1.193023 -1.313194 -1.589782e+00 1.376141e+00 -1.193792e-02 -1.992654e+00 -1.971444e+00 9.655266e-01
STEAP2\225871_at 1.372030e-01 1.417692e+00 0.573894 0.425106 1.618731e-01 9.512809e-02 1.014428e+00 0.575258 -2.500000 -9.830552e-02 ... 2.959283e-01 -6.333808e-01 0.039590 -0.342473 -2.063003e-02 2.447803e-01 -4.129984e-01 4.213754e-01 -5.859930e-01 2.143950e-01
AKAP5\207800_at -3.620299e-01 1.467300e-01 0.022813 -0.369833 -6.366712e-01 -5.551043e-02 -4.435810e-01 -0.490820 -2.500000 1.217579e+00 ... 1.125127e+00 -3.930062e-01 1.293941 1.228261 2.201931e-02 -4.610851e-01 -5.629784e-01 3.899632e-01 -4.174401e-01 1.039304e+00
C1orf57\223272_s_at 2.111401e-01 -4.749875e-01 0.036480 0.124417 1.142755e-01 8.470742e-02 -6.151017e-02 0.462233 -2.500000 8.278270e-02 ... 1.108331e-01 -6.756677e-16 0.068329 0.399466 5.782434e-01 -1.678051e-01 -8.818476e-02 -1.767640e+00 -3.915039e-01 6.837272e-01
C7orf30\226385_s_at 2.030614e-01 1.313962e-01 -0.370263 0.286871 -8.362649e-02 2.782752e-01 -1.951197e-01 0.743667 -2.500000 8.539293e-02 ... 1.817853e-01 -5.463139e-01 0.073956 0.021864 1.023310e-01 3.883427e-02 -8.576460e-02 -2.361815e-01 -4.026115e-01 6.456563e-01
BAX\208478_s_at 6.522254e-01 -9.436456e-02 0.175317 0.388198 -1.991714e+00 1.828135e+00 -1.100958e+00 0.427227 -2.500000 9.623890e-01 ... -7.993298e-01 -1.042535e+00 -0.461467 1.065001 2.012570e-01 -4.096498e-01 6.740900e-01 1.247874e+00 2.074081e-01 -2.205594e-01
MAN1C1\218918_at 2.800930e-01 2.238196e-01 0.273756 0.033026 -1.798884e-01 5.012985e-15 -5.929027e-01 1.029952 -2.500000 3.992565e-02 ... -3.852458e-01 -4.932705e-01 0.028520 0.219546 7.410833e-01 1.358944e-01 -1.385487e-01 -5.581595e-01 6.404567e-03 6.373228e-01
SGSM3\214779_s_at -1.736405e-01 -1.414195e-02 0.157717 0.864444 2.873663e-01 8.659728e-01 -8.545423e-01 1.186905 -2.500000 3.561603e-01 ... -3.425757e-01 -1.025167e+00 -0.241786 0.283231 -4.412527e-01 -2.509460e-01 3.779060e-01 3.711572e-02 3.272294e-01 3.290084e-01
LOC100508009\231516_at -7.673898e-01 1.022311e-01 -0.403244 2.475959 -1.423754e-01 7.237218e-01 1.692347e-03 -0.233173 -1.494539 2.179593e-01 ... 1.722962e+00 -7.545144e-02 1.021933 0.293152 -6.434770e-01 1.028962e+00 -1.087279e+00 1.851378e+00 7.203381e-01 2.696586e-01
ASCL3\221161_at -2.650266e-01 -1.325737e-15 -1.215097 -1.083043 -7.659126e-01 -6.959607e-01 -3.357083e-01 -0.056452 -2.324118 -1.157994e+00 ... -6.936359e-01 -6.844133e-02 1.191375 0.437159 -2.974882e-01 8.986177e-01 -2.547238e-01 -7.817435e-01 6.498341e-02 1.799959e+00
DALRD3\221934_s_at 3.619755e-02 1.309325e-01 0.605709 1.211953 -1.332750e-01 1.277186e-01 -5.967899e-01 1.513510 -2.500000 1.499669e-01 ... 6.480116e-02 -6.160821e-01 -0.041037 0.099933 -9.786916e-02 2.601851e-02 -4.701461e-01 -3.642525e-01 -1.471343e-15 5.833984e-01
DAPK2\206324_s_at -1.390897e-01 -1.290547e+00 0.091410 -0.476111 1.444616e-01 -3.197286e-02 -3.866668e-01 1.612002 -2.500000 -3.402872e-02 ... 3.369926e-01 -4.244189e-01 0.385729 0.068463 4.582797e-01 -3.322440e-01 2.358460e-01 -7.145301e-01 9.155446e-02 2.796026e-01
NUP155\206550_s_at 3.562320e-02 1.544426e-01 -0.124364 0.108441 1.064240e-01 1.284462e-01 -2.895961e-01 0.553523 -2.500000 1.393357e-02 ... 9.987133e-02 -2.454252e-01 0.117596 -0.315022 1.098486e-02 2.075482e-01 -4.582731e-01 -5.106113e-01 1.502304e-01 3.054965e-01
TMF1\215855_s_at -2.901360e-01 -4.651857e-02 0.144372 0.909981 5.628572e-01 -1.369268e-01 -8.955115e-01 0.935953 -2.500000 -1.370132e-01 ... 8.288251e-01 2.629762e-01 0.961501 -0.311875 7.434090e-01 6.298310e-01 -2.527550e-01 -5.429247e-01 -3.419514e-02 7.723467e-01
LOC100505584\216336_x_at -4.826152e-01 1.493645e-01 0.333728 0.444856 -6.005730e-01 3.951811e-01 1.598540e+00 0.604569 -2.500000 -2.678052e-01 ... -8.367125e-02 -5.951153e-01 0.174207 -0.295462 6.472981e-01 1.078961e-01 -1.850171e-02 6.805104e-01 2.702985e-01 2.198597e-01
ARID1B\1566990_x_at -8.956404e-01 -5.731021e-01 0.254280 0.525842 7.699408e-01 4.785215e-01 -1.556190e+00 0.859796 -2.500000 7.274728e-01 ... 6.735043e-01 -2.276432e+00 0.880813 0.524485 -3.145565e-01 7.334851e-01 5.185571e-01 -6.726688e-16 5.869926e-01 6.323548e-01
LOC100287105\1570171_at 3.669336e-01 -2.334089e+00 0.654090 -2.121761 6.198339e-01 4.949633e-01 -4.233259e-01 -1.688630 -2.500000 3.427960e-01 ... 6.813099e-01 -9.100254e-02 1.191427 0.239410 -5.198894e-02 1.776968e-01 1.601922e-01 4.530840e-01 3.078466e-01 6.194966e-01
UTRN\225093_at 2.110852e-01 -2.492763e-01 -0.009004 0.107533 2.733667e-01 5.203354e-02 -2.838437e-01 0.671443 -2.500000 2.746188e-02 ... 5.733829e-01 -2.718756e-01 0.359423 -0.107341 6.024932e-01 2.935127e-01 -2.982864e-01 -4.122579e-01 -3.569925e-01 7.365671e-01
HIST3H3\208572_at 4.576939e-01 -1.722807e+00 0.453559 -0.402484 1.055068e-01 1.041188e-02 -1.265198e+00 1.267053 0.456003 -1.124851e+00 ... 6.130563e-01 -1.454852e+00 -0.149284 1.013416 -6.808040e-01 2.418817e+00 6.606057e-01 -9.158784e-01 -2.269791e-01 -6.137322e-01
ZNF316\227170_at -2.507152e-01 -3.592175e-01 0.408765 1.594587 -3.635144e-02 3.424956e-01 -1.861386e+00 1.867112 -2.500000 -2.878322e-01 ... 5.983806e-01 -5.050046e-01 0.746700 0.263718 8.621742e-02 -6.545053e-01 -1.689681e-01 1.603481e-01 -7.618066e-01 3.682123e-01
ETV6\205585_at -2.991069e-01 -4.313853e-02 0.164672 0.776184 9.948931e-01 -1.696238e-01 -1.300322e-01 1.418278 -2.500000 -5.578391e-01 ... -2.082721e-01 -3.309314e-01 0.685810 0.360216 -6.797032e-02 7.947080e-02 1.413004e+00 3.015267e-01 7.651200e-01 -2.822372e-01
LOC400550\243059_at 3.272600e-01 -1.919577e+00 0.467922 0.263556 -1.140731e+00 7.985130e-01 -3.743171e-01 1.090638 -2.500000 5.887833e-01 ... -1.039014e+00 -1.392440e-01 -0.083357 1.058742 1.019857e+00 3.933854e-01 -3.014409e-01 -6.822756e-01 5.185212e-01 3.989441e-02
B3GAT2\234146_at -9.805872e-01 3.168191e-16 2.257970 -0.466146 6.546636e-01 4.193820e-01 1.345748e-01 0.456639 -0.615108 2.781852e-01 ... 1.944672e+00 1.190688e+00 1.062559 0.110046 -6.473686e-01 1.818159e+00 -2.247820e-01 -1.332021e-01 1.705207e-01 -1.003410e+00
HMGB4\230473_s_at -8.205253e-01 3.830290e-01 -0.424356 -1.008983 3.071867e-01 -8.664867e-01 1.288813e+00 -0.028566 -0.459913 1.310405e-01 ... -5.256141e-01 -4.916938e-01 1.890940 1.571216 1.680680e+00 1.689341e-01 1.654385e+00 -1.252686e+00 -3.768941e-01 4.776298e-01
CMTM6\217947_at 2.871199e-01 1.925319e-01 0.066176 0.192016 3.497992e-01 -1.374658e-01 3.960978e-02 0.357320 -2.500000 -4.072392e-01 ... 2.346540e-01 -2.945002e-01 0.515876 -0.228873 2.877550e-01 1.303690e-01 -2.318001e-01 -1.057098e-01 -3.018645e-01 3.103722e-01
BRAP\209922_at 2.883632e-01 3.656153e-01 -1.425366 0.573506 -1.313536e-01 6.251475e-01 -8.979595e-01 0.938397 -2.500000 -5.492051e-02 ... 3.275027e-02 -2.500000e+00 -1.223492 0.123448 1.318926e-01 -6.709925e-02 -1.196013e+00 -3.523880e-15 1.273002e-02 1.054824e+00
BARD1\205345_at 4.120692e-01 8.377061e-01 0.050744 0.618542 7.031704e-01 -3.162949e-02 -5.587457e-01 0.663091 -2.500000 -3.677368e-01 ... 1.196425e+00 -4.250562e-01 1.041858 -0.439212 3.288060e-01 3.428453e-01 -6.847699e-01 -6.303019e-01 -3.214081e-01 7.341888e-01
MR1\207565_s_at 5.335984e-01 -5.450360e-01 0.386406 0.285703 -5.637429e-01 -4.158945e-15 -4.111376e-01 0.770155 -2.500000 1.582400e-01 ... -5.757536e-01 -5.854919e-01 0.503805 0.092773 2.414565e-01 1.434490e-01 2.827439e-01 1.890089e-01 1.101660e-01 1.935011e-02
SORT1\212797_at 6.156220e-01 -1.648567e+00 -0.077905 0.099130 -4.616306e-01 7.210408e-01 -1.818820e-01 0.795778 -2.500000 4.109577e-01 ... -4.905258e-01 -4.931033e-01 0.246322 0.293623 7.029947e-01 3.494057e-01 -4.095347e-02 1.013247e-01 1.701813e-01 3.615812e-01
FLJ31485\243351_at 8.793840e-01 4.511735e-01 -0.101133 0.709585 1.520049e+00 1.788267e-01 5.830351e-02 0.267418 -1.713196 2.380288e-02 ... 8.141621e-01 -4.732190e-01 1.408914 -0.250313 -1.316088e+00 2.594872e-01 -1.300986e-01 -7.129734e-01 1.771152e-02 -1.846044e-02
BPI\205557_at 2.617538e-01 -8.212551e-01 1.276366 1.074141 -8.965759e-01 1.498612e+00 1.481477e+00 1.118083 -0.380919 1.550287e+00 ... -1.012536e+00 7.006553e-01 -0.088471 -1.236672 1.183149e+00 1.338136e-01 -8.761095e-02 2.194417e+00 -5.984596e-01 1.404617e+00
CECR3\1564352_at -1.160974e+00 -2.776760e-01 0.050318 -0.114216 1.531258e-01 4.567647e-01 -1.399043e+00 0.450699 -2.500000 5.172905e-01 ... 6.292538e-01 5.616309e-01 -1.321470 0.141950 2.682927e-01 6.744392e-01 1.913051e-01 7.238533e-03 -1.794781e-01 -3.650810e-15

42449 rows × 55 columns


In [44]:
design_df = DataFrame({'control': control[col_names], 'treatment': treatment[col_names]})
design_df


Out[44]:
control treatment
subject.inTrialId
GSE8581GSM213034 False False
GSE8581GSM212811 True False
GSE8581GSM213036 True False
GSE8581GSM212075 False True
GSE8581GSM211008 True False
GSE8581GSM210090 False False
GSE8581GSM212855 False False
GSE8581GSM212070 True False
GSE8581GSM212810 False True
GSE8581GSM210193 False False
GSE8581GSM210071 False True
GSE8581GSM212069 False False
GSE8581GSM212849 False False
GSE8581GSM212787 True False
GSE8581GSM210015 True False
GSE8581GSM212853 True False
GSE8581GSM212790 True False
GSE8581GSM212067 True False
GSE8581GSM211865 False False
GSE8581GSM210196 True False
GSE8581GSM212809 False False
GSE8581GSM210992 False True
GSE8581GSM211009 False True
GSE8581GSM210014 True False
GSE8581GSM213018 False False
GSE8581GSM212789 True False
GSE8581GSM210978 False False
GSE8581GSM213037 False False
GSE8581GSM212074 False True
GSE8581GSM210993 False True
GSE8581GSM210087 True False
GSE8581GSM212068 True False
GSE8581GSM210012 False False
GSE8581GSM212854 False False
GSE8581GSM212852 False True
GSE8581GSM212848 False True
GSE8581GSM211010 False False
GSE8581GSM213035 True False
GSE8581GSM210194 False True
GSE8581GSM210008 False True
GSE8581GSM212788 False True
GSE8581GSM211872 False False
GSE8581GSM212850 False False
GSE8581GSM210188 False False
GSE8581GSM210009 True False
GSE8581GSM210979 True False
GSE8581GSM213017 False False
GSE8581GSM211007 False True
GSE8581GSM210192 True False
GSE8581GSM213020 False True
GSE8581GSM210994 False False
GSE8581GSM210010 False False
GSE8581GSM210011 False True
GSE8581GSM210007 False False
GSE8581GSM213019 False False

In [45]:
%load_ext rpy2.ipython


The rpy2.ipython extension is already loaded. To reload it, use:
  %reload_ext rpy2.ipython

In [46]:
%%R -i hd_df -i design_df -o top_fit

library(limma)

design.matrix <- sapply(design_df, as.numeric)
contrast.matrix <- makeContrasts(control - treatment, levels = design.matrix)

fit <- lmFit(hd_df, design.matrix)
fit <- contrasts.fit(fit, contrast.matrix)
fit <- eBayes(fit)

contr=1
top_fit = data.frame(
    probe = rownames(fit$coefficients),
    log_fc = fit$coefficients[, contr],
    t = fit$t[, contr],
    p_value = fit$p.value[, contr],
    b = fit$lods[, contr],
    stringsAsFactors = FALSE
)

volcanoplot(fit)



In [47]:
top_probes = top_fit[top_fit.p_value < 0.001] \
.set_index(['probe']).sort(['p_value'], ascending = True)
top_probes


Out[47]:
log_fc t p_value b
probe
DNM3\1558502_s_at 1.375778 4.806308 0.000009 1.945957
CASC4\1552719_at 1.130824 4.687018 0.000014 1.650206
ZNF30\232014_at 1.224315 4.522312 0.000026 1.246204
LOC389634\1560119_at 1.204393 4.320760 0.000053 0.759543
RP9\227852_at 0.921346 4.312707 0.000055 0.740289
SLITRK6\232481_s_at 1.113788 4.223494 0.000075 0.528059
AGTPBP1\204499_at 1.030825 4.211223 0.000078 0.499022
LOC389906\1558045_a_at 1.219038 4.182897 0.000086 0.432140
GFRA2\205722_s_at 1.246157 4.171039 0.000090 0.404202
IL1F8\224230_at -1.228913 -4.110874 0.000111 0.263023
KLF8\230986_at 1.056070 4.102448 0.000114 0.243327
RGS8\234297_at 1.093250 4.083148 0.000122 0.198289
STRBP\223246_s_at 1.181890 4.008383 0.000157 0.024786
KIAA1377\236325_at 1.056305 3.958482 0.000186 -0.090131
NAPEPLD\242635_s_at 1.170465 3.953391 0.000190 -0.101814
EGFEM1P\1558411_at 1.012719 3.919803 0.000212 -0.178704
ATF7IP2\219870_at 1.012826 3.880784 0.000242 -0.267603
LOC644192\232878_at 1.116624 3.875276 0.000247 -0.280114
LOC100508546\244037_at 1.070067 3.866795 0.000254 -0.299360
RUNX1T1\216831_s_at 1.200380 3.863911 0.000256 -0.305900
ANGEL1\244685_at -1.170453 -3.844109 0.000274 -0.350735
SEMA6D\226492_at 0.998893 3.835239 0.000282 -0.370778
C9orf85\235866_at 1.025476 3.811178 0.000305 -0.425026
HHEX\215933_s_at 0.876758 3.788431 0.000329 -0.476140
CLDN15\1552845_at 1.070289 3.787280 0.000330 -0.478723
KLHL24\206551_x_at 0.886946 3.767771 0.000352 -0.522422
NHLRC3\236953_s_at 0.998617 3.766816 0.000353 -0.524558
ELMOD1\231930_at 1.035745 3.735181 0.000392 -0.595148
ARHGAP29\1558280_s_at 0.960541 3.717077 0.000416 -0.635398
MGC12488\211386_at 0.925050 3.713993 0.000420 -0.642242
GOLM1\1558248_at 1.133981 3.710102 0.000426 -0.650875
MTHFSD\230442_at 1.004106 3.708857 0.000427 -0.653636
LOC400590\1560992_at 1.137758 3.707311 0.000430 -0.657063
LOC729799\217667_at 1.135771 3.706478 0.000431 -0.658910
ZNF334\220022_at 1.084147 3.684612 0.000463 -0.707301
ZFP82\235176_at 0.920436 3.675480 0.000477 -0.727463
AK2\212172_at 1.019343 3.665957 0.000492 -0.748457
FAM173B\225668_at 1.061109 3.661258 0.000499 -0.758805
ZFP14\232911_at 0.892159 3.647028 0.000523 -0.790097
GTF2H2B\215470_at 0.959996 3.640101 0.000534 -0.805306
LOC339047\1558075_at 0.882082 3.629213 0.000554 -0.829174
CREBZF\225595_at 0.935700 3.626421 0.000559 -0.835290
LOC286272\232182_at 0.930666 3.622087 0.000566 -0.844775
LOC100507100\230187_s_at 1.060942 3.620916 0.000569 -0.847338
GATAD2B\238076_at 0.738321 3.617144 0.000576 -0.855586
ABCA8\1565778_at 1.045979 3.615368 0.000579 -0.859468
PTPRB\1560105_at 1.022702 3.613720 0.000582 -0.863069
ZNRF1\225962_at 0.900862 3.611676 0.000586 -0.867536
GPATCH8\212487_at 0.830247 3.561935 0.000687 -0.975750
TSC2\215624_at 1.116368 3.557750 0.000696 -0.984815
EDA\211127_x_at 1.070483 3.555175 0.000702 -0.990390
LOC441259\1552622_s_at 0.897353 3.551402 0.000711 -0.998555
CASP8AP2\1570001_at 1.032101 3.547817 0.000719 -1.006305
CTBP2\215377_at 1.018589 3.534816 0.000749 -1.034381
CACNA1D\243334_at 1.026184 3.529423 0.000762 -1.046007
FBXL13\1556770_a_at 0.945425 3.525738 0.000771 -1.053947
RRP9\204133_at 0.998497 3.518901 0.000788 -1.068664
ZNF521\226676_at 0.793170 3.518051 0.000790 -1.070491
ITGB1BP2\219829_at 1.042978 3.514850 0.000798 -1.077375
MAGIX\219965_s_at 1.088856 3.492431 0.000857 -1.125480
FOXJ1\205906_at 1.039986 3.492429 0.000857 -1.125482
LOC644613\1563426_a_at 1.035299 3.487102 0.000871 -1.136885
EFCAB7\226800_at 0.884007 3.483128 0.000882 -1.145385
LOC100294362\1569078_at 0.984446 3.475146 0.000905 -1.162441
GTPBP10\230684_at 1.062098 3.469254 0.000922 -1.175013
ZNF440\215892_at 0.901321 3.466112 0.000931 -1.181715
WDR52\221103_s_at 0.991380 3.461480 0.000945 -1.191585
NCRNA00182\1558515_at 0.838394 3.460203 0.000948 -1.194305
CXorf56\235442_at 0.995492 3.457467 0.000956 -1.200130
TOP3B\215781_s_at 0.946690 3.447509 0.000987 -1.221308

In [48]:
def get_gene_name(probe):
    return probe.split('\\')[0]

top_probes['gene_name'] = top_probes.index.map(get_gene_name)
top_probes


Out[48]:
log_fc t p_value b gene_name
probe
DNM3\1558502_s_at 1.375778 4.806308 0.000009 1.945957 DNM3
CASC4\1552719_at 1.130824 4.687018 0.000014 1.650206 CASC4
ZNF30\232014_at 1.224315 4.522312 0.000026 1.246204 ZNF30
LOC389634\1560119_at 1.204393 4.320760 0.000053 0.759543 LOC389634
RP9\227852_at 0.921346 4.312707 0.000055 0.740289 RP9
SLITRK6\232481_s_at 1.113788 4.223494 0.000075 0.528059 SLITRK6
AGTPBP1\204499_at 1.030825 4.211223 0.000078 0.499022 AGTPBP1
LOC389906\1558045_a_at 1.219038 4.182897 0.000086 0.432140 LOC389906
GFRA2\205722_s_at 1.246157 4.171039 0.000090 0.404202 GFRA2
IL1F8\224230_at -1.228913 -4.110874 0.000111 0.263023 IL1F8
KLF8\230986_at 1.056070 4.102448 0.000114 0.243327 KLF8
RGS8\234297_at 1.093250 4.083148 0.000122 0.198289 RGS8
STRBP\223246_s_at 1.181890 4.008383 0.000157 0.024786 STRBP
KIAA1377\236325_at 1.056305 3.958482 0.000186 -0.090131 KIAA1377
NAPEPLD\242635_s_at 1.170465 3.953391 0.000190 -0.101814 NAPEPLD
EGFEM1P\1558411_at 1.012719 3.919803 0.000212 -0.178704 EGFEM1P
ATF7IP2\219870_at 1.012826 3.880784 0.000242 -0.267603 ATF7IP2
LOC644192\232878_at 1.116624 3.875276 0.000247 -0.280114 LOC644192
LOC100508546\244037_at 1.070067 3.866795 0.000254 -0.299360 LOC100508546
RUNX1T1\216831_s_at 1.200380 3.863911 0.000256 -0.305900 RUNX1T1
ANGEL1\244685_at -1.170453 -3.844109 0.000274 -0.350735 ANGEL1
SEMA6D\226492_at 0.998893 3.835239 0.000282 -0.370778 SEMA6D
C9orf85\235866_at 1.025476 3.811178 0.000305 -0.425026 C9orf85
HHEX\215933_s_at 0.876758 3.788431 0.000329 -0.476140 HHEX
CLDN15\1552845_at 1.070289 3.787280 0.000330 -0.478723 CLDN15
KLHL24\206551_x_at 0.886946 3.767771 0.000352 -0.522422 KLHL24
NHLRC3\236953_s_at 0.998617 3.766816 0.000353 -0.524558 NHLRC3
ELMOD1\231930_at 1.035745 3.735181 0.000392 -0.595148 ELMOD1
ARHGAP29\1558280_s_at 0.960541 3.717077 0.000416 -0.635398 ARHGAP29
MGC12488\211386_at 0.925050 3.713993 0.000420 -0.642242 MGC12488
GOLM1\1558248_at 1.133981 3.710102 0.000426 -0.650875 GOLM1
MTHFSD\230442_at 1.004106 3.708857 0.000427 -0.653636 MTHFSD
LOC400590\1560992_at 1.137758 3.707311 0.000430 -0.657063 LOC400590
LOC729799\217667_at 1.135771 3.706478 0.000431 -0.658910 LOC729799
ZNF334\220022_at 1.084147 3.684612 0.000463 -0.707301 ZNF334
ZFP82\235176_at 0.920436 3.675480 0.000477 -0.727463 ZFP82
AK2\212172_at 1.019343 3.665957 0.000492 -0.748457 AK2
FAM173B\225668_at 1.061109 3.661258 0.000499 -0.758805 FAM173B
ZFP14\232911_at 0.892159 3.647028 0.000523 -0.790097 ZFP14
GTF2H2B\215470_at 0.959996 3.640101 0.000534 -0.805306 GTF2H2B
LOC339047\1558075_at 0.882082 3.629213 0.000554 -0.829174 LOC339047
CREBZF\225595_at 0.935700 3.626421 0.000559 -0.835290 CREBZF
LOC286272\232182_at 0.930666 3.622087 0.000566 -0.844775 LOC286272
LOC100507100\230187_s_at 1.060942 3.620916 0.000569 -0.847338 LOC100507100
GATAD2B\238076_at 0.738321 3.617144 0.000576 -0.855586 GATAD2B
ABCA8\1565778_at 1.045979 3.615368 0.000579 -0.859468 ABCA8
PTPRB\1560105_at 1.022702 3.613720 0.000582 -0.863069 PTPRB
ZNRF1\225962_at 0.900862 3.611676 0.000586 -0.867536 ZNRF1
GPATCH8\212487_at 0.830247 3.561935 0.000687 -0.975750 GPATCH8
TSC2\215624_at 1.116368 3.557750 0.000696 -0.984815 TSC2
EDA\211127_x_at 1.070483 3.555175 0.000702 -0.990390 EDA
LOC441259\1552622_s_at 0.897353 3.551402 0.000711 -0.998555 LOC441259
CASP8AP2\1570001_at 1.032101 3.547817 0.000719 -1.006305 CASP8AP2
CTBP2\215377_at 1.018589 3.534816 0.000749 -1.034381 CTBP2
CACNA1D\243334_at 1.026184 3.529423 0.000762 -1.046007 CACNA1D
FBXL13\1556770_a_at 0.945425 3.525738 0.000771 -1.053947 FBXL13
RRP9\204133_at 0.998497 3.518901 0.000788 -1.068664 RRP9
ZNF521\226676_at 0.793170 3.518051 0.000790 -1.070491 ZNF521
ITGB1BP2\219829_at 1.042978 3.514850 0.000798 -1.077375 ITGB1BP2
MAGIX\219965_s_at 1.088856 3.492431 0.000857 -1.125480 MAGIX
FOXJ1\205906_at 1.039986 3.492429 0.000857 -1.125482 FOXJ1
LOC644613\1563426_a_at 1.035299 3.487102 0.000871 -1.136885 LOC644613
EFCAB7\226800_at 0.884007 3.483128 0.000882 -1.145385 EFCAB7
LOC100294362\1569078_at 0.984446 3.475146 0.000905 -1.162441 LOC100294362
GTPBP10\230684_at 1.062098 3.469254 0.000922 -1.175013 GTPBP10
ZNF440\215892_at 0.901321 3.466112 0.000931 -1.181715 ZNF440
WDR52\221103_s_at 0.991380 3.461480 0.000945 -1.191585 WDR52
NCRNA00182\1558515_at 0.838394 3.460203 0.000948 -1.194305 NCRNA00182
CXorf56\235442_at 0.995492 3.457467 0.000956 -1.200130 CXorf56
TOP3B\215781_s_at 0.946690 3.447509 0.000987 -1.221308 TOP3B

In [49]:
from utils import Entrez

entrez = Entrez()

top_probes['gene_id'] = top_probes.gene_name.map(entrez.get_gene_id)
top_probes


Out[49]:
log_fc t p_value b gene_name gene_id
probe
DNM3\1558502_s_at 1.375778 4.806308 0.000009 1.945957 DNM3 26052
CASC4\1552719_at 1.130824 4.687018 0.000014 1.650206 CASC4 113201
ZNF30\232014_at 1.224315 4.522312 0.000026 1.246204 ZNF30 90075
LOC389634\1560119_at 1.204393 4.320760 0.000053 0.759543 LOC389634 389634
RP9\227852_at 0.921346 4.312707 0.000055 0.740289 RP9 6100
SLITRK6\232481_s_at 1.113788 4.223494 0.000075 0.528059 SLITRK6 84189
AGTPBP1\204499_at 1.030825 4.211223 0.000078 0.499022 AGTPBP1 23287
LOC389906\1558045_a_at 1.219038 4.182897 0.000086 0.432140 LOC389906 389906
GFRA2\205722_s_at 1.246157 4.171039 0.000090 0.404202 GFRA2 2675
IL1F8\224230_at -1.228913 -4.110874 0.000111 0.263023 IL1F8 None
KLF8\230986_at 1.056070 4.102448 0.000114 0.243327 KLF8 11279
RGS8\234297_at 1.093250 4.083148 0.000122 0.198289 RGS8 85397
STRBP\223246_s_at 1.181890 4.008383 0.000157 0.024786 STRBP 55342
KIAA1377\236325_at 1.056305 3.958482 0.000186 -0.090131 KIAA1377 57562
NAPEPLD\242635_s_at 1.170465 3.953391 0.000190 -0.101814 NAPEPLD 222236
EGFEM1P\1558411_at 1.012719 3.919803 0.000212 -0.178704 EGFEM1P 93556
ATF7IP2\219870_at 1.012826 3.880784 0.000242 -0.267603 ATF7IP2 80063
LOC644192\232878_at 1.116624 3.875276 0.000247 -0.280114 LOC644192 644192
LOC100508546\244037_at 1.070067 3.866795 0.000254 -0.299360 LOC100508546 None
RUNX1T1\216831_s_at 1.200380 3.863911 0.000256 -0.305900 RUNX1T1 862
ANGEL1\244685_at -1.170453 -3.844109 0.000274 -0.350735 ANGEL1 23357
SEMA6D\226492_at 0.998893 3.835239 0.000282 -0.370778 SEMA6D 80031
C9orf85\235866_at 1.025476 3.811178 0.000305 -0.425026 C9orf85 138241
HHEX\215933_s_at 0.876758 3.788431 0.000329 -0.476140 HHEX 3087
CLDN15\1552845_at 1.070289 3.787280 0.000330 -0.478723 CLDN15 24146
KLHL24\206551_x_at 0.886946 3.767771 0.000352 -0.522422 KLHL24 54800
NHLRC3\236953_s_at 0.998617 3.766816 0.000353 -0.524558 NHLRC3 387921
ELMOD1\231930_at 1.035745 3.735181 0.000392 -0.595148 ELMOD1 55531
ARHGAP29\1558280_s_at 0.960541 3.717077 0.000416 -0.635398 ARHGAP29 9411
MGC12488\211386_at 0.925050 3.713993 0.000420 -0.642242 MGC12488 84786
GOLM1\1558248_at 1.133981 3.710102 0.000426 -0.650875 GOLM1 51280
MTHFSD\230442_at 1.004106 3.708857 0.000427 -0.653636 MTHFSD 64779
LOC400590\1560992_at 1.137758 3.707311 0.000430 -0.657063 LOC400590 400590
LOC729799\217667_at 1.135771 3.706478 0.000431 -0.658910 LOC729799 729799
ZNF334\220022_at 1.084147 3.684612 0.000463 -0.707301 ZNF334 55713
ZFP82\235176_at 0.920436 3.675480 0.000477 -0.727463 ZFP82 284406
AK2\212172_at 1.019343 3.665957 0.000492 -0.748457 AK2 204
FAM173B\225668_at 1.061109 3.661258 0.000499 -0.758805 FAM173B 134145
ZFP14\232911_at 0.892159 3.647028 0.000523 -0.790097 ZFP14 57677
GTF2H2B\215470_at 0.959996 3.640101 0.000534 -0.805306 GTF2H2B 653238
LOC339047\1558075_at 0.882082 3.629213 0.000554 -0.829174 LOC339047 None
CREBZF\225595_at 0.935700 3.626421 0.000559 -0.835290 CREBZF 58487
LOC286272\232182_at 0.930666 3.622087 0.000566 -0.844775 LOC286272 286272
LOC100507100\230187_s_at 1.060942 3.620916 0.000569 -0.847338 LOC100507100 100507100
GATAD2B\238076_at 0.738321 3.617144 0.000576 -0.855586 GATAD2B 57459
ABCA8\1565778_at 1.045979 3.615368 0.000579 -0.859468 ABCA8 10351
PTPRB\1560105_at 1.022702 3.613720 0.000582 -0.863069 PTPRB 5787
ZNRF1\225962_at 0.900862 3.611676 0.000586 -0.867536 ZNRF1 84937
GPATCH8\212487_at 0.830247 3.561935 0.000687 -0.975750 GPATCH8 23131
TSC2\215624_at 1.116368 3.557750 0.000696 -0.984815 TSC2 7249
EDA\211127_x_at 1.070483 3.555175 0.000702 -0.990390 EDA 1896
LOC441259\1552622_s_at 0.897353 3.551402 0.000711 -0.998555 LOC441259 441259
CASP8AP2\1570001_at 1.032101 3.547817 0.000719 -1.006305 CASP8AP2 9994
CTBP2\215377_at 1.018589 3.534816 0.000749 -1.034381 CTBP2 1488
CACNA1D\243334_at 1.026184 3.529423 0.000762 -1.046007 CACNA1D 776
FBXL13\1556770_a_at 0.945425 3.525738 0.000771 -1.053947 FBXL13 222235
RRP9\204133_at 0.998497 3.518901 0.000788 -1.068664 RRP9 9136
ZNF521\226676_at 0.793170 3.518051 0.000790 -1.070491 ZNF521 25925
ITGB1BP2\219829_at 1.042978 3.514850 0.000798 -1.077375 ITGB1BP2 26548
MAGIX\219965_s_at 1.088856 3.492431 0.000857 -1.125480 MAGIX 79917
FOXJ1\205906_at 1.039986 3.492429 0.000857 -1.125482 FOXJ1 2302
LOC644613\1563426_a_at 1.035299 3.487102 0.000871 -1.136885 LOC644613 None
EFCAB7\226800_at 0.884007 3.483128 0.000882 -1.145385 EFCAB7 84455
LOC100294362\1569078_at 0.984446 3.475146 0.000905 -1.162441 LOC100294362 100294362
GTPBP10\230684_at 1.062098 3.469254 0.000922 -1.175013 GTPBP10 85865
ZNF440\215892_at 0.901321 3.466112 0.000931 -1.181715 ZNF440 126070
WDR52\221103_s_at 0.991380 3.461480 0.000945 -1.191585 WDR52 55779
NCRNA00182\1558515_at 0.838394 3.460203 0.000948 -1.194305 NCRNA00182 None
CXorf56\235442_at 0.995492 3.457467 0.000956 -1.200130 CXorf56 63932
TOP3B\215781_s_at 0.946690 3.447509 0.000987 -1.221308 TOP3B 8940

In [50]:
from utils import Kegg

kegg = Kegg(gene_ids = top_probes.gene_id)

top_probes['pathways_ids'] = top_probes.gene_id.map(kegg.get_pathways_image_links_by_gene)

from IPython.display import HTML
HTML(top_probes.to_html(escape=False))


Out[50]:
log_fc t p_value b gene_name gene_id pathways_ids
probe
DNM3\1558502_s_at 1.375778 4.806308 0.000009 1.945957 DNM3 26052 [hsa04144, hsa04721, hsa04961, hsa05100]
CASC4\1552719_at 1.130824 4.687018 0.000014 1.650206 CASC4 113201 []
ZNF30\232014_at 1.224315 4.522312 0.000026 1.246204 ZNF30 90075 []
LOC389634\1560119_at 1.204393 4.320760 0.000053 0.759543 LOC389634 389634 []
RP9\227852_at 0.921346 4.312707 0.000055 0.740289 RP9 6100 []
SLITRK6\232481_s_at 1.113788 4.223494 0.000075 0.528059 SLITRK6 84189 []
AGTPBP1\204499_at 1.030825 4.211223 0.000078 0.499022 AGTPBP1 23287 []
LOC389906\1558045_a_at 1.219038 4.182897 0.000086 0.432140 LOC389906 389906 []
GFRA2\205722_s_at 1.246157 4.171039 0.000090 0.404202 GFRA2 2675 []
IL1F8\224230_at -1.228913 -4.110874 0.000111 0.263023 IL1F8 None []
KLF8\230986_at 1.056070 4.102448 0.000114 0.243327 KLF8 11279 []
RGS8\234297_at 1.093250 4.083148 0.000122 0.198289 RGS8 85397 []
STRBP\223246_s_at 1.181890 4.008383 0.000157 0.024786 STRBP 55342 []
KIAA1377\236325_at 1.056305 3.958482 0.000186 -0.090131 KIAA1377 57562 []
NAPEPLD\242635_s_at 1.170465 3.953391 0.000190 -0.101814 NAPEPLD 222236 [hsa04723]
EGFEM1P\1558411_at 1.012719 3.919803 0.000212 -0.178704 EGFEM1P 93556 []
ATF7IP2\219870_at 1.012826 3.880784 0.000242 -0.267603 ATF7IP2 80063 []
LOC644192\232878_at 1.116624 3.875276 0.000247 -0.280114 LOC644192 644192 []
LOC100508546\244037_at 1.070067 3.866795 0.000254 -0.299360 LOC100508546 None []
RUNX1T1\216831_s_at 1.200380 3.863911 0.000256 -0.305900 RUNX1T1 862 [hsa05200, hsa05202, hsa05221]
ANGEL1\244685_at -1.170453 -3.844109 0.000274 -0.350735 ANGEL1 23357 []
SEMA6D\226492_at 0.998893 3.835239 0.000282 -0.370778 SEMA6D 80031 [hsa04360]
C9orf85\235866_at 1.025476 3.811178 0.000305 -0.425026 C9orf85 138241 []
HHEX\215933_s_at 0.876758 3.788431 0.000329 -0.476140 HHEX 3087 [hsa04950, hsa05202]
CLDN15\1552845_at 1.070289 3.787280 0.000330 -0.478723 CLDN15 24146 [hsa04514, hsa04530, hsa04670, hsa05160]
KLHL24\206551_x_at 0.886946 3.767771 0.000352 -0.522422 KLHL24 54800 []
NHLRC3\236953_s_at 0.998617 3.766816 0.000353 -0.524558 NHLRC3 387921 []
ELMOD1\231930_at 1.035745 3.735181 0.000392 -0.595148 ELMOD1 55531 []
ARHGAP29\1558280_s_at 0.960541 3.717077 0.000416 -0.635398 ARHGAP29 9411 []
MGC12488\211386_at 0.925050 3.713993 0.000420 -0.642242 MGC12488 84786 []
GOLM1\1558248_at 1.133981 3.710102 0.000426 -0.650875 GOLM1 51280 []
MTHFSD\230442_at 1.004106 3.708857 0.000427 -0.653636 MTHFSD 64779 []
LOC400590\1560992_at 1.137758 3.707311 0.000430 -0.657063 LOC400590 400590 []
LOC729799\217667_at 1.135771 3.706478 0.000431 -0.658910 LOC729799 729799 []
ZNF334\220022_at 1.084147 3.684612 0.000463 -0.707301 ZNF334 55713 []
ZFP82\235176_at 0.920436 3.675480 0.000477 -0.727463 ZFP82 284406 []
AK2\212172_at 1.019343 3.665957 0.000492 -0.748457 AK2 204 [hsa00230, hsa01100, hsa01130]
FAM173B\225668_at 1.061109 3.661258 0.000499 -0.758805 FAM173B 134145 []
ZFP14\232911_at 0.892159 3.647028 0.000523 -0.790097 ZFP14 57677 []
GTF2H2B\215470_at 0.959996 3.640101 0.000534 -0.805306 GTF2H2B 653238 []
LOC339047\1558075_at 0.882082 3.629213 0.000554 -0.829174 LOC339047 None []
CREBZF\225595_at 0.935700 3.626421 0.000559 -0.835290 CREBZF 58487 []
LOC286272\232182_at 0.930666 3.622087 0.000566 -0.844775 LOC286272 286272 []
LOC100507100\230187_s_at 1.060942 3.620916 0.000569 -0.847338 LOC100507100 100507100 []
GATAD2B\238076_at 0.738321 3.617144 0.000576 -0.855586 GATAD2B 57459 []
ABCA8\1565778_at 1.045979 3.615368 0.000579 -0.859468 ABCA8 10351 [hsa02010]
PTPRB\1560105_at 1.022702 3.613720 0.000582 -0.863069 PTPRB 5787 [hsa04520]
ZNRF1\225962_at 0.900862 3.611676 0.000586 -0.867536 ZNRF1 84937 []
GPATCH8\212487_at 0.830247 3.561935 0.000687 -0.975750 GPATCH8 23131 []
TSC2\215624_at 1.116368 3.557750 0.000696 -0.984815 TSC2 7249 [hsa04115, hsa04150, hsa04151, hsa04152, hsa04910, hsa04919, hsa05231]
EDA\211127_x_at 1.070483 3.555175 0.000702 -0.990390 EDA 1896 [hsa04060]
LOC441259\1552622_s_at 0.897353 3.551402 0.000711 -0.998555 LOC441259 441259 []
CASP8AP2\1570001_at 1.032101 3.547817 0.000719 -1.006305 CASP8AP2 9994 []
CTBP2\215377_at 1.018589 3.534816 0.000749 -1.034381 CTBP2 1488 [hsa04310, hsa04330, hsa05200, hsa05220]
CACNA1D\243334_at 1.026184 3.529423 0.000762 -1.046007 CACNA1D 776 [hsa04010, hsa04020, hsa04022, hsa04024, hsa04260, hsa04261, hsa04270, hsa04713, hsa04723, hsa04724, hsa04725, hsa04726, hsa03440, hsa03460]

In [51]:
pathways_df = DataFrame(kegg.get_all_pathways_rows(), columns=['link', 'genes'])

pathways_df = pathways_df[pathways_df.genes > 1].sort(['genes'], ascending = False)

HTML(pathways_df.to_html(escape=False))


Out[51]:
link genes
5 hsa04723 2
32 hsa05200 2
33 hsa05202 2