SS2012-StatisticalMethods.git

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
45
46
47
48
49
50
51
library(LearnBayes)
n=seq(0,100, by = 1) 
p=dpois(n,5)
plot(n,p,type='l',main='Poisson',xlab='n',ylab='f(n,Np)',col='blue')
p=dpois(n,10)
points(n,p,type='l',col='red')
p=dpois(n,50)
points(n,p,type='l',col='green')
legend('topright',legend=c('N=5','N=10','N=50'),lty=c(1,1,1),col=c('blue','red','green'))
####
c=seq(0,10,by=0.01)
e=dexp(c,1./0.5)
plot(c,e,type='l',main='Exponential Distribution',xlab='x',ylab='f(x,xi)',col=blue)
e=dexp(c,1./2)
points(c,e,type='l',col='red')
e=dexp(c,1./4)
points(c,e,type='l',col='green')
legend('topright',legend=c('0.5','2','4'),lty=c(1,1,1),col=c('blue','red','green'))
####
n2=seq(0,20, by = 1)
b=dbinom(n2,5,0.5)
plot(n2,b,type='h',main='Binomial Distribution',xlab='n',ylab='f(n;N,p)',col='black')
b=dbinom(n2,10,0.5)
points(n2+0.1,b,type='h',col='red')
b=dbinom(n2,20,0.5)
points(n2+0.2,b,type='h',col='green')
b=dbinom(n2,20,0.8)
points(n2+0.3,b,type='h',col='blue')
b=dbinom(n2,20,0.2)
points(n2+0.4,b,type='h',col='cyan')
legend('topright',legend=c('(5,0.5)','(10,0.5)','(20,0.5)','(20,0.8)','(20,0.2)'),lty=c(1,1,1,1,1),col=c('black','red','green','blue','cyan')) #To fix
####
x3=seq(0,10,by=0.01)
l=dlnorm(x3,0,0.5)
plot(x3,l,type='l',main='Log-Normal Dist.',xlab='x',ylab='f(x;mu,sigma)')
l=dlnorm(x3,0,1)
points(x3,l,type='l',col='green')
l=dlnorm(x3,2,1)
points(x3,l,type='l',col='red')
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX