コード例 #1
0
                    [['soma', '0.5', 'kap', 'gbar']],
                    [['soma', '0.5', 'na8st', 'gbar']]
                 ]
errfun = 'rms'
fitfun = 'get_v'
fitnessweights = [1]
model_dir = '../../../model/cells/dapmodel0.json'
mechanism_dir = '../../../model/channels/schmidthieber'
data_dir = '../../../data/2015_08_26b/rampIV/3.0(nA).csv'

fitter = HodgkinHuxleyFitter(variable_keys, errfun, fitfun, fitnessweights,
                 model_dir, mechanism_dir, data_dir, simulation_params={'celsius': 35})

# create cell
candidate = np.ones(len(variable_keys))  # gbars should be 1
fitter.update_cell(candidate)

# extract parameter
channel_list = get_channel_list(fitter.cell, 'soma')
ion_list = get_ionlist(channel_list)

v_exp = fitter.data.v.values
t_exp = fitter.data.t.values
i_exp = fitter.data.i.values
dt_exp = t_exp[1] - t_exp[0]

dvdt_exp = np.concatenate((np.array([(v_exp[1]-v_exp[0])/dt_exp]), np.diff(v_exp)/dt_exp))
celsius = fitter.simulation_params['celsius']

# convert units
dvdt_sc, i_inj_sc, _, _, cell_area = convert_units(fitter.cell.soma.L, fitter.cell.soma.diam,