similar to: How to put title Vertically

Hi Dieter and R community: I tried both of these three versions with ylim as suggested, none work: I am getting only single (pch = 16) not overlayed (pch =3) everytime. *vs 1* require(lattice) xyplot(Sepal.Length ~ Sepal.Width | Species , data= iris, panel= function(x, y, subscripts) { panel.xyplot(x, y, pch=16, col = "green4", ylim = c(0, 10)) panel.lmline(x, y, lty=4, col =

Dear list, I have done a scatter plot of multiple variables in the same graph, with different col and pch. I managed to do it with the following code but not know how to make a function of these so that next time if I want to do similar graph but with new variables, I dont have to copy the code and then change the old variables with the new ones but just call a function with the new

hi, in 1.9.1, the return value from ARMAacf(pacf=TRUE) is not named by lags, contrary to ?ARMAacf. the simple fix is to move names(Acf) <- down after if(pacf), with an appropriate starting lag as pacf=TRUE appears to start at lag 1 (whereas pacf=FALSE starts at lag 0). for consistency, one could argue to append 1 for lag 0 for pacf=TRUE (or start pacf=F at lag 1). however, given the

How can I make sure the residual signal, after subtracting the trend extracted through some technique, is actually trend-free ? I would greatly appreciate any suggestion about some Stationarity tests. I'd like to make sure I have got the difference between ACF and PACF right. In the following I am citing some definitions. I would appreciate your thoughts. ACF(k) estimates the correlation

I am getting the following: *** Weave Errors *** Error in driver$runcode(drobj, chunk, chunkopts) : Error in eval(expr, envir, enclos) : couldn't find function "pacf.mts" *** Source Errors *** Error in eval(expr, envir, enclos) : couldn't find function "pacf.mts" make[1]: *** [checkVignettes] Error 1 I don't really understand the new namespace mechanism,

pacf in devel seems by default to return a different number of lags than 1.7.1 for $pacf. I don't see any mention of this in the NEWS file, or any change in the documentation, so I suspect it is and error, though it may be an undocumented improvement. (Newbie question: How is the simplest way to display a function like pacf.default that is not exported from a namespace?) Paul

Dear list, according to the documentation of acf{ts} "the partial correlation coefficient is estimated by fitting autoregressive models of successively higher orders up to lag.max. " However, R seems to return the Yule-Walker estimates of the PACF by default. You can check this using c(1:10) as the series: the YW estimates are 0.7000000 and -0.1527035 for lags 1 and 2 . If the PACF

Hi! I want to plot three graphs (residuals, ACF and PACF of a model). Ideally I would use a c(2,2) disposition where the residuals plot would start at position 1,1 and span to position 1,2. Then I would plot the ACF in position 2,1 and the PACF in position 2,2. Maybe is clearer like this: -------------------------- | | | residuals | |

Dear all, I found the following behaviour with pacf() in the multivariate case, set.seed(10) x <- rnorm(1000,sd=10000) y <- rnorm(1000,sd=1) pacf(ts(cbind(x,y)),plot=FALSE,lag.max=10) Partial autocorrelations of series 'cbind(x, y)', by lag , , x x y 0.047 ( 1) 0.000 ( -1) 0.011 ( 2) 0.000 ( -2) 0.005 ( 3) 0.000 ( -3) 0.013 ( 4)

Dear list, Is there a way to reorder the xaxis using lattice. Using the following data, the x-axis is ordered as BP GH MH PF RE RP SF VT but I would like the x-axis to be ordered as PF RP BP GH VT SF RE MH. Kön Skalor Tillfälle Medelvärde 1 Kvinnor BP 1-inskrivning 36.45283 2 Kvinnor GH 1-inskrivning 38.62255 3 Kvinnor MH 1-inskrivning 62.88889 4

When I use the command for PACF, lags 5, 10, 15, and 20 are labeled. I would like to label lag 1. I would greatly appreciate if someone could tell me how to do this. Below is the command that I am using: pacf(data$R1,main="Series R1 Residuals") [[alternative HTML version deleted]]

R: Is there an easy way to get the acf and pacf for an irregular times series? That is, the acf and pacf with lag lengths that are in units of time, not observation number. Thanks, Jason Higbee Research Associate Federal Reserve Bank of St. Louis The views expressed in this email are the author's and not necessarily those of the Federal Reserve Bank of St. Louis or the Federal Reserve

Hi Folks, Here is my attempt at a simple polar plot. > pos <- seq(0,360,by=5) > tspk <- rep(c(1,0,1),c(13,47,13)) > require(plotrix) > polar.plot(tspk,pos,rp.type="s",point.symbols=17,point.col="green4") I only see half the symbols, the other half of each symbol is hidden under the circular grid. In fact if I change rp.type="r", I see the

Dear experts: I want to overlay lattice scatter plot: I do not know why the following code is not plotting subscripts ! Sorry if this question is too simple: require(lattice) xyplot(Sepal.Length ~ Sepal.Width | Species , data= iris, panel= function(x, y, subscripts) { panel.xyplot(x, y, pch=16, col = "green4", ylim = c(0, 10)) panel.lmline(x, y, lty=4, col = "green4")

Hi, I'm slowly working through Tsay's "Analysis of Financial Time Series" 3rd ed. ?I'm trying to replicate Table 2.1 on p.47, which gives PACF, AIC, and BIC for the monthly simple returns of the CRSP value-weighted index. The data: http://faculty.chicagobooth.edu/ruey.tsay/teaching/fts3/m-ibm3dx2608.txt > da <-

Hello, r-help at r-project.orgbarplot(twcons.area, beside=T, col=c("green4", "blue", "red3", "gray"), xlab="estate", ylab="number of persons", ylim=c(0, 110), legend.text=c("treated", "mix", "untreated", "NA")) produces a barplot very fine. In addition, I'd like to get the

Dear list, following a previous querry we are still stuck! As pointed out by Erin Hodges the "ts" library includes a PACF function which reports the partial correlation of population density at time t against lagged population density. However, what we are trying to calculate is the partial correlation between rate of population change, Rt=log Nt/Nt-1, against lagged population

Estimados amigos Estoy dibujando las funciones acf y pacf de una variable de una serie "zoo": > ls.str(pat="T0.5") T0.5 : 'zoo' series from 2017-11-08 23:00:00 to 2017-11-15 06:59:00   Data: num [1:9120, 1:3] 55 49.8 51 50.1 36.5 ...   Index:  POSIXct[1:9120], format: "2017-11-08 23:00:00" "2017-11-08 23:01:00" "2017-11-08

I have used acf() and pacf() in R to get the acf and pacf values at max/lag=20 but the output did not show the values associated with lag numbers. lag numbers is shown in decimals. -- Rashid Ameer View my recent publication at * http://www.emeraldinsight.com/fwd.htm?id=aob&ini=aob&doi=10.1108/17538391211282854 * Details for my works are available directly at

I'd like to stretch a plotted character vertically, to create a "sequence logo". Is there a parameter to allow stretching text() output vertically or squeeze horizontally? I know about Oliver Bembom's seqLogo library, but this generates a sequence logo plot using a separate bitmap device. I want to recreate the sequence logo *inside* an existing plot. Alternatively, is