similar to: Manual two-way demeaning of unbalanced panel data (Wansbeek/Kapteyn transformation)

I have an unexpected result with the functions lag() and diff() in the plm (panel data) package when used with transform(). These plm-specific functions are supposed to generate lags and first differences within each panel. lag() does not work properly the first time (it reproduces the same series--this is a common time series pitfall), BUT then it does work properly when it is run a second

I didn't see you got an answer posted to this question: You can't modify a pdata.frame object. Your transforms turn it back to a normal data frame and diff and lag won't work as expected. Try: Grunfeld.p <- pdata.frame(Grunfeld,c("firm","year")) tmp <- transform(Grunfeld.p, d.value = diff(Grunfeld.p$value,1)) tmp <- cbind(tmp, l.value =

Dear List, I am using the output of a ML estimation on a random effects model with first-order autocorrelation to make a further conditional test. My model is much like this (which reproduces the method on the famous Grunfeld data, for the econometricians out there it is Table 5.2 in Baltagi): library(Ecdat) library(nlme) data(Grunfeld)

Dear R users What would be the best way to approach estimating a panel regression with twoways random effects, on unbalanced data? Unfortunately, the "plm" package has no implementation of twoways random effects for unbalanced data. Currently I'm considering two approaches: - extend "plm" to cover this type of panel regression. (For the authors, cc'ed:) Would

Hi list, has anyone succeeded in using pgmm() on any dataset besides Arellano/Bond's EmplUK, as shown in the vignette? Whatever I try, I eventually get a runtime error because of a singular matrix at various points in pgmm.diff() (which gets called by pgmm()). For example, when estimating a "dynamic" version of the Grunfeld data: data(Grunfeld, package="Ecdat") grun

Taking the Grunfeld data, which is built-in in R, for example, (1)How can I construct a dataset (or dataframe) that consists of the data of all firms in 1951? (2)How can I calculate the average capital in each form over the period 1951-1954? What I can imagine is to categorize the data by firm, and then select the data between 1951 and 1954 for each firm, but how can I do it? Thanks, Miao

Dear R users I get different F-statistic results for a "within" model, when using "time" or "twoways" effects in plm() [1] and when manually specifying the time control dummies [2]. [1] vignette("plm") [2] http://cran.r-project.org/doc/contrib/Farnsworth-EconometricsInR.pdf Two examples below: library("AER") data("Grunfeld", package =

Hello Cecilia, nice hearing from you again. I must restate a couple of my old hints, though ;^) 1) please always put the authors c/c, as we are not guaranteed to browse through the r-help every day 2) please provide reproducible examples. As example(pooltest) keeps working fine, as do some other cases I tried (Grunfeld data etc.), I don't know what the problem is but evidently your data are

Dear R users, I have an unbalanced panel with an average of I=100 individuals and a total of T=1370 time intervals, i.e. T>>I. So far, I have been using the plm package. I wish to estimate a FE model like: res<-plm(x~c+v, data=pdata_frame, effect="twoways", model="within", na.action=na.omit) ?where c varies over i and t, and v represents an exogenous impact on x

Dear Ivo, dear list, (see: Message: 70 Date: Thu, 26 Mar 2009 21:39:19 +0000 From: ivowel at gmail.com Subject: [R] pgmm (Blundell-Bond) sample needed) I think I finally figured out how to replicate your supersimple GMM example with pgmm() so as to get the very same results as Stata. Having no other regressors in the formula initially drove me crazy. This was a case where simpler models are

I am working with panel data. I am using the plm package to do this. I would like to do be able to adjust for autocorrelation, as one does with glm models and correlation structures (eg corr=corARMA(q=4)) . In particular, I want to employ MA(4) error structure. Is there a way of doing this with the plm package? (Note: I do not really want to use the pggls function for various

Hi everyone and thank you for the help you could give me. My data is in a spreadsheet. The 1st column identifies the firm (with the fiscal number), the columns 2 to 11 have the variable value for 11 years. I have many variables (files like this). Each file has about 40.000 firms (rows). I transformed all the files in txt files. The data is a panel data, like this: firm revenu2007 revenue2006

Dear all I am trying to fit Variable Coefficients Models on Unbalanced Panel Data. I managed to fit such models on balanced panel data (the example from the "plm" vignette), but I failed to do so on my real, unbalanced panel data. I can reproduce the error on a modified example from the vignette: > require(plm) > data("Hedonic") > Hed <- pvcm(mv ~ crim + zn + indus

Hi, I maintain postfix mailserver with sudo enabled dovecot LDA using multiple uids. I want to enforce quotas per-unix-user (which can have multiple domains with multiple emails set up using postfixadmin), so I have set up a FS quota backend with imap_quota plugin. Unfortunately, the quota appears to work (when i disable chrooting) according to dovecot-info.log, but doesn't send anything over

Hello, I understand that the/ demean/ argument in the *ar()* function to fit an autoregressive model selects the best AR model fitted to the mean deleted observations. What is the purpose of using this demean procedure at all? Its seems silly as the post doesn't deal with R problems.... Thanks -- View this message in context:

Suppose I have longitudinal data and want to use the econometric strategy of "de-meaning" a model matrix by time. For sake of illustration 'mat' is a model matrix for 3 individuals each with 3 observations where ``1'' denotes that individual i was in group j at time t or ``0'' otherwise. mat <- matrix(c(1,1,0,0,0,0,0,0,1,0,0,0,1,1,1,0,0,0,0,0,1,0,0,0,1,1,0),

Interesting since imap_acl also makes dovecot segfault. Maybe there is a problem in the way Dovecot loads plugins? Cheers, Olivier On 16/01/2011 19:09, dovecot-request at dovecot.org wrote: > Send dovecot mailing list submissions to > dovecot at dovecot.org > > Message: 1 > Date: Sun, 16 Jan 2011 12:50:17 +0100 > From: "G.Nau" <b404_r66 at yahoo.de> >

Dear sir,    How to do bilinear time series in R?Is there any functions or packages?  thank you! -----Sincerely yours Kuangnan Fang 方匡南 敬上 department of statistics ,Economics school,Xia men University. Fujian Province (361005) China Mobile Phone：15860721915 SKYPE: ruiqwy MSN Messenger: ruiqwy@hotmail.com QQ:39863401 --- 09年3月31日，周二, r-help-request@r-project.org