Dear friends.
As promised ver 1.01 solved the former problem and now dates plot OK in
windows. But.
In the following code, "cex.axis"=0.5 works OK but
"xaxp"=c(14620,14680,4)
has seemingly no effect. The limits for xaxp was found from the default
plot and inspecting par(). I need to put the date marks at specific points.
Could I suppress them altogether and make a new xaxis corresponding to the
points ? Was xaxp not intended in this way ?
Thank you for your (so far) endless patience !
data <- read.table("vanco1.dat",h=T)
op <- par("xaxp"=c(14620,14680,4),"cex.axis"=0.5)
attach(data)
library(date)
dates <- mdy.date(month,day,year)
plot(dates,vanco,pch=c(19),xlab='',ylab="[Vancomycin]",
main="Vanco versus weekly dose")
abline(h=10,lty=2)
dm <- c(1,1,1,1,2,2,3)
dd <- c(7,17,21,31,23,25,3)
dy <- rep(2000,7)
limits <- mdy.date(dm,dd,dy)
abline(v=limits,lty=8)
x <- c(14611.37, 14620.95, 14627.94, 14634.44, 14650.85, 14663.85,
14667.58,14680.15)
mtext(c("1","2","3","4","4.5","3","2","3"),3,at=x)
par(op)
-------------- next part --------------
month day year dosis vanco
12 31 1999 1 3.8
1 3 2000 1 2.4
1 7 2000 1 3.4
1 12 2000 2 7.1
1 14 2000 2 5.3
1 17 2000 2 5.7
1 19 2000 3 7.7
1 21 2000 3 8.3
1 24 2000 4 8.7
1 26 2000 4 11.1
1 31 2000 4 12.3
2 7 2000 4.5 12.0
2 14 2000 4.5 14.9
2 18 2000 4.5 23.0
2 25 2000 3 20.1
2 28 2000 2 11.3
3 1 2000 2 13.7
3 8 2000 3 12.8
3 10 2000 3 7.4
3 13 2000 3 9.1
3 15 2000 3 12.2
3 17 2000 3 13.8
-------------- next part --------------
Troels Ring, MD
Department of Nephrology
Aalborg Hospital, Denmark
tring at mail1.stofanet.dk
> As promised ver 1.01 solved the former problem and now dates plot OK in > windows. But. > In the following code, "cex.axis"=0.5 works OK but "xaxp"=c(14620,14680,4) > has seemingly no effect. The limits for xaxp was found from the default > plot and inspecting par(). I need to put the date marks at specific points. > Could I suppress them altogether and make a new xaxis corresponding to the > points ? Was xaxp not intended in this way ? > Thank you for your (so far) endless patience ! > > data <- read.table("vanco1.dat",h=T) > op <- par("xaxp"=c(14620,14680,4),"cex.axis"=0.5) > attach(data) > library(date) > dates <- mdy.date(month,day,year) > plot(dates,vanco,pch=c(19),xlab='',ylab="[Vancomycin]", > main="Vanco versus weekly dose") > abline(h=10,lty=2) > dm <- c(1,1,1,1,2,2,3) > dd <- c(7,17,21,31,23,25,3) > dy <- rep(2000,7) > limits <- mdy.date(dm,dd,dy) > abline(v=limits,lty=8) > x <- c(14611.37, 14620.95, 14627.94, 14634.44, 14650.85, 14663.85, > 14667.58,14680.15) > mtext(c("1","2","3","4","4.5","3","2","3"),3,at=x) > par(op)I wanted to tell you to use something like plot(dates, vanco, axes=FALSE) axis(.....) box() as a workaround, but plot(dates,vanco,axes=FALSE) draws the "dates"-axis! Bug??? Something like plot(dates, vanco, xaxt="n") axis(1, temp <- seq(14620, 14680, 4), as.character(as.date(temp)), cex.axis=0.5) should work. Regards, Uwe Ligges -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._