data <- read.csv('H:/Collaboration/Clinic/Mary_Alice_Nading/BoxPlotData.csv') attach(data) GROUP <- as.factor(GROUP) # par(mfrow=c(1,2)) boxplot(NUCLEARE ~ GROUP, subset=(GROUP==1 | GROUP==2 | GROUP==3 | GROUP==4 | GROUP==5 | GROUP==6), xlab = 'GROUP', ylab = "Nuclear score", ylim = c(-1,6), boxwex = 0.25, at=1:6-0.2, col = "yellow") boxplot(NUCLEARD ~ GROUP, subset=(GROUP==1 | GROUP==2 | GROUP==3 | GROUP==4 | GROUP==5 | GROUP==6), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:6+0.2, col = "orange") legend(2, 6, c("Nuclear Epidermal", "Nuclear Dermal"), fill = c("yellow", "orange")) title('Box Plot of Nuclear score') boxplot(CYTOPLAS ~ GROUP, subset=(GROUP==1 | GROUP==2 | GROUP==3 | GROUP==4 | GROUP==5 | GROUP==6), xlab = 'GROUP', ylab = "Cytoplasmic score", ylim = c(-1,6), boxwex = 0.25, at=1:6-0.2, col = "yellow") boxplot(CYTOPL_A ~ GROUP, subset=(GROUP==1 | GROUP==2 | GROUP==3 | GROUP==4 | GROUP==5 | GROUP==6), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:6+0.2, col = "orange") legend(2, 6, c("Cytoplasmic Epidermal", "Cytoplasmic Dermal"), fill = c("yellow", "orange")) title('Box Plot of Cytoplasmic score') boxplot(NUCLEARE ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==10 | GROUP==11), xlab = 'GROUP', ylab = "Nuclear score", ylim = c(-1,6), boxwex = 0.25, at=1:4-0.2, col = "yellow") boxplot(NUCLEARD ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==10 | GROUP==11), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:4+0.2, col = "orange") legend(2, 6, c("Nuclear Epidermal", "Nuclear Dermal"), fill = c("yellow", "orange")) title('Box Plot of Nuclear score') boxplot(CYTOPLAS ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==10 | GROUP==11), xlab = 'GROUP', ylab = "Cytoplasmic score", ylim = c(-1,6), boxwex = 0.25, at=1:4-0.2, col = "yellow") boxplot(CYTOPL_A ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==10 | GROUP==11), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:4+0.2, col = "orange") legend(2, 6, c("Cytoplasmic Epidermal", "Cytoplasmic Dermal"), fill = c("yellow", "orange")) title('Box Plot of Cytoplasmic score') boxplot(NUCLEARE ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==6 | GROUP==7 | GROUP==8 | GROUP==9), xlab = 'GROUP', ylab = "Nuclear score", ylim = c(-1,6), boxwex = 0.25, at=1:6-0.2, col = "yellow") boxplot(NUCLEARD ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==6 | GROUP==7 | GROUP==8 | GROUP==9), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:6+0.2, col = "orange") legend(2, 6, c("Nuclear Epidermal", "Nuclear Dermal"), fill = c("yellow", "orange")) title('Box Plot of Nuclear score') boxplot(CYTOPLAS ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==6 | GROUP==7 | GROUP==8 | GROUP==9), xlab = 'GROUP', ylab = "Cytoplasmic score", ylim = c(-1,6), boxwex = 0.25, at=1:6-0.2, col = "yellow") boxplot(CYTOPL_A ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==6 | GROUP==7 | GROUP==8 | GROUP==9), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:6+0.2, col = "orange") legend(2, 6, c("Cytoplasmic Epidermal", "Cytoplasmic Dermal"), fill = c("yellow", "orange")) title('Box Plot of Cytoplasmic score') boxplot(NUCLEARE ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==13 | GROUP==14), xlab = 'GROUP', ylab = "Nuclear score", ylim = c(-1,6), boxwex = 0.25, at=1:4-0.2, col = "yellow") boxplot(NUCLEARD ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==13 | GROUP==14), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:4+0.2, col = "orange") legend(2, 6, c("Nuclear Epidermal", "Nuclear Dermal"), fill = c("yellow", "orange")) title('Box Plot of Nuclear score') boxplot(CYTOPLAS ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==13 | GROUP==14), xlab = 'GROUP', ylab = "Cytoplasmic score", ylim = c(-1,6), boxwex = 0.25, at=1:4-0.2, col = "yellow") boxplot(CYTOPL_A ~ GROUP, subset=(GROUP==4 | GROUP==5 | GROUP==13 | GROUP==14), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:4+0.2, col = "orange") legend(2, 6, c("Cytoplasmic Epidermal", "Cytoplasmic Dermal"), fill = c("yellow", "orange")) title('Box Plot of Cytoplasmic score') boxplot(NUCLEARE ~ GROUP, subset=(GROUP==1 | GROUP==2 | GROUP==3 | GROUP==4 | GROUP==5 | GROUP==6 | GROUP==10 | GROUP==11 | GROUP==12), xlab = 'GROUP', ylab = "Nuclear score", ylim = c(-1,6), boxwex = 0.25, at=1:9-0.2, col = "yellow") boxplot(NUCLEARD ~ GROUP, subset=(GROUP==1 | GROUP==2 | GROUP==3 | GROUP==4 | GROUP==5 | GROUP==6 | GROUP==10 | GROUP==11 | GROUP==12), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:9+0.2, col = "orange") legend(2, 6, c("Nuclear Epidermal", "Nuclear Dermal"), fill = c("yellow", "orange")) title('Box Plot of Nuclear score') boxplot(CYTOPLAS ~ GROUP, subset=(GROUP==1 | GROUP==2 | GROUP==3 | GROUP==4 | GROUP==5 | GROUP==6 | GROUP==10 | GROUP==11 | GROUP==12), xlab = 'GROUP', ylab = "Cytoplasmic score", ylim = c(-1,6), boxwex = 0.25, at=1:9-0.2, col = "yellow") boxplot(CYTOPL_A ~ GROUP, subset=(GROUP==1 | GROUP==2 | GROUP==3 | GROUP==4 | GROUP==5 | GROUP==6 | GROUP==10 | GROUP==11 | GROUP==12), ylim = c(-1,6), boxwex = 0.25, add=TRUE, at=1:9+0.2, col = "orange") legend(2, 6, c("Cytoplasmic Epidermal", "Cytoplasmic Dermal"), fill = c("yellow", "orange")) title('Box Plot of Cytoplasmic score')