matlab error

??? Index exceeds matrix dimensions.

Error in ==> motionEstTSS at 154 imwrite((Best_Predicted_Frame_mov(1,c).cdata),filename,'jpg');

Any ideas?

Line 153: filename = strcat('Best_Predicted_Frame',num2str(c),'.jpg');

Line 154: imwrite((Best_Predicted_Frame_mov(1,c).cdata),filename,'jpg');

[Yoav Livneh]

What is the variable c? It appears the array Best_Predicted_Frame_mov is smaller than the value of c.

Answer 1

[Fangjun Jiang]

The following will cause the same error. Check the size of your variable Best_Predicted_Frame_mov and the value of variable c.

a=rand(10,1);
b=a(11,1);

Answer 2

[Thandava Moorthy]

Error using fz_costo_new (line 19)

Error due to multiple causes.

Error in lsqnonlin (line 206)

initVals.F = feval (funfcn {3}, xCurrent, varargin {:});

Error in parameterestim_sana (line 108)

[bestpar, resnorm, residual, exitflag, output, lmd, jacobian] = lsqnonlin (@fz_costo_new, startpar, [0 0 0],

[5 5 5], options); - Show complete stack trace

Caused by:

Error using fz_costo_new (line 19)

Invalid setting in 'MCBS_project_than / Constant2' for parameter 'Value'. - Show complete stack trace

Error using fz_costo_new (line 19)

Error evaluating parameter 'Value' in 'MCBS_project_than / Constant2' - Show complete stack trace

Error using fz_costo_new (line 19)

Unrecognized function or variable 'K1'. - Show complete stack trace

Error using fz_costo_new (line 19)

Variable 'K1' has been deleted from base workspace.

Suggested Actions:

• Undo the deletion. - Fix

• Load a file into base workspace. - Fix

• Create a new variable. - Fix

- Show complete stack trace

Error using fz_costo_new (line 19)

Invalid setting in 'MCBS_project_than / Constant3' for parameter 'Value'. - Show complete stack trace

Error using fz_costo_new (line 19)

Error evaluating parameter 'Value' in 'MCBS_project_than / Constant3' - Show complete stack trace

Error using fz_costo_new (line 19)

Unrecognized function or variable 'K2'. - Show complete stack trace

Error using fz_costo_new (line 19)

Variable 'K2' has been deleted from base workspace.

Suggested Actions:

• Undo the deletion. - Fix

• Load a file into base workspace. - Fix

• Create a new variable. - Fix

- Show complete stack trace

Error using fz_costo_new (line 19)

Invalid setting in 'MCBS_project_than / Ra_INC_b // V' for parameter 'Value'. - Show complete stack trace

Error using fz_costo_new (line 19)

Error evaluating parameter 'Value' in 'MCBS_project_than / Ra_INC_b // V' - Show complete stack trace

Error using fz_costo_new (line 19)

Unrecognized function or variable 'RaINCb'. - Show complete stack trace

Error using fz_costo_new (line 19)

Variable 'RaINCb' has been deleted from base workspace.

Suggested Actions:

• Undo the deletion. - Fix

• Load a file into base workspace. - Fix

• Create a new variable. - Fix

- Show complete stack trace

Failure in initial objective function evaluation. LSQNONLIN cannot continue.

211 rethrow (userFcn_ME)

K >>

this the error im facing right now

my parameter estimation file is below,

clear

clear global

clc

close all

%format long

global SRex Km Ke V IL6b RaIL6 Tv t_start t_end dataIL6GLP1 t

% -------------------------------------------------------------------------

% NB: GLYCEMIA has to be expressed in (mmol/l

% INSULINEMIA has to be expressed in (mU/l)

% -------------------------------------------------------------------------

scelta=menu('Select case','HC1 (Muscelli)','HC2 (Nauck)');

switch scelta

case 1

% -------------------------------------------------------------------------

% OGTT MUSCELLI DATA NORMO

% E. Muscelli, A. Mari, A. Natali, B.D. Astiarraga, S. Camastra, S. Frascerra, J.J. Holst, E. Ferrannini, Impact of incretin hormones on ?-cell function in subjects with normal or impaired glucose tolerance, Am. J. Physiol. Endocrinol. Metab. 291 (2006) E1144-E1150.

t_start=1;

t=[0 10 20 30 40 60 90 120 150 180];

dataGLP1out=[-0.0611404034260573;0.153761026536856;0.381887441294696;0.577225771343123;0.785935420734692;1.02732340392451;1.28878665794757;1.46511064263731;1.60857337534354;1.75203610804977;1.79782967573179;1.72437887843715;1.63132839759097;1.49215509632929;1.40577198333615;1.31946203761713;1.21963443800674;1.14002844377781;1.05353557987351;0.927440928858545;0.867507785455261;0.788267627596867]';

dataIL6out=[-0.0611404034260573;0.153761026536856;0.381887441294696;0.577225771343123;0.785935420734692;1.02732340392451;1.28878665794757;1.46511064263731;1.60857337534354;1.75203610804977;1.79782967573179;1.72437887843715;1.63132839759097;1.49215509632929;1.40577198333615;1.31946203761713;1.21963443800674;1.14002844377781;1.05353557987351;0.927440928858545;0.867507785455261;0.788267627596867]'*(1/6); %conversion from pmol/l to mU/l

% GIP+GLP1

% incretinaOGTT=[145.716 270.0798 339.5178 329.7588 322.815 318.4284 272.0028 280.0878 284.4702 254.0166]';

SRex = 0.050;

Km= 0.009;

Ke = 0.059;

V = 14;

IL6b =1.7;

RaIL6 = Ke*IL6b*V;

Tv =67;

%t_start= 1;

t_end = 155;

%RaINCb = 217;

%K1= 1;

%K2=0.1;

%t_end=180;

% ----------------------------------------------------------------

case 2

%OGTT NAUCK DATA HEALTHY

% M.A. Nauck, A. El-Ouaghlidi, B. Gabrys, K. Hücking, J.J Holst, C.F Deacon, B. Gallwitz, W.E. Schmidt, J.J. Meier, Secretion of incretin hormones (GIP and GLP-1) and incretin effect after oral glucose in first-degree relatives of patients with type 2 diabetes, Regul Pept. 122 (2004) 209-217.

% t=[0 15 30 45 60 75 90 120 150 180 210 240]';

% glicemiaOGTT=[5.54 8.04 11 11.54 10.64 8.73 7.27 5.73 5.14 4.64 4.63 4.68]';

%insulinemiaOGTT=[45.45 193.9 478.78 654.5 706.06 472.73 318.18 169.7 109.1 60.61 36.36 26.67]'*(1/6); %conversion from pmol/l to mU/l

%GIP+GLP1

%incretinaOGTT=[92.2602 298.347 418.8954 423.4032 431.433 402.3756 370.2276 289.7586 221.4258 178.7952 142.5594 104.3616]';

%Gb = 5.5;

%Ib = 7.6;

%Incb=90;

% k7 ESTIMATED FROM IIVG

%k7=0.158;

%tfinal=240;

% -------------------------------------------------------------------------

end

% -------------------------------------------------------------------------

dataIL6GLP1 = [dataGLP1out dataIL6out];

%FIXED PARAMETERS SIMULATOR

%D = 417; %quantity in mmol equivalent to 75g

%kabs=0.22;

%f=0.90;

% HGB_b=0.76875; % hepatic balance

V=14;

%M_ipog=0.03;

%alfa=1.0;

%p=2;

% K1=(p/(p+1))*HGB_b/((G_b)*V);

% K2=(1/(p+1))*HGB_b/(I_b*V);

%K3=0.0718; % different from zero only when G>10

%K4=0.717; % different from zero only when G>10

%K6=0.1;

%K9=0.1;

%Ra_INC_b=K6*INC_b*V;

%K5 = 6.91;

% ------------------------------------

% initial values for the parameters to be estimated

%K_0=0.013;

K1_0 = 1;

K2_0 = 0.1;

RaINCb_0=217;

% BETA_0=1;

% K5_0 = 25;

startpar = [K1_0 K2_0 RaINCb_0];

% -------------------------------------

options=optimset('Display','Iter','TolFun',1E-20,'TolX',1E-7,'MaxFunEvals',3000,'MaxIter',500)';

[bestpar, resnorm, residual, exitflag, output, lmd, jacobian] = lsqnonlin(@fz_costo_new, startpar, [0 0 0], [5 5 5], options);

%bestpar, resnorm, residual, exitflag, output, lmd, jacobian

%K_0=bestpar(1)

K1_0=bestpar(1);

K2_0=bestpar(2);

RaINCb_0=bestpar(3);

% BETA_0=bestpar(4)

% K5_0=bestpar(5)

% beta=-K7*(G_b^1.3)-K8*INC_b+K9*I_b;

% kgri = Kmax;

% E error covariance matrix

[M, n] = size (residual);

E = inv (jacobian 'jacobian *);

% -------------------------------------

% standard deviation of the parameters

% Std_K = sqrt (S (1,1));

std_K1 = sqrt (S (1,1));

std_K2 = sqrt (E (2,2));

std_RaINCb = sqrt (S (3,3));

% std_BETA = sqrt (E (4,4));

% std_K5 = sqrt (E (5.5));

% -------------------------------------

% percent coefficient of variation of the estimated parameters

CV_K% = (std_K / K) * 100

CV_K1 = (std_Gama / Gama) * 100;

CV_K2 = (std_K8 / K8) * 100;

CV_RaINCb = (std_M / M) * 100;

% CV_BETA = (std_BETA / BETA) * 100

% CV_K5 = (std_K5 / K5) * 100

% -------------------------------------

Figures (1)

sim ('MCBS_project_than', t)% replace 'Nuovo_OGTT_Ela' with the name of your Simulink model

plot (t, GLP1_data_VICT, 'ro')

xlabel ('time (min)')

ylabel ('GLP1 (mmol / L)')

hold on

plot (tout, GLP1out, 'k')

figures (2)

plot (t, IL6_data_VICT, 'ro')

xlabel ('time (min)')

ylabel ('insulin (mU / L)')

hold on

plot (tout, IL6mout, 'k')

% figures (3)

% plot (t, incretinaOGTT, 'ro')

% xlabel ('time (min)')

% ylabel ('incretin (ng / L)')

% hold on

% plot (tout, INCout, 'k')

RMSEogtt = sqrt (resnorm / (3 * (length (GLP1_data_VICT))));

Here is my fz_costo file is below,

function matrix_scarti = fz_costo_new (parameters)

global SRex Km Ke V IL6b RaIL6 Tv t_start t_end RaINCb K1 K2 dataIL6GLP1 t

K1 = parameters (1);

RaINCb = parameters (2);

K2 = parameters (3);

% = SREX parameters (4);

% K5 = parameters (5);

[A, b] = size (dataIL6GLP1);

FSD = [0.015 * ones (a, 1) 0.04 * ones (a, 1) 0.055 * ones (a, 1)];

% Kgri = Kmax;

sum = 0;

sim ('MCBS_project_than.slx', t)% replace 'Nuovo_OGTT_Ela' with the name of your Simulink model

GLP1out = [145; 144.983612661427; 144.901703747925; 144.492852854599; 142.465839448341; 132.749472474337; 108.241140537250; 92.4628797737528; 92.4628797737525; 92.4628797737519; 68.2279536519477; 50.7472523011548; 33.8415846009488; 21.2687903843176; 13.9937064557888; 10.9994796749988; 11.4180931459149; 14.5415293407342; 20.0193104618235; 27.7738825859954; 36.3345009549500; 41.2479481722555; 46.0212817210750; 52.2773763959316 ; 59.6286225521930; 64.1960138400612; 68.6209578472389; 74.0524490606560; 80.0781769637733; 84.0553183987940; 87.9069951033905; 92.5303362329276; 97.6227906453387; 101.109030739800; 104.483442279190; 108.446350178841; 112.749123768602; 116.491781335296; 119.409022293921; 122.035845387021; 125.005902690712; 128.471252013136; 131.842494913110; 134.516072867313; 136.685604332842; 138.910242173325; 141.561107134264; 144.450000817259; 146.439495297789 ; 148.364989754182; 150.624951232535; 153.081271789160; 153.445087161108; 153.445087161109; 153.445087161110; 153.822978968231; 153.921610327064; 153.592396612185; 152.718148326212; 151.227750654245; 149.064818839237; 146.122637858476; 142.186503719168; 136.910373048788; 130.549331288480; 126.518596337588; 122.615350598220; 117.913182108577; 112.782992639757; 109.320903978364; 105.968704803433; 101.982898623359; 97.6112882541214; 94.5740347954464; 91.6348560242062; 88.2125133162925; 84.5167634657075; 81.2900817399232; 78.7524321119111; 76.4596643837533; 73.8826534798798; 70.8969835206899; 67.9974252983936; 65.6838074788047; 63.7934853451588; 61.8567033201319; 59.5628704131568; 57.0761073888142; 54.9351938536194; 53.2769056118132; 51.7765564555988; 50.0632992309416; 48.0583232652025; 46.1210455667978; 44.5974632964747; 43.3612361717790; 42.0830121359225; 40.5512222632452; 38.8857669344694; 37.7524180467645; 36.6555526168083; 35.3567616263242; 33.9361517534880; 32.7022963788562; 31.7476465401317; 30.8939737425392; 29.9254512372861; 28.7866295194129; 27.6742771941943; 26.7947847680259; 26.0863879093472; 25.3627708468259; 24.4961377658798; 23.5451718016503; 22.8920459490068; 22.2603053623720; 21.5187156121648; 20.7138600097065; 20.0138015765728; 19.4662088378172; 18.9738397252296; 18.4180031299590; 17.7697457073534; 17.1397822461560; 16.6386419603784; 16.2311208381292; 15.8145064543026; 15.3185709285879; 14.7786937143971; 14.4052210112797; 14.0442157938336; 13.6212145853906; 13.1624546925866; 12.7629878790706; 12.4491078253367; 12.1672717959376; 11.8499870258959; 11.4804001133382; 11.1214581254571; 10.8347704525427; 10.6010354554832; 10.3629268993920; 10.1338337719011];6742771941943; 26.7947847680259; 26.0863879093472; 25.3627708468259; 24.4961377658798; 23.5451718016503; 22.8920459490068; 22.2603053623720; 21.5187156121648; 20.7138600097065; 20.0138015765728; 19.4662088378172; 18.9738397252296; 18.4180031299590; 17.7697457073534; 17.1397822461560; 16.6386419603784; 16.2311208381292; 15.8145064543026; 15.3185709285879; 14.7786937143971; 14.4052210112797; 14.0442157938336; 13.6212145853906; 13.1624546925866; 12.7629878790706; 12.4491078253367; 12.1672717959376; 11.8499870258959; 11.4804001133382; 11.1214581254571; 10.8347704525427; 10.6010354554832; 10.3629268993920; 10.1338337719011];6742771941943; 26.7947847680259; 26.0863879093472; 25.3627708468259; 24.4961377658798; 23.5451718016503; 22.8920459490068; 22.2603053623720; 21.5187156121648; 20.7138600097065; 20.0138015765728; 19.4662088378172; 18.9738397252296; 18.4180031299590; 17.7697457073534; 17.1397822461560; 16.6386419603784; 16.2311208381292; 15.8145064543026; 15.3185709285879; 14.7786937143971; 14.4052210112797; 14.0442157938336; 13.6212145853906; 13.1624546925866; 12.7629878790706; 12.4491078253367; 12.1672717959376; 11.8499870258959; 11.4804001133382; 11.1214581254571; 10.8347704525427; 10.6010354554832; 10.3629268993920; 10.1338337719011];6386419603784; 16.2311208381292; 15.8145064543026; 15.3185709285879; 14.7786937143971; 14.4052210112797; 14.0442157938336; 13.6212145853906; 13.1624546925866; 12.7629878790706; 12.4491078253367; 12.1672717959376; 11.8499870258959; 11.4804001133382; 11.1214581254571; 10.8347704525427; 10.6010354554832; 10.3629268993920; 10.1338337719011];6386419603784; 16.2311208381292; 15.8145064543026; 15.3185709285879; 14.7786937143971; 14.4052210112797; 14.0442157938336; 13.6212145853906; 13.1624546925866; 12.7629878790706; 12.4491078253367; 12.1672717959376; 11.8499870258959; 11.4804001133382; 11.1214581254571; 10.8347704525427; 10.6010354554832; 10.3629268993920; 10.1338337719011];

IL6mout = [0; 1.26182921468215e-09; 4.54104709955408e-08; 1.21015748531682e-06; 3.03879713374870e-05; 0.000730547588521312; 0.00713914690745592; 0.0155321449711194; 0.0155321449711196; 0.0155321449711200; 0.544992357138957; 1.76922838426363, 4.13961280420290, 7.62604889160560; 11.9673090440306; 17.2718386116334; 23.7583415934445; 31.8961621135888 ; 42.4796557306135; 56.3668716061309; 71.5109446323153; 80.1783859350137; 88.5982459233266; 99.6348742443182; 112.608101921523; 120.660872052794; 128.464858317559; 138.046626201386; 148.681061256891; 155.692816064709; 162.485543392943; 170.641678902605; 179.629214479080; 185.775583471871; 191.726515843495; 198.717657594361; 206.311779463682; 212.918080962869; 218.058867232653; 222.689321369683; 227.928134893007; 234.043026551546; 239.996030246296 ; 244.710945553444; 248.532559203805; 252.455304211680; 257.132054604697; 262.232498490425; 265.738546517624; 269.133569405437; 273.120697981357; 277.457638628103; 278.092113952785; 278.092113952786; 278.092113952788; 277.975694874590; 276.730626803023; 274.501243978538; 271.577690132246; 267.948809019530; 263.517567186271; 258.011346719300; 250.941367887279; 241.599276013405; 230.367745163149; 223.259647753706; 216.375530174427; 208.080166232962; 199.025708741683; 192.922357710239; 187.010694981398; 179.979185757652; 172.263438992590; 166.908886363226; 161.725666321843; 155.688081394602; 149.164919164879; 143.469133255423; 138.997588939066; 134.956280534399; 130.410787251421; 125.142069648992; 120.021439889762; 115.941473130762; 112.612121817787; 109.197110607724; 105.150148161534; 100.759157297821; 96.9782907659850; 94.0587577497344; 91.4153822095034; 88.3935137867370; 84.8545826569043; 81.4307608461797; 78.7448626123270; 76.5697123868758; 74.3163874834087; 71.6136292972194; 68.6710789863812; 66.6755409773384; 64.7422736741467; 62.4506602837038; 59.9405540868332; 57.7593088648997; 56.0808773173657; 54.5784487852616; 52.8703803331654; 50.8595112285663; 48.8908731892597; 47.3405820333084; 46.0967346076938; 44.8218727470769; 43.2925516849798; 41.6105897289461; 40.4618646981877; 39.3490121604440; 38.0402800343805; 36.6164853658321; 35.3772929081137; 34.4165779935424; 33.5513169754757; 32.5711937282662; 31.4256586071914; 30.3082563029243; 29.4254960893719; 28.7120913843049; 27.9787086018549; 27.1032591778311; 26.1464857769986; 25.4910894332737; 24.8558231104652; 24.1090910748338; 23.2958490012782; 22.5869008647253; 22.0383964087143; 21.5445211139147; 20.9852109998943; 20.3312697679168; 19.6919903201571; 19.1874652824177; 18.7806070933237; 18.3621174755572; 17.9586839105913];0808773173657; 54.5784487852616; 52.8703803331654; 50.8595112285663; 48.8908731892597; 47.3405820333084; 46.0967346076938; 44.8218727470769; 43.2925516849798; 41.6105897289461; 40.4618646981877; 39.3490121604440; 38.0402800343805; 36.6164853658321; 35.3772929081137; 34.4165779935424; 33.5513169754757; 32.5711937282662; 31.4256586071914; 30.3082563029243; 29.4254960893719; 28.7120913843049; 27.9787086018549; 27.1032591778311; 26.1464857769986; 25.4910894332737; 24.8558231104652; 24.1090910748338; 23.2958490012782; 22.5869008647253; 22.0383964087143; 21.5445211139147; 20.9852109998943; 20.3312697679168; 19.6919903201571; 19.1874652824177; 18.7806070933237; 18.3621174755572; 17.9586839105913];0808773173657; 54.5784487852616; 52.8703803331654; 50.8595112285663; 48.8908731892597; 47.3405820333084; 46.0967346076938; 44.8218727470769; 43.2925516849798; 41.6105897289461; 40.4618646981877; 39.3490121604440; 38.0402800343805; 36.6164853658321; 35.3772929081137; 34.4165779935424; 33.5513169754757; 32.5711937282662; 31.4256586071914; 30.3082563029243; 29.4254960893719; 28.7120913843049; 27.9787086018549; 27.1032591778311; 26.1464857769986; 25.4910894332737; 24.8558231104652; 24.1090910748338; 23.2958490012782; 22.5869008647253; 22.0383964087143; 21.5445211139147; 20.9852109998943; 20.3312697679168; 19.6919903201571; 19.1874652824177; 18.7806070933237; 18.3621174755572; 17.9586839105913];0402800343805; 36.6164853658321; 35.3772929081137; 34.4165779935424; 33.5513169754757; 32.5711937282662; 31.4256586071914; 30.3082563029243; 29.4254960893719; 28.7120913843049; 27.9787086018549; 27.1032591778311; 26.1464857769986; 25.4910894332737; 24.8558231104652; 24.1090910748338; 23.2958490012782; 22.5869008647253; 22.0383964087143; 21.5445211139147; 20.9852109998943; 20.3312697679168; 19.6919903201571; 19.1874652824177; 18.7806070933237; 18.3621174755572; 17.9586839105913];0402800343805; 36.6164853658321; 35.3772929081137; 34.4165779935424; 33.5513169754757; 32.5711937282662; 31.4256586071914; 30.3082563029243; 29.4254960893719; 28.7120913843049; 27.9787086018549; 27.1032591778311; 26.1464857769986; 25.4910894332737; 24.8558231104652; 24.1090910748338; 23.2958490012782; 22.5869008647253; 22.0383964087143; 21.5445211139147; 20.9852109998943; 20.3312697679168; 19.6919903201571; 19.1874652824177; 18.7806070933237; 18.3621174755572; 17.9586839105913];1874652824177; 18.7806070933237; 18.3621174755572; 17.9586839105913];1874652824177; 18.7806070933237; 18.3621174755572; 17.9586839105913];

IL6pout = [1.70000000000000, 1.70000000000000, 1.70000000000021, 1.70000000002828, 1.70000000357964, 1.70000043431404, 1.70001436407676, 1.70004884124751, 1.70004884124751, 1.70004884124751, 1.70126650489711, 1.70777906939199, 1.73136092720512, 1.78958581350041, 1.89997802773520, 2.09084904859125, 2.40482030777192, 2.91569509089136, 3.75547034978718, 5.11726398527806, 6.88114145523957, 8.00029997908417, 9.15150365623831; 10.7430544990301 ; 12.7141375640905; 13.9837479454439; 15.2420076945532; 16.8190122331484; 18.6038247438260; 19.7979809608084; 20.9653960899742; 22.3788614704863; 23.9487366055987; 25.0292483789033; 26.0794301524095; 27.3174584532564; 28.6665733550989; 29.8433547984939; 30.7615334991771; 31.5895479814835; 32.5271558835455; 33.6225897330541; 34.6897330545602; 35.5360715864890; 36.2227719605825; 36.9275700961055; 37.7679404773858; 38.6844446715731; 39.3151607886974 ; 39.9258602779983; 40.6429798550784; 41.4228451795985; 41.5376144592738; 41.5376144592741; 41.5376144592744; 41.6760571105431; 41.8009542226496; 41.9283720945476; 42.0344990249708; 42.1061599995427; 42.1192042398629; 42.0330032417251; 41.7723804727791; 41.2016570732203; 40.2289709290159; 39.4782801429815; 38.6675740467923; 37.5973188186864; 36.3317585958622; 35.4301585666142; 34.5258386260988; 33.4162456426874; 32.1631576144867; 31.2746545769948; 30.4025676197217; 29.3740710831817; 28.2500900729190; 27.2595531134709; 26.4762619272024; 25.7652110882443; 24.9625650634213; 24.0288766822081; 23.1187347931318; 22.3912980212714; 21.7963603272821; 21.1856518285508; 20.4612491041929; 19.6746660649582; 18.9967995597651; 18.4722660146339; 17.9972970531794; 17.4543937166049; 16.8185653863078; 16.2035819572520; 15.7204496827217; 15.3287698109098; 14.9233250159095; 14.4371593418835; 13.9081330656940; 13.5487505525276; 13.2007298011943; 12.7883865834215; 12.3370076159192; 11.9448448222340; 11.6422909980412; 11.3715886677929; 11.0641285489764; 10.7023663174178; 10.3485790593916; 10.0694424540814; 9.84506792891107, 9.61546387090487, 9.34024229389131, 9.03787412867614, 8.83081877588084, 8.63037722540569, 8.39485627527745, 8.13891801546056, 7.91622956715545, 7.74285644622853, 7.58683148080522, 7.41037830603540, 7.20435343093618, 7.00374410312760; 6.84474119208449, 6.71586412534209, 6.58372574859290, 6.42619766816235, 6.25435524614080, 6.13609411085890, 6.02161472408488, 5.88725045576106, 5.74120613714676, 5.61396013539640, 5.51478895634015, 5.42561226099964, 5.32490413133828, 5.20736390913955, 5.09281273595135, 5.00189689142804, 4.92819903296630, 4.85274049588297, 4.78006407344437];0641285489764; 10.7023663174178; 10.3485790593916; 10.0694424540814; 9.84506792891107, 9.61546387090487, 9.34024229389131, 9.03787412867614, 8.83081877588084, 8.63037722540569, 8.39485627527745, 8.13891801546056, 7.91622956715545, 7.74285644622853, 7.58683148080522, 7.41037830603540, 7.20435343093618, 7.00374410312760, 6.84474119208449, 6.71586412534209, 6.58372574859290, 6.42619766816235, 6.25435524614080, 6.13609411085890, 6.02161472408488; 5.88725045576106, 5.74120613714676, 5.61396013539640, 5.51478895634015, 5.42561226099964, 5.32490413133828, 5.20736390913955, 5.09281273595135, 5.00189689142804, 4.92819903296630, 4.85274049588297, 4.78006407344437];0641285489764; 10.7023663174178; 10.3485790593916; 10.0694424540814; 9.84506792891107, 9.61546387090487, 9.34024229389131, 9.03787412867614, 8.83081877588084, 8.63037722540569, 8.39485627527745, 8.13891801546056, 7.91622956715545, 7.74285644622853, 7.58683148080522, 7.41037830603540, 7.20435343093618, 7.00374410312760, 6.84474119208449, 6.71586412534209, 6.58372574859290, 6.42619766816235, 6.25435524614080, 6.13609411085890, 6.02161472408488; 5.88725045576106, 5.74120613714676, 5.61396013539640, 5.51478895634015, 5.42561226099964, 5.32490413133828, 5.20736390913955, 5.09281273595135, 5.00189689142804, 4.92819903296630, 4.85274049588297, 4.78006407344437];58683148080522; 7.41037830603540, 7.20435343093618, 7.00374410312760, 6.84474119208449, 6.71586412534209, 6.58372574859290, 6.42619766816235, 6.25435524614080, 6.13609411085890, 6.02161472408488, 5.88725045576106, 5.74120613714676, 5.61396013539640, 5.51478895634015, 5.42561226099964, 5.32490413133828, 5.20736390913955, 5.09281273595135, 5.00189689142804, 4.92819903296630, 4.85274049588297, 4.78006407344437];58683148080522; 7.41037830603540, 7.20435343093618, 7.00374410312760, 6.84474119208449, 6.71586412534209, 6.58372574859290, 6.42619766816235, 6.25435524614080, 6.13609411085890, 6.02161472408488, 5.88725045576106, 5.74120613714676, 5.61396013539640, 5.51478895634015, 5.42561226099964, 5.32490413133828, 5.20736390913955, 5.09281273595135, 5.00189689142804, 4.92819903296630, 4.85274049588297, 4.78006407344437];

model_out = [GLP1out IL6mout IL6pout];

matrix = zeros (a, 2);

for j = 1: 2

for i = 1: a

v (i, j) = (FSD (i, j) ^ 2) * (abs (dataIL6GLP1 (i, j))) ^ 2;

matrix (i, j) = (dataIL6GLP1 (i, j) - model_out (i, j)) / sqrt (v (i, j));

end

end

wastes_matrix = matrix; % WEIGHTED RESIDUALS!