Ejemplo n.º 1
0
def search_one_asymmetry(**kwargs):
    sess = PINGDataSession()  # username/password stored in an env variable for me.
    sess.login()

    return sess.regress_multistep(X='MRI_cort_area_ctx_supramarginal_AI',
                                  Y='MRI_cort_area_ctx_rostralmiddlefrontal_AI',
                                  **kwargs)
Ejemplo n.º 2
0
def search_all_pairwise(plot=True, **kwargs):
    """For each pair of variables, look for a significant regression slope."""

    sess = PINGDataSession()  # username/password stored in an env variable for me.
    sess.login()

    all_data = compute_all_asymmetries(prefix=['MRI_cort_area', 'MRI_cort_thick',
                                               'MRI_subcort_vol', 'DTI_fiber_vol'])
    results = []

    keys = list(set(all_data.keys()) - set(('SubjID',)))
    for ki in range(len(keys)):
        key1 = keys[ki]
        if skip_key(key1):
            continue

        for ii in range(ki + 1, len(keys)):
            key2 = keys[ii]
            if skip_key(key2) or skip_pairing(key1, key2):
                continue

            try:
                result = sess.regress_multistep(key1, key2, plot=plot, **kwargs)
                results.append(result)
            except Exception as e:
                print("Exception: %s" % str(e))  # force print
            else:
                if plot:
                    if result is None:
                        plt.close()
                        plt.close()
                    else:
                        plt.show()