symulacje

Menu główne

Start
Kody podstawowe
- Informacja
- - - - - - - - - -
- Rozdział 2
  - 2.1
  - 2.2
  - 2.3
- Rozdział 3
  - 3.1
  - 3.2
  - 3.3
  - 3.4
  - 3.5
  - 3.6
- Rozdział 4
  - 4.1
  - 4.2
  - 4.3
  - 4.4
  - 4.5
  - 4.6
  - 4.7
  - 4.8
  - 4.9
  - 4.10
  - 4.11
  - 4.12
  - 4.13
  - 4.14
  - 4.15
  - 4.16
  - 4.17
  - 4.18
- Rozdział 5
  - 5.1
  - 5.2
  - 5.3
  - 5.4
  - 5.5
  - 5.6
  - 5.7
  - 5.8
  - 5.9
  - 5.10
  - 5.11
  - 5.12
  - 5.13
  - 5.14
  - 5.15
  - 5.16
- Rozdział 6
  - 6.1
  - 6.2
  - 6.3
  - 6.4
  - 6.5
  - 6.6
  - 6.7
  - 6.8
  - 6.9
  - 6.10
  - 6.11
  - 6.12
  - 6.13
  - 6.14
  - 6.15
  - 6.16
  - 6.17
- Rozdział 7
  - 7.1
  - 7.2
  - 7.3
  - 7.4
  - 7.5
  - 7.6
- Rozdział 8
  - 8.1
  - 8.2
  - 8.3
  - 8.4
  - 8.5
  - 8.6
  - 8.7
  - 8.8
  - 8.9
  - 8.10
  - 8.11
  - 8.12
  - 8.13
  - 8.14
  - 8.15
Kody rozszerzone
- Informacja
- - - - - - - - - -
- Rozdział 2
  - 2.1
  - 2.2
  - 2.3
- Rozdział 3
  - 3.1
  - 3.2
  - 3.3
  - 3.4
  - 3.5
  - 3.6
- Rozdział 4
  - 4.1
  - 4.2
  - 4.3
  - 4.4
  - 4.5
  - 4.6
  - 4.7
  - 4.8
  - 4.9
  - 4.10
  - 4.11
  - 4.12
  - 4.13
  - 4.14
  - 4.15
  - 4.16
  - 4.17
  - 4.18
- Rozdział 5
  - 5.1
  - 5.2
  - 5.3
  - 5.4
  - 5.5
  - 5.6
  - 5.7
  - 5.8
  - 5.9
  - 5.10
  - 5.11
  - 5.12
  - 5.13
  - 5.14
  - 5.15
  - 5.16
- Rozdział 6
  - 6.1
  - 6.2
  - 6.3
  - 6.4
  - 6.5
  - 6.6
  - 6.7
  - 6.8
  - 6.9
  - 6.10
  - 6.11
  - 6.12
  - 6.13
  - 6.14
  - 6.15
  - 6.16
  - 6.17
- Rozdział 7
  - 7.1
  - 7.2
  - 7.3
  - 7.4
  - 7.5
  - 7.6
- Rozdział 8
  - 8.1
  - 8.2
  - 8.3
  - 8.4
  - 8.5
  - 8.6
  - 8.7
  - 8.8
  - 8.9
  - 8.10
  - 8.11
  - 8.12
  - 8.13
  - 8.14
  - 8.15
Linki
- Odnośniki do stron www
Download
- Download

6.12

Kody podstawowe > Rozdział 6

Generowanie wartości losowych z rozkładu chi-kwadrat z wykorzystaniem algorytmu Metropolisa-Hastingsa

n=1000
x=c()
x[1]=2
a=1
for (i in 1:n)
{

s=2*a*runif(1)-a
x[i+1]=x[i]+s
while (dchisq(x[i+1],3)/dchisq(x[i],3)<runif(1)) { s=2*a*runif(1)-a; x[i+1]=x[i]+s}

}
par(mfrow=c(1,2))
plot(x,type="l")
hist(x)