comparison shm_csr.r @ 2:e85fec274cde draft

Uploaded

1:faae21ba5c63

	pc = ggplot(genesForPlot, aes(x = factor(1), y=Freq, fill=Gene))
	pc = pc + geom_bar(width = 1, stat = "identity") + scale_fill_manual(labels=genesForPlot$label, values=c("IGA1" = "lightblue1", "IGA2" = "blue4"))
	pc = pc + coord_polar(theta="y") + scale_y_continuous(breaks=NULL)
	pc = pc + theme(panel.background = element_rect(fill = "white", colour="black"), text = element_text(size=13, colour="black"))
	pc = pc + xlab(" ") + ylab(" ") + ggtitle(paste("IGA subclasses", "( n =", sum(genesForPlot$Freq), ")"))
	write.table(genesForPlot, "IGA.txt", sep="\t",quote=F,row.names=F,col.names=T)

	png(filename="IGA.png")
	print(pc)
	dev.off()
}

	pc = ggplot(genesForPlot, aes(x = factor(1), y=Freq, fill=Gene))
	pc = pc + geom_bar(width = 1, stat = "identity") + scale_fill_manual(labels=genesForPlot$label, values=c("IGG1" = "olivedrab3", "IGG2" = "red", "IGG3" = "gold", "IGG4" = "darkred"))
	pc = pc + coord_polar(theta="y") + scale_y_continuous(breaks=NULL)
	pc = pc + theme(panel.background = element_rect(fill = "white", colour="black"), text = element_text(size=13, colour="black"))
	pc = pc + xlab(" ") + ylab(" ") + ggtitle(paste("IGG subclasses", "( n =", sum(genesForPlot$Freq), ")"))
	write.table(genesForPlot, "IGG.txt", sep="\t",quote=F,row.names=F,col.names=T)

	png(filename="IGG.png")
	print(pc)
	dev.off()
}

2:e85fec274cde

	pc = ggplot(genesForPlot, aes(x = factor(1), y=Freq, fill=Gene))
	pc = pc + geom_bar(width = 1, stat = "identity") + scale_fill_manual(labels=genesForPlot$label, values=c("IGA1" = "lightblue1", "IGA2" = "blue4"))
	pc = pc + coord_polar(theta="y") + scale_y_continuous(breaks=NULL)
	pc = pc + theme(panel.background = element_rect(fill = "white", colour="black"), text = element_text(size=13, colour="black"))
	pc = pc + xlab(" ") + ylab(" ") + ggtitle(paste("IGA subclasses", "( n =", sum(genesForPlot$Freq), ")"))
	write.table(genesForPlot, "IGA_pie.txt", sep="\t",quote=F,row.names=F,col.names=T)

	png(filename="IGA.png")
	print(pc)
	dev.off()
}

	pc = ggplot(genesForPlot, aes(x = factor(1), y=Freq, fill=Gene))
	pc = pc + geom_bar(width = 1, stat = "identity") + scale_fill_manual(labels=genesForPlot$label, values=c("IGG1" = "olivedrab3", "IGG2" = "red", "IGG3" = "gold", "IGG4" = "darkred"))
	pc = pc + coord_polar(theta="y") + scale_y_continuous(breaks=NULL)
	pc = pc + theme(panel.background = element_rect(fill = "white", colour="black"), text = element_text(size=13, colour="black"))
	pc = pc + xlab(" ") + ylab(" ") + ggtitle(paste("IGG subclasses", "( n =", sum(genesForPlot$Freq), ")"))
	write.table(genesForPlot, "IGG_pie.txt", sep="\t",quote=F,row.names=F,col.names=T)

	png(filename="IGG.png")
	print(pc)
	dev.off()
}