R/HRQoLplot.R

In idaejin/PROreg: Patient Reported Outcomes Regression Analysis

HRQoLplot <- function(data,legend=FALSE,title="Short Form-36 Health Survey",dimlabel.cex=NULL,legend.cex=1,linewidth=3,title.cex=1,lty=1){

  number <- 0
  number.cex <- rep(1,8)

  maxmin <- data.frame(
    PF=c(20,0),
    RP=c(4,0),
    BP=c(9,0),
    GH=c(20,0),
    VT=c(20,0),
    SF=c(8,0),
    RE=c(3,0),
    MH=c(13,0))

  colnames(maxmin) <- c("PF (20)","RP (4)","BP (9)","GH (20)","VT (20)","SF (8)","RE (3)","MH (13)")

  colnames(data) <- c("PF (20)","RP (4)","BP (9)","GH (20)","VT (20)","SF (8)","RE (3)","MH (13)")


  dat <- rbind(maxmin,data)

  radarchart(dat, axistype=number, seg=4, pty=32, plty=lty,plwd=linewidth, na.itp=FALSE,cglcol="black",
             title=title,pcol=brewer.pal(8,"Set1"),vlcex=dimlabel.cex,cex.main=title.cex,axislabcol="black",
             calcex=number.cex)

  # Legend
  if (legend==TRUE){
    legend("topright", legend=c(rownames(dat[-c(1,2),])),text.col=brewer.pal(8,"Set1"),bty="n",cex=legend.cex,
           lty=lty,lwd=legend.cex,col=brewer.pal(8,"Set1"))
  }

}