Guest User

Untitled

a guest
Jan 17th, 2019
87
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.55 KB | None | 0 0
  1. set.seed(1)
  2. df <- data.frame(a=rnorm(10),b=rnorm(10))
  3. summarydist <- function(x) {
  4. y1 <- summary(x)
  5. y2 <- IQR(x)
  6. names(y2) <- "IQR"
  7. require(moments)
  8. y3 <- skewness(x)
  9. names(y3) <- "Skewness"
  10. y4 <- kurtosis(x)
  11. names(y4) <- "kurtosis"
  12. c(y1,y2,y3,y4)
  13.  
  14. }
  15. sapply(df,summarydist)
  16.  
  17. # a b
  18. #Min. -0.8356000 -2.2150000
  19. #1st Qu. -0.5462000 -0.0377500
  20. #Median 0.2566000 0.4919000
  21. #Mean 0.1322000 0.2488000
  22. #3rd Qu. 0.5537000 0.9132000
  23. #Max. 1.5950000 1.5120000
  24. #IQR 1.0998807 0.9509302
  25. #Skewness 0.2961938 -1.1871418
  26. #kurtosis 2.2752871 3.8598299
  27.  
  28. # library(fBasics)
  29. basicStats(df)
  30. a b
  31. nobs 10.000000 10.000000
  32. NAs 0.000000 0.000000
  33. Minimum -0.835629 -2.214700
  34. Maximum 1.595281 1.511781
  35. 1. Quartile -0.546187 -0.037748
  36. 3. Quartile 0.553693 0.913182
  37. Mean 0.132203 0.248845
  38. Median 0.256576 0.491872
  39. Sum 1.322028 2.488450
  40. SE Mean 0.246843 0.338210
  41. LCL Mean -0.426195 -0.516240
  42. UCL Mean 0.690600 1.013930
  43. Variance 0.609314 1.143862
  44. Stdev 0.780586 1.069515
  45. Skewness 0.252895 -1.013599
  46. Kurtosis -1.157017 0.126462
  47.  
  48. DF <- basicStats(df)[c(3:8,15:16),]
  49. rbind(DF, IQR=DF[4,]-DF[3,])
  50. a b
  51. Minimum -0.835629 -2.214700
  52. Maximum 1.595281 1.511781
  53. 1. Quartile -0.546187 -0.037748
  54. 3. Quartile 0.553693 0.913182
  55. Mean 0.132203 0.248845
  56. Median 0.256576 0.491872
  57. Skewness 0.252895 -1.013599
  58. Kurtosis -1.157017 0.126462
  59. IQR 1.099880 0.950930
Add Comment
Please, Sign In to add comment