Genetic Interaction

Box-plot for each attribute cut in 5 bins

gi <- read.csv("D:/0-DOUTORADO/GeneticInteraction/calculos/wholeModel/R/WholeProcessado2.csv")

labels <- names(gi[4:length(gi) - 1])

labels

##  [1] "max_degInt" "max_betInt" "max_degppi" "max_degreg" "max_degmet"
##  [6] "max_regin"  "max_regout" "max_metin"  "max_metout" "max_betppi"
## [11] "max_betreg" "max_betmet" "min_degInt" "min_betInt" "min_degppi"
## [16] "min_degreg" "min_degmet" "min_regin"  "min_regout" "min_metin" 
## [21] "min_metout" "min_betppi" "min_betreg" "min_betmet" "CNInt"     
## [26] "CNppi"      "CNreg"      "CNmet"      "FSWInt"     "FSWppi"    
## [31] "FSWreg"     "FSWmet"     "JCI"        "JCppi"      "JCreg"     
## [36] "JCmet"      "shortest"

for (name in labels) {
    lm <- lm(gi$score ~ gi[[name]])

    cut <- cut(gi[[name]], breaks = 5)
    boxplot(gi$score ~ cut, xlab = name, ylab = "S-score", varwidth = TRUE)
    abline(lm, lwd = 2, col = "red")
}

plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 plot of chunk unnamed-chunk-1 

## Error: 'breaks' are not unique

plot of chunk unnamed-chunk-1

breaks <- c(-Inf,unique(quantile(gi[[name]], na.rm=T)),Inf)

Box-plot for each attribute cut in 20 bins

gi <- read.csv("D:/0-DOUTORADO/GeneticInteraction/calculos/wholeModel/R/WholeProcessado2.csv")

labels <- names(gi[4:length(gi) - 1])

labels

##  [1] "max_degInt" "max_betInt" "max_degppi" "max_degreg" "max_degmet"
##  [6] "max_regin"  "max_regout" "max_metin"  "max_metout" "max_betppi"
## [11] "max_betreg" "max_betmet" "min_degInt" "min_betInt" "min_degppi"
## [16] "min_degreg" "min_degmet" "min_regin"  "min_regout" "min_metin" 
## [21] "min_metout" "min_betppi" "min_betreg" "min_betmet" "CNInt"     
## [26] "CNppi"      "CNreg"      "CNmet"      "FSWInt"     "FSWppi"    
## [31] "FSWreg"     "FSWmet"     "JCI"        "JCppi"      "JCreg"     
## [36] "JCmet"      "shortest"

for (name in labels) {
    lm <- lm(gi$score ~ gi[[name]])
    breaks <- c(-Inf, unique(quantile(gi[[name]], na.rm = T)), Inf)
    cut <- cut(gi[[name]], breaks = 20)
    boxplot(gi$score ~ cut, xlab = name, ylab = "S-score")
    abline(lm, lwd = 2, col = "red")
}

plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 plot of chunk unnamed-chunk-2 

## Error: 'breaks' are not unique

plot of chunk unnamed-chunk-2