similar to: trouble switching to 'plm' from 'xtabond' and Stata

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

Hello list, I'm following the paper (http://www.jstatsoft.org/v27/i02/paper) on how to use "plm" to run panel regressions, and am having trouble with what I believe should be something very basic. When I run the command (p.9 in the paper): R> dynformula(emp~wage+capital,log=list(capital=FALSE,TRUE),lag=list(emp=2,c(2,3)),diff=list(FALSE,capital=TRUE)) I see: emp ~ wage +

Hello list, I'm following the paper (http://www.jstatsoft.org/v27/i02/paper) on how to use "plm" to run panel regressions, and am having trouble with what I believe should be something very basic. When I run the command (p.9 in the paper): R> dynformula(emp~wage+capital,log=list(capital=FALSE,TRUE),lag=list(emp=2,c(2,3)),diff=list(FALSE,capital=TRUE)) I see: emp ~ wage +

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

Hello All! I calculate a variogram using the function variog (package geoR) afterwards I use variofit to fit a spherical model (see code below). Now I just changed the units of the variable (in this case MPa to kPa just a factor of 1000). If I do so, I get a different fit and therefore different ranges etc. Why? The semi-variance is of course 6 orders of magnitude higher but the values

Dear R Experts--- Sorry for all the questions yesterday and today. I am trying to use Yves Croissant's pgmm function in the plm package with Blundell-Bond moments. I have read the Blundell-Bond paper, and want to run the simplest model first, d[i,t] = a*d[i,t-1] + fixed[i] + u[i,t] . no third conditioning variables yet. the full set of moment conditions recommended for system-GMM,

Hi all, I am attempting to use latin hypercube sampling to sample different variable functions in a series of simultaneous differential equations. There is very little code online about lhs or clhs, so from different other help threads I have seen, it seems I need to create a probability density function for each variable function, and then use latin hypercube sampling on this pdf. So far, I

Hello, I am using {plm} to estimate panel models. I want to estimate a model that includes fixed effects for time and individual, but has a random individual effect for the coefficient on the independent variable. That is, I would like to estimate the model: Y_it = a_i + a_t + B_i * X_it + e_it Where i denotes individuals, t denotes time, X is my independent variable, and B (beta) is the

Hello, I've been trying to fit Arrellano-Bond model with pgmm but I am getting very strange errors. I've looked around and found no reference to them. I've specified the model in dozens of different ways, and each seems to give me a new kind of error. This leads me to believe this has to do with the way the data is specified, but I can't see anything thats wrong with. My

When using the PLM package (version 1.2-8), I encounter the probem that calling the FGLS estimator evokes strange behavior, when choosing the "random" effects model. After calling the PGGLS function to estimate FGLS, PLM gives me a warning, stating that the "random" model has been replaced with the "pooling" model. I would, however, really like to estimate the random

dear giovanni--- thanks for answering on r-help to me as well as privately. I very much appreciate your responding. I read the plm vignette. I don't have the book, so I can't consult it. :-(. I am going to post this message now (rather than just email it privately), because other amateurs may have similar questions in the future, and find this message and your answers via google.

Hi, I'm trying to fit a fixed effect (LSDV) panelmodel with R. I have a dataset with y as dependent, x1&x2 as indeps, t as time index and i as an id-variable for each individual. There are three observations for each individual (t=1, t=2, t=3). I want to try a simple regression, but with individual intercepts: ------------------------------------------------- # reading in some data ...

Dear, I built the generalized method of moments model to estimate the sales rank in the bookstore using plm package in R. The equation is: data1.gmm <- pgmm(dynformula(lnsales_rank ~ ln_price + avg_ham_rate + avg_spam_rate + num_of_ham+ num_of_spam + ship_code2 +ship_code3 +ship_code4+ ship_code5+ ship_code6 + ship_ code7, lag = list(0, 0, 0, 0,0,0,0,0,0,0,0,0), log =FALSE), data=data,

It seems to only works for zoo or ts objects? I tried to run it for xts objects, and error message occurs. Once I coerce it to be an zoo object (by as.zoo), it works. Error message: Error in model.frame.default(formula = dynformula(PIh - PI ~ L(X, 0:i) + : variable lengths differ (found for 'L(X, 0:i)') In addition: Warning messages: 1: In zoo(coredata(x), order.by = index(x), ...)

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

Dear Limin, might be just about anything. Could you please provide a reproducible example? Best, Giovanni ----------------- Original message ---------------------- Message: 51 Date: Mon, 17 May 2010 10:36:03 +0800 (CST) From: ??? <dlmsos at 163.com> To: r-help at r-project.org Subject: [R] pgmm function Message-ID: <b2cba0.35fc.128a41e3684.Coremail.dlmsos at 163.com> Content-Type:

Hi, I want to estimate a dynamic paneldata model with the following code, but unfortenately I received the error message below. form<-PB~Activity+Solvency+Cap_Int

Hello, So I am afraid I am having a recurring problem that I just can't figure out. I am using the plm package to conduct a panel analysis - although I am not sure if the problem is arising as a result of the plm package or something more general. I am trying to run a fixed effects model with effects over time and individual. The model has various lags, and the problem is that these lags do

Hi, all! Anyone know an easy way to specify the following model. Panel dataset, with stock through time, by firm. I want to run a model of y on a bunch of explanatory variables, and one dummy for each firm, which is 1 for observations that come from firm i, and 0 everywhere else. I have over 200 firms (and a factor variable that contains a firm identifier). Any easy way of going about

Dear Sir/ Madam, I am writing about the panel data for my bachelor degree. I would really appreciate if You could help dealing with R functions. I am trying to estimate the panel data lm model with plm function. When i include 3dummy variables into the regression it dont appear in the sumarry of the model, but when i estimate a simple lm model it appears. Why is it so? What should i do to