#CD45+
View(treated)
ggplot(treated, aes(x=Mice, y=CD45 )) +
  geom_point(size=2, shape=23)


#CD3+
ggplot(data=treated, aes(x=Mice , y=CD3 ))+
  geom_point() + xlab('Mice')


#CD4+
ggplot(data=treated, aes(x=Mice , y=CD4 ))+
  geom_point() + xlab('Mice')


#CD4+/FOXP3+
ggplot(data=treated, aes(x=Mice , y= CD4FOXP3 ))+
  geom_point() + xlab('Mice')


#CD8
ggplot(data=treated, aes(x=Mice , y=CD8 ))+
  geom_point() + xlab('Mice')

#CD45+
View(untreated)
ggplot(untreated, aes(x=Mice, y=CD45 )) +
  geom_point(size=2, shape=23)


#CD3+
ggplot(data=untreated, aes(x=Mice , y=CD3 ))+
  geom_point() + xlab('Mice')


#CD4+
ggplot(data=untreated, aes(x=Mice , y=CD4 ))+
  geom_point() + xlab('Mice')


#CD4+/FOXP3+
ggplot(data=untreated, aes(x=Mice , y= CD4FOXP3 ))+
  geom_point() + xlab('Mice')


#CD8
ggplot(data=untreated, aes(x=Mice , y=CD8 ))+
  geom_point() + xlab('Mice')

View(untreated_and_treated)
#CD45+
ggplot(data=untreated_and_treated, aes(x=Mice , y=CD45 , color= Treatment))+
  geom_point() + xlab('Mice')


#CD3+
ggplot(data=untreated_and_treated, aes(x=Mice , y=CD3 , color= Treatment))+
  geom_point() + xlab('Mice')


#CD4+
ggplot(data=untreated_and_treated, aes(x=Mice , y=CD4 , color= Treatment))+
  geom_point() + xlab('Mice')


#CD4+/FOXP3+
ggplot(data=untreated_and_treated, aes(x=Mice , y=CD4FOXP3 , color= Treatment))+
  geom_point() + xlab('Mice')


#CD8
ggplot(data=untreated_and_treated, aes(x=Mice , y=CD8 , color= Treatment))+
  geom_point() + xlab('Mice')

#box and whisker plot without 0.25 and 0.75 lines
View(untreated_and_treated)
ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD45", 
            desc_stat = "mean_sd")


ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD3", 
            desc_stat = "mean_sd")


ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD4", 
            desc_stat = "mean_sd")


ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD4FOXP3", 
            desc_stat = "mean_sd")


ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD8", 
            desc_stat = "mean_sd")

#box and whisker plot with 0.25 and 0.75 lines
ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD45", 
            desc_stat = "mean_sd",
            error.plot = "errorbar",            
            add = "mean")
`fun.y` is deprecated. Use `fun` instead.`fun.ymin` is deprecated. Use `fun.min` instead.`fun.ymax` is deprecated. Use `fun.max` instead.

ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD3", 
            desc_stat = "mean_sd",
            error.plot = "errorbar",            
            add = "mean")
`fun.y` is deprecated. Use `fun` instead.`fun.ymin` is deprecated. Use `fun.min` instead.`fun.ymax` is deprecated. Use `fun.max` instead.

ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD4", 
            desc_stat = "mean_sd",
            error.plot = "errorbar",            
            add = "mean")
`fun.y` is deprecated. Use `fun` instead.`fun.ymin` is deprecated. Use `fun.min` instead.`fun.ymax` is deprecated. Use `fun.max` instead.

ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD4FOXP3", 
            desc_stat = "mean_sd",
            error.plot = "errorbar",            
            add = "mean")
`fun.y` is deprecated. Use `fun` instead.`fun.ymin` is deprecated. Use `fun.min` instead.`fun.ymax` is deprecated. Use `fun.max` instead.

ggerrorplot(untreated_and_treated, x = "Treatment", y = "CD8", 
            desc_stat = "mean_sd",
            error.plot = "errorbar",            
            add = "mean")
`fun.y` is deprecated. Use `fun` instead.`fun.ymin` is deprecated. Use `fun.min` instead.`fun.ymax` is deprecated. Use `fun.max` instead.

#box and whisker plot  colored box 
ggplot(untreated_and_treated, aes(x=as.factor(Treatment), y=CD45)) + 
  geom_boxplot(fill="slateblue", alpha=0.2) + 
  xlab("Treatment")


ggplot(untreated_and_treated, aes(x=as.factor(Treatment), y=CD4)) + 
  geom_boxplot(fill="slateblue", alpha=0.2) + 
  xlab("Treatment")


ggplot(untreated_and_treated, aes(x=as.factor(Treatment), y=CD3)) + 
  geom_boxplot(fill="slateblue", alpha=0.2) + 
  xlab("Treatment")


ggplot(untreated_and_treated, aes(x=as.factor(Treatment), y=CD4FOXP3)) + 
  geom_boxplot(fill="slateblue", alpha=0.2) + 
  xlab("Treatment")


ggplot(untreated_and_treated, aes(x=as.factor(Treatment), y=CD8)) + 
  geom_boxplot(fill="slateblue", alpha=0.2) + 
  xlab("Treatment")

#box plot with individual points 
There were 35 warnings (use warnings() to see them)
a <- ggplot(untreated_and_treated, aes(Treatment, CD45))
a + geom_boxplot(outlier.shape = NA) + geom_jitter(width = 0.2)


b <- ggplot(untreated_and_treated, aes(Treatment, CD4))
b + geom_boxplot(outlier.shape = NA) + geom_jitter(width = 0.2)


c <- ggplot(untreated_and_treated, aes(Treatment, CD3))
c + geom_boxplot(outlier.shape = NA) + geom_jitter(width = 0.2)


d <- ggplot(untreated_and_treated, aes(Treatment, CD4FOXP3))
d + geom_boxplot(outlier.shape = NA) + geom_jitter(width = 0.2)


e <- ggplot(untreated_and_treated, aes(Treatment, CD8))
e + geom_boxplot(outlier.shape = NA) + geom_jitter(width = 0.2)

LS0tCnRpdGxlOiAic2NhdHRlcnBsb3QiCm91dHB1dDogaHRtbF9ub3RlYm9vawotLS0KCgotLS0KdGl0bGU6ICJSIE5vdGVib29rIgpvdXRwdXQ6CiAgaHRtbF9kb2N1bWVudDoKICAgIGRmX3ByaW50OiBwYWdlZAotLS0KCgpgYGB7cn0KI0NENDUrClZpZXcodHJlYXRlZCkKZ2dwbG90KHRyZWF0ZWQsIGFlcyh4PU1pY2UsIHk9Q0Q0NSApKSArCiAgZ2VvbV9wb2ludChzaXplPTIsIHNoYXBlPTIzKQoKI0NEMysKZ2dwbG90KGRhdGE9dHJlYXRlZCwgYWVzKHg9TWljZSAsIHk9Q0QzICkpKwogIGdlb21fcG9pbnQoKSArIHhsYWIoJ01pY2UnKQoKI0NENCsKZ2dwbG90KGRhdGE9dHJlYXRlZCwgYWVzKHg9TWljZSAsIHk9Q0Q0ICkpKwogIGdlb21fcG9pbnQoKSArIHhsYWIoJ01pY2UnKQoKI0NENCsvRk9YUDMrCmdncGxvdChkYXRhPXRyZWF0ZWQsIGFlcyh4PU1pY2UgLCB5PSBDRDRGT1hQMyApKSsKICBnZW9tX3BvaW50KCkgKyB4bGFiKCdNaWNlJykKCiNDRDgKZ2dwbG90KGRhdGE9dHJlYXRlZCwgYWVzKHg9TWljZSAsIHk9Q0Q4ICkpKwogIGdlb21fcG9pbnQoKSArIHhsYWIoJ01pY2UnKQoKYGBgCgoKYGBge3J9CiNDRDQ1KwpWaWV3KHVudHJlYXRlZCkKZ2dwbG90KHVudHJlYXRlZCwgYWVzKHg9TWljZSwgeT1DRDQ1ICkpICsKICBnZW9tX3BvaW50KHNpemU9Miwgc2hhcGU9MjMpCgojQ0QzKwpnZ3Bsb3QoZGF0YT11bnRyZWF0ZWQsIGFlcyh4PU1pY2UgLCB5PUNEMyApKSsKICBnZW9tX3BvaW50KCkgKyB4bGFiKCdNaWNlJykKCiNDRDQrCmdncGxvdChkYXRhPXVudHJlYXRlZCwgYWVzKHg9TWljZSAsIHk9Q0Q0ICkpKwogIGdlb21fcG9pbnQoKSArIHhsYWIoJ01pY2UnKQoKI0NENCsvRk9YUDMrCmdncGxvdChkYXRhPXVudHJlYXRlZCwgYWVzKHg9TWljZSAsIHk9IENENEZPWFAzICkpKwogIGdlb21fcG9pbnQoKSArIHhsYWIoJ01pY2UnKQoKI0NEOApnZ3Bsb3QoZGF0YT11bnRyZWF0ZWQsIGFlcyh4PU1pY2UgLCB5PUNEOCApKSsKICBnZW9tX3BvaW50KCkgKyB4bGFiKCdNaWNlJykKYGBgCgoKYGBge3J9ClZpZXcodW50cmVhdGVkX2FuZF90cmVhdGVkKQojQ0Q0NSsKZ2dwbG90KGRhdGE9dW50cmVhdGVkX2FuZF90cmVhdGVkLCBhZXMoeD1NaWNlICwgeT1DRDQ1ICwgY29sb3I9IFRyZWF0bWVudCkpKwogIGdlb21fcG9pbnQoKSArIHhsYWIoJ01pY2UnKQoKI0NEMysKZ2dwbG90KGRhdGE9dW50cmVhdGVkX2FuZF90cmVhdGVkLCBhZXMoeD1NaWNlICwgeT1DRDMgLCBjb2xvcj0gVHJlYXRtZW50KSkrCiAgZ2VvbV9wb2ludCgpICsgeGxhYignTWljZScpCgojQ0Q0KwpnZ3Bsb3QoZGF0YT11bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIGFlcyh4PU1pY2UgLCB5PUNENCAsIGNvbG9yPSBUcmVhdG1lbnQpKSsKICBnZW9tX3BvaW50KCkgKyB4bGFiKCdNaWNlJykKCiNDRDQrL0ZPWFAzKwpnZ3Bsb3QoZGF0YT11bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIGFlcyh4PU1pY2UgLCB5PUNENEZPWFAzICwgY29sb3I9IFRyZWF0bWVudCkpKwogIGdlb21fcG9pbnQoKSArIHhsYWIoJ01pY2UnKQoKI0NEOApnZ3Bsb3QoZGF0YT11bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIGFlcyh4PU1pY2UgLCB5PUNEOCAsIGNvbG9yPSBUcmVhdG1lbnQpKSsKICBnZW9tX3BvaW50KCkgKyB4bGFiKCdNaWNlJykKYGBgCgpgYGB7cn0KI2JveCBhbmQgd2hpc2tlciBwbG90IHdpdGhvdXQgMC4yNSBhbmQgMC43NSBsaW5lcwpWaWV3KHVudHJlYXRlZF9hbmRfdHJlYXRlZCkKZ2dlcnJvcnBsb3QodW50cmVhdGVkX2FuZF90cmVhdGVkLCB4ID0gIlRyZWF0bWVudCIsIHkgPSAiQ0Q0NSIsIAogICAgICAgICAgICBkZXNjX3N0YXQgPSAibWVhbl9zZCIpCgpnZ2Vycm9ycGxvdCh1bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIHggPSAiVHJlYXRtZW50IiwgeSA9ICJDRDMiLCAKICAgICAgICAgICAgZGVzY19zdGF0ID0gIm1lYW5fc2QiKQoKZ2dlcnJvcnBsb3QodW50cmVhdGVkX2FuZF90cmVhdGVkLCB4ID0gIlRyZWF0bWVudCIsIHkgPSAiQ0Q0IiwgCiAgICAgICAgICAgIGRlc2Nfc3RhdCA9ICJtZWFuX3NkIikKCmdnZXJyb3JwbG90KHVudHJlYXRlZF9hbmRfdHJlYXRlZCwgeCA9ICJUcmVhdG1lbnQiLCB5ID0gIkNENEZPWFAzIiwgCiAgICAgICAgICAgIGRlc2Nfc3RhdCA9ICJtZWFuX3NkIikKCmdnZXJyb3JwbG90KHVudHJlYXRlZF9hbmRfdHJlYXRlZCwgeCA9ICJUcmVhdG1lbnQiLCB5ID0gIkNEOCIsIAogICAgICAgICAgICBkZXNjX3N0YXQgPSAibWVhbl9zZCIpCgpgYGAKCmBgYHtyfQojYm94IGFuZCB3aGlza2VyIHBsb3Qgd2l0aCAwLjI1IGFuZCAwLjc1IGxpbmVzCmdnZXJyb3JwbG90KHVudHJlYXRlZF9hbmRfdHJlYXRlZCwgeCA9ICJUcmVhdG1lbnQiLCB5ID0gIkNENDUiLCAKICAgICAgICAgICAgZGVzY19zdGF0ID0gIm1lYW5fc2QiLAogICAgICAgICAgICBlcnJvci5wbG90ID0gImVycm9yYmFyIiwgICAgICAgICAgICAKICAgICAgICAgICAgYWRkID0gIm1lYW4iKQoKZ2dlcnJvcnBsb3QodW50cmVhdGVkX2FuZF90cmVhdGVkLCB4ID0gIlRyZWF0bWVudCIsIHkgPSAiQ0QzIiwgCiAgICAgICAgICAgIGRlc2Nfc3RhdCA9ICJtZWFuX3NkIiwKICAgICAgICAgICAgZXJyb3IucGxvdCA9ICJlcnJvcmJhciIsICAgICAgICAgICAgCiAgICAgICAgICAgIGFkZCA9ICJtZWFuIikKCmdnZXJyb3JwbG90KHVudHJlYXRlZF9hbmRfdHJlYXRlZCwgeCA9ICJUcmVhdG1lbnQiLCB5ID0gIkNENCIsIAogICAgICAgICAgICBkZXNjX3N0YXQgPSAibWVhbl9zZCIsCiAgICAgICAgICAgIGVycm9yLnBsb3QgPSAiZXJyb3JiYXIiLCAgICAgICAgICAgIAogICAgICAgICAgICBhZGQgPSAibWVhbiIpCgpnZ2Vycm9ycGxvdCh1bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIHggPSAiVHJlYXRtZW50IiwgeSA9ICJDRDRGT1hQMyIsIAogICAgICAgICAgICBkZXNjX3N0YXQgPSAibWVhbl9zZCIsCiAgICAgICAgICAgIGVycm9yLnBsb3QgPSAiZXJyb3JiYXIiLCAgICAgICAgICAgIAogICAgICAgICAgICBhZGQgPSAibWVhbiIpCgpnZ2Vycm9ycGxvdCh1bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIHggPSAiVHJlYXRtZW50IiwgeSA9ICJDRDgiLCAKICAgICAgICAgICAgZGVzY19zdGF0ID0gIm1lYW5fc2QiLAogICAgICAgICAgICBlcnJvci5wbG90ID0gImVycm9yYmFyIiwgICAgICAgICAgICAKICAgICAgICAgICAgYWRkID0gIm1lYW4iKQoKYGBgCgpgYGB7cn0KI2JveCBhbmQgd2hpc2tlciBwbG90ICBjb2xvcmVkIGJveCAKZ2dwbG90KHVudHJlYXRlZF9hbmRfdHJlYXRlZCwgYWVzKHg9YXMuZmFjdG9yKFRyZWF0bWVudCksIHk9Q0Q0NSkpICsgCiAgZ2VvbV9ib3hwbG90KGZpbGw9InNsYXRlYmx1ZSIsIGFscGhhPTAuMikgKyAKICB4bGFiKCJUcmVhdG1lbnQiKQoKZ2dwbG90KHVudHJlYXRlZF9hbmRfdHJlYXRlZCwgYWVzKHg9YXMuZmFjdG9yKFRyZWF0bWVudCksIHk9Q0Q0KSkgKyAKICBnZW9tX2JveHBsb3QoZmlsbD0ic2xhdGVibHVlIiwgYWxwaGE9MC4yKSArIAogIHhsYWIoIlRyZWF0bWVudCIpCgpnZ3Bsb3QodW50cmVhdGVkX2FuZF90cmVhdGVkLCBhZXMoeD1hcy5mYWN0b3IoVHJlYXRtZW50KSwgeT1DRDMpKSArIAogIGdlb21fYm94cGxvdChmaWxsPSJzbGF0ZWJsdWUiLCBhbHBoYT0wLjIpICsgCiAgeGxhYigiVHJlYXRtZW50IikKCmdncGxvdCh1bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIGFlcyh4PWFzLmZhY3RvcihUcmVhdG1lbnQpLCB5PUNENEZPWFAzKSkgKyAKICBnZW9tX2JveHBsb3QoZmlsbD0ic2xhdGVibHVlIiwgYWxwaGE9MC4yKSArIAogIHhsYWIoIlRyZWF0bWVudCIpCgpnZ3Bsb3QodW50cmVhdGVkX2FuZF90cmVhdGVkLCBhZXMoeD1hcy5mYWN0b3IoVHJlYXRtZW50KSwgeT1DRDgpKSArIAogIGdlb21fYm94cGxvdChmaWxsPSJzbGF0ZWJsdWUiLCBhbHBoYT0wLjIpICsgCiAgeGxhYigiVHJlYXRtZW50IikKYGBgCgpgYGB7cn0KI2JveCBwbG90IHdpdGggaW5kaXZpZHVhbCBwb2ludHMgCmEgPC0gZ2dwbG90KHVudHJlYXRlZF9hbmRfdHJlYXRlZCwgYWVzKFRyZWF0bWVudCwgQ0Q0NSkpCmEgKyBnZW9tX2JveHBsb3Qob3V0bGllci5zaGFwZSA9IE5BKSArIGdlb21faml0dGVyKHdpZHRoID0gMC4yKQoKYiA8LSBnZ3Bsb3QodW50cmVhdGVkX2FuZF90cmVhdGVkLCBhZXMoVHJlYXRtZW50LCBDRDQpKQpiICsgZ2VvbV9ib3hwbG90KG91dGxpZXIuc2hhcGUgPSBOQSkgKyBnZW9tX2ppdHRlcih3aWR0aCA9IDAuMikKCmMgPC0gZ2dwbG90KHVudHJlYXRlZF9hbmRfdHJlYXRlZCwgYWVzKFRyZWF0bWVudCwgQ0QzKSkKYyArIGdlb21fYm94cGxvdChvdXRsaWVyLnNoYXBlID0gTkEpICsgZ2VvbV9qaXR0ZXIod2lkdGggPSAwLjIpCgpkIDwtIGdncGxvdCh1bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIGFlcyhUcmVhdG1lbnQsIENENEZPWFAzKSkKZCArIGdlb21fYm94cGxvdChvdXRsaWVyLnNoYXBlID0gTkEpICsgZ2VvbV9qaXR0ZXIod2lkdGggPSAwLjIpCgplIDwtIGdncGxvdCh1bnRyZWF0ZWRfYW5kX3RyZWF0ZWQsIGFlcyhUcmVhdG1lbnQsIENEOCkpCmUgKyBnZW9tX2JveHBsb3Qob3V0bGllci5zaGFwZSA9IE5BKSArIGdlb21faml0dGVyKHdpZHRoID0gMC4yKQpgYGAKCg==