plot with the axes at xlim and ylim

I'd like to make a plot with the axes drawn at the values of xlim and
ylim. The default plot draws the axes slightly outside these values. I
have been experimenting with the par(mgp) setting, but specifying
par(mgp=c(2,0,-1)) gives me

"invalid value specified for graphics parameter "mgp".".

A more complicated mathod that almost seems to work is doing the plot
with "axes=F" and then adding the axes with "axis" statements:

xlim <- c(0,10)
ylim <- c(0,11)
plot(c(1:9),c(1:9),axes=F,xlim=xlim,ylim=ylim)
axis(1,pos=ylim[1])
axis(2,pos=xlim[1])
axis(3,pos=ylim[2],labels=F)
axis(4,pos=xlim[2],labels=F)

The problem here is that the vertical axes are not long enough.

Is there an easier way to get what I want?

I'd like to make a plot with the axes drawn at the values of xlim and
ylim. The default plot draws the axes slightly outside these values. I
have been experimenting with the par(mgp) setting, but specifying
par(mgp=c(2,0,-1)) gives me

"invalid value specified for graphics parameter "mgp".".

A more complicated mathod that almost seems to work is doing the plot
with "axes=F" and then adding the axes with "axis" statements:

xlim <- c(0,10)
ylim <- c(0,11)
plot(c(1:9),c(1:9),axes=F,xlim=xlim,ylim=ylim)
axis(1,pos=ylim[1])
axis(2,pos=xlim[1])
axis(3,pos=ylim[2],labels=F)
axis(4,pos=xlim[2],labels=F)

The problem here is that the vertical axes are not long enough.

Is there an easier way to get what I want?

Thanks,
Arnold Kester

--
Met vriendelijke groet,
Arnold Kester

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.

r-help mailing list -- Read

Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch

On 19 Mar 2002, Arnold Kester wrote:

I'd like to make a plot with the axes drawn at the values of xlim and
ylim. The default plot draws the axes slightly outside these values. I
have been experimenting with the par(mgp) setting, but specifying

Is there an easier way to get what I want?

par(xaxs="i", yaxs="i")


-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

i am using R 1.5 on solaris 2.8 and redhat 7.2 (what i mention
below happens for all of these and on R-1.4.1 on solaris 2.8)

i am fitting an addtive model (y=roweffect+columneffct+error) to the
following data using rlm

   8.631371  8.219719
  12.407817 12.373610
  10.558580 10.284926
  11.525872 11.504441
  12.644998 12.583985
  12.772225 12.629980
  12.482169 12.271868
  12.039199 12.061671
  11.892378 11.688070
 11.725565 11.655078
 11.251226 11.383293
  9.199893  7.860003
 12.275562 12.240147
 13.638981 13.381291
 13.205941 13.060880
 11.056920 10.707565

if the above matrix is x, when i do the following:
  m <- dim(x)[1]
  n <- dim(x)[2]
  rows<- as.factor(rep(1:m,n))
  cols <- as.factor(rep(1:n,rep(m,n)))

  z  <-rlm(as.vector(x)~cols+rows,psi=psi.bisquare)

i get an NA in the estimate z$coef[13]

this relates to the 12th row  effect which appears to be an
outlier. is this how one interprets an NA  here?

one consequence is that i cant get standard errors
through summary.rlm(z). is there a "work around" this?