|
|
@ -235,13 +235,37 @@ else if (typediff=='gauss') {
|
|
|
|
|
|
|
|
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
for
|
|
|
|
nbrsm <- c()
|
|
|
|
|
|
|
|
for (d in 1:2) {
|
|
|
|
|
|
|
|
if (!is.null(get(paste0('RSM',d)))) {
|
|
|
|
|
|
|
|
nbrsm[d] <- length(get(paste0('RSM',d)))
|
|
|
|
|
|
|
|
for (i in seq(2,nbrsm[d]+1)) {
|
|
|
|
|
|
|
|
rsm[[d]][i] <- get(paste0("RSM",d))[i-1] # RSM doit etre une liste de vecteurs
|
|
|
|
|
|
|
|
if (!is.null(THR) & rsm[[d]][i]<0) {
|
|
|
|
|
|
|
|
stop('Error 198: With the THR option, the numbers defined in RSM1 and RSM2 must nonnegative')
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
if (!is.null(DIF) & is.null(PCM)) {
|
|
|
|
|
|
|
|
PCM <- as.matrix(DIF)
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
matdiff <- as.matrix(rep(0,nbitems))
|
|
|
|
|
|
|
|
matdisc <- as.matrix(rep(0,nbitems))
|
|
|
|
|
|
|
|
matpmin <- as.matrix(rep(0,nbitems))
|
|
|
|
|
|
|
|
matpmax <- as.matrix(rep(0,nbitems))
|
|
|
|
|
|
|
|
matacc <- as.matrix(rep(0,nbitems))
|
|
|
|
|
|
|
|
matmu <- as.matrix(rep(0,length(DIM)))
|
|
|
|
|
|
|
|
matmu <- as.matrix(rep(0,(length(DIM)+1)*length(DIM)/2 ))
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
if (nbdisc !=0) {
|
|
|
|
|
|
|
|
matdisc <- as.matrix(DISC)
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
else {
|
|
|
|
|
|
|
|
matdisc <- rep(1,nbitems)
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|