R/Normal_Distr_Test_graph.R

In yikeshu0611/NormalDistrTest: To realize normal distrbution tests

#' To plot graphs of transformed data
#'
#' @param x your data
#'
#' @return several plots
#' @export
#'
#' @examples
#' x<-1:100
#' Normal_Distr_Test_graph(x)
Normal_Distr_Test_graph<-function(x){
  x=na.omit(x)
  addd=0
  if (length(x) >2){
    xt1=x
    addd=addd+1
    xt2=suppressWarnings(log(x))
    if ((!(NaN %in% xt2)) & (!(-Inf %in% xt2))& (!(Inf %in% xt2))){addd=addd+1}
    xt3=suppressWarnings(x^2)
    if ((!(NaN %in% xt3)) & (!(-Inf %in% xt3))& (!(Inf %in% xt3))){addd=addd+1}
    xt4=suppressWarnings(x^3)
    if ((!(NaN %in% xt4)) & (!(-Inf %in% xt4))& (!(Inf %in% xt4))){addd=addd+1}
    xt5=suppressWarnings(sqrt(x))
    if ((!(NaN %in% xt5)) & (!(-Inf %in% xt5))& (!(Inf %in% xt5))){addd=addd+1}
    xt6=suppressWarnings(1/x)
    if ((!(NaN %in% xt6)) & (!(-Inf %in% xt6))& (!(Inf %in% xt6))){addd=addd+1}
    xt7=suppressWarnings(1/(x^2))
    if ((!(NaN %in% xt7)) & (!(-Inf %in% xt7))& (!(Inf %in% xt7))){addd=addd+1}
    xt8=suppressWarnings(1/(x^3))
    if ((!(NaN %in% xt8)) & (!(-Inf %in% xt8))& (!(Inf %in% xt8))){addd=addd+1}
    xt9=suppressWarnings(1/sqrt(x))
    if ((!(NaN %in% xt9)) & (!(-Inf %in% xt9))& (!(Inf %in% xt9))){addd=addd+1}
    xt10=suppressWarnings(asin(sqrt(x)))
    if ((!(NaN %in% xt10)) & (!(-Inf %in% xt10))& (!(Inf %in% xt10))){addd=addd+1}
    if (addd==1){
      hist(xt1, freq = F,main ="identity")
      points(density(xt1),cex=0.5)
      lines(density(xt1))
    }else if((addd>=2) & (addd<=3)){
      par(mfrow=c(1,3))
      hist(xt1, freq = F,main ="identity")
      points(density(xt1),cex=0.5)
      lines(density(xt1))
      if ((!(NaN %in% xt2)) &(!(-Inf %in% xt2))& (!(Inf %in% xt2))){
        hist(xt2, freq = F,main="log(x)")
        points(density(xt2),cex=0.5)
        lines(density(xt2))
      }
      if ((!(NaN %in% xt3)) & (!(-Inf %in% xt3))& (!(Inf %in% xt3))){
        hist(xt3, freq = F,main="x^2")
        points(density(xt3),cex=0.5)
        lines(density(xt3))
      }
      if ((!(NaN %in% xt4)) &  (!(-Inf %in% xt4))& (!(Inf %in% xt4))){
        hist(xt4, freq = F,main="x^3")
        points(density(xt4),cex=0.5)
        lines(density(xt4))
      }
      if ((!(NaN %in% xt5)) & (!(-Inf %in% xt5))& (!(Inf %in% xt5))){
        hist(xt5, freq = F,main="sqrt(x)")
        points(density(xt5),cex=0.5)
        lines(density(xt5))
      }
      if ((!(NaN %in% xt6)) & (!(-Inf %in% xt6))& (!(Inf %in% xt6))){
        hist(xt6, freq = F,main="1/x")
        points(density(xt6),cex=0.5)
        lines(density(xt6))
      }
      if ((!(NaN %in% xt7)) & (!(-Inf %in% xt7))& (!(Inf %in% xt7))){
        hist(xt7, freq = F,main="1/(x^2)")
        points(density(xt7),cex=0.5)
        lines(density(xt7))
      }
      if ((!(NaN %in% xt8)) & (!(-Inf %in% xt8))& (!(Inf %in% xt8))){
        hist(xt8, freq = F,main="1/(x^3)")
        points(density(xt8),cex=0.5)
        lines(density(xt8))
      }
      if ((!(NaN %in% xt9)) & (!(-Inf %in% xt9))& (!(Inf %in% xt9))){
        hist(xt9, freq = F,main="1/sqrt(x)")
        points(density(xt9),cex=0.5)
        lines(density(xt9))
      }
      if ((!(NaN %in% xt10)) & (!(-Inf %in% xt10))& (!(Inf %in% xt10))){
        hist(xt10, freq = F,main="asin(sqrt(x))")
        points(density(xt10),cex=0.5)
        lines(density(xt10))
      }
    }else if((addd>=4) & (addd<=6)){
      par(mfrow=c(2,3))
      hist(xt1, freq = F,main ="identity")
      points(density(xt1),cex=0.5)
      lines(density(xt1))
      if ((!(NaN %in% xt2)) &(!(-Inf %in% xt2))& (!(Inf %in% xt2))){
        hist(xt2, freq = F,main="log(x)")
        points(density(xt2),cex=0.5)
        lines(density(xt2))
      }
      if ((!(NaN %in% xt3)) & (!(-Inf %in% xt3))& (!(Inf %in% xt3))){
        hist(xt3, freq = F,main="x^2")
        points(density(xt3),cex=0.5)
        lines(density(xt3))
      }
      if ((!(NaN %in% xt4)) &  (!(-Inf %in% xt4))& (!(Inf %in% xt4))){
        hist(xt4, freq = F,main="x^3")
        points(density(xt4),cex=0.5)
        lines(density(xt4))
      }
      if ((!(NaN %in% xt5)) & (!(-Inf %in% xt5))& (!(Inf %in% xt5))){
        hist(xt5, freq = F,main="sqrt(x)")
        points(density(xt5),cex=0.5)
        lines(density(xt5))
      }
      if ((!(NaN %in% xt6)) & (!(-Inf %in% xt6))& (!(Inf %in% xt6))){
        hist(xt6, freq = F,main="1/x")
        points(density(xt6),cex=0.5)
        lines(density(xt6))
      }
      if ((!(NaN %in% xt7)) & (!(-Inf %in% xt7))& (!(Inf %in% xt7))){
        hist(xt7, freq = F,main="1/(x^2)")
        points(density(xt7),cex=0.5)
        lines(density(xt7))
      }
      if ((!(NaN %in% xt8)) & (!(-Inf %in% xt8))& (!(Inf %in% xt8))){
        hist(xt8, freq = F,main="1/(x^3)")
        points(density(xt8),cex=0.5)
        lines(density(xt8))
      }
      if ((!(NaN %in% xt9)) & (!(-Inf %in% xt9))& (!(Inf %in% xt9))){
        hist(xt9, freq = F,main="1/sqrt(x)")
        points(density(xt9),cex=0.5)
        lines(density(xt9))
      }
      if ((!(NaN %in% xt10)) & (!(-Inf %in% xt10))& (!(Inf %in% xt10))){
        hist(xt10, freq = F,main="asin(sqrt(x))")
        points(density(xt10),cex=0.5)
        lines(density(xt10))
      }
    }else if((addd>=7) & (addd<=9)){
      par(mfrow=c(3,3))
      hist(xt1, freq = F,main ="identity")
      points(density(xt1),cex=0.5)
      lines(density(xt1))
      if ((!(NaN %in% xt2)) &(!(-Inf %in% xt2))& (!(Inf %in% xt2))){
        hist(xt2, freq = F,main="log(x)")
        points(density(xt2),cex=0.5)
        lines(density(xt2))
      }
      if ((!(NaN %in% xt3)) & (!(-Inf %in% xt3))& (!(Inf %in% xt3))){
        hist(xt3, freq = F,main="x^2")
        points(density(xt3),cex=0.5)
        lines(density(xt3))
      }
      if ((!(NaN %in% xt4)) &  (!(-Inf %in% xt4))& (!(Inf %in% xt4))){
        hist(xt4, freq = F,main="x^3")
        points(density(xt4),cex=0.5)
        lines(density(xt4))
      }
      if ((!(NaN %in% xt5)) & (!(-Inf %in% xt5))& (!(Inf %in% xt5))){
        hist(xt5, freq = F,main="sqrt(x)")
        points(density(xt5),cex=0.5)
        lines(density(xt5))
      }
      if ((!(NaN %in% xt6)) & (!(-Inf %in% xt6))& (!(Inf %in% xt6))){
        hist(xt6, freq = F,main="1/x")
        points(density(xt6),cex=0.5)
        lines(density(xt6))
      }
      if ((!(NaN %in% xt7)) & (!(-Inf %in% xt7))& (!(Inf %in% xt7))){
        hist(xt7, freq = F,main="1/(x^2)")
        points(density(xt7),cex=0.5)
        lines(density(xt7))
      }
      if ((!(NaN %in% xt8)) & (!(-Inf %in% xt8))& (!(Inf %in% xt8))){
        hist(xt8, freq = F,main="1/(x^3)")
        points(density(xt8),cex=0.5)
        lines(density(xt8))
      }
      if ((!(NaN %in% xt9)) & (!(-Inf %in% xt9))& (!(Inf %in% xt9))){
        hist(xt9, freq = F,main="1/sqrt(x)")
        points(density(xt9),cex=0.5)
        lines(density(xt9))
      }
      if ((!(NaN %in% xt10)) & (!(-Inf %in% xt10))& (!(Inf %in% xt10))){
        hist(xt10, freq = F,main="asin(sqrt(x))")
        points(density(xt10),cex=0.5)
        lines(density(xt10))
      }
    }else if(addd==10){
      par(mfrow=c(4,3))
      hist(xt1, freq = F,main ="identity")
      points(density(xt1),cex=0.5)
      lines(density(xt1))
      if ((!(NaN %in% xt2)) &(!(-Inf %in% xt2))& (!(Inf %in% xt2))){
        hist(xt2, freq = F,main="log(x)")
        points(density(xt2),cex=0.5)
        lines(density(xt2))
      }
      if ((!(NaN %in% xt3)) & (!(-Inf %in% xt3))& (!(Inf %in% xt3))){
        hist(xt3, freq = F,main="x^2")
        points(density(xt3),cex=0.5)
        lines(density(xt3))
      }
      if ((!(NaN %in% xt4)) &  (!(-Inf %in% xt4))& (!(Inf %in% xt4))){
        hist(xt4, freq = F,main="x^3")
        points(density(xt4),cex=0.5)
        lines(density(xt4))
      }
      if ((!(NaN %in% xt5)) & (!(-Inf %in% xt5))& (!(Inf %in% xt5))){
        hist(xt5, freq = F,main="sqrt(x)")
        points(density(xt5),cex=0.5)
        lines(density(xt5))
      }
      if ((!(NaN %in% xt6)) & (!(-Inf %in% xt6))& (!(Inf %in% xt6))){
        hist(xt6, freq = F,main="1/x")
        points(density(xt6),cex=0.5)
        lines(density(xt6))
      }
      if ((!(NaN %in% xt7)) & (!(-Inf %in% xt7))& (!(Inf %in% xt7))){
        hist(xt7, freq = F,main="1/(x^2)")
        points(density(xt7),cex=0.5)
        lines(density(xt7))
      }
      if ((!(NaN %in% xt8)) & (!(-Inf %in% xt8))& (!(Inf %in% xt8))){
        hist(xt8, freq = F,main="1/(x^3)")
        points(density(xt8),cex=0.5)
        lines(density(xt8))
      }
      if ((!(NaN %in% xt9)) & (!(-Inf %in% xt9))& (!(Inf %in% xt9))){
        hist(xt9, freq = F,main="1/sqrt(x)")
        points(density(xt9),cex=0.5)
        lines(density(xt9))
      }
      if ((!(NaN %in% xt10)) & (!(-Inf %in% xt10))& (!(Inf %in% xt10))){
        hist(xt10, freq = F,main="asin(sqrt(x))")
        points(density(xt10),cex=0.5)
        lines(density(xt10))
      }
    }
  }else{
  stop("Sample size is less than 2")
  }
  par(mfrow=c(1,1))
}