Same data for the same code, why do I get different numbers of DEGs. Hello everyone, I encountered a problem when using the samr package to find differential expressed genes. The same data for the same code, why do I get a different number of DEGs? For example, I got 108 DEGs in one analysis, but I got 368 DEGs when I repeated the following code. My code is as follows:

x=read.csv("expr.csv",sep = ',',row.names = 1,header = T)
#View(head(x))
Pheno<-read.csv('target.csv',sep = ',',header = T,row.names = 1)
R<-rownames(pheno)
x=x[,R]
y=pheno$group
Data=list(x=as.matrix(x),y=y,geneid=rownames(x),genenames=as.character(rownames(x)),logged2=TRUE)
Samr.obj<-samr(data,resp.type="Two class unpaired", nperms=1000)
Delta.table <- samr.compute.delta.table(samr.obj, min.foldchange=1.2,nvals=100)
Siggenes.table <- samr.compute.siggenes.table(samr.obj, del=0, data, delta.table,all.genes=TRUE)
a <- siggenes.table$genes.up; # all up regulated genes
b <- siggenes.table$genes.lo; ​​# all down regulated genes
c <- rbind(a,b)
#View(head(c))
Lo <- c[as.numeric(c[,8])<20,]#FDR<20%
For (i in 1:nrow(lo))
{
 Tp <- as.numeric(as.vector(as.matrix(lo[i,1])))-1;
 Lo[i,3] <- as.character(as.vector(as.matrix(x[tp,1])));
}
Write.csv(lo,"DEGs.csv")