【悬赏30论坛币】序列模式(arulesSequences) 有计算lift、confidence功能吗？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

看到序列模式里面仅仅计算了support，请问如何计算lift、confidence?
目前程序如下:

1. library(arulesSequences)
   
2. data.seq <- as(data.frame(item=data$late_class),"transactions")
   
3. transactionInfo(data.seq)$sequenceID<-data$user_id
   
4. transactionInfo(data.seq)$eventID<-data$rank
   
5. head(transactionInfo(data.seq),20)
   
6. results<-cspade(data.seq,parameter=list(support=0,maxlen=4,maxsize=NULL),control=list(verbose=TRUE))
   
7. results<-sort(result,by="support")
   
8. inspect(results[1:10])

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：confidence Sequences sequence arules Rules library control support result

扫了一遍这个包的pdf，没发现有计算confidence、lift功能，所以自己写了一个

1. pp <- results@data@p
   
2. ii <- results@data@i
   
3. conf <- vector("numeric",length=length(pp)-1)
   
4. lift <- vector("numeric",length=length(pp)-1)
   
5. support <- results@quality$support
   
6. names(conf) <- names(lift) <- names(support) <- rep('WAIT',times=length(conf))
   
7. for(i in 1:(length(pp)-1)){
   
8.         index <- (pp[i]+1):pp[i+1]
   
9.         names(conf)[i] <- names(lift)[i] <- names(support)[i] <- paste(ii[index],collapse=' ')
   
10.         if(length(index)>1){
    
11.                 temp <- support[which(names(support)==paste(ii[index][1:(length(index)-1)],collapse=' '))]
    
12.                 conf[i] <- support[i]/temp
    
13.                 lift[i] <- conf[i]/support[which(names(support)==ii[index][length(index)])]
    
14.         }
    
15. }
    
16. final <- cbind(support,conf,lift)[which((lift>2&conf>0.3)|lift<0.5&conf>0.3),]

复制代码