Re: How to obtain standard errors of local regression coefficients in spgwr::ggwr()?
Does ?spgwr::ggwr state that any such values are returned? It says:
"a SpatialPointsDataFrame (may be gridded) or SpatialPolygonsDataFrame
object (see package "sp") with fit.points, weights, GWR coefficient
estimates, dispersion if a "quasi"-family is used, and the residuals of
type "type" in its "data" slot."
Nothing else is returned. You may like to try the GWmodel package, which
may offer some further output, but this is a largely unexplored area, and
in any case inference is not encouraged for GW output as the technique is
Hope this clarifies.
PS. Please post plain text, not HTML-formatted.
On Sun, 6 Jun 2021, mr via R-sig-Geo wrote:
> Dear list,
> I am using spgwr::ggwr() to fit generalized geographically weighted regression with Poisson model and log-link function. The results provide local coefficient estimates, but i am missing how to get their standard errors (or t statistics) to compute pseudo p-values.
> Below is a toy example using SpatialEpi::NYleukemia dataset:
> # --------
> ## Load data
> population <- NYleukemia$data$population
> cases <- ceiling(NYleukemia$data$cases * 100)
> centroids <- latlong2grid(NYleukemia$geo[, 2:3])
> # data frame
> nyleuk <- data.frame(centroids, cases, population)
> # set coordinates as vector
> coordny <- cbind(centroids[,1],centroids[,2])
> # set a kernel bandwidth
> bw <- 0.5
> # fit ggwr()
> m_pois <- ggwr(cases ~ offset(log(population)),
> data = nyleuk, gweight = gwr.Gauss,
> adapt = bw, family = poisson(link="log"),
> type="working", coords = coordny)
> # returns spatial point with coefficients
> # but no standard errors :(
> # -------
> Is there any way to get standard errors of the local coefficients?
> Thank you, Manuel
> Sent from Mail for Windows 10
> [[alternative HTML version deleted]]
> R-sig-Geo mailing list
> [hidden email] > https://stat.ethz.ch/mailman/listinfo/r-sig-geo >