R/yf.R

In ternvis: Visualisation, Verification and Calibration of Ternary Probabilistic Forecasts

yf <-
function(p=cbind(1,1,1)/3,M=tsetup()$M32) {
   n <- nrow(p)
   yout <- NA*p[,1] # default for fault finding
   for (i in 1:n){
    if(!is.na(p[i,2]))
                  yout[i] <-  (M %*% p[i,])[2]
   }
   cbind(yout)
}