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 = []
