#Programs to solve Example 1.1 to Example 1.5

#Example 1.1
q1=1 - pt(2.05, 14)
q2=qt(0.253, 8)
ans=cbind(q1, q2)
cname=c('q1','q2')
print(ans, name.cols=cname)

###############################################################

#Example 1.2
q1=pchisq(11.5, 9)
q2=qchisq(0.085, 13)
ans=cbind(q1, q2)
cname=c('q1','q2')
print(ans, name.cols=cname)

###############################################################

#Example 1.3
q1=1 - pf(5.87, 5, 20)
q2=qf(0.95, 4, 10)
ans=cbind(q1, q2)
cname=c('q1','q2')
print(ans, name.cols=cname)

###############################################################

#Example 1.4
a=dbinom(0, 5, 0.8)
b=dbinom(2, 5, 0.8)
c=1-pbinom(3, 5, 0.8)
d=pbinom(4, 5, 0.8)
e=pbinom(4, 5, 0.8) - pbinom(2, 5, 0.8)
f=1-pbinom(1, 5, 0.8)
ans=cbind(a, b, c, d, e, f)
cname=c('a','b', 'c', 'd', 'e', 'f')
print(ans, name.cols=cname)

###############################################################

#Example 1.5
a=dpois(4, 3.5)
b=ppois(1, 3.5)
c=1-ppois(2, 3.5)
ans=cbind(a, b, c)
cname=c('a','b', 'c')
print(ans, name.cols=cname)