[R] Query: how to modify the plot of acf
Matthias Braeunig
mb.atelier at web.de
Fri Aug 18 20:56:36 CEST 2006
Hi Stefano,
the manual tells us that we can access components of an acf object
directly by acf.obj[.], but assignment ]<- does not work this way.
One way of doing what you want is to assign NA to x$acf[x$lag==0] like so:
x <- acf(runif(100))
x$acf[1] <- NA
plot(x)
But I suppose what you actually want is to have a reasonable ylim to be
in effect when plotting your acf. I have struggled with this myself and
found the following solution:
In plot.acf the formula for the white noise confidence interval is
something along the lines of
acf.clim <- function(x=stop("your acf"), n=x$n.used, CL=.95)
{
cl = qnorm((1+CL)/2)/sqrt(n)
attr(cl,"CL") <- CL
return( cl)
}
Using this function you can plot
x <- acf(runif(100),plot=F)
plot(x, ylim=2*acf.clim(x)*c(-1,1)) # adjust magnification as needed
As for your second question we first prevent plotting of x-axis and then
add our custom axis:
x <- acf(runif(1000),100,plot=F)
plot(x, ylim=2*acf.clim(x)*c(-1,1),xaxt='n') # note option 'xaxt'!
axis(1,at=12*(0:8))
Hope that helps,
Matthias
NB: It would be a good idea to include the acf.clim(.) inside of the acf
object... Now it is somewhat hidden inside of plot.acf
x$clim <- acf.clim(x)
str(x)
Stefano Sofia wrote:
> I need to modify the graph of the autocorrelation. I tried to do it through plot.acf but with no success.
>
> 1. I would like to get rid of the lag zero
> 2. I would like to have numbers on the x-axis only at lags 12, 24, 36, 48, 60, ...
>
> Could anybody help me in this?
>
> Any help will be appreciated
> Thank you for your attention
> Stefano
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
>
--
Journalist: "Mr. Gandhi, what do you think of Western civilization?"
Gandhi: "I think it would be a good idea."
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
~ Matthias Braeunig, Dipl.Phys.
~ University Medical Center Freiburg
~ Institute of Environmental Medicine and Hospital Epidemiology
~ Department of Evaluation Research in Complementary Medicine
~ Hugstetter Str. 55, D - 79106 Freiburg, Germany
~ Phone: +49 (0)761 270-8306, Fax: -8343
~ Web: http://kompmed.uniklinik-freiburg.de
~ Email: matthias.braeunig at uniklinik-freiburg.de
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
More information about the R-help
mailing list