WebSVN – bdplot – Diff – /bdgraphics/R/4_bdplot.R


### bdplot.R
#### This is a collection of functions that must bpre reachable by the user This is called 
 
 
 
## strategy missing for xytype
## Could we test if there is an active graphics device? And if not, run bdopen()?
bdplot <- function(x, y=NULL, z=NULL,
                   ## The rest is irrelevant for the user
                   cex.plot=bdp()$cex.plot,
                   xaxt="n",
                   xlab=bdp()$xlab, ylab=bdp()$ylab,
                   xlim=NULL,
                   ylim=NULL,
                   col=bdp()$col,main=NULL,
                   ## only used in case of plot mode
                   type=bdp()$type,
                   ## Only used in case of histogram plotting
                   border=bdp()$border,
                   ...)
{
  prm <- bdp()
 
 
  par(cex.main=prm$cex.main, cex.lab=prm$cex.lab, cex.axis=prm$cex.axis)
  par(tcl=prm$tcl)
  par(lwd=prm$lwd)
  par(bty=prm$bty)
 
  
  ##Set margins
  mar <- prm$mar.nolab
  if( !is.na(prm$xlab) ){ mar[1] <- prm$mar.lab[1] }
  if( !is.na(prm$ylab) ){ mar[2] <- prm$mar.lab[2] }
  if( any(!is.na(prm$toplab),!is.null(main)) ) { mar[3] <- prm$mar.lab[3] }
  if( !is.na(prm$rightlab) ){ mar[4] <- prm$mar.lab[4] }
  par(mar=mar)
 
  par(mgp=prm$mgp.global)
  
  ## Do the plot
  if(!is.null(y)){
    ##       cat("The Bacher/Delff Plotting System (R), (C), TM operating in xy-plotting mode.\n")
 
    if(is.null(xlim))
      xlim <- range(x,bdp()$xlim,na.rm=TRUE)
    if(is.null(ylim))
      ylim <- range(y,bdp()$ylim,na.rm=TRUE)
    
    plot(x, y,
         type=type,
         xaxt=xaxt, yaxt="n",
         cex=cex.plot,
         xlab=xlab,
         ylab=ylab,
         xlim=xlim,
         ylim=ylim,
         col=col,
         main=main,
         ...)
  } else if (is.numeric(x)){
    if(prm$method=="barplot"){
      ## this does not provide the full x,y functionality of barplot
      barplot(height=x,
              ## maybe these two are wrong/stupid?
              cex.axis=prm$cex.lab,
              cex.names=prm$cex.lab,
              col=bdp()$hcol,
              border=border,
              xlab=xlab,
              ylab=ylab,
              main=main,
              ...)
      prm$draw.xaxis <- FALSE
      prm$draw.yaxis <- FALSE
      
    } else {
 
      if(is.null(xlim))
        xlim <- range(c(1:length(x),bdp()$xlim),na.rm=TRUE)
 
      if(is.null(ylim))
        ylim <- range(x,bdp()$ylim,na.rm=TRUE)
      
      plot(x,
           type=type,
           xaxt="n", yaxt="n",
           cex=cex.plot,
           xlab=xlab,
           ylab=ylab,
           xlim=xlim,
           ylim=ylim,
           col=col,
           main=main,
           ...)
    }
  } else if( class(x) == "acf"){
    ## Not cleaned up/checked
    plot(x, xlab="", ylab="", xaxt="n", yaxt="n", cex=prm$cex.plot, xlim=prm$xlim, ylim=prm$ylim,col=col,...)
  }  else if( class(x) == "histogram"){
    ## Notice: the color of the bars can ONLY be set width bdp(hcol="color").
    plot(x,
         xlab=xlab, ylab=ylab,
         main=prm$main,
         xaxt="n",yaxt="n",col=bdp()$hcol,
         border=border,
         ...)
  } else if( class(x) == "trellis"){
    ## Not cleaned up/checked
    ## very experimental
    cat("The Bacher/Delff Plotting System (R), (C), TM operating in lattice mode.\n Remeber that labels must written in the lattice object.\n")
    ##lattice.options(layout.widths = prm$lattice.width,
    ##                layout.heights = prm$lattice.height)
    trellis.par.set(prm$myLatticeSettings()) 
    plot(x,col=col,...)
    ## with lattice/trellis, the axis drawing doesn't work.
    prm$draw.xaxis <- FALSE
    prm$draw.yaxis <- FALSE
  } else if (class(x)=="princomp") {
    plot(x,
         main=prm$main,
         ...)
  } else if (prm$method=="image.plot"){
    ### This could be checked in the beginning by is.null(z)
    cat("The Bacher/Delff Plotting System (R), (C), TM operating in image plotting mode.\nx and y axis can not be configured.\n")
    image.plot(x=x,y=y,z=z,
               xlab=xlab,
               ylab=ylab,
               ...)
  } 
  
###Grid stuff
  ## This really hould be done before adding the rest of the plot contents.
  ## We could make default values for grid.v og grid.h. Man maa kunne lave noget kvalificeret ud fra range og noget heltalsdivision
  
 
  if( prm$grid ){
    grid(col=prm$grid.col)
  }
  ## vertical lines
  if( !is.na(prm$grid.v) ){
    if(prm$grid.v=="Def"){
      grid(nx=NULL,ny=NA,col=prm$grid.col)
    } else
    {
      abline(v=prm$grid.v, lty="dotted", col=prm$grid.col)
    }
  }
 
  if( !is.na(prm$grid.h) ){
    if(prm$grid.h=="Def"){
      grid(nx=NA,ny=NULL,col=prm$grid.col)
    }else{
      abline(h=prm$grid.h, lty="dotted", col=prm$grid.col)}
    }
 
### Axis stuff This part should be improved. A function that draws
### axis should be run for all four axis. And it should be possible to
### supply whatever vector to base it on (especially relevant for
### taxis and raxis)
  
  ## this is because the use of mgp gives a ridiculous warning when too small
  options(warn=-1)
  if(prm$draw.xaxis){
    ##    mgp.old <- par()$mgp
    ##    par(mgp) <- mgp.xaxis
    ## is the abscissa a time object?
    if( "POSIXt"%in%class(x[1])){
      axis.POSIXct(1, x,mgp=prm$mgp.xaxis, lwd=prm$lwd)
    } else {
      axis(1, mgp=prm$mgp.xaxis, lwd=prm$lwd)
    }
 
    
  }
 
  if(prm$draw.yaxis){ axis(2, y,mgp=prm$mgp.yaxis, lwd=prm$lwd) }
  ## top axis
  if(prm$draw.taxis){ axis(3, mgp=prm$mgp.taxis, lwd=prm$lwd) }
  ## axis to the right
  if(prm$draw.raxis){ axis(4, mgp=prm$mgp.raxis, lwd=prm$lwd) }
 
  ## switch back on warnings
  options(warn=0)
 
 
  ##Title stuff
  ## Hvorfor??
  scale <- 1
  ##  if( prm$type=="hist" & !is.na(prm$xlab) )
#  if( !is.na(prm$xlab)  ){ mtext(prm$xlab, line=prm$xlabLine, side=1, cex=prm$cex.lab/scale) }
  ##  if( prm$type=="hist" &!is.na(prm$ylab) )
#  if( !is.na(prm$ylab)  ){ mtext(prm$ylab, line=prm$ylabLine, side=2, cex=prm$cex.lab/scale) }
  if( !is.na(prm$toplab) ){ mtext(prm$toplab, side=3, line=0.25, cex=prm$cex.lab/scale) }
  if( !is.na(prm$rightlab) ){
    mtext(prm$rightlab, side=4, line=0.75, cex=prm$cex.lab/scale,mgp=prm$mgp.raxis)
##    mtext(rightlab, side=4, line=0.5, cex=bdp()$cex.lab,mgp=bdp()$mgp.raxis)
  }
}
 
 
 
 

Rev 29	Rev 31
1	`### bdplot.R`	1	`### bdplot.R`
2	`#### This is a collection of functions that must bpre reachable by the user This is called`	2	`#### This is a collection of functions that must bpre reachable by the user This is called`
3		3
4		4
5		5
6	`## strategy missing for xytype`	6	`## strategy missing for xytype`
7	`## Could we test if there is an active graphics device? And if not, run bdopen()?`	7	`## Could we test if there is an active graphics device? And if not, run bdopen()?`
8	`bdplot <- function(x, y=NULL, z=NULL,`	8	`bdplot <- function(x, y=NULL, z=NULL,`
9	`## The rest is irrelevant for the user`	9	`## The rest is irrelevant for the user`
10	`cex.plot=bdp()$cex.plot,`	10	`cex.plot=bdp()$cex.plot,`
11	`xaxt="n",`	11	`xaxt="n",`
12	`xlab=bdp()$xlab, ylab=bdp()$ylab,`	12	`xlab=bdp()$xlab, ylab=bdp()$ylab,`
13	`xlim=NULL,`	13	`xlim=NULL,`
14	`ylim=NULL,`	14	`ylim=NULL,`
15	`col=bdp()$col,main=NULL,`	15	`col=bdp()$col,main=NULL,`
16	`## only used in case of plot mode`	16	`## only used in case of plot mode`
17	`type=bdp()$type,`	17	`type=bdp()$type,`
18	`## Only used in case of histogram plotting`	18	`## Only used in case of histogram plotting`
19	`border=bdp()$border,`	19	`border=bdp()$border,`
20	`...)`	20	`...)`
21	`{`	21	`{`
22	`prm <- bdp()`	22	`prm <- bdp()`
23		23
24		24
25	`par(cex.main=prm$cex.main, cex.lab=prm$cex.lab, cex.axis=prm$cex.axis)`	25	`par(cex.main=prm$cex.main, cex.lab=prm$cex.lab, cex.axis=prm$cex.axis)`
26	`par(tcl=prm$tcl)`	26	`par(tcl=prm$tcl)`
27	`par(lwd=prm$lwd)`	27	`par(lwd=prm$lwd)`
28	`par(bty=prm$bty)`	28	`par(bty=prm$bty)`
29		29
30		30
31	`##Set margins`	31	`##Set margins`
32	`mar <- prm$mar.nolab`	32	`mar <- prm$mar.nolab`
33	`if( !is.na(prm$xlab) ){ mar[1] <- prm$mar.lab[1] }`	33	`if( !is.na(prm$xlab) ){ mar[1] <- prm$mar.lab[1] }`
34	`if( !is.na(prm$ylab) ){ mar[2] <- prm$mar.lab[2] }`	34	`if( !is.na(prm$ylab) ){ mar[2] <- prm$mar.lab[2] }`
35	`if( any(!is.na(prm$toplab),!is.null(main)) ) { mar[3] <- prm$mar.lab[3] }`	35	`if( any(!is.na(prm$toplab),!is.null(main)) ) { mar[3] <- prm$mar.lab[3] }`
36	`if( !is.na(prm$rightlab) ){ mar[4] <- prm$mar.lab[4] }`	36	`if( !is.na(prm$rightlab) ){ mar[4] <- prm$mar.lab[4] }`
37	`par(mar=mar)`	37	`par(mar=mar)`
38		38
39	`par(mgp=prm$mgp.global)`	39	`par(mgp=prm$mgp.global)`
40		40
41	`## Do the plot`	41	`## Do the plot`
42	`if(!is.null(y)){`	42	`if(!is.null(y)){`
43	`## cat("The Bacher/Delff Plotting System (R), (C), TM operating in xy-plotting mode.\n")`	43	`## cat("The Bacher/Delff Plotting System (R), (C), TM operating in xy-plotting mode.\n")`
44		44
45	`if(is.null(xlim))`	45	`if(is.null(xlim))`
46	`xlim <- range(x,bdp()$xlim,na.rm=TRUE)`	46	`xlim <- range(x,bdp()$xlim,na.rm=TRUE)`
47	`if(is.null(ylim))`	47	`if(is.null(ylim))`
48	`ylim <- range(y,bdp()$ylim,na.rm=TRUE)`	48	`ylim <- range(y,bdp()$ylim,na.rm=TRUE)`
49		49
50	`plot(x, y,`	50	`plot(x, y,`
51	`type=type,`	51	`type=type,`
52	`xaxt=xaxt, yaxt="n",`	52	`xaxt=xaxt, yaxt="n",`
53	`cex=cex.plot,`	53	`cex=cex.plot,`
54	`xlab=xlab,`	54	`xlab=xlab,`
55	`ylab=ylab,`	55	`ylab=ylab,`
56	`xlim=xlim,`	56	`xlim=xlim,`
57	`ylim=ylim,`	57	`ylim=ylim,`
58	`col=col,`	58	`col=col,`
59	`main=main,`	59	`main=main,`
60	`...)`	60	`...)`
61	`} else if (is.numeric(x)){`	61	`} else if (is.numeric(x)){`
62	`if(prm$method=="barplot"){`	62	`if(prm$method=="barplot"){`
63	`## this does not provide the full x,y functionality of barplot`	63	`## this does not provide the full x,y functionality of barplot`
64	`barplot(height=x,`	64	`barplot(height=x,`
65	`## maybe these two are wrong/stupid?`	65	`## maybe these two are wrong/stupid?`
66	`cex.axis=prm$cex.lab,`	66	`cex.axis=prm$cex.lab,`
67	`cex.names=prm$cex.lab,`	67	`cex.names=prm$cex.lab,`
68	`col=bdp()$hcol,`	68	`col=bdp()$hcol,`
69	`border=border,`	69	`border=border,`
70	`xlab=xlab,`	70	`xlab=xlab,`
71	`ylab=ylab,`	71	`ylab=ylab,`
72	`main=main,`	72	`main=main,`
73	`...)`	73	`...)`
74	`prm$draw.xaxis <- FALSE`	74	`prm$draw.xaxis <- FALSE`
75	`prm$draw.yaxis <- FALSE`	75	`prm$draw.yaxis <- FALSE`
76		76
77	`} else {`	77	`} else {`
78		78
79	`if(is.null(xlim))`	79	`if(is.null(xlim))`
80	`xlim <- range(c(1:length(x),bdp()$xlim),na.rm=TRUE)`	80	`xlim <- range(c(1:length(x),bdp()$xlim),na.rm=TRUE)`
81		81
82	`if(is.null(ylim))`	82	`if(is.null(ylim))`
83	`ylim <- range(x,bdp()$ylim,na.rm=TRUE)`	83	`ylim <- range(x,bdp()$ylim,na.rm=TRUE)`
84		84
85	`plot(x,`	85	`plot(x,`
86	`type=type,`	86	`type=type,`
87	`xaxt="n", yaxt="n",`	87	`xaxt="n", yaxt="n",`
88	`cex=cex.plot,`	88	`cex=cex.plot,`
89	`xlab=xlab,`	89	`xlab=xlab,`
90	`ylab=ylab,`	90	`ylab=ylab,`
91	`xlim=xlim,`	91	`xlim=xlim,`
92	`ylim=ylim,`	92	`ylim=ylim,`
93	`col=col,`	93	`col=col,`
94	`main=main,`	94	`main=main,`
95	`...)`	95	`...)`
96	`}`	96	`}`
97	`} else if( class(x) == "acf"){`	97	`} else if( class(x) == "acf"){`
98	`## Not cleaned up/checked`	98	`## Not cleaned up/checked`
99	`plot(x, xlab="", ylab="", xaxt="n", yaxt="n", cex=prm$cex.plot, xlim=prm$xlim, ylim=prm$ylim,col=col,...)`	99	`plot(x, xlab="", ylab="", xaxt="n", yaxt="n", cex=prm$cex.plot, xlim=prm$xlim, ylim=prm$ylim,col=col,...)`
100	`} else if( class(x) == "histogram"){`	100	`} else if( class(x) == "histogram"){`
101	`## Notice: the color of the bars can ONLY be set width bdp(hcol="color").`	101	`## Notice: the color of the bars can ONLY be set width bdp(hcol="color").`
102	`plot(x,`	102	`plot(x,`
103	`xlab=xlab, ylab=ylab,`	103	`xlab=xlab, ylab=ylab,`
104	`main=prm$main,`	104	`main=prm$main,`
105	`xaxt="n",yaxt="n",col=bdp()$hcol,`	105	`xaxt="n",yaxt="n",col=bdp()$hcol,`
106	`border=border,`	106	`border=border,`
107	`...)`	107	`...)`
108	`} else if( class(x) == "trellis"){`	108	`} else if( class(x) == "trellis"){`
109	`## Not cleaned up/checked`	109	`## Not cleaned up/checked`
110	`## very experimental`	110	`## very experimental`
111	`cat("The Bacher/Delff Plotting System (R), (C), TM operating in lattice mode.\n Remeber that labels must written in the lattice object.\n")`	111	`cat("The Bacher/Delff Plotting System (R), (C), TM operating in lattice mode.\n Remeber that labels must written in the lattice object.\n")`
112	`##lattice.options(layout.widths = prm$lattice.width,`	112	`##lattice.options(layout.widths = prm$lattice.width,`
113	`## layout.heights = prm$lattice.height)`	113	`## layout.heights = prm$lattice.height)`
114	`trellis.par.set(prm$myLatticeSettings())`	114	`trellis.par.set(prm$myLatticeSettings())`
115	`plot(x,col=col,...)`	115	`plot(x,col=col,...)`
116	`## with lattice/trellis, the axis drawing doesn't work.`	116	`## with lattice/trellis, the axis drawing doesn't work.`
117	`prm$draw.xaxis <- FALSE`	117	`prm$draw.xaxis <- FALSE`
118	`prm$draw.yaxis <- FALSE`	118	`prm$draw.yaxis <- FALSE`
119	`} else if (class(x)=="princomp") {`	119	`} else if (class(x)=="princomp") {`
120	`plot(x,`	120	`plot(x,`
121	`main=prm$main,`	121	`main=prm$main,`
122	`...)`	122	`...)`
123	`} else if (prm$method=="image.plot"){`	123	`} else if (prm$method=="image.plot"){`
124	`### This could be checked in the beginning by is.null(z)`	124	`### This could be checked in the beginning by is.null(z)`
125	`cat("The Bacher/Delff Plotting System (R), (C), TM operating in image plotting mode.\nx and y axis can not be configured.\n")`	125	`cat("The Bacher/Delff Plotting System (R), (C), TM operating in image plotting mode.\nx and y axis can not be configured.\n")`
126	`image.plot(x=x,y=y,z=z,`	126	`image.plot(x=x,y=y,z=z,`
127	`xlab=xlab,`	127	`xlab=xlab,`
128	`ylab=ylab,`	128	`ylab=ylab,`
129	`...)`	129	`...)`
130	`}`	130	`}`
131		131
-		132	`###Grid stuff`
-		133	`## This really hould be done before adding the rest of the plot contents.`
132	`## We could make default values for grid.v og grid.h. Man maa kunne lave noget kvalificeret ud fra range og noget heltalsdivision`	134	`## We could make default values for grid.v og grid.h. Man maa kunne lave noget kvalificeret ud fra range og noget heltalsdivision`
133		135
134	`##Grid stuff`	-
135	`if( prm$grid ){`	136	`if( prm$grid ){`
136	`grid(col=prm$grid.col)`	137	`grid(col=prm$grid.col)`
137	`}`	138	`}`
138	`## vertical lines`	139	`## vertical lines`
139	`if( !is.na(prm$grid.v) ){`	140	`if( !is.na(prm$grid.v) ){`
140	`if(prm$grid.v=="Def"){`	141	`if(prm$grid.v=="Def"){`
141	`grid(nx=NULL,ny=NA,col=prm$grid.col)`	142	`grid(nx=NULL,ny=NA,col=prm$grid.col)`
142	`} else`	143	`} else`
143	`{`	144	`{`
144	`abline(v=prm$grid.v, lty="dotted", col=prm$grid.col)`	145	`abline(v=prm$grid.v, lty="dotted", col=prm$grid.col)`
145	`}`	146	`}`
146	`}`	147	`}`
147		148
148	`if( !is.na(prm$grid.h) ){`	149	`if( !is.na(prm$grid.h) ){`
149	`if(prm$grid.h=="Def"){`	150	`if(prm$grid.h=="Def"){`
150	`grid(nx=NA,ny=NULL,col=prm$grid.col)`	151	`grid(nx=NA,ny=NULL,col=prm$grid.col)`
151	`}else{`	152	`}else{`
152	`abline(h=prm$grid.h, lty="dotted", col=prm$grid.col)}`	153	`abline(h=prm$grid.h, lty="dotted", col=prm$grid.col)}`
153	`}`	154	`}`
154		155
-		156	`### Axis stuff This part should be improved. A function that draws`
-		157	`### axis should be run for all four axis. And it should be possible to`
-		158	`### supply whatever vector to base it on (especially relevant for`
155	`## Axis stuff`	159	`### taxis and raxis)`
-		160
-		161	`## this is because the use of mgp gives a ridiculous warning when too small`
-		162	`options(warn=-1)`
156	`if(prm$draw.xaxis){`	163	`if(prm$draw.xaxis){`
157	`## mgp.old <- par()$mgp`	164	`## mgp.old <- par()$mgp`
158	`## par(mgp) <- mgp.xaxis`	165	`## par(mgp) <- mgp.xaxis`
159	`options(warn=-1)`	166	`## is the abscissa a time object?`
160	`if( "POSIXt"%in%class(x[1])){`	167	`if( "POSIXt"%in%class(x[1])){`
161	`axis.POSIXct(1, x,mgp=prm$mgp.xaxis, lwd=prm$lwd)`	168	`axis.POSIXct(1, x,mgp=prm$mgp.xaxis, lwd=prm$lwd)`
162	`} else {`	169	`} else {`
163	`axis(1, mgp=prm$mgp.xaxis, lwd=prm$lwd)`	170	`axis(1, mgp=prm$mgp.xaxis, lwd=prm$lwd)`
164	`}`	171	`}`
165	`options(warn=0)`	-
166		172
167	`}`	173	`}`
-		174
168	`if(prm$draw.yaxis){ axis(2, mgp=prm$mgp.yaxis, lwd=prm$lwd) }`	175	`if(prm$draw.yaxis){ axis(2, y,mgp=prm$mgp.yaxis, lwd=prm$lwd) }`
-		176	`## top axis`
-		177	`if(prm$draw.taxis){ axis(3, mgp=prm$mgp.taxis, lwd=prm$lwd) }`
169	`## axis to the right`	178	`## axis to the right`
170	`if(prm$draw.raxis){ axis(4, mgp=prm$mgp.raxis, lwd=prm$lwd) }`	179	`if(prm$draw.raxis){ axis(4, mgp=prm$mgp.raxis, lwd=prm$lwd) }`
-		180
-		181	`## switch back on warnings`
-		182	`options(warn=0)`
-		183
-		184
171	`##Title stuff`	185	`##Title stuff`
172	`## Hvorfor??`	186	`## Hvorfor??`
173	`scale <- 1`	187	`scale <- 1`
174	`## if( prm$type=="hist" & !is.na(prm$xlab) )`	188	`## if( prm$type=="hist" & !is.na(prm$xlab) )`
175	`# if( !is.na(prm$xlab) ){ mtext(prm$xlab, line=prm$xlabLine, side=1, cex=prm$cex.lab/scale) }`	189	`# if( !is.na(prm$xlab) ){ mtext(prm$xlab, line=prm$xlabLine, side=1, cex=prm$cex.lab/scale) }`
176	`## if( prm$type=="hist" &!is.na(prm$ylab) )`	190	`## if( prm$type=="hist" &!is.na(prm$ylab) )`
177	`# if( !is.na(prm$ylab) ){ mtext(prm$ylab, line=prm$ylabLine, side=2, cex=prm$cex.lab/scale) }`	191	`# if( !is.na(prm$ylab) ){ mtext(prm$ylab, line=prm$ylabLine, side=2, cex=prm$cex.lab/scale) }`
178	`if( !is.na(prm$toplab) ){ mtext(prm$toplab, side=3, line=0.25, cex=prm$cex.lab/scale) }`	192	`if( !is.na(prm$toplab) ){ mtext(prm$toplab, side=3, line=0.25, cex=prm$cex.lab/scale) }`
179	`if( !is.na(prm$rightlab) ){`	193	`if( !is.na(prm$rightlab) ){`
180	`mtext(prm$rightlab, side=4, line=0.75, cex=prm$cex.lab/scale,mgp=prm$mgp.raxis)`	194	`mtext(prm$rightlab, side=4, line=0.75, cex=prm$cex.lab/scale,mgp=prm$mgp.raxis)`
181	`## mtext(rightlab, side=4, line=0.5, cex=bdp()$cex.lab,mgp=bdp()$mgp.raxis)`	195	`## mtext(rightlab, side=4, line=0.5, cex=bdp()$cex.lab,mgp=bdp()$mgp.raxis)`
182	`}`	196	`}`
183	`}`	197	`}`
184		198
185		199
186		200
187		201

Subversion Repositories bdplot

(root)/bdgraphics/R/4_bdplot.R – Rev 29 → 31