wsf[^̌vZ_fO\wKfƂā\xRR[hCStanR[hłB
(2018/09/24 )

͂߂-------------------------------------------------------------------

EŐV̏CɂĂ
  https://kkatahira.github.io/cmbd-book/index.html
@QƂĂBȉ̏͂̃TCg甲̂łB

ERR[hRStudioœƂz肵Ă܂B

EKvȃpbP[WCXg[ĂȂꍇ́CCXg[ĎgpĂBႦ΁ftidyversefgpKvƂ́C 
> install.packages('tidyverse') 
RStudiõR}hCőłĂ΃CXg[邱Ƃł܂B tidyverséCقƂǂ̃R[hŃOt̍쐬ȂǂɎgp̂Ȃ̂ŁC ̃R}hŃCXg[ĂĂB܂C gridExtra}ׂۂɂ悭g܂B

E2͂̃t@C̓fBNg (tH_) gchapter2hɓĂCƂ悤ɁC{̏͂ƂɃfBNgĂ܂B

EfBNgɂ鑼̃t@CǂݍނƂ̂ŁC{Iɂ͂̃fBNg\̂܂܁CƃfBNg̃t@ĈfBNgڂēĂBƃfBNg݊JĂt@ĈfBNgɂɂ́CRStudiõj[Session -> Set Working Directory -> To Source File LocationI܂B

Eɑ2͂Ƒ4͂ł́Cǂ̃R[hł邾̃t@CɎ߂悤ɂ܂B {̓V~[Vɂf[^̐ƁCp[^Cʂ̕`Ȃǂ͕悢Ǝv̂łCw҂ɂƂĂ͑璷łĂ̃R[hɂ܂Ƃ܂ĂƂ肩₷ƍlłB̂߁CvO~Oł͖{]܂ȂƂRs[&y[XgĂӏ܂BB


2-------------------------------------------------------------------

RW_simulation.R (}2.2̍쐬Ɏgp)
RWf̃V~[V܂B

plot_softmax.R (}2.4̍쐬Ɏgp)
\tg}bNX֐vbg܂B

qlearning_simulation.R (}2.5̍쐬Ɏgp)
QwKf̃V~[VɂIf[^𐶐C̃f[^ɑ΂ŖސɂQwKf̃p[^܂B

qlearning_plot_ll_function.R (}2.6̍쐬Ɏgp)
}2.5Ɠf[^𐶐Cɑ΂QwKf̑ΐޓx֐vbg܂ (R[h̑Oqlearning_simulation.RƓł)B

qlearning_model_comparison.R (}2.7̍쐬Ɏgp)
}2.5Ɠf[^𐶐Cɑ΂QwKfCWSLSfC_f̍ŖސsCfr܂B

4-------------------------------------------------------------------

4͈ȍ~̃R[hɂẮCMCMC@pꍇ́CRStañCXg[KvɂȂ܂BRStañCXg[@͗Ⴆ΂QlɂĂB

Bernoulli_bayes.R (}4.1̍쐬Ɏgp)
xk[Cs (RCgX̗) ̃xCYB

qlearning_simulation_Bayes_single.R (}4.2̍쐬Ɏgp)
QwKf̃V~[Vɂl̑If[^𐶐C̃f[^xCYCMAP܂ (r̂ߍŖސs)B
StanR[hqlearning_single_subject.stanĂяo܂B

qlearning_groupdata.R (}4.5̍쐬Ɏgp)
QwKf̃V~[Vɂ40l̑If[^𐶐C KwxCY@CŒʕ́ClxMAPC lx̍Ŗސs܂B
StanR[hłqlearning_group.stanĂяo܂B

5-------------------------------------------------------------------

5͂́Cf[^̐R[hƗCf[^data.frameƂāCcsvt@CɂɊi[܂B ܂CŖސMAP̂߂̍œKɂoptim̑ւɁC萫\̂悢RsolnppbP[Wsolnpp܂B

KvɉĈȉ̍ƂsĂB
EƃfBNgɁfdatafƂf[^i[pfBNgB

EƃfBNgɁfresultsfƂf[^i[pfBNgB

E> install.packages('Rsolnp')@RsolnppbP[WCXg[B

[f[^]
generate_group_data_fqlearning.R
f[^Ɏg{R[hłB̃R[hsƁC40l̃V~[Vɂf[^C./data/ ȉcsvt@Co͂܂B generate_data_fqlearning_single_subject.R ͂̃R[h̐ݒ蕔ύXāC1l̃f[^𐶐悤ɂ̂łB

[p[^CfI]
model_selection_single_subject.R (\5.1̍쐬Ɏgp)
l̎Q҂̃f[^̃fĨfBgenerate_group_data_fqlearning.Rō쐬f[^ǂݍ݁CẽfI܂B
f֐i[model_functions_FQ.R, ŖސEMAPɗp֐`parameter_fit_functions.RɉC 6StanR[h ( smodel_qlearning_single.stan, smodel_fqlearning_single.stan, smodel_dfqlearning_single.stan, smodel_qlearning_single_WBIC.stan, smodel_fqlearning_single_WBIC.stan, smodel_dfqlearning_single_WBIC.stan) ǂ݂݂܂BStanR[h̍ŏ3͒ʏ̃xCY肨WAIČvZpC㔼3WBIČvẐ߂̂̂łBWBICvZꍇt@C𕪂Ȃ@̂łCł͂ꂼ̃R[h̗₷dĕʂ̃t@Cɂ܂Bmodel_selection_group.Rł́Ĉ悤ɂ邱Ƃłꂼ̃fStanR[hɂ܂Ƃ߂Ă܂B
܂C3̃f̓lXg֌Wɂ (WIQwKfCFwKf͂ꂼDF-QwKf̓`ɂȂĂ) ߁CstanR[h͈ɂ܂Ƃ߂邱Ƃ\łBmodel_functions_FQ.RlɊȗ邱Ƃ\łB

model_selection_group.R (}5.4̍쐬Ɏgp)
Wcf[^̃fĨfBgenerate_group_data_fqlearning.Rō쐬f[^ǂݍ݁CẽfI܂B
ŖސCMAP̂parameter_fit_functions_group.Rǂݍ݂܂B
StanR[hCsmodel_qlearning_group.stan, smodel_fqlearning_group.stan, smodel_dfqlearning_group.stanǂݍ݂܂B
_ʃfI@sɂMatlabSPMKvɂȂ܂B Rł̎\Ǝv̂łC̓_͉߂ĒǋL܂B

[̑]

plot_log_posterior_fqlearning.R (}5.2̍쐬Ɏgp)
Ӗޓx̐̂߂̑ΐmxvbgB

fqlearning_laplace_approx.R (}5.3̍쐬Ɏgp)
txŒ肵F-QwKfpvXߎ̐B

6-------------------------------------------------------------------

y6-1 (QԔr)z

[f[^]
generate_data_qlearning_group_comparison.R (}6.2̍쐬Ɏgp)
6-1̃f[^ɎgR[hłB̃R[hsƁC40l̃V~[Vɂf[^C./data/ ȉcsvt@Co͂܂B

[ftBbeBO]
qlearning_group_comparison_ML_MAP.R (6.2.1̌vZɎgp)
6-2̃f[^ɑ΂C6.2.1 (p[^̓_lpQԔr), 6.2.2 (Œʃf̃fIɂQԔr)ŏЉ͂s܂B
model_functions_Q_mulitiple_group.R, parameter_fit_functions_multiple_group.Rǂݍ݂܂B

qlearning_group_comparison_HB.R (6.2.3̌vZ; }6.3̍쐬Ɏgp)
6-2̃f[^KwxCY@ŉ͂܂B
smodel_qlearning_multiple_group.stang܂B

y6-2 (f[^Ƃ̑֕)z

[f[^]
generate_data_qlearning_trait_random_slope.R (}6.5̍쐬Ɏgp)
6-2̃f[^ɎgR[hłB

[ftBbeBO]
qlearning_trait_ML_MAP.R
6-2̃f[^ŖސCMAPŉ͂܂B
model_functions_Q_mulitiple_group.R, parameter_fit_functions_multiple_group.Rǂݍ݂܂B

qlearning_trait_HB.R
6-2̃f[^KwxCY@ŉ͂܂B
StanR[hCsmodel_qlearning_trait.stan, smodel_qlearning_trait_random_slope.stang܂B
