R/misc.R

In gdepalma/SuscTesting:

theme_dbets <- function () {
  theme_fivethirtyeight() +
    theme(axis.line = element_line(colour = "black"),
          panel.grid.major = element_line(color='gray90'),
          panel.grid.minor = element_blank(),
          axis.text=element_text(size=11),
          axis.title=element_text(size=11),
          plot.title=element_text(size=15,hjust = 0))
}

parse_data=function(a1){

  if(any(names(a1)!=c("MIC","DIA"))==TRUE) stop("File not formatted correctly.")

  xcens=rep(0,nrow(a1))
  ycens=rep(0,nrow(a1))
  xobs=rep(NA,nrow(a1))
  yobs=rep(NA,nrow(a1))

  for(i in 1:nrow(a1)){
    if(grepl('>',a1[i,1])==TRUE){
      xobs[i]=substr(a1[i,1],2,nchar(a1[i,1]))
      xcens[i]=1
    }else if(grepl('<',a1[i,1])==TRUE){
      xobs[i]=substr(a1[i,1],2,nchar(a1[i,1]))
      xcens[i]=-1
    }else
      xobs[i]=a1[i,1]

    if(grepl('>',a1[i,2])==TRUE){
      yobs[i]=substr(a1[i,2],2,nchar(a1[i,2]))
      ycens[i]=1
    }else if(grepl('<',a1[i,2])==TRUE){
      yobs[i]=substr(a1[i,2],2,nchar(a1[i,2]))
      ycens[i]=-1
    }else
      yobs[i]=a1[i,2]
  }

  xobs=as.numeric(xobs)
  yobs=as.numeric(yobs)

  return(list(MIC=xobs,DIA=yobs,xcens=xcens,ycens=ycens))

}