示例#1
0
            print "Caught exception when trying to add values from  %s to statistics collector" % filename
            continue


#Serialization
if pickle_location != '':
    for aminoacid in aminoacids:
        statistics_collector_from_pdb.restype_av_scores[aminoacid].pickle_res_type_average_scores(pickle_location+aminoacid+'.txt')


#deserialize archived files
if archive_listfile != '' and pdb_listfile == '':
    if add_arch == True:
        for archive in archive_list:
            f = file(archive, 'rb')
            statistics_collector_from_pdb.add_archived_data( cPickle.load(f) )
            f.close()
            statistics_collector_from_archive = statistics_collector_from_pdb
    elif add_arch ==False:
        for archive in archive_list:
            f = file(archive, 'rb')
            statistics_collector_from_archive.add_archived_data( cPickle.load(f) )
            f.close()



#Deserialization of one particular file (faster)
#aminoacids = ['GLU']
# = '/home/rahel/uni/wise1314/python/pdbstats/pickle11k/GLU.txt'
#statistics_collector_from_archive.add_archived_data(cPickle.load(file(archive, 'rb')))
for filename in all_file_names:
    pe_instance = PoseEnergies()
    pe_instance.loadFile(filename)
    statistics_collector_from_pdb.add_pose_energies(pe_instance)


for aminoacid in aminoacids:
    statistics_collector_from_pdb.restype_av_scores[aminoacid].pickle_res_type_average_scores('pickle-test/'+aminoacid+'.txt')


for aminoacid in aminoacids:
    filename = 'pickle-test/'+ aminoacid + '.txt'
    f = file(filename, 'rb')
    loaded_object = cPickle.load(f)
    statistics_collector_from_archive.add_archived_data( loaded_object)
    #archived_res_type_average_scores[aminoacid] = loaded_object
    f.close()

if not statistics_collector_from_pdb.calculate_averages_and_stddevs('GLU', 'faketerm1') == statistics_collector_from_archive.calculate_averages_and_stddevs('GLU', 'faketerm1'):
	error3_seen = True
elif not statistics_collector_from_pdb.calculate_averages_and_stddevs('ASP', 'faketerm1') == statistics_collector_from_archive.calculate_averages_and_stddevs('ASP', 'faketerm1'):
	error3_seen = True
elif not statistics_collector_from_pdb.calculate_averages_and_stddevs('MET', 'faketerm2') == statistics_collector_from_archive.calculate_averages_and_stddevs('MET', 'faketerm2'):
	error3_seen = True
#elif not statistics_collector_from_pdb.calculate_averages_and_stddevs('ALA', 'faketerm4') == statistics_collector_from_archive.calculate_averages_and_stddevs('ALA', 'faketerm4'):
#	error3_seen = True