return Durbin factor objects

Author: rsbivand

R/ML_models.R

@@ -154,6 +154,7 @@ errorsarlm <- function(formula, data = list(), listw, na.action, weights=NULL,
                     attr(dvars, "f") <- Durbin
                     attr(dvars, "inds") <- inds
                     attr(dvars, "zero_fill") <- zero_fill
+                    attr(dvars, "formula_durbin_factors") <- formula_durbin_factors
                 }
 		x <- cbind(x, WX)
 		m <- NCOL(x)
@@ -507,6 +508,8 @@ errorsarlm <- function(formula, data = list(), listw, na.action, weights=NULL,
                 hf_calls=get("hf_calls", envir=env), intern_classic=iC,
                 pWinternal=pWinternal, weights=weights, emixedImps=emixedImps,
                 dvars=dvars), class=c("Sarlm"))
+        attr(ret, "have_factor_preds") <- have_factor_preds
+
         rm(env)
         GC <- gc()
 	if (zero.policy) {
@@ -680,6 +683,7 @@ lagsarlm <- function(formula, data = list(), listw,
                     attr(dvars, "f") <- Durbin
                     attr(dvars, "inds") <- inds
                     attr(dvars, "zero_fill") <- zero_fill
+                    attr(dvars, "formula_durbin_factors") <- formula_durbin_factors
                 }
 		x <- cbind(x, WX)
 		m <- NCOL(x)
@@ -887,6 +891,7 @@ lagsarlm <- function(formula, data = list(), listw,
                 f_calls=get("f_calls", envir=env),
                 hf_calls=get("hf_calls", envir=env), intern_classic=iC),
                 class=c("Sarlm"))
+        attr(ret, "have_factor_preds") <- have_factor_preds
         rm(env)
         GC <- gc()
 	if (zero.policy) {
@@ -1051,6 +1056,7 @@ sacsarlm <- function(formula, data = list(), listw, listw2=NULL, na.action,
                     attr(dvars, "f") <- Durbin
                     attr(dvars, "inds") <- inds
                     attr(dvars, "zero_fill") <- zero_fill
+                    attr(dvars, "formula_durbin_factors") <- formula_durbin_factors
                 }
 		x <- cbind(x, WX)
 		m <- NCOL(x)
@@ -1329,6 +1335,7 @@ sacsarlm <- function(formula, data = list(), listw, listw2=NULL, na.action,
             optimHess=FALSE, insert=FALSE, interval1=interval1,
             interval2=interval2, timings=do.call("rbind", timings)[, c(1, 3)]),
             class=c("Sarlm"))
+        attr(ret, "have_factor_preds") <- have_factor_preds
         rm(env)
         GC <- gc()
         if (is.null(llprof)) ret$llprof <- llprof

R/SLX_WX.R

@@ -100,6 +100,7 @@ lmSLX <- function(formula, data = list(), listw, na.action, weights=NULL, Durbin
             attr(dvars, "f") <- Durbin
             attr(dvars, "inds") <- inds
             attr(dvars, "zero_fill") <- zero_fill
+            attr(dvars, "formula_durbin_factors") <- formula_durbin_factors
         }
 	x <- cbind(x, WX)
 	rm(WX)
@@ -214,6 +215,7 @@ lmSLX <- function(formula, data = list(), listw, na.action, weights=NULL, Durbin
         tms <- as.character(lm.model$terms)
         attr(lm.model, "SLX_call") <- paste(y_name, " ", tms[1], " ", tms[3],
             ", data = ", d_name, ", listw = ", listw_name, sep="")
+        attr(lm.model, "have_factor_preds") <- have_factor_preds
         class(lm.model) <- c("SlX", class(lm.model))
         lm.model
 }

R/spBreg.R

@@ -130,6 +130,7 @@ spBreg_lag <- function(formula, data = list(), listw, na.action, Durbin, type,
             attr(dvars, "f") <- Durbin
             attr(dvars, "inds") <- inds
             attr(dvars, "zero_fill") <- zero_fill
+            attr(dvars, "formula_durbin_factors") <- formula_durbin_factors
         }
 	x <- cbind(x, WX)
 	m <- NCOL(x)
@@ -386,6 +387,7 @@ spBreg_lag <- function(formula, data = list(), listw, na.action, Durbin, type,
     attr(res, "acc_rate") <- acc_rate
     attr(res, "dvars") <- dvars
     attr(res, "MH") <- priors$rhoMH
+    attr(res, "have_factor_preds") <- have_factor_preds
     class(res) <- c("MCMC_sar_G", class(res))
     res
 
@@ -613,6 +615,7 @@ spBreg_err <- function(formula, data = list(), listw, na.action, Durbin, etype,
             attr(dvars, "f") <- Durbin
             attr(dvars, "inds") <- inds
             attr(dvars, "zero_fill") <- zero_fill
+            attr(dvars, "formula_durbin_factors") <- formula_durbin_factors
         }
 	x <- cbind(x, WX)
         xcolnames <- colnames(x)
@@ -948,6 +951,7 @@ spBreg_err <- function(formula, data = list(), listw, na.action, Durbin, etype,
     attr(res, "n") <- n
     attr(res, "k") <- k
     attr(res, "MH") <- priors$lambdaMH
+    attr(res, "have_factor_preds") <- have_factor_preds
     class(res) <- c("MCMC_sem_G", class(res))
     res
 
@@ -1112,6 +1116,7 @@ spBreg_sac <- function(formula, data = list(), listw, listw2=NULL, na.action,
             attr(dvars, "f") <- Durbin
             attr(dvars, "inds") <- inds
             attr(dvars, "zero_fill") <- zero_fill
+            attr(dvars, "formula_durbin_factors") <- formula_durbin_factors
         }
 	x <- cbind(x, WX)
 	m <- NCOL(x)
@@ -1383,6 +1388,7 @@ spBreg_sac <- function(formula, data = list(), listw, listw2=NULL, na.action,
     attr(res, "acc_rate1") <- acc_rate1
     attr(res, "acc_rate2") <- acc_rate2
     attr(res, "dvars") <- dvars
+    attr(res, "have_factor_preds") <- have_factor_preds
     class(res) <- c("MCMC_sac_G", class(res))
     res