File: odds_probs.R

package info (click to toggle)
r-cran-bayesfactor 0.9.12-4.7%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 1,492 kB
  • sloc: cpp: 1,555; sh: 16; makefile: 7
file content (44 lines) | stat: -rw-r--r-- 1,439 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
## ----echo=FALSE,message=FALSE,results='hide'----------------------------------
options(markdown.HTML.stylesheet = 'extra/manual.css')
library(knitr)
options(digits=3)
require(graphics)
set.seed(2)

## ----message=FALSE,results='hide',echo=FALSE----------------------------------
library(BayesFactor)
options(BFprogress = FALSE)
bfversion = BFInfo()
session = sessionInfo()[[1]]
rversion = paste(session$version.string," on ",session$platform,sep="")

## -----------------------------------------------------------------------------
data(puzzles)
bf = anovaBF(RT ~ shape*color + ID, whichRandom = "ID", data = puzzles)
bf

## -----------------------------------------------------------------------------
prior.odds = newPriorOdds(bf, type = "equal")
prior.odds

## -----------------------------------------------------------------------------
priorOdds(prior.odds) <- c(4,3,2,1)
prior.odds

## -----------------------------------------------------------------------------
post.odds = prior.odds * bf
post.odds

## -----------------------------------------------------------------------------
post.prob = as.BFprobability(post.odds)
post.prob

## -----------------------------------------------------------------------------
post.prob / .5

## -----------------------------------------------------------------------------
post.prob[1:3]

## -----------------------------------------------------------------------------
post.prob[1:3] / 1