similar to: corrigendum on fixed effects and R2 in within models

Hi R-helpers, I have two questions I hope you could help me with them: In the plm package how can I calculate the R2 within, R2 between and R2 overall? Is there any special reason to not display these values? When using first differences do I need to have some special care with dummies (both year dummies and industry dummies)? (A friend who works with Stata told me that there is

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

Hello again Max. A correction to my response from yesterday. Things were better than they seemed. I thought it over, checked Arellano's panel book and Driscoll and Kraay (Rev. Econ. Stud. 1998) and finally realized that vcovSCC does what you want: in fact, despite being born primarily for dealing with cross-sectional correlation, 'SCC' standard errors are robust to "both

Hi, Try this: which(duplicated(res10Percent)) # [1] 117 125 157 189 213 235 267 275 278 293 301 327 331 335 339 367 369 371 379 #[20] 413 415 417 441 459 461 477 479 505 res10PercentSub1<-subset(res10Percent[which(duplicated(res10Percent)),],dummy==1)? #most of the duplicated are dummy==1 res10PercentSub0<-subset(res10Percent[which(duplicated(res10Percent)),],dummy==0)

Thank you all!!! The problem was the decimal symbol! My data was saved in a txt file, so I?ve introduced the dec="," in ?read.table? and it worked. What I?ve done was coeficientes<-read.table("coeficientes.txt",sep="\t",header=T,dec=",") Then, subset worked fine coeficientesWanted<-subset(coeficientes,b1>0) Thanks again, Cec?lia Carmo

Hello, Could anybody offer any advice about implementing a Kohonen self organizing map for categorical data? Specifically I am wondering if there are any pre-existent packages that can deal with categorical data and/or how one would compare the input vector of categoricals with the self organizing map nodes. Thanks in advance. George Chen ----- Original Message ----- From: r-help-request at

--- the forwarded message follows --- -------------- next part -------------- An embedded message was scrubbed... From: "Cecilia Carmo" <cecilia.carmo at ua.pt> Subject: Re: [R] subset dataframe/list Date: Mon, 01 Jun 2009 21:33:15 +0100 Size: 3657 URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090601/921f7638/attachment-0002.mht>

Hi, Try this: final3New<-read.table(file="real_data_cecilia.txt",sep="\t") dim(final3New) #[1] 5369??? 5 #Inside the split within split, dummy==1 for the first row.? For lists that have many rows, I selected the row with dummy==0 (from the rest) using the #condition that the absolute difference between the dimensions of those rows and the first row dimension was minimum

Hi everyone: I need to count the number of banks of each firm in my data. The firm is identified by the fiscal number. The banks of each firm appears like this: Firm Banks 500600700 Citybank 500600700 CGD 500600700 BES 500600800 Citybank 500600800 Bank1 500600900 CGD I want to obtain the following dataframe: Firm

Hi R-helpers, I want to determine the coefficients of the following regression for several subsets, and I want to save it in a dataframe: The data is in ?regaccdis?, ?regaccdis$caedois? is the column that defines the subsets and the function I have runned is coef(plm(ff,data=regaccdis,na.action=na.omit,model="pooling",subset=(regaccdis$caedois==i))) I?ve created a dataframe named

Hi R-helpers I want a function that performs the opposite of lag() with panel data. I have transformed my data before with pdata.frame(mydata, index=c("groupindex", “timeindex")) And then I’ve done lag(mydata, -1) but it doesn’t work. The error message was: Error in rep(1, ak) : invalid ''times'' argument Thank you in advance, Cecília Carmo

Hi everyone. I have two questions. I’ve found some other questions and answers similar to these but they didn’t solve my problem. I’m working with a panel of firm/years observations (see my reproducible example). I’m using the plm package. My panel not only is unbalanced but also have some gap’s in years. #reproducible example

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

Hi everyone, Could anyone tell me what means the follow error message Error in xj[i] : invalid subscript type 'closure' It happens when I run the function plm, like this: >ff<-totaccz~lactivoz+varvolnegz >ss<-plm(ff,data=regaccdis,na.omit) Error in xj[i] : invalid subscript type 'closure' > coef(ss) (Intercept) lactivoz varvolnegz 0.02571212 6.94227541

Hi All, I am doing econometric modeling of panel data (fixed effects). We currently use Eviews to do this, but I have discovered a bug in Eviews 7 and am exploring the use of R to build panel data models / longitudinal data models. I looked at the plm package but do not see how I can incorporate AR terms in the model using the plm package. I have an Eviews model with two AR terms, AR(1) and

Hello to all on the list, I'm trying to estimate a fixed effects model from a large (unbalanced) panel data set. I have no problems when using only an individual effect or only a time effect, but I get an error message when I try for a "twoways" effect. Here is some of the code: paneldata27 is the entire panel data set: > dim(paneldata27) [1] 1178831 8 >

Hi everyone, I’m working with a panel of firm/years observations. My panel not only is unbalanced but also have some gap’s in years. For example, firm 1 has 1999, 2000, 2001, 2004, 2005, firm 2 has 2000, 2001, 2003, 2005, and so on. I’m using the plm package and what I’m asking is how can I handle with this gap’s ? Thank you very much, Cecília Carmo Universidade de Aveiro

Hi everyone! I?m running the pooltest in plm package, like this pooltest(cstfin12~lmaccdiscrz+lcobjur+lliq+lcollateral+ldimensao, data = dados3, model = "within") But I got the following error: Error in FUN(X[[1L]], ...) : insufficient number of observations My data is an unbalanced panel with 20907 observations (6971 individuals and years 2001 to 2007). This is not enough?

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 R-helpers, First of all, sorry for those who have (eventually) already received that request. The mail has been bumped several times, so I am not sure the list has received it... and I need help (if you have time)! ;-) I have a very simple question and I really hope that someone could help me I would like to estimate a simple fixed effect regression model with clustered standard errors by