% Huntington 2001
% PND11 male and female rats administered malathion for 11 days (5, 50, or 150 mg/kg)

% mean BW for males and females approximately 25 g


rat_parameters_postnatal 	        % rat parameters
TSTOP=528;      			  % Simulation duration
BW0=0.025;   	 			   % kg
AGE0=11.0; 



%11 days of dosing
DT = xlsRead("Huntington_2001_dosing.xls", "dosing", "A1:A23");  % Dosing Times (hr)
DORAL=xlsRead("Huntington_2001_dosing.xls", "dosing", "D1:D23"); % Dose in ug/kg



prepare @clear @all

load @format=ascii @separator=comma @ignoreheader @file=Huntington_2001_repeatdose_1a.csv Huntington_1
tdata1=Huntington_1

load @format=ascii @separator=comma @ignoreheader @file=Huntington_2001_repeatdose_2a.csv Huntington_2
tdata2=Huntington_2

load @format=ascii @separator=comma @ignoreheader @file=Huntington_2001_repeatdose_3a.csv Huntington_3
tdata3=Huntington_3


start@nocallback

%5 mg/kg
%plot(_time, _blbes, "-", tdata1(:,1), tdata1(:,2),tdata1(:,1), tdata1(:,5), "*", 'Huntingdon_plasma_repeatdose_1.aps')
%plot(_time, _rbcce, "-", tdata1(:,1), tdata1(:,3), tdata1(:,1), tdata1(:,6),"*", 'Huntingdon_RBC_repeatdose_1.aps')
%plot(_time, _bce, "-", tdata1(:,1), tdata1(:,4), tdata1(:,1), tdata1(:,7), "*", 'Huntingdon_brain_repeatdose_1.aps')


%50 mg/kg
%plot(_time, _blbes, "-", tdata2(:,1), tdata2(:,2),tdata2(:,1), tdata2(:,5), "*", 'Huntingdon_plasma_repeatdose_2.aps')
%plot(_time, _rbcce, "-", tdata2(:,1), tdata2(:,3),tdata2(:,1), tdata2(:,6), "*", 'Huntingdon_RBC_repeatdose_2.aps')
%plot(_time, _bce, "-", tdata2(:,1), tdata2(:,4), tdata2(:,1), tdata2(:,7),"*", 'Huntingdon_brain_repeatdose_2.aps')


%150 mg/kg
plot(_time, _blbes, "-", tdata3(:,1), tdata3(:,2),tdata3(:,1), tdata3(:,5), "*", 'Huntingdon_plasma_repeatdose_3.aps')
plot(_time, _rbcce, "-", tdata3(:,1), tdata3(:,3),tdata3(:,1), tdata3(:,6), "*", 'Huntingdon_RBC_repeatdose_3.aps')
plot(_time, _bce, "-", tdata3(:,1), tdata3(:,4),tdata3(:,1), tdata3(:,7), "*", 'Huntingdon_brain_repeatdose_3.aps')



