Log in Page Discussion History Go to the site toolbox

540Project:Leonrcode

From BluWiki

Data files for Cho

#ARIS#

#Plot 3 graphs on one plot

#how much do we scale the text
clab <- 1.5
cmain <- 2.3
clegend <- 1.5

par(mfrow=c(1,3))

#Cell cycle data
kim<-c(0.1612334, 0.2219552, 0.3117359 ,0.2899309, 0.2549579, 0.3304386, 0.2734528, 0.2818993,0.2977723, 0.3501287 ,0.3837607, 0.3989134, 0.4427269,
 0.3498900, 0.4376838, 0.4682092, 0.4731219)
plot(cex.lab=clab, cex.main=cmain,1:17,kim,xaxp=c(1,17,2),ylim=c(0,1),xlim=c(1,17),ylab="ARI (measure of agreeement)",xlab="Number of Principal Components",main="Standardized Cell Cycle Data",pch=15)
lines(1:17,kim,pch=16)

#Xhang's points
points(1:5,c( 0.3690727, 0.4484942, 0.4701616, 0.4860360, 0.4860360),col="red",pch=16)
lines(1:5,c( 0.3690727, 0.4484942, 0.4701616, 0.4860360, 0.4860360),col="red",pch=16)

#Leon's points 
shortMKL<-c(0.040748943, 0.048012299 ,0.027792279, 0.011486479, 0.013285407,
0.013892259, 0.016708111, 0.026103777, 0.014887695, 0.008729655,
0.015223416, 0.017379713 ,0.013759798, 0.024321858, 0.023029511,)

longMKL<-c(0.4536869,0.4390780,0.4678440,0.4200678)
fileNum<-c(6,9,12,15)

points(2:16, shortMKL, col="blue",pch=17)
lines(2:16, shortMKL, col="blue",pch=17)

points(fileNum, longMKL, col="purple",pch=8)
lines(fileNum, longMKL, col="purple",pch=8)
legend(cex=clegend, locator(1), c("PCA", "fPCA", "Short MKL Representation", "Long MKL Representation"),col=c("black",,"red","blue","purple"),pch=c(15,16,17,8))


#Plot Less Noisy
kim<-> c(0.5610905, 0.9227669, 1.0000000 ,1.0000000,
1.0000000, 1.0000000 ,1.0000000 ,1.0000000,
1.0000000, 1.0000000, 1.0000000, 1.0000000,
1.0000000, 1.0000000, 1.0000000, 1.0000000,
1.0000000, 1.0000000, 1.0000000, 1.0000000)

numKimTimePoints<-length(kim)
numTimePoints<-20

plot(cex.lab=clab, cex.main=cmain,1:numKimTimePoints,kim,ylim=c(0,1),xlim=c(1,17),ylab="ARI (measure of agreeement)",xlab="Number of Principal Components",main="Less Noisy Simulation, 3 Clusters",pch=15)
lines(1:numKimTimePoints,kim,pch=15)

#Xhang's points
points(1:5,c(0.9899832, 1.0000000 ,1.0000000 ,1.0000000 ,1.0000000),col="red",pch=16)
lines(1:5,c(0.9899832, 1.0000000 ,1.0000000, 1.0000000 ,1.0000000),col="red",pch=16)

#Leon's points 
shortMKL<-c(0.01324889, 0.01290301, 0.01075094 ,0.01220738, 0.01390845, 0.01170633, 0.01594497, 0.01157196, 0.01297928, 0.01150743, 0.02118876, 0.02118876,
0.02118876, 0.02118876, 0.01743758, 0.00000000 ,0.00000000, 0.00000000)
points(2:(numTimePoints-1), shortMKL, col="blue",pch=17)
lines(2:(numTimePoints-1), shortMKL, col="blue",pch=17)

longMKL<- c(0.980068,  1.000000,1.000000, 1.000000)
fileNum<-c(6,9,12,15)

points(fileNum, longMKL, col="purple",pch=8)
lines(fileNum, longMKL, col="purple",pch=8)
legend(cex=clegend,locator(1), c("PCA", "fPCA", "Short MKL Representation", "Long MKL Representation"),col=c("black",,"red","blue","purple"),pch=c(15,16,17,8))


#Plot Ultra Noisy

kim<-c( 0.1500926, 0.1374338, 0.1419731, 0.4291463, 0.5220689, 0.5372731, 0.5988467, 0.6308278, 0.6535010, 0.6597730, 0.6525915, 0.6816172, 0.6997159,
 0.6756447, 0.6971203, 0.7251389, 0.8015260, 0.7931116, 0.7839299, 0.8080807)

numKimTimePoints<-length(kim)
numTimePoints<-20

plot(cex.lab=clab, cex.main=cmain,1:numKimTimePoints,kim,ylim=c(0,1),xlim=c(1,numTimePoints),ylab="ARI (measure of agreeement)",xlab="Number of Principal Components",main="Noisy Simulation, 3 Clusters",pch=15)
lines(1:numTimePoints,kim,pch=16)

#Xhang's points
points(1:5,c(0.4407172, 0.9314342 ,0.9314342 ,0.9314342, 0.9314342),col="red",pch=16)
lines(1:5,c(0.4407172, 0.9314342 ,0.9314342 ,0.9314342, 0.9314342),col="red",pch=16)

#Leon's points 

shortMKL<-c(0.026792891,  0.012161651 , 0.025596719,  0.006320796 , 0.003200292,0.003765788,  0.013165935 , 0.008788605,  0.002590373 , 0.022641987,
0.004125476,  0.006155053,  0.029017783, -0.001625406,  0.001739478, 0.000000000,  0.000000000,  0.000000000)

longMKL<- c(0.4536869, 0.4390780, 0.4678440, 0.4200678)
fileNum<-c(6,9,12,15)

points(2:(numTimePoints-1), shortMKL, col="blue",pch=17)
lines(2:(numTimePoints-1), shortMKL, col="blue",pch=17)

points(fileNum, longMKL, col="purple",pch=8)
lines(fileNum, longMKL, col="purple",pch=8)

legend(cex=clegend,locator(1), c("PCA", "fPCA", "Short MKL Representation", "Long MKL Representation"),col=c("black",,"red","blue","purple"),pch=c(15,16,17,8))




Site Toolbox:

Personal tools
GNU Free Documentation License 1.2
This page was last modified on 12 April 2006, at 04:08.
Disclaimers - About BluWiki