R Console Page 1
> #########8TH RECORD - BADM##########
>
> #Analyze the sample data with one sample and two sample test for th
e odinal and scale data
> #odinal test- 1) one sample( wilcoxon signed rank test)
> #odinal test- 1) two sample( mann whitney u test)
> #scale test- 1) one sample( one sample t test )
> #scale test- 1) two sample( two sample t test )
>
> #TEST-01
> mtcars
mpg cyl disp hp drat wt qsec vs am gear ca
rb
Mazda RX4 21.0 6 160.0 110 3.90 2.620 16.46 0 1 4
4
Mazda RX4 Wag 21.0 6 160.0 110 3.90 2.875 17.02 0 1 4
4
Datsun 710 22.8 4 108.0 93 3.85 2.320 18.61 1 1 4
1
Hornet 4 Drive 21.4 6 258.0 110 3.08 3.215 19.44 1 0 3
1
Hornet Sportabout 18.7 8 360.0 175 3.15 3.440 17.02 0 0 3
2
Valiant 18.1 6 225.0 105 2.76 3.460 20.22 1 0 3
1
Duster 360 14.3 8 360.0 245 3.21 3.570 15.84 0 0 3
4
Merc 240D 24.4 4 146.7 62 3.69 3.190 20.00 1 0 4
2
Merc 230 22.8 4 140.8 95 3.92 3.150 22.90 1 0 4
2
Merc 280 19.2 6 167.6 123 3.92 3.440 18.30 1 0 4
4
Merc 280C 17.8 6 167.6 123 3.92 3.440 18.90 1 0 4
4
Merc 450SE 16.4 8 275.8 180 3.07 4.070 17.40 0 0 3
3
Merc 450SL 17.3 8 275.8 180 3.07 3.730 17.60 0 0 3
3
Merc 450SLC 15.2 8 275.8 180 3.07 3.780 18.00 0 0 3
3
Cadillac Fleetwood 10.4 8 472.0 205 2.93 5.250 17.98 0 0 3
4
Lincoln Continental 10.4 8 460.0 215 3.00 5.424 17.82 0 0 3
4
Chrysler Imperial 14.7 8 440.0 230 3.23 5.345 17.42 0 0 3
4
Fiat 128 32.4 4 78.7 66 4.08 2.200 19.47 1 1 4
R Console Page 2
1
Honda Civic 30.4 4 75.7 52 4.93 1.615 18.52 1 1 4
2
Toyota Corolla 33.9 4 71.1 65 4.22 1.835 19.90 1 1 4
1
Toyota Corona 21.5 4 120.1 97 3.70 2.465 20.01 1 0 3
1
Dodge Challenger 15.5 8 318.0 150 2.76 3.520 16.87 0 0 3
2
AMC Javelin 15.2 8 304.0 150 3.15 3.435 17.30 0 0 3
2
Camaro Z28 13.3 8 350.0 245 3.73 3.840 15.41 0 0 3
4
Pontiac Firebird 19.2 8 400.0 175 3.08 3.845 17.05 0 0 3
2
Fiat X1-9 27.3 4 79.0 66 4.08 1.935 18.90 1 1 4
1
Porsche 914-2 26.0 4 120.3 91 4.43 2.140 16.70 0 1 5
2
Lotus Europa 30.4 4 95.1 113 3.77 1.513 16.90 1 1 5
2
Ford Pantera L 15.8 8 351.0 264 4.22 3.170 14.50 0 1 5
4
Ferrari Dino 19.7 6 145.0 175 3.62 2.770 15.50 0 1 5
6
Maserati Bora 15.0 8 301.0 335 3.54 3.570 14.60 0 1 5
8
Volvo 142E 21.4 4 121.0 109 4.11 2.780 18.60 1 1 4
2
> attach(mtcars)
The following objects are masked from mtcars (pos = 3):
am, carb, cyl, disp, drat, gear, hp, mpg, qsec, vs, wt
> mpg
[1] 21.0 21.0 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 17.8 16.4 17.3
15.2 10.4 10.4 14.7 32.4 30.4 33.9 21.5
[22] 15.5 15.2 13.3 19.2 27.3 26.0 30.4 15.8 19.7 15.0 21.4
> #HO= Data has normality
> #H1= Data has No normality
> [Link](mpg)
Shapiro-Wilk normality test
data: mpg
W = 0.94756, p-value = 0.1229
> #hence the p value is more than 0.05 , Accept H0 the data is norma
R Console Page 3
l
> #Hence the data is Normal we use the one sample t test , one sam
ple data
> mpg
[1] 21.0 21.0 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 17.8 16.4 17.3
15.2 10.4 10.4 14.7 32.4 30.4 33.9 21.5
[22] 15.5 15.2 13.3 19.2 27.3 26.0 30.4 15.8 19.7 15.0 21.4
> #HO: >=22
> #H1: !=22
> [Link](mpg,mu=22,alternative="[Link]")
One Sample t-test
data: mpg
t = -1.7921, df = 31, p-value = 0.08288
alternative hypothesis: true mean is not equal to 22
95 percent confidence interval:
17.91768 22.26357
sample estimates:
mean of x
20.09062
> #Hence the p value > 0.05 , Accept H0
>
>
> #TEST-02
> hp
[1] 110 110 93 110 175 105 245 62 95 123 123 180 180 180 205 215
230 66 52 65 97 150 150 245 175 66
[27] 91 113 264 175 335 109
> #testing normality
> #H0: THE DATA IS NORMAL
> #H1: THE DATA IS NOT NORMAL
> [Link](hp)
Shapiro-Wilk normality test
data: hp
W = 0.93342, p-value = 0.04881
> #hence the p value is less than o.o5 ,Accept the H1 , the data has
no normality
> #hence the data is not normal we use the wilcoxon signed rank test
( one sample data)
> hp
[1] 110 110 93 110 175 105 245 62 95 123 123 180 180 180 205 215
230 66 52 65 97 150 150 245 175 66
[27] 91 113 264 175 335 109
R Console Page 4
> #H0; hp<200
> #H1; hp !=200
> [Link](hp,mu=12,alternative="[Link]")
Wilcoxon signed rank test with continuity correction
data: hp
V = 528, p-value = 8.27e-07
alternative hypothesis: true location is not equal to 12
Warning message:
In [Link](hp, mu = 12, alternative = "[Link]") :
cannot compute exact p-value with ties
> #Hence the p value lesser than 0.05 , Accept H0
>
>
>
> #TEST-03
> vs
[1] 0 0 1 1 0 1 0 1 1 1 1 0 0 0 0 0 0 1 1 1 1 0 0 0 0 1 0 1 0 0 0 1
> mpg
[1] 21.0 21.0 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 17.8 16.4 17.3
15.2 10.4 10.4 14.7 32.4 30.4 33.9 21.5
[22] 15.5 15.2 13.3 19.2 27.3 26.0 30.4 15.8 19.7 15.0 21.4
> #testing the normality
> #H0: THE DATA IS NORMAL
> #H1: THE DATA IS NOT NORMAL
> [Link](mpg)
Shapiro-Wilk normality test
data: mpg
W = 0.94756, p-value = 0.1229
> #hence the p value greater than 0.05, ACCEPT H0 , mpg data is norm
al
> #we need to use the scale test of two sample t test(two sample data
)
> #H0; HAS impact of vs on mpg
> #H1; NO impact of vs on mpg
> vs
[1] 0 0 1 1 0 1 0 1 1 1 1 0 0 0 0 0 0 1 1 1 1 0 0 0 0 1 0 1 0 0 0 1
> mpg
[1] 21.0 21.0 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 17.8 16.4 17.3
15.2 10.4 10.4 14.7 32.4 30.4 33.9 21.5
[22] 15.5 15.2 13.3 19.2 27.3 26.0 30.4 15.8 19.7 15.0 21.4
> [Link](mpg~vs,alternative="[Link]")
R Console Page 5
Welch Two Sample t-test
data: mpg by vs
t = -4.6671, df = 22.716, p-value = 0.0001098
alternative hypothesis: true difference in means between group 0 and
group 1 is not equal to 0
95 percent confidence interval:
-11.462508 -4.418445
sample estimates:
mean in group 0 mean in group 1
16.61667 24.55714
> #Hence the p-value<0.05 , ACCEPT H1
>
>
>
> #TEST-04
> am
[1] 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 1 1 1 1 1 1 1
> disp
[1] 160.0 160.0 108.0 258.0 360.0 225.0 360.0 146.7 140.8 167.6 167.
6 275.8 275.8 275.8 472.0 460.0 440.0
[18] 78.7 75.7 71.1 120.1 318.0 304.0 350.0 400.0 79.0 120.3 95.
1 351.0 145.0 301.0 121.0
> #testing the normality
> #H0: THE DATA IS NORMAL
> #H1: THE DATA IS NOT NORMAL
> [Link](disp)
Shapiro-Wilk normality test
data: disp
W = 0.92001, p-value = 0.02081
> #Hence p-value lesser than 0,05 , accept H1
> #need to use the odinal test of mann whitney u test
> #H0; AM HAS NO IMPACT ON DISP
> #H1; AM HAS IMPACT ON DISP
> [Link](disp,am,alternative="[Link]")
Wilcoxon rank sum test with continuity correction
data: disp and am
W = 1024, p-value = 2.75e-12
alternative hypothesis: true location shift is not equal to 0
Warning message:
In [Link](disp, am, alternative = "[Link]") :
R Console Page 6
cannot compute exact p-value with ties
> #hence p value lesser than 0.05 , ACCEPT H
>
>
>
>
> #############THANK---YOU################
