SHARE
TWEET

CY Leung ratings, 201403

a guest Mar 14th, 2014 522 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. # Date file was retrieved on March 14, 2014, from
  2. # http://hkupop.hku.hk/chinese/popexpress/ce2012/dataset/tp1403013.sav
  3.  
  4. ############################
  5. # To load in the data file #
  6. ############################
  7.  
  8. # Install the "foreign" package, to read in SPSS files
  9. install.packages("foreign")
  10.  
  11. # Load in the "foreign" package
  12. library("foreign")
  13.  
  14. # Read in the data file
  15. cy201403 <- read.spss("tp1403013.sav", to.data.frame=T)
  16.  
  17. # Check if file has been read in correctly
  18. head(cy201403)
  19.  
  20.  
  21. ######################
  22. # Basic descriptions #
  23. ######################
  24.  
  25. # Number of subjects = 1017
  26. dim(cy201403)[1]
  27.  
  28. # Number of raters = 998
  29. sum(table(cy201403$CE_rating))
  30.  
  31.  
  32.  
  33. ##################
  34. # On raw ratings #
  35. ##################
  36.  
  37. # A quick scan of the distribution of the raw ratings
  38. table(cy201403$CE_rating)
  39.  
  40. # Mean of raw ratings = 47.4008
  41. mean(cy201403$CE_rating, na.rm=T)
  42.  
  43. # Median of raw ratings = 50
  44. median(cy201403$CE_rating, na.rm=T)
  45.  
  46. # Histogram of the raw ratings
  47. hist(cy201403$CE_rating, main="CY Leung's rating (raw), \nMarch 2013",
  48.         xlab="Raw Scores", ylab="Frequency", right=F, axes=F,
  49.         breaks=seq(0, 100, 10), col="lightgray")
  50. axis(1, at=seq(0, 100, 10))
  51. axis(2, at=seq(0, 300, 50))
  52.  
  53. # Indicate mean on plot
  54. abline(v=mean(cy201403$CE_rating, na.rm=T), col="red", lwd=4)
  55.  
  56. # Indicate median on plot
  57. abline(v=median(cy201403$CE_rating, na.rm=T), col="blue", lwd=4)
  58.  
  59. # Add labels
  60. text(mean(cy201403$CE_rating, na.rm=T), 200, "Mean", col="red", lwd=4, pos=2)
  61. text(median(cy201403$CE_rating, na.rm=T), 200, "Median", col="blue", lwd=4, pos=4)
  62.  
  63.  
  64. ######################
  65. # On weighed ratings #
  66. ######################
  67.  
  68. # Calculate weighed ratings
  69. cy201403$weighedRating <- cy201403$CE_rating * cy201403$weight
  70.  
  71. # Mean of raw ratings = 47.48411
  72. mean(cy201403$weighedRating, na.rm=T)
  73.  
  74. # Median of raw ratings = 41.9557
  75. median(cy201403$weighedRating, na.rm=T)
  76.  
  77. # Histrogram of the weighed ratings
  78. hist(cy201403$weighedRating, main="CY Leung's rating (weighed), \nMarch 2013",
  79.         xlab="Weighed Scores", ylab="Frequency", right=F, axes=F,
  80.         breaks=seq(0, 500, 10), col="lightgray")
  81. axis(1, at=seq(0, 500, 50))
  82. axis(2, at=seq(0, 140, 20))
  83.  
  84. # Indicate mean on plot
  85. abline(v=mean(cy201403$weighedRating, na.rm=T), col="red", lwd=4)
  86.  
  87. # Indicate median on plot
  88. abline(v=median(cy201403$weighedRating, na.rm=T), col="blue", lwd=4)
  89.  
  90. # Add labels
  91. text(mean(cy201403$weighedRating, na.rm=T), 100, "Mean", col="red", lwd=4, pos=4)
  92. text(median(cy201403$weighedRating, na.rm=T), 100, "Median", col="blue", lwd=4, pos=2)
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top