RLAB 5-10

ARKADEEP MUKHERJEE
17BEC1001

LAB5
BINOMIAL DISTRIBUTION

Code:

dbinom(2,20,0.1)
pbinom(1,20,0.1,lower.tail = F)
sum(dbinom(1:3,20,0.10))
x=0:20
y=dbinom(x,size=20,prob=0.10)
plot(x,y,type="h")
dbinom(0:20,20,0.10)

dbinom(0:20,20,0.10)
z=pbinom(0:20,20,0.10)

data.frame(x,y,z)

Graph:
POISSON DISTRIBUTION

QUES 1:

CODE:

x=0:100 p=dpois(x,4)
mean=sum(x*p)
var=sum(x^2*p)-(mean)^2
mean
var

QUES 2:

CODE:

x=0:8
p=dbinom(x,8,0.3)
plot(x,p,type="h",main="Pmf for binomial distribution",col=2)
LAB 6
Arkadeep Mukherjee
17bec1001

x=seq(-10,10,by=0.1)
y=dnorm(x,mean=2.5,sd=0.5)
plot(x,y,type="l")
y1=pnorm(x,mean=2.5,sd=0.5)
plot(x,y1)
x=seq(-3,3,length=200)
y=dnorm(x,mean=0,sd=1)
plot(x,y,type="l")
x=seq(-3,0,length=100)
y=dnorm(x,0,1)
polygon(c(-3,x,0),c(0,y,0),col="red")
pnorm(0,mean=0,sd=1)

x=seq(-3,3,length=200)
y=dnorm(x,mean=0,sd=1)
plot(x,y,type="l")
x=seq(1,2,length=100)
y=dnorm(x,0,1)
polygon(c(1,x,2),c(0,y,0),col="red")
pnorm(2,mean=0,sd=1)-pnorm(1,mean=0,sd=1)
0.1359051

x=seq(-3,3,length=200)
y=dnorm(x,mean=0,sd=1)
plot(x,y,type="l")
qnorm(0.40,mean=0,sd=1)
x=seq(-3,-0.2533,length=100)
y=dnorm(x,mean=0,sd=1)
polygon(c(-3,x,-0.2533),c(0,y,0),col="blue")
text(-1,0.2,"0.40")

-0.2533471
LAB 7
Arkadeep Mukherjee
17BEC1001

1. xbar=9900
muo=10000
sigma=120
n=30
z=(xbar-muo)/(sigma/sqrt(n))
z
alpha=0.05
zalpha=qnorm(alpha)
zalpha
if(z<zalpha) print("reject Ho") else print("accept ho")

xbar=9900
> muo=10000
> sigma=120
> n=30
> z=(xbar-muo)/(sigma/sqrt(n))
> z
[1] -4.564355
> alpha=0.05
> zalpha=qnorm(alpha)
> zalpha
[1] -1.644854
> if(z<zalpha) print("reject Ho") else print("accept ho")
[1] "reject Ho"

2. xbar=2.1
muo=2
sigma=0.25
n=35
z=(xbar-muo)/(sigma/sqrt(n))
z
alpha=0.05
zalpha=qnorm(alpha)
zalpha
if(z>zalpha) print("reject Ho") else print("accept ho")

xbar=2.1
> muo=2
> sigma=0.25
> n=35
> z=(xbar-muo)/(sigma/sqrt(n))
> z
[1] 2.366432
> alpha=0.05
> zalpha=qnorm(alpha)
> zalpha
[1] -1.644854
> if(z>zalpha) print("reject Ho") else print("accept ho")
[1] "reject Ho"

3. xbar=14.6
muo=15.4
sigma=2.5
n=35
z=(xbar-muo)/(sigma/sqrt(n))
z
alpha=0.05
zalpha=qnorm(alpha/2)
zalpha
if(abs(z)<abs(zalpha)) print("reject Ho") else print("accept ho")

xbar=14.6
> muo=15.4
> sigma=2.5
> n=35
> z=(xbar-muo)/(sigma/sqrt(n))
> z
[1] -1.893146
> alpha=0.05
> zalpha=qnorm(alpha/2)
> zalpha
[1] -1.959964
> if(abs(z)<abs(zalpha)) print("reject Ho") else print("accept ho")
[1] "reject Ho"

>

4. p=85/148
P=60/100
n=148
Q=1-P
z=(p-P)/sqrt(P*Q/n)
z
alpha=0.05
zalpha=qnorm(alpha)
if(z<zalpha) print("reject Ho") else print("accept Ho")
p=85/148
> P=60/100
> n=148
> Q=1-P
> z=(p-P)/sqrt(P*Q/n)
> z
[1] -0.6375983
> alpha=0.05
> zalpha=qnorm(alpha)
> if(z<zalpha) print("reject Ho") else print("accept Ho")
[1] "accept Ho"
 >

5. p=30/214
P=12/100
n=214
Q=1-P
z=(p-P)/sqrt(P*Q/n)
z
alpha=0.05
zalpha=qnorm(1-alpha)
if(z<zalpha) print("reject Ho") else print("accept Ho")
p=30/214
> P=12/100
> n=214
> Q=1-P
> z=(p-P)/sqrt(P*Q/n)
> z
[1] 0.908751
> alpha=0.05
> zalpha=qnorm(1-alpha)
> if(z<zalpha) print("reject Ho") else print("accept Ho")
[1] "reject Ho"

>

6. p1=0.78
p2=0.75
n1=150
n2=120
alpha=0.05
P=(n1*p1+n2*p2)/(n1+n2)
Q=1-P
SE=sqrt((P*Q/n1)+(P*Q/n2))
zcal=(p1-p2)/SE
zcal
ztab=qnorm(1-alpha)
ztab
if(z<ztab) print("accept Ho") else print("reject Ho")

p1=0.78
> p2=0.75
> n1=150
> n2=120
> alpha=0.05
> P=(n1*p1+n2*p2)/(n1+n2)
> Q=1-P
> SE=sqrt((P*Q/n1)+(P*Q/n2))
> zcal=(p1-p2)/SE
> zcal
[1] 0.5791405
> ztab=qnorm(1-alpha)
> ztab
[1] 1.644854
> if(z<ztab) print("accept Ho") else print("reject Ho")
[1] "accept Ho"
7.p1=45/80
p2=56/103
n1=80
n2=103
alpha=0.01
P=(n1*p1+n2*p2)/(n1+n2)
Q=1-P
SE=sqrt((P*Q/n1)+(P*Q/n2))
zcal=(p1-p2)/SE
zcal
ztab=qnorm(1-alpha)
ztab
if(z<ztab) print("accept Ho") else print("reject Ho")

p1=45/80
> p2=56/103
> n1=80
> n2=103
> alpha=0.01
> P=(n1*p1+n2*p2)/(n1+n2)
> Q=1-P
> SE=sqrt((P*Q/n1)+(P*Q/n2))
> zcal=(p1-p2)/SE
> zcal
[1] 0.2538201
> ztab=qnorm(1-alpha)
> ztab
[1] 2.326348
> if(z<ztab) print("accept Ho") else print("reject Ho")
[1] "accept Ho"
LAB 8
Arkadeep Mukherjee
17BEC1001
1.
x=c(0.593,0.142,0.529,0.691,0.231,0.793,0.519,0.392,0.418)
xbar=mean(x)
alpha=0.05
mu=0.3
sd=sqrt(var(x))
n=length(x)
s=sd*sqrt(n/(n-1))
t=(xbar-mu)/(s/sqrt(n))
tv=qt(1-alpha,df=n-1)
if(t<tv)print("Accept H0") else print("Reject H0")

OUTPUT:
Reject H0

2.
x1=c(175,168,168,190,156,181,182,175,174,179)
x2=c(185,169,173,173,188,186,175,174,179,180)
alpha=0.05
n1=length(x1)
n2=length(x2)
x1bar=mean(x1)
x2bar=mean(x2)
sd1=sqrt(var(x1))
sd2=sqrt(var(x2))
s=sqrt(n1*(sd1)^2+n2*(sd2)^2)/(n1+n2-2)
t=abs(x1bar-x2bar)/s*(sqrt((1/n1)+(1/n2)))
tv=qt(1-(alpha/2),n1+n2-2)
if(t<tv)print("Accept H0") else print("Reject H0")
OUTPUT: Accept H0
3.
A=c(14.1,10.1,14.7,13.7,14.0)
B=c(14.0,14.5,13.7,12.7,14.1)
alpha=0.1
n1=length(A)
n2=length(B)
S1square=(var(A)*n1)/(n1-1)
S2square=(var(B)*n2)/(n2-1)
F=S1square/S2square
TV=qf(1-alpha/2,n1-1,n2-1)
if(F<TV)print("Accept H0") else print("Reject H0")
OUTPUT:
Reject H0

LAB 9
Arkadeep Mukherjee
17BEC1001
CHI SQUARE TEST
1.
m=2
n=3
data<-matrix(c(35,42,61,48,51,68),ncol=n,byrow=T)
alpha=0.05
v1=chisq.test(data)
v2=qchisq(1-alpha,(m-1)*(n-1))
if(v1[1]<v2)print("Accept H0")else print("Reject H0")
OUTPUT:
Accept H0

2.
n=5
alpha=0.05
N=320
P<-0.5
x=c(0:n)
obf<-c(14,56,110,88,40,12)
exf<-(dbinom(x,n,P)*N)
chisq<-sum((obf-exf)^2/exf)
cv=chisq
tv=qchisq(1-alpha,n-1)
if(cv<=tv)print("Accept H0") else print("Reject H0")
Output:
Accept H0
LAB 10
Arkadeep Mukherjee
17 BEC1001
ANOVA
1. group1=c(551,457,450,731,499,632)
group2=c(595,580,508,583,633,517)
group3=c(639,615,511,573,648,677)
group4=c(417,449,517,438,415,555)
group5=c(563,631,522,613,656,679)
group=data.frame(cbind(group1,group2,group3,group4,group5))
summary(group)
stgr=stack(group)
crd=aov(values~ind,data=stgr)
summary(crd)
boxplot(group,ylab="Aggregates",main="Absorption of Moisture")

output:
> group1=c(551,457,450,731,499,632)
> group2=c(595,580,508,583,633,517)
> group3=c(639,615,511,573,648,677)
> group4=c(417,449,517,438,415,555)
> group5=c(563,631,522,613,656,679)
> group=data.frame(cbind(group1,group2,group3,group4,group5))
> summary(group)
group1 group2 group3 group4 group5
Min. :450.0 Min. :508.0 Min. :511.0 Min. :415.0 Min. :522.
0
1st Qu.:467.5 1st Qu.:532.8 1st Qu.:583.5 1st Qu.:422.2 1st Qu.:575.
5
Median :525.0 Median :581.5 Median :627.0 Median :443.5 Median :622.
0
Mean :553.3 Mean :569.3 Mean :610.5 Mean :465.2 Mean :610.
7
3rd Qu.:611.8 3rd Qu.:592.0 3rd Qu.:645.8 3rd Qu.:500.0 3rd Qu.:649.
8
Max. :731.0 Max. :633.0 Max. :677.0 Max. :555.0 Max. :679.
0
> stgr=stack(group)
> crd=aov(values~ind,data=stgr)
> summary(crd)
Df Sum Sq Mean Sq F value Pr(>F)
ind 4 85356 21339 4.302 0.00875 **
Residuals 25 124020 4961
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> boxplot(group,ylab="Aggregates",main="Absorption of Moisture")
2. Oper1 Oper2 Oper3 Oper4 Oper5 Oper6
M1 42.5 39.3 39.6 39.9 42.9 43.6
M2 39.8 40.1 40.5 42.3 42.5 43.1
M3 40.2 40.5 41.3 43.4 44.9 45.1
M4 41.3 42.2 43.5 44.2 45.9 42.3

data=read.table(file.choose(),header=TRUE)
time=c(t(as.matrix(data)))
f=c("Oper1","Oper2","Oper3","Oper4","Oper5")
g=c("M1","M2","M3","M4")
k=ncol(data)
n=nrow(data)
Operators=gl(k,1,n*k,factor(f))
Machines=gl(n,k,n*k,factor(g))
anova=aov(time~Machines+Operators)
summary(anova)
interaction.plot(Operators,Machines,time)
par(mfrow=c(1,2))
plot(time~Machines+Operators,main="Product time")