开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 人工智能/神经网络 > 查看源码
c_rev.gp
资源名称：apriori(c).zip　[点击查看]
上传用户：lengbin
上传日期：2010-03-31
资源大小：121k
文件大小：2k
源码类别：

人工智能/神经网络

开发平台：

C/C++

c_rev.gp：源码内容
  1. #-----------------------------------------------------------------------
  2. # File    : c_rev.gp
  3. # Contents: visualization of the dependence between posterior confidence
  4. #           and reversed confidence, i.e. the confidence of the reversed
  5. #           rule antecedent <- consequent
  6. # Author  : Christian Borgelt
  7. # History : 23.09.1998 file created
  8. #-----------------------------------------------------------------------
  9. # set terminal postscript
  11. set view 50,330,1
  12. set xrange [0.001:0.999]
  13. set yrange [0.001:0.999]
  14. set zrange [0:1]
  15. set isosamples 20
  16. set hidden3d
  18. set xlabel "c_post"
  19. set ylabel "c_rev"
  21. min(a,b) = a < b ? a : b
  23. diff(x,y,p) = (1 -p -(y*p)/x +y*p >= 0) ? abs(x -p) : 0
  25. quot(x,y,p) = (1 -p -(y*p)/x +y*p >= 0) ? 1 -min(x/p,p/x) : 0
  27. i(x,y,p)    = (y*p*log((y*p)/(p*(y*p)/x)) 
  28.             + (p-y*p)*log((p-y*p)/(p*(1-(y*p)/x))) 
  29.             + ((y*p)/x-y*p)*log(((y*p)/x-y*p)/((1-p)*(y*p)/x)) 
  30.             + (1-p-(y*p)/x+y*p)*log((1-p-(y*p)/x+y*p) 
  31.                                   /((1-p)*(1-(y*p)/x)))) /log(2)
  33. info(x,y,p) = (1 -p -(y*p)/x +y*p >= 0) ? i(x,y,p) : 0
  35. chi2(x,y,p) = (1 -p -(y*p)/x +y*p >= 0) 
  36.             ? (p*((y*p)/x)-y*p)**2 /(p*(1-p)*((y*p)/x)*(1-(y*p)/x)) : 0
  38. set zlabel "d_diff"
  39. set title "prior confidence 0.2"
  40. splot diff(x,y,0.2)
  41. pause -1 "Hit return to continue"
  43. set title "prior confidence 0.3"
  44. splot diff(x,y,0.3)
  45. pause -1 "Hit return to continue"
  47. set title "prior confidence 0.4"
  48. splot diff(x,y,0.4)
  49. pause -1 "Hit return to continue"
  51. set zlabel "d_quot"
  52. set title "prior confidence 0.2"
  53. splot quot(x,y,0.2)
  54. pause -1 "Hit return to continue"
  56. set title "prior confidence 0.3"
  57. splot quot(x,y,0.3)
  58. pause -1 "Hit return to continue"
  60. set title "prior confidence 0.4"
  61. splot quot(x,y,0.4)
  62. pause -1 "Hit return to continue"
  64. set zlabel "d_info"
  65. set title "prior confidence 0.2"
  66. splot info(x,y,0.2)
  67. pause -1 "Hit return to continue"
  69. set title "prior confidence 0.3"
  70. splot info(x,y,0.3)
  71. pause -1 "Hit return to continue"
  73. set title "prior confidence 0.4"
  74. splot info(x,y,0.4)
  75. pause -1 "Hit return to continue"
  77. set zlabel "d_chi^2"
  78. set title "prior confidence 0.2"
  79. splot chi2(x,y,0.2)
  80. pause -1 "Hit return to continue"
  82. set title "prior confidence 0.3"
  83. splot chi2(x,y,0.3)
  84. pause -1 "Hit return to continue"
  86. set title "prior confidence 0.4"
  87. splot chi2(x,y,0.4)
  88. pause -1 "Hit return to continue"