?? rats_demo.m.svn-base
字號:
%
% 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');
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -