1_Descriptive_Statistic .ipynb

Kernel: R (R-Project)

tabl <- read.csv('mouseWeight1.txt')
names(tabl)

1. 'дієта'
2. 'вага'

tabl[10:15,]
tabl[60:65,]

	дієта	вага
10	1	30.43728
11	1	20.40348
12	1	21.79013
13	1	21.06970
14	1	20.63963
15	1	22.32389

	дієта	вага
60	1	22.41765
61	2	24.30747
62	2	20.26251
63	2	23.27339
64	2	23.19406
65	2	25.05985

w <- tabl$вага
w1 <- tabl$вага[tabl$дієта==1]
w2 <- tabl$вага[tabl$дієта==2]

w1

1. 23.3422098507307
2. 17.9799685009413
3. 21.0958631701005
4. 23.0860782352691
5. 25.2555575952414
6. 25.5626329113725
7. 20.4903074657829
8. 20.4375390180754
9. 21.3565206547252
10. 30.4372815769861
11. 20.4034793738554
12. 21.7901253581698
13. 21.069695292415
14. 20.6396335259849
15. 22.3238908251481
16. 19.315251143612
17. 21.6018886364804
18. 22.622448783711
19. 24.6530579563722
20. 25.5407804183289
21. 19.8954263144445
22. 21.690702054744
23. 27.5805707351191
24. 24.5366788510043
25. 17.9780968569553
26. 21.4650001239849
27. 23.3823133586978
28. 22.5366471358651
29. 22.0435927078604
30. 20.6943305669254
31. 28.6239804743051
32. 28.2741461981801
33. 26.4463308422104
34. 23.4101460034366
35. 27.288417968264
36. 27.981335524665
37. 22.4123963276562
38. 22.1951833030105
39. 23.5344445322319
40. 24.9956932400151
41. 27.2979150047325
42. 27.6978888081471
43. 22.7110502328542
44. 18.8777842560575
45. 20.3694377217311
46. 25.8371928736841
47. 20.4800586156757
48. 29.6750422102166
49. 19.9592798545229
50. 28.9190950928138
51. 21.3585871729976
52. 20.7542275583147
53. 24.8046960844302
54. 22.6225243061017
55. 22.4804454082181
56. 22.8167671985793
57. 22.9884999254735
58. 28.1601356396926
59. 27.3264035244652
60. 22.4176518802258

mn1 <- mean(w1)
md1 <- median(w1)
rng1 <- diff(range(w1))
sd1 <- sd(w1)
qnt1 <- quantile(w1,c(0.25,0.75))

mn2 <- mean(w2)
md2 <- median(w2)
rng2 <- diff(range(w2))
sd2 <- sd(w2)
qnt2 <- quantile(w2,c(0.25,0.75))

cat( '            контроль | дієта\n')
cat( 'середнє     ', mn1, mn2 ,'\n')
cat( 'медіана     ', md1,md2 ,'\n')
cat( 'діапазон    ', rng1, rng2,'\n')
cat( 'станд відх  ', sd1, sd2 ,'\n')
cat( '1ша квартиль', qnt1[1], qnt2[1],'\n')
cat( '3тя квартиль', qnt1[2], qnt2[2],'\n')

            контроль | дієта
середнє      23.39207 24.61242 
медіана      22.62249 24.15584 
діапазон     12.45918 14.08372 
станд відх   3.077543 2.876705 
1ша квартиль 21.08932 22.91824 
3тя квартиль 25.54624 26.83143 

xl=range( w )

#windows()
par(mfrow=c(2,1))

hist(w1,col=rgb(1,0,0,0.5), breaks=30, xlim=xl, freq = T, main='Control', xlab='Weight, g', ylab='Frequency')
hist(w2,col=rgb(1,0,0,0.5), breaks=30, xlim=xl, freq = T, main='Diet', xlab='Weight, g', ylab='Frequency')

#abline(v=c(mean(w1),median(w1)) ,col=c('red','green') )
#abline(v=c(mean(w2),median(w2)) ,col=c('red','green') )

#windows()

hist(w1,col=rgb(1,0,0,0.5), breaks=30, xlim=xl, freq = T, main='Control vs Diet', xlab='Weight, g', ylab='Frequency')
hist(w2,col=rgb(0,1,0,0.5), breaks=30, xlim=xl, freq = T, add=T)
legend( "topright", c('Control','Diet'), fill=c(rgb(1,0,0,0.5),rgb(0,1,0,0.5)) )

boxplot(list(Control=w1,Diet=w2),   col=c('red','green'),
        main='Control vs Diet' , ylab='Weight, g')

boxplot( вага~дієта,tabl, col=c('red','green') , main='Control vs Diet', ylab='Weight, g', names=c('Control','Diet'))