latex.table <- function(coef.mat,se.mat,digits=3,table.command=TRUE) {   
  nc <- ncol(coef.mat)
  coef.mat <- round(coef.mat,3)
  se.mat <- round(se.mat,3)
  text.now <- NULL
  if (table.command) {
    text.now <- c(text.now,"\\begin{table}\n")
    text.now <- c(text.now,"\\centering\n")
  }
  text.now <- c(text.now,"\\begin{tabular}[R]{",rep("c",nc+1),"}\n")
  text.now <- c(text.now,"\\hline\n")
  for (j in 1:ncol(coef.mat)) {
    text.now <- c(text.now," & ", colnames(coef.mat)[j])
  }
  text.now <- c(text.now,"\\\\\n\\hline\n")
  for (i in 1:nrow(coef.mat)) {
    ##print coef estimates
    text.now <- c(text.now,rownames(coef.mat)[i])
    for (j in 1:ncol(coef.mat)) {
      if (is.na(coef.mat[i,j])) {
        text.now <- c(text.now," & ")
      } else {
        text.now <- c(text.now," & ", coef.mat[i, j])
      }
    }
    text.now <- c(text.now,"\\\\\n")
    ##  print SEs
    for (j in 1:ncol(coef.mat)) {
      if (is.na(se.mat[i,j])) {
        text.now <- c(text.now," & ")
      } else {
        text.now <- c(text.now," & ", sprintf("(%s)", se.mat[i,j]))
      }
    }
    text.now <- c(text.now,"\\\\[1mm]\n")
  }
  text.now <- c(text.now,"\\\\\n")
  text.now <- c(text.now,"\\hline")
  text.now <- c(text.now,"\n")
  text.now <- c(text.now,"\\end{tabular}\n")
  if (table.command) text.now <- c(text.now,"\\end{table}\n")  
  paste(text.now,collapse="")
}