if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")

BiocManager::install("mixtools")
BiocManager::install("bcellViper")
BiocManager::install("viper")

library(viper)
library(mixtools)
library(bcellViper)
library(pheatmap)
library(RColorBrewer)
library (dplyr)
library(Seurat)
library(patchwork)
library(ggplot2)
library(sctransform)
library(viridis)

sdat<-seu@assaysseurat_clusters

summary(rowSums(sdat))

summary(rowSums(abs(sdat)))
plot(hist(rowSums(sdat)))
plot(hist(rowSums(abs(sdat))))

comparision with cluster 5;

for (i in 1:11){

print (i-1)

id1<-which(cluster==i-1)
id2<-which(cluster==1-i)

reg<-regul_HT4d_ref_unpr

signature12<-rowTtest(sdat[,id1],sdat[,id2])
signature<-signature12
signature <- (qnorm(signaturestatistic))[, 1]

nullmodel <- ttestNull(sdat[,id1],sdat[,id2], per = 100,
repos = TRUE, verbose = FALSE)

mrs <- msviper(signature, reg, nullmodel, verbose = FALSE)

c<-paste(i, "i", sep="")
out<-paste(c, "pdf", sep=".")

out.rda<-paste(c, "rda", sep=".")

pdf(file=out)
summary(mrs)
plot(mrs, cex = .7)

mrs <- ledge(mrs)

dev.off()
}

########################