[Splm-commits] r41 - pkg

[noreply at r-forge.r-project.org]

Author: gpiras
Date: 2009-10-13 16:08:59 +0200 (Tue, 13 Oct 2009)
New Revision: 41

Removed:
   pkg/felag.R
Log:
eliminate all functions erroneously  uploaded

Deleted: pkg/felag.R
===================================================================
--- pkg/felag.R	2009-10-13 14:08:44 UTC (rev 40)
+++ pkg/felag.R	2009-10-13 14:08:59 UTC (rev 41)
@@ -1,55 +0,0 @@
-`felag` <-
-function(y,x,wy,ysms,xsms,ytms, xtms, wytms, wysms, beta,sige,yt,xt,N,T,NT,k,effects,method, rho,listw,inde){
-		mx<-apply(x,2,mean)
-		intercept <- mean(y)- mean(wy)*rho -  mx%*%beta
-if (effects=="spfe"){
-	res.sfe <- as.matrix(ysms) - as.matrix(wysms) *rho - xsms %*% as.matrix(beta) - as.numeric(intercept)
-	xhat <- x %*% as.matrix(beta) + rep(res.sfe,T) + as.numeric(intercept)
-	res.t.sfe <- res.sfe / sqrt(sige / T* rep(1,N) + diag(as.numeric(sige) * xsms%*% crossprod(xt)%*% t(xsms) ))
-	res.t.con <- as.numeric(intercept) / sqrt(as.numeric(sige) / NT + as.numeric(sige) * t(as.matrix(mx)) %*% 	crossprod(xt) %*% as.matrix(mx))
-	N.vars <- k + N
-	res.e <- y - xhat - rho* wy
-FE.out<-list(res.sfe=res.sfe, res.t.sfe=res.t.sfe, intercept=intercept, 	res.t.con=res.t.con,xhat=xhat,N.vars=N.vars,res.e=res.e)
-	}
-if (effects== "tpfe")	{
-	res.tfe <- as.matrix(ytms) - as.matrix(wytms)* rho - xtms %*% as.matrix(beta) - as.numeric(intercept)
-	xhat <- x %*% as.matrix(beta) + rep(res.tfe,each=N) + as.numeric(intercept)
-	res.t.tfe <- res.tfe / sqrt(sige / N* rep(1,T) + diag(as.numeric(sige) * xtms%*% crossprod(xt)%*% t(xtms) ))
-	res.t.con <- as.numeric(intercept) / sqrt(as.numeric(sige) / NT + as.numeric(sige) * t(as.matrix(mx)) %*% 	crossprod(xt) %*% as.matrix(mx))
-	N.vars <- k + T
-	res.e <- y - xhat - rho* wy
-FE.out<-list(res.tfe=res.tfe, res.t.tfe=res.t.tfe, intercept=intercept, 	res.t.con=res.t.con,xhat=xhat,N.vars=N.vars,res.e=res.e)
-		}
-if (effects== "sptpfe"){
-	res.sfe <- as.matrix(ysms) - as.matrix(wysms) * rho - xsms %*% as.matrix(beta) - as.numeric(intercept)
-	res.tfe <- as.matrix(ytms) - as.matrix(wytms) * rho - xtms %*% as.matrix(beta) - as.numeric(intercept)
-	res.t.sfe <- res.sfe / sqrt(sige / T* rep(1,N) + diag(as.numeric(sige) * xsms%*% crossprod(xt)%*% t(xsms) ))
-	res.t.tfe <- res.tfe / sqrt(sige / N* rep(1,T) + diag(as.numeric(sige) * xtms%*% crossprod(xt)%*% t(xtms) ))
-	res.t.con <- as.numeric(intercept) / sqrt(as.numeric(sige) / NT + as.numeric(sige) * t(as.matrix(mx)) %*% crossprod(xt) %*% as.matrix(mx))
-	xhat<- x %*% as.matrix(beta) + rep(res.sfe,T) + rep(res.tfe,each=N) + as.numeric(intercept)
-	N.vars <- k + N + T - 1
-	res.e <- y - xhat - rho* wy
-FE.out<-list(res.tfe=res.tfe, res.t.tfe=res.t.tfe, res.sfe=res.sfe, res.t.sfe=res.t.sfe, intercept=intercept, res.t.con=res.t.con,xhat=xhat,N.vars=N.vars,res.e=res.e)
-		}
-if (effects=="pooled") {
-	xhat <-   x %*% as.matrix(beta)
-	res.e <- y - xhat - rho* wy
-	FE.out<-list(xhat=xhat,N.vars=k,res.e=res.e)
-	}
-if (method=="eigen"){
-	IrWi<-invIrW(listw,rho)
-	xtb <- xt %*% beta
-	yhat <- unlist(tapply(xhat,inde, function(u) IrWi %*% u))
-	ywhat <- unlist(tapply(xtb,inde, function(u) IrWi %*% u))
-	r1 <- as.matrix(yt - mean(yt))
-	r2 <- as.matrix(ywhat - mean(ywhat))
-	r1r2 <- crossprod(r1,r2)
-	r1r1 <- crossprod(r1)
-	r2r2 <- crossprod(r2)
-	res.corr <- as.numeric(r1r2^2) / (as.numeric(r1r1)*as.numeric(r2r2))
-}
-else res.corr <- NULL 
-FE.out <- list(FE.out, res.corr=res.corr)
-FE.out
-	}
-