
%
% rats example from WinBUGS
%

dataStruct = struct( ...
    'x', [8.0, 15.0, 22.0, 29.0, 36.0], ...
    'xbar', 22, ...
    'N', 30, ...
    'T', 5, ...
    'Y', [ 151, 199, 246, 283, 320; ...
	   145, 199, 249, 293, 354; ...
	   147, 214, 263, 312, 328; ...
	   155, 200, 237, 272, 297; ...
	   135, 188, 230, 280, 323; ...
	   159, 210, 252, 298, 331; ...
	   141, 189, 231, 275, 305; ...
	   159, 201, 248, 297, 338; ...
	   177, 236, 285, 350, 376; ...
	   134, 182, 220, 260, 296; ...
	   160, 208, 261, 313, 352; ...
	   143, 188, 220, 273, 314; ...
	   154, 200, 244, 289, 325; ...
	   171, 221, 270, 326, 358; ...
	   163, 216, 242, 281, 312; ...
	   160, 207, 248, 288, 324; ...
	   142, 187, 234, 280, 316; ...
	   156, 203, 243, 283, 317; ...
	   157, 212, 259, 307, 336; ...
	   152, 203, 246, 286, 321; ...
	   154, 205, 253, 298, 334; ...
	   139, 190, 225, 267, 302; ...
	   146, 191, 229, 272, 302; ...
	   157, 211, 250, 285, 323; ...
	   132, 185, 237, 286, 331; ...
	   160, 207, 257, 303, 345; ...
	   169, 216, 261, 295, 333; ...
	   157, 205, 248, 289, 316; ...
	   137, 180, 219, 258, 291; ...
	   153, 200, 244, 286, 324 ] );    

init0 = struct( ...
    'alpha', [ 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, ...
    	       250, 250, 250, 250, 250, 250, 250, 250, 250, 250, 250, ...
	       250, 250, 250, 250, 250, 250, 250, 250 ], ...
    'beta', [ 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, ...
    	      6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6 ], ...
    'alpha_c', 150, ...
    'beta_c', 10, ...
    'tau_c', 1, ...
    'alpha_tau', 1, ...
    'beta_tau', 1 );

[samples, stats, structArray] = matbugs(dataStruct, ...
		fullfile(pwd, 'rats_model.txt'), ...
		'init', init0, ...
                'nChains', 1, ...
		'view', 1, 'nburnin', 1000, 'nsamples', 500, ...
		'thin', 10, ...
		'monitorParams', {'alpha_c', 'alpha_tau', 'alpha0', ...
	       			'beta_c', 'beta_tau'}, ...
                'Bugdir', 'C:/Program Files/WinBUGS14');

