In [ ]:

    

import pandas as pd
import numpy as np


    

In [ ]:

    

manual_data01_45 = pd.read_excel('manual_count/manual_dataset01_Kr-78_4,5min_incid.xls')
manual_data01_85 = pd.read_excel('manual_count/manual_dataset01_Kr-78_8,5min_incid.xls')
manual_data02 = pd.read_excel('manual_count/manual_dataset02.xls')

auto_data01_45 = pd.read_csv('auto_count/auto_dataset01_Kr-78_4,5min_incid.txt')
auto_data01_85 = pd.read_csv('auto_count/auto_dataset01_Kr-78_8,5min_incid.txt')
auto_data02 = pd.read_csv('auto_count/auto_dataset02.txt')


    

In [ ]:

    

auto = np.asarray(list(auto_data01_45[(auto_data01_45['initial_radius'] == 20) &
                                      (auto_data01_45['delta_radius'] == 2)].auto_count) +
                  list(auto_data01_85[(auto_data01_85['initial_radius'] == 5) &
                                      (auto_data01_85['delta_radius'] == 4)].auto_count) +
                  list(auto_data02[(auto_data02['initial_radius'] == 10) &
                                      (auto_data02['delta_radius'] == 8)].auto_count))


    

In [ ]:

    

manual = np.asarray(list(manual_data01_45.manual_count) + 
                    list(manual_data01_85.manual_count) +
                    list(manual_data02.manual_count))


    

In [ ]:

    

smp_mean = np.mean(auto / manual)
smp_stdmean = np.std(auto / manual) / np.sqrt(len(auto))


    

In [ ]:

    

print('mu = ', smp_mean, ', sigma = ', smp_stdmean)


    

In [ ]: