Thuis
Contacten

    Hoofdpagina


Docent: drs. Rob Flohr

Dovnload 2 Mb.

Docent: drs. Rob Flohr



Pagina2/25
Datum05.12.2018
Grootte2 Mb.

Dovnload 2 Mb.
1   2   3   4   5   6   7   8   9   ...   25

> library(Bolstad)

> pi=c(0.4,0.5,0.6)

> pi.prior=c(0.30,0.40,0.30)

> results=binodp(5,20,pi,pi.prior)

Conditional distribution of x given pi and n:


0 1 2 3 4 5 6 7 8 9 10

0.4 0 5e-04 0.0031 0.0123 0.0350 0.0746 0.1244 0.1659 0.1797 0.1597 0.1171

0.5 0 0e+00 0.0002 0.0011 0.0046 0.0148 0.0370 0.0739 0.1201 0.1602 0.1762

0.6 0 0e+00 0.0000 0.0000 0.0003 0.0013 0.0049 0.0146 0.0355 0.0710 0.1171

11 12 13 14 15 16 17 18 19 20

0.4 0.0710 0.0355 0.0146 0.0049 0.0013 0.0003 0.0000 0.0000 0e+00 0

0.5 0.1602 0.1201 0.0739 0.0370 0.0148 0.0046 0.0011 0.0002 0e+00 0

0.6 0.1597 0.1797 0.1659 0.1244 0.0746 0.0350 0.0123 0.0031 5e-04 0


Joint distribution:
0 1 2 3 4 5 6 7 8 9 10

[1,] 0 1e-04 9e-04 0.0037 0.0105 0.0224 0.0373 0.0498 0.0539 0.0479 0.0351

[2,] 0 0e+00 1e-04 0.0004 0.0018 0.0059 0.0148 0.0296 0.0481 0.0641 0.0705

[3,] 0 0e+00 0e+00 0.0000 0.0001 0.0004 0.0015 0.0044 0.0106 0.0213 0.0351

11 12 13 14 15 16 17 18 19 20

[1,] 0.0213 0.0106 0.0044 0.0015 0.0004 0.0001 0.0000 0e+00 0e+00 0

[2,] 0.0641 0.0481 0.0296 0.0148 0.0059 0.0018 0.0004 1e-04 0e+00 0

[3,] 0.0479 0.0539 0.0498 0.0373 0.0224 0.0105 0.0037 9e-04 1e-04 0


Marginal distribution of x:
0 1 2 3 4 5 6 7 8 9 10

[1,] 0 2e-04 0.001 0.0042 0.0124 0.0287 0.0536 0.0837 0.1126 0.1333 0.1408

11 12 13 14 15 16 17 18 19 20

[1,] 0.1333 0.1126 0.0837 0.0536 0.0287 0.0124 0.0042 0.001 2e-04 0



Prior Likelihood Posterior

0.4 0.3 0.0223941059 0.78037052

0.5 0.4 0.0059143066 0.20609666

0.6 0.3 0.0003883481 0.01353282


> #Continue prior: beta verdeling (10,10)

> x=seq(0,1)

> curve(dbeta(x,10,10))



> binobp(5,20,10,10)

Posterior Mean : 0.375

Posterior Variance : 0.0057165

Posterior Std. Deviation : 0.0756073


Prob. Quantile

------ ---------

0.005 0.1961166

0.01 0.2108954

0.025 0.2336393

0.05 0.2541426

0.5 0.3728969

0.95 0.5030548

0.975 0.5282049

0.99 0.5572392

0.995 0.5768111

>

1   2   3   4   5   6   7   8   9   ...   25

  • > x=seq(0,1) > curve(dbeta(x,10,10))
  • > binobp(5,20,10,10)

  • Dovnload 2 Mb.