我需要计算2个项目的支持。对于给定的数据,我有一个条件概率矩阵,将一个项目放在购物篮中的第一个项目的概率如下:
c("Hammer", "Nails", "Screws", "Screwdriver", "Wrench", "Level", "Drill", "Brush")
p<-c(0.02, 0.2, 0.2, 0.2, 0.05, 0.01, 0.02, 0.3)
所以对于1-item支持我很简单
t(CondProbs)%*%p
您能给我以下几点的提示吗?
通过说2个项目的支持,我的意思是计算每对产品的出现机会P(X,Y)。我有1个项目支持:
sup1=c(0.15,0.49,0.64,0.57,0.19,0.09,0.12,0.58)
借助此功能,我需要2个项目的支持。我找到了一个示例,但是我可以确定,我应该在不使用 for 函数的情况下进行计算。
>for (j in 1:n){
+ if(sup1[j]>=supthresh){
+ for (k in j:n){
+ if (sup1[k]>=supthresh){
+ sup2[j,k]<-data[,j]%*%data[,k]
+ sup2[k,j]<-sup2[j,k] } } } }
> sup2<-sup2/nobs
> sup2
p.s。 supthresh 是规则的阈值,这里我们有supthresh<-0.2