Scribd
Upload
50 views12 pages

Experiment - 1 Mat2001: Name Suyash Yadav Reg. 19BME0320

1. The document describes an experiment conducted by Suyash Yadav to analyze the distribution of data. Various statistical measures like mean, median, mode, variance, standard deviation, and coefficient of variation are calculated. 2. Additional statistical analysis is performed to determine skewness, kurtosis, and whether the distribution is platykurtic or not based on the beta and gamma coefficients. The distribution is found to be negatively skewed and platykurtic. 3. The key results are: Mean = 208.9831, Median = 209.7826, Mode = 216.1538, Variance = 601.928, Standard Deviation = 24.53423, Coefficient of

Uploaded by

SUYASH yadav
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
50 views12 pages

Experiment - 1 Mat2001: Name Suyash Yadav Reg. 19BME0320

1. The document describes an experiment conducted by Suyash Yadav to analyze the distribution of data. Various statistical measures like mean, median, mode, variance, standard deviation, and coefficient of variation are calculated. 2. Additional statistical analysis is performed to determine skewness, kurtosis, and whether the distribution is platykurtic or not based on the beta and gamma coefficients. The distribution is found to be negatively skewed and platykurtic. 3. The key results are: Mean = 208.9831, Median = 209.7826, Mode = 216.1538, Variance = 601.928, Standard Deviation = 24.53423, Coefficient of

Uploaded by

SUYASH yadav
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

Experiment -1

Mat2001

Name= SUYASH YADAV


Reg.=19BME0320
R CODE

> mid=seq(175,245,10)

> mid

[1] 175 185 195 205 215 225 235 245

> f=c(52,68,85,92,100,95,70,28)

> fr.distr=data.frame(mid,f)

> fr.distr
mid f
1 175 52 2 185
68 3 195 85 4
205 92

5 215 100

6 225 95

7 235 70

8 245 28

> mean=(sum(mid*f))/sum(f)

> mean

[1] 208.9831

> midx=seq(175,245,10)

> frequency=c(52,68,85,92,100,95,70,28) > fr.dist<-data.frame(midx,frequency)

> fr.dist

midx frequency

1 175 52 2 185
68 3 195 85 4
205 92

5 215 100

6 225 95 7
235 70 8
245 28
> cl=cumsum(frequency)

> cl

[1] 52 120 205 297 397 492 562 590

> n=sum(frequency)

>n

[1] 590

> ml=min(which(cl>=n/2))

> ml

[1] 4

> h=10

>h

[1] 10

> f=frequency[ml]
>f

[1] 92

> c=cl[ml-1]

>c

[1] 205

> l=mid[ml]-h/2

>l

[1] 200

> median=l+(((n/2)-c)/f)*h

> median
[1] 209.7826

> m=which(frequency==max(frequency))

>m

[1] 5

> fm=frequency[m]

> fm

[1] 100

> f1=frequency[m-1]

> f2=frequency[m+1]

> f1

[1] 92

> f2

[1] 95

> l=midx[m]-h/2
>l

[1] 210

> mode=l+((fm-f1)/(2*fm-f1-f2))*h

> mode

[1] 216.1538

> x=var(f)

>x

[1] 601.928

> sd=sqrt(var(f))
> sd

[1] 24.53423

> cv=100*(sd/mean)

> cv

[1] 11.73982

OUTPUT
MEAN=208.9831

MEDIAN=209.7826

MODE=216.1538

VARIANCE=601.928

STANDARD DEVIATION=24.53423

COEFFICIENT OF VARIATION=11.73982
R CODE
> mid=seq(175,245,10)

> mid

[1] 175 185 195 205 215 225 235 245

> f=c(52,68,85,92,100,95,70,28)

> fr.distr=data.frame(mid,f)

> fr.distr mid f

1 175 52
2 185 68
3 195 85
4 205 92
5 215 100

6 225 95
7 235 70
8 245 28
> y=rep(mid,f)

>y

[1] 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175
175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175 175
175 175 175 175 175 175

[46] 175 175 175 175 175 175 175 185 185 185 185 185 185 185 185 185 185 185 185
185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185
185 185 185 185 185 185

[91] 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185 185
185 185 185 185 185 185 185 185 185 185 185 195 195 195 195 195 195 195 195 195
195 195 195 195 195 195

[136] 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195
195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195
195 195 195 195 195 195

[181] 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195 195
195 195 195 195 195 195 205 205 205 205 205 205 205 205 205 205 205 205 205 205
205 205 205 205 205 205

[226] 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205
205 205 205 205
205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205
205 205
[271] 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205 205
205 205 205 205 205 205 205 205 215 215 215 215 215 215 215 215 215 215 215 215
215 215 215 215 215 215

[316] 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215
215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215
215 215 215 215 215 215
[361] 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215
215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 215 225 225
225 225 225 225 225 225

[406] 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225
225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225
225 225 225 225 225 225

[451] 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225
225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225 225
225 225 225 235 235 235

[496] 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235
235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235
235 235 235 235 235 235

[541] 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235
235 235 235 245 245 245 245 245 245 245 245 245 245 245 245 245 245 245 245 245
245 245 245 245 245 245

[586] 245 245 245 245 245

> n=length(y)

>n

[1] 590

> mean(y)

[1] 208.9831

> m3=sum((y-mean)^3)/n

> m3

[1] -447.3203

> m2=sum((y-mean)^2)/n

> m2
[1] 388.0336

> beta1=((m3)^2/(m2)^3

+)

> beta1=((m3)^2/(m2)^3)

> beta1
[1] 0.003424753

> m4=sum((y-mean)^4)/n

> m4

[1] 306260.9

> beta2=(m4/(m2)^2)

> beta2

[1] 2.034009

> gam2=beta2-3

> gam2

[1] -0.9659908

OUTPUT
Skewness=0.003424753

Beta coefficient of kurtosis (2)= 2.034009

)
Fisher’s Gamma coefficient of kurtosis ( 2  )= ] -0.9659908

  0 negatively skewed distribution

2<3  Curve is a platykurtic curve

 < 0  platykurtic curve

You might also like