% Quantification of blot data done with KF's procedure; last poin is the background!!!!

time = [0 10 20 30 45 60 90 105 2*60 2.5*60 3*60 4*60 6*60 8*60 10*60 12*60 15*60 18*60 24*60];

% STAT 1 protein nuclear concentration - based on Shaofei blots
stat1_n_data = [0.0517 0.0748 0.1072 0.4866 0.1902 0.3916 0.2654 0.1868 0.2440 0.3383 0.3055 0.3561 0.5573 0.5945 0.7978 0.9731 0.8738 0.8535 1.0000 0.004];

% STAT 1 protein cytoplasmic concentration - based on Shaofei blots
%stat1c_data = [0.3451 0.3907 0.5272 0.3140 0.4404 0.4229 0.4216 0.4772 0.3508 0.5472 0.5984 0.6481 0.6893 0.8575 0.8999 1.0000 0.9680 0.9969 0.9557 ???];

stat1c_data = [0.2989 0.3445 0.4330 0.3055 0.4201 0.4355 0.3780 0.4513 0.3325 0.5036 0.5803 0.5963 0.6376 0.8557 0.7842 0.8766 0.9131 0.9474 1.0000 0.0041];
new_data = [60.5216   74.6588   89.4471   49.8039 75.4824   68.7647   64.9765   75.5059 63.7608   90.1020  103.0000  110.3686 119.8980  146.3176  140.5255  145.8196 147.5961  149.0902  148.3294    8.7529];
%phosphorylated STAT1 in cytoplasm

%[0.3521 0.6714 0.8989 0.8863 0.9519 1.0000 0.9655 0.9811 0.8175 0.660 0.9085 0.8225 0.8167 0.9382 0.7489 0.8067 0.7710 0.7534 0.6874 0.2527 ???];

stat1pc_data = [0.3516 0.6503 0.8946 0.8763 0.9316 1.0000 0.9568 0.9562 0.8197 0.6531 0.8844 0.8158 0.7982 0.8991 0.7460 0.7898 0.7247 0.7220 0.6828 0.4188];

%phosphorylated STAT1 in nucleus

stat1pn_data = [0.2889 0.6257 0.6253 1.0000 0.8326 0.9500 0.7531 0.7023 0.7444 0.8395 0.8750 0.8983 0.8694 0.7886 0.7574 0.7688 0.6630 0.7187 0.4934 0.2950];

% phosphorylated STAT2 in cytoplasm
stat2pc_time = [0 10 20 30 60 2*60 5*60];
stat2pc_data = [0.6420 1.0000 0.8370 0.8769 0.9315 0.8303 0.9773 0.6425];

% total STAT2 in cytoplasm
stat2c_data = [1.0000    0.9067    0.8106    0.7353    0.6927    0.7052    0.8491    0.7965    0.3983];

% phosphorylated STAT2 in nucleus
stat2pn_data = [0.4584    0.9540    1.0000    0.8651    0.7955    0.7284    0.7161    0.4646    0.4584];

% total STAT2 in nucleus
stat2n_data = [0.7607    0.9021    0.8813    0.8631    0.8134    0.7975    0.8566    1.0000    0.4186];

%Homodimer of phosphorylated STAT1 (time points without 60 min!!!)
time_hom = [0 10 20 30 45 90 105 2*60 2.5*60 3*60 4*60 6*60 8*60 10*60 12*60 15*60 18*60 24*60];
s1_s1_data =[0.1270 0.4135 0.8067 0.7401 1.0000 0.7379 0.6262 0.4703 0.4073 0.3297  0.3008 0.3416 0.2040 0.1487 0.0988 0.1934 0.2838 0.1478 0.1280];

%Heterodimer of phosphorylated STAT1 and STAT2
time_het = [0 10 20 30 60 2*60 5*60 10*60];
s1_s2_data =[0.4442  1.0000  0.7818  0.7462  0.7129  0.6823  0.7501  0.5029  0.4442];

% STAT1|STAT2|IRF9 complex
isgf3_data = [0.4912 0.9177 0.9061 0.7784 0.8686 0.8832 1.0000 0.9196 0.4912];

%Heterodimer of phosphorylated STAT1 and STAT3

s1_s3_data = [0.1618 0.4043 0.7182 0.7068 1.0000 0.7173 0.6606 0.4591 0.4183 0.3413 0.2967 0.3406 0.1413 0.1258 0.0886 0.2363 0.3092 0.1723 0.1294];

%Total IRF1 in nucleus
irf1n_data = [0.2705 0.3066 0.3097 0.3600 0.2687 0.4751 0.8199 0.9917 1.0000 0.9926 0.7386 0.4790 0.6178 0.5879 0.3633 0.3435 0.5199 0.3934 0.3206 0.1276];

%Total IRF1 in cytoplasm??
irf1c_data = [0.6416 0.5885 0.5904 0.6298 0.6616 0.5999 0.4332 0.3974 0.4279 0.6291 1.0000 0.9868 0.9364 0.8392 0.8285 0.8480 0.8085 0.8509 0.7278 0.4120];


%data quantified in Galveston


% time 0'	10'	20'	30'	45'	60'	90'	105'	120'	2.5h	3h	4h	6h	8h	12h	14h	16h	18h	24h

irf1_cyt_g = [1.87 1.31 1.9 2.38 2.13	1.94 0.85 0.76 1.69 5.03 8.84 9.95 6.74	5.51 4.19 4.55 5.04	4.97 3.14];
irf1_nuc_g = [12.12	8.02 10.84 12.29 10.41 19.47 38.7 51.96 53.94 61.09 36.5 17.17 22.49 16.59 9.69 8.69 18.07 9.53 7.47];


stat1p_c_g=[0.05 12.95 26.63 25.91 26.08 30.11 28.92 28.98 20.94 21.04 26.43 25.53 26.35 38.38 24.12 34.36 22.03 23.1 16.88];
stat1_c_g= [11.01 13.28 16.53 12.21 14.89 15.67 15.38 16.34 14.39 18.76 19.97 22.41 23.77 39.68 39.62 47.55 42.11 48.31 58];

stat1p_n_g = [0.39 9.44 8.72 22.71 15.02 23.78 15.2 14.19 15.85 19.13 17.43 15.05 26.42 14.86 13.25 19.19 15.89 25.24 16.18];
stat1_n_g = [3.02 4.25 4.67 11.22 6.22 10.67 8.41 7.4 8.25 10.25 9.35 9.48 12.01 17.05 23.06 31.07 22.2 22.56 30.95];




lysis = [0.5704    0.4987    0.6444    0.8108    0.9034    1.1153    1.2065    1.2467    0.4987];
time_lysis = [0 1 2 5 10 15 20 25];