wmat <- function(file,data){ #inverse of rmat...write a .mat
                             #file for Ox to read

          write(dim(data),file=file)
	  write(t(data),file=file,ncol=ncol(data),append=T)
          }