R : Copyright 2004, The R Foundation for Statistical Computing
Version 1.9.0  (2004-04-12), ISBN 3-900051-00-3

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.

[Previously saved workspace restored]

> invisible(options(echo = TRUE))
> 
> library(gregmisc);

Attaching package 'gregmisc':


	The following object(s) are masked from package:stats :

	 lowess 

> 
> #fig 1
> 
> A <- matrix(scan("epoch/meta-summary1.data"),nrow=4);
Read 72 items
> 
> bitmap(file="new-summary1.png", 
+        type="png16m", 
+        width=10,
+        height=4, 
+        res=150);
> 
> for (i in 1:length(A)) {
+  A[i] <- A[i] + 1;
+ }
> 
> par (mai=c(.4,.6,.2,0), xpd=NA);
> 
> barplot2(A, col=c("cyan", "green", "blue", "black"), horiz=F,  
+          ylab = "privilege",
+          xlab = "epochs (one bar is a column of metaworld data)",
+          log="y", 
+          cex.axis=2, cex.lab=3, beside=T);
> 
> dev.off();
null device 
          1 
> 
> 
> 
> #fig2a
> 
> A <- matrix(scan("epoch/meta-summary2a.data"), nrow=2); 
Read 36 items
> 
> bitmap(file="new-summary2a.png", 
+        type="png16m", 
+        width=5,
+        height=2, 
+        res=150);
> 
> par (mai=c(.4,.6,.2,0), xpd=NA);
> 
> barplot(A, col=c("red","black"),
+          horiz=F, space=1,
+          ylab = "abundance",
+          xlab = "epochs", 
+          cex.axis=2, cex.lab=3);
> 
> 
> dev.off();
null device 
          1 
> 
> 
> 
> #fig 2b
> 
> A <- matrix(scan("epoch/meta-summary2b.data"), nrow=2); 
Read 36 items
> 
> bitmap(file="new-summary2b.png", 
+        type="png16m", 
+        width=5,
+        height=2, 
+        res=150);
> 
> par (mai=c(.4,.6,.2,0), xpd=NA);
> 
> barplot(A, col=c("yellow","black"),
+          horiz=F, space=1,
+          ylab = "abundance", xlab = "epochs", 
+          cex.axis=2, cex.lab=3);
> 
> dev.off();
null device 
          1 
> 
> 
> 
> #fig 2c 
> 
> bitmap(file="new-summary2c.png", 
+        type="png16m", 
+        width=5,
+        height=2, 
+        res=150);
> 
> par (mai=c(.4,.6,.2,0), xpd=NA);
> 
> A <- matrix(scan("epoch/meta-summary2c.data"), nrow=2); 
Read 36 items
> 
> barplot(A, col=c("purple","black"),
+          horiz=F, space=1,
+          ylab = "abundance", xlab = "epochs", 
+          cex.axis=2, cex.lab=3);
> 
> 
> dev.off();
null device 
          1 
> 
> 
> 
> #fig 2d
> 
> A <- matrix(scan("epoch/meta-summary2d.data"), nrow=2); 
Read 36 items
> 
> bitmap(file="new-summary2d.png", 
+        type="png16m", 
+        width=5,
+        height=2, 
+        res=150);
> 
> par (mai=c(.4,.6,.2,0), xpd=NA);
> 
> barplot(A, col=c("gray","black"),
+          horiz=F, space=1,
+          ylab = "abundance", xlab = "epochs", 
+          cex.axis=2, cex.lab=3);
> 
> dev.off();
null device 
          1 
> 
> 
> proc.time()
[1] 1.7031250 0.1015625 3.9296875 0.9531250 0.1640625
>