Я хотел бы сгруппировать данные два на два, зеленый (fociPositiveYFPPoleta) / красный (fociPositiveCapns1) в разное время 1, 3, 6, 9, 24 и 48 часов. Это возможно?
#Poleta and Capns1 foci per h
fociPositiveYFPPoleta <- c(percentagepositivegfppoleta1hUV15, percentagepositivegfppoleta3hUV15, percentagepositivegfppoleta6hUV15, percentagepositivegfppoleta9hUV15, percentagepositivegfppoleta24hUV15, percentagepositivegfppoleta48hUV15)
fociPositiveCapns1 <- c(percentagepositivecapns11hUV15, percentagepositivecapns13hUV15, percentagepositivecapns16hUV15, percentagepositivecapns19hUV15, percentagepositivecapns124hUV15, percentagepositivecapns148hUV15)
allsamples3 <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12)
rbind (allsamples3, fociPositiveYFPPoleta, fociPositiveCapns1)
#Barplot
dev.new(height=8,width=12)
par(mar=c(10.1, 4.1, 4.1, 2.1))
barplot(fociPositiveYFPPoleta, fociPositiveCapns1, type="0", col=c("#009966", "#990000"), main = "%tage of cells with > 30 Capns1 Foci/Cell",
xlab = "",
ylab = "%tage of cells with > 30 Capns1 foci/cell",ylim = c(0,100))
text(1:length(fociPositiveYFPPoleta, fociPositiveCapns1), -3, labels = c("1h", "1h", "3h", "3h", "6h", "6h", "9h", "9h", "24h", "24h", "48h", "48h"), adj=0.8, srt=0, xpd=TRUE)
legend(x=0,y=100,c("UV15", "UV15"),cex=1,col=c("#009966","#990000"),pch="x")
dev.print(pdf, "30CapnspercentOTUV15Barplot")