disp(matlabpath); D:\LeSage\data;D:\LeSage\C_code;D:\LeSage\distrib\C_files;D:\LeSage\var_bvar;D:\LeSage\util;D:\LeSage\spatial;D:\LeSage\regress;D:\LeSage\optimize;D:\LeSage\graphs;D:\LeSage\gibbs;D:\LeSage\distrib;D:\LeSage\diagn;D:\LeSage\coint;D:\MATLABR11\WORKSPACE\Neural_Nets\Synth_tr;D:\MATLABR11\WORKSPACE\Neural_Nets\Cushing;D:\MATLABR11\WORKSPACE\sys502;D:\MATLABR11\WORKSPACE\sys302;D:\MATLABR11\WORKSPACE\Phil_Sim;D:\MATLABR11\toolbox\matlab\general;D:\MATLABR11\toolbox\matlab\ops;D:\MATLABR11\toolbox\matlab\lang;D:\MATLABR11\toolbox\matlab\elmat;D:\MATLABR11\toolbox\matlab\elfun;D:\MATLABR11\toolbox\matlab\specfun;D:\MATLABR11\toolbox\matlab\matfun;D:\MATLABR11\toolbox\matlab\datafun;D:\MATLABR11\toolbox\matlab\polyfun;D:\MATLABR11\toolbox\matlab\funfun;D:\MATLABR11\toolbox\matlab\sparfun;D:\MATLABR11\toolbox\matlab\graph2d;D:\MATLABR11\toolbox\matlab\graph3d;D:\MATLABR11\toolbox\matlab\specgraph;D:\MATLABR11\toolbox\matlab\graphics;D:\MATLABR11\toolbox\matlab\uitools;D:\MATLABR11\toolbox\matlab\strfun;D:\MATLABR11\toolbox\matlab\iofun;D:\MATLABR11\toolbox\matlab\timefun;D:\MATLABR11\toolbox\matlab\datatypes;D:\MATLABR11\toolbox\matlab\winfun;D:\MATLABR11\toolbox\matlab\demos;D:\MATLABR11\toolbox\compiler;D:\MATLABR11\toolbox\finance\finance;D:\MATLABR11\toolbox\finance\calendar;D:\MATLABR11\toolbox\finance\findemos;D:\MATLABR11\toolbox\fuzzy\fuzzy;D:\MATLABR11\toolbox\fuzzy\fuzdemos;D:\MATLABR11\toolbox\stats;D:\MATLABR11\toolbox\nnet\nnet;D:\MATLABR11\toolbox\nnet\nndemos;D:\MATLABR11\toolbox\nnet\nnutils;D:\MATLABR11\toolbox\nnet\nnobsolete;D:\MATLABR11\toolbox\signal\signal;D:\MATLABR11\toolbox\signal\siggui;D:\MATLABR11\toolbox\signal\sigdemos;D:\MATLABR11\toolbox\splines;D:\MATLABR11\toolbox\optim;D:\MATLABR11\toolbox\ident;D:\MATLABR11\toolbox\tour;D:\MATLABR11\work;D:\MATLABR11\toolbox\local feval('whos') Name Size Bytes Class K5_5 89530x2 1432480 double array W5 917x917 6727112 double array ans 3x3 72 double array dat_1 918x3 22032 double array dat_2 918x3 22032 double array Grand total is 1025466 elements using 8203728 bytes disp(cd); d:\matlabr11\workspace\sys502 clear disp(feature('ArrayEditing')) 1 feval('whos') Name Size Bytes Class data 26x7 1456 double array Grand total is 182 elements using 1456 bytes y = data(:,1); feval('whos') Name Size Bytes Class data 26x7 1456 double array y 26x1 208 double array Grand total is 208 elements using 1664 bytes x = data(:,2:7); feval('whos') Name Size Bytes Class data 26x7 1456 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 364 elements using 2912 bytes X = [ones(size(x,1),1),x]; feval('whos') Name Size Bytes Class X 26x7 1456 double array data 26x7 1456 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 546 elements using 4368 bytes edit dbstack ; disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m mdbstatus 'd:\matlabr11\workspace\sys502\poisson_reg.m' dbstack ; disp(which('poisson')); d:\matlabr11\workspace\sys502\poisson.m clear 'd:\matlabr11\workspace\sys502\poisson.m' [beta,mu,dev,df,se]=poisson(y,X); feval('whos') Name Size Bytes Class X 26x7 1456 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array mu 26x1 208 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 588 elements using 4704 bytes beta beta = 3.3903 -0.068311 0.00033078 0.21243 0.0090362 0.092624 0.13154 feval('whos') Name Size Bytes Class X 26x7 1456 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array mu 26x1 208 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 588 elements using 4704 bytes output output = 3.3903 1.2262e-181 -0.068311 0.016492 0.00033078 1.851e-005 0.21243 0.00048597 0.0090362 0.011437 0.092624 0.017286 0.13154 0.036575 feval('whos') Name Size Bytes Class X 26x7 1456 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 602 elements using 4816 bytes feval('whos') Name Size Bytes Class X 26x7 1456 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 602 elements using 4816 bytes disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbstop at 5 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' Poisson_reg dbstack > In d:\matlabr11\workspace\sys502\Poisson_reg.m at line 5 feval('whos') Name Size Bytes Class X 26x7 1456 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 602 elements using 4816 bytes dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 609 elements using 4872 bytes b b = 20.792 -3.354 0.019833 11.37 0.46189 4.8289 6.5231 feval('whos') Name Size Bytes Class X 26x7 1456 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 609 elements using 4872 bytes dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 610 elements using 4880 bytes dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 610 elements using 4880 bytes eval( 'if exist(''diag'',''var''), mauifunc(diag), end', ''); eval( 'if exist(''X'',''var''), mauifunc(X), end', ''); Columns 1 through 6 1 2 967 0 39 0 1 2 815 0 33 2 1 3 900 1 35 1 1 3 1007 1 24 2 1 3 1100 0 50 1.5 1 4 897 0 25 1 1 5 1400 0 30 1 1 3 2261 0 29 2 1 4 1290 1 33 1.5 1 4 2104 0 40 1 1 8 2240 1 50 2 1 2 641 0 25 0 1 3 862 0 25 0 1 4 1043 0 30 0 1 4 1325 0 50 0 1 2 782 1 25 0 1 3 1126 1 30 0 1 4 1226 0 37 2 1 2 929 0 30 1 1 4 1137 0 25 0 1 3 743 0 25 0 1 3 596 0 50 0 1 2 803 0 27 0 1 2 696 0 30 1 1 2 691 0 30 2 1 3 1023 1 30 1 Column 7 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 feval('whos') Name Size Bytes Class X 26x7 1456 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 610 elements using 4880 bytes eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); inv(X'*diag(mu)*X) Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 8.581100e-017. > In d:\matlabr11\workspace\sys502\Poisson_reg.m at line 11 ans = Columns 1 through 6 2.7931e-027 5.2784e-028 -6.0748e-031 -9.0232e-028 -8.3315e-029 -2.9347e-028 5.2784e-028 1.0139e-028 -1.1322e-031 -1.7404e-028 -1.5967e-029 -5.6499e-029 -6.0748e-031 -1.1322e-031 1.582e-034 2.1305e-031 1.697e-032 4.866e-032 -9.0232e-028 -1.7404e-028 2.1305e-031 3.1569e-028 2.6619e-029 8.5401e-029 -8.3315e-029 -1.5967e-029 1.697e-032 2.6619e-029 2.551e-030 9.4341e-030 -2.9347e-028 -5.6499e-029 4.866e-032 8.5401e-029 9.4341e-030 3.9677e-029 -3.561e-028 -6.4754e-029 8.8964e-032 1.1705e-028 9.9325e-030 3.0588e-029 Column 7 -3.5611e-028 -6.4754e-029 8.8964e-032 1.1705e-028 9.9325e-030 3.0588e-029 8.2623e-029 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes dbstep Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 8.581100e-017. > In d:\matlabr11\workspace\sys502\Poisson_reg.m at line 11 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 1.8573e+022 8.921e+023 2.9358e+026 1.9052e+027 2.0424e+027 1.1003e+018 8.3243e+021 1.4002e+034 1.0415e+029 2.7969e+031 6.7997e+035 3.0587e+019 8.5595e+019 1.091e+021 3.0108e+027 4.343e+025 1.4032e+028 1.6313e+028 1.165e+025 1.0269e+018 1.1872e+016 6.6573e+019 2.8129e+018 1.6845e+020 1.908e+022 3.3433e+026 b b = 19.792 -3.354 0.019833 11.37 0.46189 4.8289 6.5231 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes eval( 'if exist(''diag'',''var''), mauifunc(diag), end', ''); eval( 'if exist(''diag'',''var''), mauifunc(diag), end', ''); eval( 'if exist(''REGRESSION'',''var''), mauifunc(REGRESSION), end', ''); eval( 'if exist(''for'',''var''), mauifunc(for), end', ''); eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''lp'',''var''), mauifunc(lp), end', ''); eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 8.1583 eval( 'if exist(''devold'',''var''), mauifunc(devold), end', ''); eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 1.8573e+022 8.921e+023 2.9358e+026 1.9052e+027 2.0424e+027 1.1003e+018 8.3243e+021 1.4002e+034 1.0415e+029 2.7969e+031 6.7997e+035 3.0587e+019 8.5595e+019 1.091e+021 3.0108e+027 4.343e+025 1.4032e+028 1.6313e+028 1.165e+025 1.0269e+018 1.1872e+016 6.6573e+019 2.8129e+018 1.6845e+020 1.908e+022 3.3433e+026 eval( 'if exist(''exp'',''var''), mauifunc(exp), end', ''); eval( 'if exist(''abs'',''var''), mauifunc(abs), end', ''); eval( 'if exist(''devold'',''var''), mauifunc(devold), end', ''); eval( 'if exist(''reweighted'',''var''), mauifunc(reweighted), end', ''); eval( 'if exist(''ny'',''var''), mauifunc(ny), end', ''); eval( 'if exist(''freedom'',''var''), mauifunc(freedom), end', ''); disp(which('poisson')); d:\matlabr11\workspace\sys502\poisson.m dbstop at 39 in 'd:\matlabr11\workspace\sys502\poisson.m' eval( 'if exist(''values'',''var''), mauifunc(values), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 1.8573e+022 8.921e+023 2.9358e+026 1.9052e+027 2.0424e+027 1.1003e+018 8.3243e+021 1.4002e+034 1.0415e+029 2.7969e+031 6.7997e+035 3.0587e+019 8.5595e+019 1.091e+021 3.0108e+027 4.343e+025 1.4032e+028 1.6313e+028 1.165e+025 1.0269e+018 1.1872e+016 6.6573e+019 2.8129e+018 1.6845e+020 1.908e+022 3.3433e+026 eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''reweighted'',''var''), mauifunc(reweighted), end', ''); eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 8.1583 eval( 'if exist(''yadj'',''var''), mauifunc(yadj), end', ''); eval( 'if exist(''disp'',''var''), mauifunc(disp), end', ''); eval( 'if exist(''my'',''var''), mauifunc(my), end', ''); eval( 'if exist(''se'',''var''), mauifunc(se), end', ''); 0.11797 0.028489 7.7246e-005 0.060897 0.0035729 0.038908 0.062922 dbquit » feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbclear at 5 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' [beta,mu,dev,df,se]=poisson(y,x); dbstack > In d:\matlabr11\workspace\sys502\poisson.m at line 39 feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x6 1248 double array y 26x1 208 double array yadj 26x1 208 double array Grand total is 267 elements using 2136 bytes mu mu = 53.5 55.5 56.5 58.5 64.5 44.5 49.5 70.5 72.5 82.5 85.5 45.5 47.5 49.5 56.5 60.5 62.5 64.5 66.5 35.5 38.5 43.5 46.5 46.5 50.5 65.5 feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x6 1248 double array y 26x1 208 double array yadj 26x1 208 double array Grand total is 267 elements using 2136 bytes eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''reweighted'',''var''), mauifunc(reweighted), end', ''); eval( 'if exist(''reweighted'',''var''), mauifunc(reweighted), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); dbstep feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x6 1248 double array y 26x1 208 double array yadj 26x1 208 double array Grand total is 267 elements using 2136 bytes dbstep feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x6 1248 double array y 26x1 208 double array yadj 26x1 208 double array Grand total is 267 elements using 2136 bytes eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 53.5 55.5 56.5 58.5 64.5 44.5 49.5 70.5 72.5 82.5 85.5 45.5 47.5 49.5 56.5 60.5 62.5 64.5 66.5 35.5 38.5 43.5 46.5 46.5 50.5 65.5 dev dev = 0.12056 feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x6 1248 double array y 26x1 208 double array yadj 26x1 208 double array Grand total is 267 elements using 2136 bytes eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); dbstep feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x6 1248 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 293 elements using 2344 bytes eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); z z = 3.9703 4.0074 4.0254 4.0605 4.1589 3.7843 3.8919 4.2485 4.2767 4.4067 4.4427 3.8067 3.8502 3.8919 4.0254 4.0944 4.1272 4.1589 4.1897 3.5554 3.6377 3.7613 3.8287 3.8287 3.9121 4.1744 feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x6 1248 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 293 elements using 2344 bytes eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''nx'',''var''), mauifunc(nx), end', ''); 6 size(mu*ones(1,nx)) ans = 26 6 feval('whos') Name Size Bytes Class ans 1x2 16 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x6 1248 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 295 elements using 2360 bytes eval( 'if exist(''y'',''var''), mauifunc(y), end', ''); 53 55 56 58 64 44 49 70 72 82 85 45 47 49 56 60 62 64 66 35 38 43 46 46 50 65 eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); dbquit » feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes disp(which('poisson')); d:\matlabr11\workspace\sys502\poisson.m dbstop at 42 in 'd:\matlabr11\workspace\sys502\poisson.m' disp(which('poisson')); d:\matlabr11\workspace\sys502\poisson.m clear 'd:\matlabr11\workspace\sys502\poisson.m' disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m clear 'd:\matlabr11\workspace\sys502\poisson_reg.m' disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbstop at 5 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' Poisson_reg dbstack > In d:\matlabr11\workspace\sys502\Poisson_reg.m at line 5 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes b b = 1 1 1 1 1 1 1 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); eval( 'if exist(''X'',''var''), mauifunc(X), end', ''); Columns 1 through 6 1 2 967 0 39 0 1 2 815 0 33 2 1 3 900 1 35 1 1 3 1007 1 24 2 1 3 1100 0 50 1.5 1 4 897 0 25 1 1 5 1400 0 30 1 1 3 2261 0 29 2 1 4 1290 1 33 1.5 1 4 2104 0 40 1 1 8 2240 1 50 2 1 2 641 0 25 0 1 3 862 0 25 0 1 4 1043 0 30 0 1 4 1325 0 50 0 1 2 782 1 25 0 1 3 1126 1 30 0 1 4 1226 0 37 2 1 2 929 0 30 1 1 4 1137 0 25 0 1 3 743 0 25 0 1 3 596 0 50 0 1 2 803 0 27 0 1 2 696 0 30 1 1 2 691 0 30 2 1 3 1023 1 30 1 Column 7 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 eval( 'if exist(''X'',''var''), mauifunc(X), end', ''); Columns 1 through 6 1 2 967 0 39 0 1 2 815 0 33 2 1 3 900 1 35 1 1 3 1007 1 24 2 1 3 1100 0 50 1.5 1 4 897 0 25 1 1 5 1400 0 30 1 1 3 2261 0 29 2 1 4 1290 1 33 1.5 1 4 2104 0 40 1 1 8 2240 1 50 2 1 2 641 0 25 0 1 3 862 0 25 0 1 4 1043 0 30 0 1 4 1325 0 50 0 1 2 782 1 25 0 1 3 1126 1 30 0 1 4 1226 0 37 2 1 2 929 0 30 1 1 4 1137 0 25 0 1 3 743 0 25 0 1 3 596 0 50 0 1 2 803 0 27 0 1 2 696 0 30 1 1 2 691 0 30 2 1 3 1023 1 30 1 Column 7 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 1.8573e+022 8.921e+023 2.9358e+026 1.9052e+027 2.0424e+027 1.1003e+018 8.3243e+021 1.4002e+034 1.0415e+029 2.7969e+031 6.7997e+035 3.0587e+019 8.5595e+019 1.091e+021 3.0108e+027 4.343e+025 1.4032e+028 1.6313e+028 1.165e+025 1.0269e+018 1.1872e+016 6.6573e+019 2.8129e+018 1.6845e+020 1.908e+022 3.3433e+026 dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes dbstep Warning: Matrix is singular to working precision. > In d:\matlabr11\workspace\sys502\Poisson_reg.m at line 11 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes b b = NaN NaN NaN NaN NaN NaN NaN feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes mu mu = Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf 9.4908e+290 Inf Inf Inf Inf Inf Inf Inf Inf Inf 1.9562e+282 Inf Inf Inf Inf feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf 9.4908e+290 Inf Inf Inf Inf Inf Inf Inf Inf Inf 1.9562e+282 Inf Inf Inf Inf dbquit » feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes output output = 3.3903 1.2262e-181 -0.068311 0.016492 0.00033078 1.851e-005 0.21243 0.00048597 0.0090362 0.011437 0.092624 0.017286 0.13154 0.036575 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes b = output(:,1) + .1 b = 3.4903 0.031689 0.10033 0.31243 0.10904 0.19262 0.23154 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbstop at 13 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbclear at 13 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbclear at 5 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m clear 'd:\matlabr11\workspace\sys502\poisson_reg.m' disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbstop at 8 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' Poisson_reg dbstack > In d:\matlabr11\workspace\sys502\Poisson_reg.m at line 8 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf 9.4908e+290 Inf Inf Inf Inf Inf Inf Inf Inf Inf 1.9562e+282 Inf Inf Inf Inf dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes mu mu = 3.3521e+045 6.1009e+038 4.4628e+042 7.4937e+046 9.5649e+051 8.3836e+041 1.2338e+064 4.1347e+101 4.0179e+059 1.6845e+095 7.9496e+101 5.7291e+030 2.5206e+040 3.4576e+048 5.9355e+061 1.0904e+037 1.8935e+052 1.0267e+057 4.242e+043 1.983e+052 1.3054e+035 7.8406e+029 6.4727e+037 2.3685e+033 1.7389e+033 5.9202e+047 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes eval( 'if exist(''sqrt'',''var''), mauifunc(sqrt), end', ''); eval( 'if exist(''nx'',''var''), mauifunc(nx), end', ''); eval( 'if exist(''disp'',''var''), mauifunc(disp), end', ''); eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 3.3521e+045 6.1009e+038 4.4628e+042 7.4937e+046 9.5649e+051 8.3836e+041 1.2338e+064 4.1347e+101 4.0179e+059 1.6845e+095 7.9496e+101 5.7291e+030 2.5206e+040 3.4576e+048 5.9355e+061 1.0904e+037 1.8935e+052 1.0267e+057 4.242e+043 1.983e+052 1.3054e+035 7.8406e+029 6.4727e+037 2.3685e+033 1.7389e+033 5.9202e+047 eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); disp(which('poisson')); d:\matlabr11\workspace\sys502\poisson.m dbstop at 42 in 'd:\matlabr11\workspace\sys502\poisson.m' feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 3.3521e+045 6.1009e+038 4.4628e+042 7.4937e+046 9.5649e+051 8.3836e+041 1.2338e+064 4.1347e+101 4.0179e+059 1.6845e+095 7.9496e+101 5.7291e+030 2.5206e+040 3.4576e+048 5.9355e+061 1.0904e+037 1.8935e+052 1.0267e+057 4.242e+043 1.983e+052 1.3054e+035 7.8406e+029 6.4727e+037 2.3685e+033 1.7389e+033 5.9202e+047 dbquit » feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbclear at 8 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' poisson(y,X) dbstack > In d:\matlabr11\workspace\sys502\poisson.m at line 42 feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 319 elements using 2552 bytes eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 0.12056 eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 0.12056 eval( 'if exist(''disp'',''var''), mauifunc(disp), end', ''); eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 53.5 55.5 56.5 58.5 64.5 44.5 49.5 70.5 72.5 82.5 85.5 45.5 47.5 49.5 56.5 60.5 62.5 64.5 66.5 35.5 38.5 43.5 46.5 46.5 50.5 65.5 eval( 'if exist(''lp'',''var''), mauifunc(lp), end', ''); 3.9797 4.0164 4.0342 4.069 4.1667 3.7955 3.902 4.2556 4.2836 4.4128 4.4485 3.8177 3.8607 3.902 4.0342 4.1026 4.1352 4.1667 4.1972 3.5695 3.6507 3.7728 3.8395 3.8395 3.922 4.1821 eval( 'if exist(''log'',''var''), mauifunc(log), end', ''); eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 0.12056 eval( 'if exist(''reweighted'',''var''), mauifunc(reweighted), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''least'',''var''), mauifunc(least), end', ''); eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 53.5 55.5 56.5 58.5 64.5 44.5 49.5 70.5 72.5 82.5 85.5 45.5 47.5 49.5 56.5 60.5 62.5 64.5 66.5 35.5 38.5 43.5 46.5 46.5 50.5 65.5 eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); dbstep feval('whos') Name Size Bytes Class beta 7x1 56 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 326 elements using 2608 bytes eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''x'',''var''), mauifunc(x), end', ''); Columns 1 through 6 1 2 967 0 39 0 1 2 815 0 33 2 1 3 900 1 35 1 1 3 1007 1 24 2 1 3 1100 0 50 1.5 1 4 897 0 25 1 1 5 1400 0 30 1 1 3 2261 0 29 2 1 4 1290 1 33 1.5 1 4 2104 0 40 1 1 8 2240 1 50 2 1 2 641 0 25 0 1 3 862 0 25 0 1 4 1043 0 30 0 1 4 1325 0 50 0 1 2 782 1 25 0 1 3 1126 1 30 0 1 4 1226 0 37 2 1 2 929 0 30 1 1 4 1137 0 25 0 1 3 743 0 25 0 1 3 596 0 50 0 1 2 803 0 27 0 1 2 696 0 30 1 1 2 691 0 30 2 1 3 1023 1 30 1 Column 7 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 beta beta = 3.3924 -0.069604 0.00033732 0.21159 0.0089995 0.091533 0.1309 feval('whos') Name Size Bytes Class beta 7x1 56 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 326 elements using 2608 bytes eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''squares'',''var''), mauifunc(squares), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''least'',''var''), mauifunc(least), end', ''); dbstep feval('whos') Name Size Bytes Class beta 7x1 56 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 326 elements using 2608 bytes dbstep dbstep dbstep dbstep feval('whos') Name Size Bytes Class beta 7x1 56 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 326 elements using 2608 bytes eval( 'if exist(''ones'',''var''), mauifunc(ones), end', ''); mu mu = 50.927 55.047 60.662 62.424 62.925 41.809 48.336 80.669 66.356 71.898 84.422 45.848 46.075 47.786 62.918 59.411 65.099 65.009 57.917 41.369 38.832 46.277 43.254 46.971 51.386 60.449 feval('whos') Name Size Bytes Class beta 7x1 56 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 326 elements using 2608 bytes eval( 'if exist(''z'',''var''), mauifunc(z), end', ''); 3.9703 4.0074 4.0254 4.0605 4.1589 3.7843 3.8919 4.2485 4.2767 4.4067 4.4427 3.8067 3.8502 3.8919 4.0254 4.0944 4.1272 4.1589 4.1897 3.5554 3.6377 3.7613 3.8287 3.8287 3.9121 4.1744 eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 50.927 55.047 60.662 62.424 62.925 41.809 48.336 80.669 66.356 71.898 84.422 45.848 46.075 47.786 62.918 59.411 65.099 65.009 57.917 41.369 38.832 46.277 43.254 46.971 51.386 60.449 eval( 'if exist(''nx'',''var''), mauifunc(nx), end', ''); 7 eval( 'if exist(''nargin'',''var''), mauifunc(nargin), end', ''); eval( 'if exist(''if'',''var''), mauifunc(if), end', ''); eval( 'if exist(''lp'',''var''), mauifunc(lp), end', ''); 3.9304 4.0082 4.1053 4.1339 4.1419 3.7331 3.8782 4.3904 4.195 4.2753 4.4358 3.8253 3.8303 3.8667 4.1418 4.0845 4.1759 4.1745 4.059 3.7225 3.6592 3.8346 3.7671 3.8495 3.9394 4.1018 eval( 'if exist(''beta'',''var''), mauifunc(beta), end', ''); 3.3924 -0.069604 0.00033732 0.21159 0.0089995 0.091533 0.1309 eval( 'if exist(''z'',''var''), mauifunc(z), end', ''); 3.9703 4.0074 4.0254 4.0605 4.1589 3.7843 3.8919 4.2485 4.2767 4.4067 4.4427 3.8067 3.8502 3.8919 4.0254 4.0944 4.1272 4.1589 4.1897 3.5554 3.6377 3.7613 3.8287 3.8287 3.9121 4.1744 eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 0.12056 eval( 'if exist(''devold'',''var''), mauifunc(devold), end', ''); 1000000 eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 50.927 55.047 60.662 62.424 62.925 41.809 48.336 80.669 66.356 71.898 84.422 45.848 46.075 47.786 62.918 59.411 65.099 65.009 57.917 41.369 38.832 46.277 43.254 46.971 51.386 60.449 eval( 'if exist(''exp'',''var''), mauifunc(exp), end', ''); eval( 'if exist(''offset'',''var''), mauifunc(offset), end', ''); 0 eval( 'if exist(''exp'',''var''), mauifunc(exp), end', ''); eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 50.927 55.047 60.662 62.424 62.925 41.809 48.336 80.669 66.356 71.898 84.422 45.848 46.075 47.786 62.918 59.411 65.099 65.009 57.917 41.369 38.832 46.277 43.254 46.971 51.386 60.449 eval( 'if exist(''exp'',''var''), mauifunc(exp), end', ''); eval( 'if exist(''offset'',''var''), mauifunc(offset), end', ''); 0 eval( 'if exist(''lp'',''var''), mauifunc(lp), end', ''); 3.9304 4.0082 4.1053 4.1339 4.1419 3.7331 3.8782 4.3904 4.195 4.2753 4.4358 3.8253 3.8303 3.8667 4.1418 4.0845 4.1759 4.1745 4.059 3.7225 3.6592 3.8346 3.7671 3.8495 3.9394 4.1018 eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 50.927 55.047 60.662 62.424 62.925 41.809 48.336 80.669 66.356 71.898 84.422 45.848 46.075 47.786 62.918 59.411 65.099 65.009 57.917 41.369 38.832 46.277 43.254 46.971 51.386 60.449 eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 50.927 55.047 60.662 62.424 62.925 41.809 48.336 80.669 66.356 71.898 84.422 45.848 46.075 47.786 62.918 59.411 65.099 65.009 57.917 41.369 38.832 46.277 43.254 46.971 51.386 60.449 eval( 'if exist(''yadj'',''var''), mauifunc(yadj), end', ''); 53 55 56 58 64 44 49 70 72 82 85 45 47 49 56 60 62 64 66 35 38 43 46 46 50 65 eval( 'if exist(''log'',''var''), mauifunc(log), end', ''); eval( 'if exist(''yadj'',''var''), mauifunc(yadj), end', ''); 53 55 56 58 64 44 49 70 72 82 85 45 47 49 56 60 62 64 66 35 38 43 46 46 50 65 eval( 'if exist(''nx'',''var''), mauifunc(nx), end', ''); 7 eval( 'if exist(''nx'',''var''), mauifunc(nx), end', ''); 7 eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); eval( 'if exist(''nargout'',''var''), mauifunc(nargout), end', ''); eval( 'if exist(''if'',''var''), mauifunc(if), end', ''); eval( 'if exist(''lp'',''var''), mauifunc(lp), end', ''); 3.9304 4.0082 4.1053 4.1339 4.1419 3.7331 3.8782 4.3904 4.195 4.2753 4.4358 3.8253 3.8303 3.8667 4.1418 4.0845 4.1759 4.1745 4.059 3.7225 3.6592 3.8346 3.7671 3.8495 3.9394 4.1018 eval( 'if exist(''beta'',''var''), mauifunc(beta), end', ''); 3.3924 -0.069604 0.00033732 0.21159 0.0089995 0.091533 0.1309 eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 0.12056 eval( 'if exist(''devold'',''var''), mauifunc(devold), end', ''); 1000000 eval( 'if exist(''abs'',''var''), mauifunc(abs), end', ''); eval( 'if exist(''nargin'',''var''), mauifunc(nargin), end', ''); eval( 'if exist(''Iteratively'',''var''), mauifunc(Iteratively), end', ''); eval( 'if exist(''size'',''var''), mauifunc(size), end', ''); eval( 'if exist(''y'',''var''), mauifunc(y), end', ''); 53 55 56 58 64 44 49 70 72 82 85 45 47 49 56 60 62 64 66 35 38 43 46 46 50 65 dbquit » feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes disp(which('poisson')); d:\matlabr11\workspace\sys502\poisson.m dbclear at 42 in 'd:\matlabr11\workspace\sys502\poisson.m' b b = 3.4903 0.031689 0.10033 0.31243 0.10904 0.19262 0.23154 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 7x7 392 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 659 elements using 5272 bytes Xb ??? Undefined function or variable 'Xb'. X*b ans = 104.83 89.307 98.204 107.93 119.69 96.532 147.58 233.98 137.24 219.27 234.63 70.823 93.028 111.76 142.24 85.282 120.37 131.27 100.46 120.42 80.857 68.834 87.063 76.848 76.539 110 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes disp(which('poisson')); d:\matlabr11\workspace\sys502\poisson.m dbstop at 42 in 'd:\matlabr11\workspace\sys502\poisson.m' poisson(y,X) dbstack > In d:\matlabr11\workspace\sys502\poisson.m at line 42 feval('whos') Name Size Bytes Class dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 319 elements using 2552 bytes eval( 'if exist(''sqrt'',''var''), mauifunc(sqrt), end', ''); eval( 'if exist(''disp'',''var''), mauifunc(disp), end', ''); eval( 'if exist(''ones'',''var''), mauifunc(ones), end', ''); eval( 'if exist(''lp'',''var''), mauifunc(lp), end', ''); 3.9797 4.0164 4.0342 4.069 4.1667 3.7955 3.902 4.2556 4.2836 4.4128 4.4485 3.8177 3.8607 3.902 4.0342 4.1026 4.1352 4.1667 4.1972 3.5695 3.6507 3.7728 3.8395 3.8395 3.922 4.1821 eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 0.12056 eval( 'if exist(''reweighted'',''var''), mauifunc(reweighted), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); dbstep feval('whos') Name Size Bytes Class beta 7x1 56 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 326 elements using 2608 bytes eval( 'if exist(''disp'',''var''), mauifunc(disp), end', ''); eval( 'if exist(''nx'',''var''), mauifunc(nx), end', ''); 7 dbstep feval('whos') Name Size Bytes Class beta 7x1 56 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 326 elements using 2608 bytes eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); lp lp = 3.9304 4.0082 4.1053 4.1339 4.1419 3.7331 3.8782 4.3904 4.195 4.2753 4.4358 3.8253 3.8303 3.8667 4.1418 4.0845 4.1759 4.1745 4.059 3.7225 3.6592 3.8346 3.7671 3.8495 3.9394 4.1018 feval('whos') Name Size Bytes Class beta 7x1 56 double array dev 1x1 8 double array devold 1x1 8 double array lp 26x1 208 double array mu 26x1 208 double array mx 1x1 8 double array my 1x1 8 double array nx 1x1 8 double array ny 1x1 8 double array offset 1x1 8 double array x 26x7 1456 double array y 26x1 208 double array yadj 26x1 208 double array z 26x1 208 double array Grand total is 326 elements using 2608 bytes eval( 'if exist(''ones'',''var''), mauifunc(ones), end', ''); eval( 'if exist(''sqrt'',''var''), mauifunc(sqrt), end', ''); eval( 'if exist(''sqrt'',''var''), mauifunc(sqrt), end', ''); eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); eval( 'if exist(''end'',''var''), mauifunc(end), end', ''); eval( 'if exist(''my'',''var''), mauifunc(my), end', ''); 26 eval( 'if exist(''nx'',''var''), mauifunc(nx), end', ''); 7 eval( 'if exist(''diag'',''var''), mauifunc(diag), end', ''); eval( 'if exist(''inv'',''var''), mauifunc(inv), end', ''); eval( 'if exist(''nx'',''var''), mauifunc(nx), end', ''); 7 eval( 'if exist(''lp'',''var''), mauifunc(lp), end', ''); 3.9304 4.0082 4.1053 4.1339 4.1419 3.7331 3.8782 4.3904 4.195 4.2753 4.4358 3.8253 3.8303 3.8667 4.1418 4.0845 4.1759 4.1745 4.059 3.7225 3.6592 3.8346 3.7671 3.8495 3.9394 4.1018 eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 53.5 55.5 56.5 58.5 64.5 44.5 49.5 70.5 72.5 82.5 85.5 45.5 47.5 49.5 56.5 60.5 62.5 64.5 66.5 35.5 38.5 43.5 46.5 46.5 50.5 65.5 eval( 'if exist(''sum'',''var''), mauifunc(sum), end', ''); eval( 'if exist(''dev'',''var''), mauifunc(dev), end', ''); 0.12056 eval( 'if exist(''offset'',''var''), mauifunc(offset), end', ''); 0 eval( 'if exist(''devold'',''var''), mauifunc(devold), end', ''); 1000000 eval( 'if exist(''nargin'',''var''), mauifunc(nargin), end', ''); eval( 'if exist(''abs'',''var''), mauifunc(abs), end', ''); eval( 'if exist(''exp'',''var''), mauifunc(exp), end', ''); eval( 'if exist(''devold'',''var''), mauifunc(devold), end', ''); 1000000 eval( 'if exist(''x'',''var''), mauifunc(x), end', ''); Columns 1 through 6 1 2 967 0 39 0 1 2 815 0 33 2 1 3 900 1 35 1 1 3 1007 1 24 2 1 3 1100 0 50 1.5 1 4 897 0 25 1 1 5 1400 0 30 1 1 3 2261 0 29 2 1 4 1290 1 33 1.5 1 4 2104 0 40 1 1 8 2240 1 50 2 1 2 641 0 25 0 1 3 862 0 25 0 1 4 1043 0 30 0 1 4 1325 0 50 0 1 2 782 1 25 0 1 3 1126 1 30 0 1 4 1226 0 37 2 1 2 929 0 30 1 1 4 1137 0 25 0 1 3 743 0 25 0 1 3 596 0 50 0 1 2 803 0 27 0 1 2 696 0 30 1 1 2 691 0 30 2 1 3 1023 1 30 1 Column 7 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 eval( 'if exist(''nargin'',''var''), mauifunc(nargin), end', ''); eval( 'if exist(''beta'',''var''), mauifunc(beta), end', ''); 3.3924 -0.069604 0.00033732 0.21159 0.0089995 0.091533 0.1309 eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 53.5 55.5 56.5 58.5 64.5 44.5 49.5 70.5 72.5 82.5 85.5 45.5 47.5 49.5 56.5 60.5 62.5 64.5 66.5 35.5 38.5 43.5 46.5 46.5 50.5 65.5 eval( 'if exist(''sum'',''var''), mauifunc(sum), end', ''); eval( 'if exist(''offset'',''var''), mauifunc(offset), end', ''); 0 eval( 'if exist(''devold'',''var''), mauifunc(devold), end', ''); 1000000 eval( 'if exist(''devold'',''var''), mauifunc(devold), end', ''); 1000000 eval( 'if exist(''Iteratively'',''var''), mauifunc(Iteratively), end', ''); eval( 'if exist(''offset'',''var''), mauifunc(offset), end', ''); 0 dbquit » feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes disp(which('poisson')); d:\matlabr11\workspace\sys502\poisson.m dbclear at 42 in 'd:\matlabr11\workspace\sys502\poisson.m' b = output(:,1) b = 3.3903 -0.068311 0.00033078 0.21243 0.0090362 0.092624 0.13154 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbstop at 8 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' Poisson_reg dbstack > In d:\matlabr11\workspace\sys502\Poisson_reg.m at line 8 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes eval( 'if exist(''y'',''var''), mauifunc(y), end', ''); 53 55 56 58 64 44 49 70 72 82 85 45 47 49 56 60 62 64 66 35 38 43 46 46 50 65 mu mu = 50.7 54.964 60.604 62.363 62.799 41.773 48.207 79.88 66.24 71.311 84.285 45.747 45.966 47.686 62.718 59.275 64.9 64.954 57.754 41.224 38.745 46.261 43.088 46.879 51.344 60.332 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes eval( 'if exist(''X'',''var''), mauifunc(X), end', ''); Columns 1 through 6 1 2 967 0 39 0 1 2 815 0 33 2 1 3 900 1 35 1 1 3 1007 1 24 2 1 3 1100 0 50 1.5 1 4 897 0 25 1 1 5 1400 0 30 1 1 3 2261 0 29 2 1 4 1290 1 33 1.5 1 4 2104 0 40 1 1 8 2240 1 50 2 1 2 641 0 25 0 1 3 862 0 25 0 1 4 1043 0 30 0 1 4 1325 0 50 0 1 2 782 1 25 0 1 3 1126 1 30 0 1 4 1226 0 37 2 1 2 929 0 30 1 1 4 1137 0 25 0 1 3 743 0 25 0 1 3 596 0 50 0 1 2 803 0 27 0 1 2 696 0 30 1 1 2 691 0 30 2 1 3 1023 1 30 1 Column 7 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 eval( 'if exist(''X'',''var''), mauifunc(X), end', ''); Columns 1 through 6 1 2 967 0 39 0 1 2 815 0 33 2 1 3 900 1 35 1 1 3 1007 1 24 2 1 3 1100 0 50 1.5 1 4 897 0 25 1 1 5 1400 0 30 1 1 3 2261 0 29 2 1 4 1290 1 33 1.5 1 4 2104 0 40 1 1 8 2240 1 50 2 1 2 641 0 25 0 1 3 862 0 25 0 1 4 1043 0 30 0 1 4 1325 0 50 0 1 2 782 1 25 0 1 3 1126 1 30 0 1 4 1226 0 37 2 1 2 929 0 30 1 1 4 1137 0 25 0 1 3 743 0 25 0 1 3 596 0 50 0 1 2 803 0 27 0 1 2 696 0 30 1 1 2 691 0 30 2 1 3 1023 1 30 1 Column 7 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes eval( 'if exist(''mu'',''var''), mauifunc(mu), end', ''); 50.7 54.964 60.604 62.363 62.799 41.773 48.207 79.88 66.24 71.311 84.285 45.747 45.966 47.686 62.718 59.275 64.9 64.954 57.754 41.224 38.745 46.261 43.088 46.879 51.344 60.332 dbstep feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes mu mu = 50.7 54.964 60.604 62.363 62.799 41.773 48.207 79.88 66.24 71.311 84.285 45.747 45.966 47.686 62.718 59.275 64.9 64.954 57.754 41.224 38.745 46.261 43.088 46.879 51.344 60.332 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbcont » feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes disp(which('poisson_reg')); d:\matlabr11\workspace\sys502\Poisson_reg.m dbclear at 8 in 'd:\matlabr11\workspace\sys502\poisson_reg.m' Poisson_reg feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes b b = 3.3903 -0.068311 0.00033078 0.21243 0.0090362 0.092624 0.13154 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes b = b +.01 b = 3.4003 -0.058311 0.010331 0.22243 0.019036 0.10262 0.14154 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes Poisson_reg feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes b b = -6.3431 -0.080301 0.010196 0.28161 0.021319 0.12695 0.084328 feval('whos') Name Size Bytes Class X 26x7 1456 double array ans 26x1 208 double array b 7x1 56 double array beta 7x1 56 double array data 26x7 1456 double array dev 1x1 8 double array df 1x1 8 double array i 1x1 8 double array mu 26x1 208 double array output 7x2 112 double array se 7x1 56 double array x 26x6 1248 double array y 26x1 208 double array Grand total is 636 elements using 5088 bytes dbquit ??? Error using ==> dbquit Debugging commands are only allowed in debug mode.