This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
def get_jacobian(net, x, noutputs): | |
x = x.squeeze() | |
n = x.size()[0] | |
x = x.repeat(noutputs, 1) | |
x.requires_grad_(True) | |
y = net(x) | |
y.backward(torch.eye(noutputs)) | |
return x.grad.data |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
from scipy import stats | |
def pearsonr_ci(x,y,alpha=0.05): | |
''' calculate Pearson correlation along with the confidence interval using scipy and numpy | |
Parameters | |
---------- | |
x, y : iterable object such as a list or np.array | |
Input for correlation calculation |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import unittest | |
import numpy as np | |
import numpy.linalg as linalg | |
def makeUnit(x): | |
"""Normalize entire input to norm 1. Not what you want for 2D arrays!""" | |
return x / linalg.norm(x) |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
from __future__ import division | |
import numpy as np | |
import pandas as pd | |
import random | |
def sample(data): | |
sample = [random.choice(data) for _ in xrange(len(data))] | |
return sample | |
def bootstrap_t_test(treatment, control, nboot = 1000, direction = "less"): |