Shared1_Descriptive_Statistic .ipynbOpen in CoCalc
Author: Artem Mischenko
Views : 14
In [1]:
tabl <- read.csv('mouseWeight1.txt') names(tabl)
  1. 'дієта'
  2. 'вага'
In [2]:
tabl[10:15,] tabl[60:65,]
дієтавага
101 30.43728
111 20.40348
121 21.79013
131 21.06970
141 20.63963
151 22.32389
дієтавага
601 22.41765
612 24.30747
622 20.26251
632 23.27339
642 23.19406
652 25.05985
In [3]:
w <- tabl$вага w1 <- tabl$вага[tabl$дієта==1] w2 <- tabl$вага[tabl$дієта==2]
In [4]:
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
In [5]:
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))
In [6]:
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
In [7]:
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') )
In [8]:
#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)) )
In [9]:
boxplot(list(Control=w1,Diet=w2), col=c('red','green'), main='Control vs Diet' , ylab='Weight, g')
In [10]:
boxplot( вага~дієта,tabl, col=c('red','green') , main='Control vs Diet', ylab='Weight, g', names=c('Control','Diet'))
In [ ]: