1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
# set up plot to disk as per my special requirements# pointsize: the default pointsize of plotted text. The allowed range is# '[6, 48]': values outside that range are reset to 12.# this works for word (emf is better than wmf)# ..later.. I like to reduce the size to 80% in MS word - it makes the axis# lables and other writing a bit smaller (I spuppose I could use cex=.8 or somthing)myemf <- function(filename, wdt=15, hgt = 12, pdir=plot_dir, has_title = F){ if(!exists('pdir')) plot_dir <- plotDir() filename <- paste(pdir, '/', filename, '.emf', sep='') win.metafile(file=filename, width=wdt/2.54, height=hgt/2.54) cat('Plotting to', filename, '\n') my_mai <- c(.925,.795,.05,.2) # c(bottom, left, top, right) if(has_title) my_mai[3] <- my_mai[3] + .15 # add space for main par(mai = my_mai)}###############################From: McGehee, Robert <Robert.McGehee_at_geodecapital.com>#Date: Sat 14 Jan 2006 - 01:27:08 EST# creates plot axis labels in scientific notationsciNotation <- function(x, digits = 1) { if (length(x) > 1) { return(append(sciNotation(x[1]), sciNotation(x[-1]))) } if (!x) return(0) exponent <- floor(log10(x)) base <- round(x/10^exponent,digits) as.expression(substitute(base %*% 10^exponent, list(base = base, exponent = exponent)))}## these functions are called from my_emf######################################################################## yymmdd date formatyymmdd <- function(thisDate = Sys.Date()){ return(format(thisDate, '%y%m%d'))}######################################################################## plotDir <- function(pdir = 'plot', ...){ plot_dir <- paste(pdir, '/', yymmdd(...), sep='') dir.create(plot_dir, rec=T) return(plot_dir)}
|
