Correlation

#Library set library(“readxl”) library(FSA) library(DescTools) require(reshape2) library(ggplot2) library(dplyr) library(outliers) library(ggpubr) library(emmeans) library(car)

#Open data data <- read_excel(“D:/Research/Result/Drug treatment project/Data/NO DRUG/2019-Project 1 DATA -Vish-Thong.xlsx”,sheet = “Whole”) names(data) data\(Group=as.factor(data\)Group) data\(Genotype=as.factor(data\)Genotype) data\(GroupxGenotype=as.factor(data\)GroupxGenotype) data\(EXP=as.factor(data\)EXP) data\(REG=as.factor(data\)REG)

#data handling data=select(data,D2S,D2L,NR3,D1R) data <- as.data.frame(data) class(data) data1=melt(data,id= “D1R”,measure.vars=c(“NR3”, “D2S”, “D2L”))

head(data1)

sp <- ggscatter(data1, x = “D1R”, y = “value”, color = “variable”, palette = “jco”, add = “reg.line”, conf.int = TRUE) sp + stat_cor(aes(color = variable))