@@ -94,10 +94,19 @@ lm.RStests <- function(model, listw, zero.policy=attr(listw, "zero.policy"), tes
94
94
warning(" Out-of-range p-value: reconsider test arguments" )
95
95
names(p.value ) <- " "
96
96
method <- " Rao's score (a.k.a Lagrange multiplier) diagnostics for spatial dependence"
97
- data.name <- paste(" \n " , paste(strwrap(paste(" model: " ,
98
- gsub(" [ ]+" , " " , paste(deparse(model $ call ),
99
- sep = " " , collapse = " " )))), collapse = " \n " ),
100
- " \n weights: " , listw_name , " \n " , sep = " " )
97
+ fmla <- gsub(" [ ]+" , " " , paste(deparse(model $ call ), sep = " " ,
98
+ collapse = " " ))
99
+ if (GNM ) {
100
+ if (! is.null(attr(model , " SLX_call" ))) {
101
+ fmla <- paste(" lm(" , attr(model , " SLX_call" ), " )" , sep = " " )
102
+ } else {
103
+ fmla <- gsub(" [ ]+" , " " , paste(deparse(model $ terms ),
104
+ sep = " " , collapse = " " ))
105
+ }
106
+ }
107
+ data.name <- paste(" \n " , paste(strwrap(paste(" model: " , fmla )),
108
+ collapse = " \n " ),
109
+ " \n test weights: " , listw_name , " \n " , sep = " " )
101
110
tres [[i ]] <- list (statistic = statistic , parameter = parameter ,
102
111
p.value = p.value , method = method , data.name = data.name )
103
112
class(tres [[i ]]) <- " htest"
0 commit comments