similar to: Residuals with Elliot-Rothenberg-Stock Unit Root Test

Dear all, Does anyone know if the PLM package (to run Panel Data Analysis) accepts quarterly data? The package vignette and documentation only use annual data -and the only time index available seems to work for years. José Mr José Luis Iparraguirre Senior Research Economist Economic Research Institute of Northern Ireland 2 -14 East Bridge Street Belfast BT1 3NQ Northern

Does anyone know how to read SPSS .por (ie Portable) files? The foreign package only deals with SPSS .sav files and not with those with a "por" extension. Thanks, José Mr José Luis Iparraguirre D'Elia Senior Research Economist Economic Research Institute of Northern Ireland Floral Buildings 2 - 14 East Bridge Street Belfast BT1 3NQ Northern Ireland United Kingdom

I am using the function prop.test (base package), which returns a list with class "htest". All I want to do is to assign one of its values to a variable, but I do not want R to print the results and added warning message whenever I invoke the function. How can I prevent R from printing on using prop.test? Regards, Jos? Luis Mr Jos? Luis Iparraguirre Senior Research Economist Economic

I could not find a package to run panel data unit root tests in R (there's a STATA module, though - PANELUNIT- and routines for Sarno and Taylor's MADF test, Levin-Lin-Chu's test, Im-Pesaran-Shin's Test, etc, to do this). Hence, I am toying with the idea of having a go at writing up one for R, but would like to know whether anyone is already working on this. José Mr

Dear R users, Consider the first two columns of a data frame like this: > z[,1:2] x y 1 1 1 2 2 2 3 3 3 4 1 4 Imagine that y represents the times that the value x happens in a population. But z is not exactly a frequency table, because in z we have x=1 twice. So, the x=1 in the first line and the x=1 in the fourth are not the same, differing according to a third variable in

Dear all, Say I have the following dataset: > DF x y z [1] 1 1 1 [2] 2 2 2 [3] 3 3 NA [4] 4 NA 4 [5] NA 5 5 And I want to omit all the rows which have NA, but only in columns X and Y, so that I get: x y z 1 1 1 2 2 2 3 3 NA If I use na.omit(DF), I would delete the row for which z=NA, obtaining thus x y z 1 1 1 2 2 2

I am having difficulty installing the package "lattice". First I tried downloading it from the CRAN site (in the normal way) : the message said it worked but when I typed library I got an error message ("there is no package"). Second I tried installing it from the local zip. I have reproduced the result below. utils:::menuInstallLocal() package 'lattice'

Hi R-users,   I have this code here: library(numDeriv)   fprime <- function(z) { alp  <- 2.0165;   rho  <- 0.868;   # simplified expressions   a      <- alp-0.5   c1     <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)   c2     <- sqrt(rho)/(1-rho)   t1     <- exp(-z/(1-rho))   t2     <- (z/(2*c2))^a   bes1   <- besselI(z*c2,a)   t1bes1 <- t1*bes1   c1*t1bes1*t2 }   ## Newton

Hi, I have a Data Set x and I want to check with a Kolmogorow-Smirnow-Test, if x comes from a Subbotin Distribution, whose density function is: function(x,location,scale,tail) # Exponential power (=Subbotin) { const<- 2*scale*tail^(1/tail) *gamma(1+1/tail) z<- (x-location)/scale exp(-1/tail*abs(z)^tail)/const } How can i do this? Thank you --

R-helpers I have 120 small Excel sheets to read and I am using library(xlsReadWrite): one example below. I had hoped to read sheets by looping over a list of numbers in their name (eg Book1.xls, Book2.xls, etc). I thought I had seen examples which used eval-parse-paste in this way. However, I have not been able to get it to work.. 1. is this a feasible approach? 2. if not

Dear All, It's been several times I observed a very strange behaviour with R 1.2.0 under Windows NT 4.0 (with sp5): accuracies are lost in many cases. Things are back to normal when I exit R and restart it. Here is an example of my last encounter with this: > bt1 <- c(7.3, 8.1, 8.4, 7, 6, 7.2, 6.7, 8.3, 9, 7.4, 8.6, 6.2, 6.9, 7.6, 5.3, 5.6) > bt1 [1] 7 8 8 7 6 7 6 8 9 7 8 6 6 7 5 5

Hi r-users,   I hope somebody can help me to understand the error message.  Here is my code; ## Newton iteration newton_gam <- function(z) { n   <- length(z)   r   <- runif(n)   tol <- 1E-6   cdf <- vector(length=n, mode="numeric")   fprime <- vector(length=n, mode="numeric")   f   <- vector(length=n, mode="numeric")     for (i in 1:1000)   {

Saben si el estimador VARHAC de Den Haan y Levin (1994) ha sido implementado en R? Solamente encontre rutinas para RATS y GAUSS en la pagina de Den Haan. Saludos, Jose Jose Luis Iparraguirre [[alternative HTML version deleted]]

Hi I'm trying to measure the effect of a policy intervention (Box and Tiao, 1975). This query has to do with the coding of the model rather than with the particulars of my dataset, so I'm not providing the actual dataset (or a simulated one) in this case, apart from some general description. The time series are of length n=34 (annual observations between 1977 and 2010). The policy

I'm analizing the Argentina stock market (merv) I download the data from yahoo library(tseries) Argentina <- get.hist.quote(instrument="^MERV","1996-10-08","2003-11-03", quote="Close") merv <- na.remove(log(Argentina)) I made the Augmented Dickey-Fuller test to analyse if merv have unit root: adf.test(merv,k=13) Dickey-Fuller = -1.4645,

Hi, Imagine I have the following data frame: > a <- c(1,NA,3) > b <- c(2,NA,NA) > c <- data.frame(cbind(a,b)) > c a b 1 1 2 2 NA NA 3 3 NA I want to delete the second row. If I use na.omit, that would also affect the third row. I tried to use a loop and an ifelse clause with is.na to get R identify that row in which all records are missing, as opposed to the first

Hola a todos. Estoy atascado en un sitio que no me he encontrado antes, debe ser fácil y rápido, a ver si alguien me puede ayudar. Puesto fácil, tengo: > aa <- matrix(11:19,3,3) > aa [,1] [,2] [,3] [1,] 11 14 17 [2,] 12 15 18 [3,] 13 16 19 Quiero la fila cuyo primer elemento es 12 y el segundo 15. > aa[aa[,1:2] == (c(12,15)),] Error: (subscript) logical

I'm analizing the Argentina stock market (merv) I download the data from yahoo library(tseries) Argentina <- get.hist.quote(instrument="^MERV","1996-10-08","2003-11-03", quote="Close") merv <- na.remove(log(Argentina)) I made the Augmented Dickey-Fuller test to analyse if merv have unit root: adf.test(merv,k=13) Dickey-Fuller = -1.4645,

Dear R-list member, an update of package 'urca' has been uploaded to CRAN (Mirror: Austria). In the updated release unit root and cointegration tests encountered in applied econometric analysis are implemented. The package is written in 'pure' R and utilises S4 classes. In particular, the Johansen procedure with likelihood ratio tests for the inclusion of a linear trend,

Hello all, I tried to extract the p-value from adf.test in tseries; however, I got the error message such as > ht=adf.test(list.var$aa) > ht$p-value Error in ht$p - value : non-numeric argument to binary operator > ht Augmented Dickey-Fuller Test data: list.var$aa Dickey-Fuller = -2.3147, Lag order = 4, p-value = 0.4461 alternative hypothesis: stationary > ht$data [1]