arem.gp
上传用户:lengbin
上传日期:2010-03-31
资源大小:121k
文件大小:2k
开发平台:

C/C++

  1. #-----------------------------------------------------------------------
  2. # File    : arem.gp
  3. # Contents: visualization of the selection behaviour
  4. #           of the additional rule evaluation measures
  5. # Author  : Christian Borgelt
  6. # History : 22.09.1998 file created
  7. #-----------------------------------------------------------------------
  8. # set terminal postscript
  9. set view 50,325,1
  10. set hidden3d
  11. set isosamples 20,20
  12. set xrange [0.001:0.999]
  13. set yrange [0.001:0.999]
  14. set zrange [0:1]
  15. min(a,b) = a < b ? a : b
  16. diff(x,y) = abs(x -y)
  17. quot(x,y) = 1 -min(x/y, y/x)
  18. info(x,y,s) = (x >= y*s) && (1-x >= s *(1-y)) 
  19.             ? (s*y*log(y/x) +(x-s*y)*log((x-s*y)/(x*(1-s))) 
  20.             +  s*(1-y)*log((1-y)/(1-x)) 
  21.             + (1-x-s*(1-y))*log((1-x-s*(1-y))/((1-x)*(1-s)))) /log(2) : 0
  22. chi2(x,y,s) = (x >= y*s) && (1-x >= s *(1-y)) 
  23.             ? (x*s -y*s)**2 /(x*(1-x)*s*(1-s)) : 0
  24. set title "d_diff (independent of antecedent support)"
  25. set xlabel "c_prior"
  26. set ylabel "c_post"
  27. set zlabel "d_diff"
  28. splot diff(x,y)
  29. pause -1 "Hit return to continue"
  30. set title "d_quot (independent of antecedent support)"
  31. set zlabel "d_quot"
  32. splot quot(x,y)
  33. pause -1 "Hit return to continue"
  34. set zlabel "d_info"
  35. set title "d_info with antecedent support 0.2"
  36. splot info(x,y,0.2)
  37. pause -1 "Hit return to continue"
  38. set title "d_info with antecedent support 0.3"
  39. splot info(x,y,0.3)
  40. pause -1 "Hit return to continue"
  41. set title "d_info with antecedent support 0.4"
  42. splot info(x,y,0.4)
  43. pause -1 "Hit return to continue"
  44. set zlabel "d_chi^2"
  45. set title "d_chi^2 with antecedent support 0.2"
  46. splot chi2(x,y,0.2)
  47. pause -1 "Hit return to continue"
  48. set title "d_chi^2 with antecedent support 0.3"
  49. splot chi2(x,y,0.3)
  50. pause -1 "Hit return to continue"
  51. set title "d_chi^2 with antecedent support 0.4"
  52. splot chi2(x,y,0.4)
  53. pause -1 "Hit return to continue"