Commit 775043e5 authored by Hao-Ting Wang's avatar Hao-Ting Wang

add correlation value to the heatmap

parent 224c3b12
......@@ -218,8 +218,8 @@
}
],
"source": [
"# find the best variates among the 6 selected\n",
"# by calculating the FWE-corrected p-value based on FDR\n",
"# find the best mode\n",
"# by calculating the FWE-corrected p-value based on false discovery rate\n",
"permute = 1000\n",
"alpha = 0.05\n",
"\n",
......@@ -232,10 +232,8 @@
" perm_model = SCCA(C=C, penalty=['l1', 'l1'], n_component=n_mod, verbose=False)\n",
" perm_model.fit(X_clean, cur_y)\n",
" perm_d = perm_model._cancorr\n",
"\n",
" permD[i, :] = perm_d\n",
"\n",
"\n",
"# calculate the FWE-corrected p value\n",
"p_j = (1 + np.sum(d < np.repeat(permD[1:, 0:1], n_mod, axis=1),0)) / float(permute)\n",
"sig = p_j < alpha\n",
......@@ -282,5 +280,5 @@
}
},
"nbformat": 4,
"nbformat_minor": 2
"nbformat_minor": 4
}
This diff is collapsed.
This diff is collapsed.
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment