R/plot_Indv.R

In sarahodell/magicsim: What the Package Does (One Line, Title Case)

#' plot_Indv ---------------------------------------------------------------
#' Plot all chromosomes of a single Indv object
#' @param breaktable A break table from function make_breaktable
#' @param het Whether or not the ChromPair is heterozygous. Default is TRUE
#'
#' @return a dataframe with the parental donor and location of breakpoints
#' @export
#'

plot_Indv <- function(breaktable,het=T){
  if(het==F){
    breaktable=breaktable[breaktable$hom=="H1",]
  }
  chr_f=sprintf("%.0f_%s",breaktable$chr,breaktable$hom)
  breaktable$chr_f=factor(chr_f,levels=unique(chr_f))
  p <- ggplot(breaktable,aes(x=start/1e6,y=1,color=donor)) +
    geom_segment(aes(xend=end/1e6,color=donor,yend=1),lineend="butt",size=10) +
    facet_grid(chr_f ~ .) +
    ggtitle("Chromosome Breakpoints") +
    xlab("Position (Mb)")  +
    theme(axis.text.y=element_blank(),axis.title.y=element_blank(),axis.ticks.y=element_blank())
  print(p)
}