Nothing
`KappaRepet` <-
function(Obs, Fit, TSS=FALSE)
{
if(sum(Obs)==0) stop("\n The observed data only contains 0")
tab <- as.data.frame(matrix(0, nrow=101, ncol=2)) ### il faut prÈciser que: "nrow = 101" sinon le data frame se remplit de "NA" au fur et ‡ mesure plutÙt que de "0" (et si il y a des NA cela pose problËme plus loin).
if(length(unique(Fit))==1){
Misc<-table(as.vector(Fit) >= as.numeric(unique(Fit)), Obs) ### Robin modified here...(avant Áa plantait l‡)
if(TSS!=TRUE) a <- KappaStat(Misc)
else a <- TSS.Stat(Misc)
TP <- Misc[4]
TN <- Misc[1]
ca0 <- (TN * 100)/sum(Misc[,1])
ca1 <- (TP * 100)/sum(Misc[,2])
if(is.na(ca0)) ca0<-0
if(is.na(ca1)) ca1<-0
if(TSS!=TRUE) invisible(list(Kappa=0, CutOff=unique(Fit), TP=TP, se=ca1, TN=TN, sp=ca0))
else invisible(list(TSS=0, CutOff=unique(Fit), TP=TP, se=ca1, TN=TN, sp=ca0))
}
else{
Quant <- quantile(Fit)
for(j in 0:100){
Seuil <- Quant[1] + (j*((Quant[5] - Quant[1])/100))
Misc<-table(Fit >= Seuil, Obs)
if(TSS!=TRUE) a <- KappaStat(Misc) else a <- TSS.Stat(Misc)
if(!is.na(a)) if(a > 0) {tab[j+1, 1] <- Seuil; tab[j+1, 2] <- a}
rm(Misc, Seuil)
}
t <- max(tab[,2],na.rm=TRUE)
seuil <- tab[tab[,2]==t,1] ### Note: Ici il se peut qu'on aie plus de 1 seuil...dans ce cas le plus bas est gardÈ.
if(t > 0) {
Misc<-table(Fit >= seuil[1], Obs)
TP <- Misc[4]
TN <- Misc[1]
ca0 <- (TN * 100)/sum(Misc[,1])
ca1 <- (TP * 100)/sum(Misc[,2])
if(TSS!=TRUE) invisible(list(Kappa = t, CutOff = seuil[1], TP = TP, se = ca1, TN = TN, sp = ca0))
else invisible(list(TSS = t, CutOff = seuil[1], TP = TP, se = ca1, TN = TN, sp = ca0))
}
else {
if(TSS!=TRUE) invisible(list(Kappa = 0, CutOff = 0, TP = 0, se = 0, TN = 0, sp = 0))
else invisible(list(TSS = 0, CutOff = 0, TP = 0, se = 0, TN = 0, sp = 0))
}
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.