mat.or.vec <- function(nr,nc) 
        if(nc==1) numeric(nr) else matrix(0,nr,nc)
