filter_type = 'FIR' taps = 25 window = 'blackmanharris' transd = True mains = 50 dirList=os.listdir(data_path) for fname in dirList: file = data_path+fname f = h5py.File(file) group_name = f.attrs['group_name'] number_in_group = f.attrs['number_in_group'] species = f.attrs['species'] location = f.attrs['location'] subject = session.query(db.Subject).\ filter_by(species=species, group_name=group_name, number_in_group=number_in_group).first() if not subject: subject = db.Subject(species=species, group_name=group_name, number_in_group=number_in_group) session.add(subject) session.commit() print file conditions = [(v,t,e,s,rem) for v in visits for t in tasks for e in eyes for s in sensors for rem in remicas] for visit, task_type, eye, sensor_type, rem in conditions: base = str(visit)+'/'+task_type+'/'+eye+'/'+sensor_type+'/'+rem base_filtered = base+'/filter_'+filter_type+'_'+str(taps)+'_'+window #If this particular set of conditions doesn't exist for this subject, just continue to the next set of conditions try: f[base_filtered] except KeyError:
#robjects.r("data<-read.table('%s')"%(base+'LangleyRtable')) #robjects.r("attach(data)") for d in dependents: print d if d == 'Has_Parent': independents = ['Age', 'Gender', 'Relationship_with_Parent', 'Heard_Through_Medium'] else: independents = ['Age', 'Gender', 'Relationship_with_Parent', 'Heard_Through_Medium', 'Same_Age_as_Parent',\ 'Same_City_as_Parent', 'Same_Country_as_Parent', 'Same_Gender_as_Parent', 'Same_Relationship_to_Parent_as_They_Had_to_Their_Parent',\ 'Heard_Through_Same_Medium_as_Parent', 'Has_Parent'] for i in independents: print i ax = plt.subplot(1,1,1) handles = {} factors = session.query(db.LangleyParticipant).values(getattr(db.LangleyParticipant, i)) factors = unique([q for q in factors]).flatten() print factors factors = [q for q in factors if (q!=None and q!='unknown' and q!='0')] print factors n_factors = len(factors) if n_factors<2: continue fs = [] f2s = [] Ds = [] p_krusks = [] p_KSs = [] skews = [] medians = []