<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"Microsoft YaHei UI";
panose-1:2 11 5 3 2 2 4 2 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"\@Microsoft YaHei UI";}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.DefaultFontHxMailStyle
{mso-style-name:"Default Font HxMail Style";
font-family:"Calibri",sans-serif;
color:windowtext;
font-weight:normal;
font-style:normal;
text-decoration:none none;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style></head><body lang=PT link=blue vlink="#954F72" style='word-wrap:break-word'><div class=WordSection1><p class=MsoNormal><span lang=EN-GB>Ok, thank you for your comments.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB>Regards, Manuel<span class=DefaultFontHxMailStyle><o:p></o:p></span></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><o:p> </o:p></span></p><p class=MsoNormal>Sent from <a href="https://go.microsoft.com/fwlink/?LinkId=550986">Mail</a> for Windows 10</p><p class=MsoNormal><span class=DefaultFontHxMailStyle><o:p> </o:p></span></p><div style='mso-element:para-border-div;border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal style='border:none;padding:0cm'><b>From: </b><a href="mailto:binbinlu@whu.edu.cn">binbinlu@whu.edu.cn</a><br><b>Sent: </b>4 de julho de 2021 09:24<br><b>To: </b><a href="mailto:manuel.ribeiro@tecnico.ulisboa.pt">mr</a>; <a href="mailto:r-sig-geo@r-project.org">r-sig-geo</a><br><b>Subject: </b>Re: [R-sig-Geo] Gwmodel an spgwr return too different local coefficients?</p></div><p class=MsoNormal><span class=DefaultFontHxMailStyle><o:p> </o:p></span></p><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>Hi,<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>In GWmodel, we didn't process the offset term, and that's why the results are very different from each other.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>Moreover, the bandwidths for them are also different, but this won't cause a big difference.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'><o:p> </o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>Cheers,<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>Binbin<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'><o:p> </o:p></span></p></div><p class=MsoNormal><img border=0 width=210 height=1 style='width:2.1875in;height:.0104in' id="Horizontal_x0020_Line_x0020_1" src="cid:image001.png@01D771AE.2B5A8600"><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'><o:p></o:p></span></p><div><div style='margin-left:7.5pt;margin-top:7.5pt;margin-right:7.5pt;margin-bottom:7.5pt'><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Verdana",sans-serif;color:black'>binbinlu@whu.edu.cn<o:p></o:p></span></p></div></div></div><blockquote style='margin-left:6.0pt;margin-top:5.0pt;margin-bottom:5.0pt'><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'><div><div><p class=MsoNormal style='background:#EFEFEF'><b><span style='font-size:9.0pt;font-family:"Tahoma",sans-serif;color:black'>From:</span></b><span style='font-size:9.0pt;font-family:"Tahoma",sans-serif;color:black'> <a href="mailto:r-sig-geo@r-project.org">mr via R-sig-Geo</a><o:p></o:p></span></p></div><div><p class=MsoNormal style='background:#EFEFEF'><b><span style='font-size:9.0pt;font-family:"Tahoma",sans-serif;color:black'>Date:</span></b><span style='font-size:9.0pt;font-family:"Tahoma",sans-serif;color:black'> 2021-07-04 01:11<o:p></o:p></span></p></div><div><p class=MsoNormal style='background:#EFEFEF'><b><span style='font-size:9.0pt;font-family:"Tahoma",sans-serif;color:black'>To:</span></b><span style='font-size:9.0pt;font-family:"Tahoma",sans-serif;color:black'> <a href="mailto:r-sig-geo@r-project.org">r-sig-geo@r-project.org</a><o:p></o:p></span></p></div><div><p class=MsoNormal style='background:#EFEFEF'><b><span style='font-size:9.0pt;font-family:"Tahoma",sans-serif;color:black'>Subject:</span></b><span style='font-size:9.0pt;font-family:"Tahoma",sans-serif;color:black'> [R-sig-Geo] Gwmodel an spgwr return too different local coefficients?<o:p></o:p></span></p></div></div></div><div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>Dear all, </span><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>I am using spgwr::ggwr() and GWmodel::ggwr.basic to fit a generalized geographically weighted regression with Poisson model and log-link function.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>Local coefficients between packages are very different, but i couldn't find out why. If any body can point me what am i missing, i apreciate a lot.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>Below is an example using spData::nc.sids. I would expect local coefficents will be around the observed mean, but this doesnt occur using the GWmodel::ggwr.basic function:<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>library (spData)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>library(GWmodel)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>library(spgwr)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>library(sp)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>data(nc.sids)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>nc.sids$SID79 <- nc.sids$SID79 + 1 <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'># formula of glm poisson model<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>f1 <- SID79 ~ offset(log(BIR79))<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'># GWmodel<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>sp::coordinates(nc.sids) <- ~ x + y<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>dM <- gw.dist(dp.locat = coordinates(nc.sids))<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>bw_gw <- bw.ggwr(formula(f1), nc.sids, family ="poisson", approach="AIC",<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> kernel="gaussian", adaptive = T, p = 2, theta = 0, longlat=F, dMat = dM)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>m_gw <- ggwr.basic(formula(f1), data = nc.sids, <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> regression.points = nc.sids, bw = bw_gw, <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> family = "poisson", kernel = "gaussian", <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> adaptive = T, cv = F, tol = 1e-5, <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> maxiter = 1000, dMat = dM)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'># spgwr<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>dfsids <- as.data.frame(nc.sids)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>xycoord <- cbind(dfsids$x,dfsids$y)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>bw_sp <- ggwr.sel(formula(f1), data = dfsids , family = poisson(link = "log"), <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> gweight = gwr.Gauss, adapt = T, coords = xycoord, <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> tol = 1e-5 , verbose = T) <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>m_sp <- ggwr(formula(f1), data = dfsids, gweight = gwr.Gauss,<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> adapt = bw_sp, family = poisson(link="log"), <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> type="working", coords = xycoord) <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'># compare results<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>obs <- mean(log(nc.sids$SID79 / nc.sids$BIR79))<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>gw <- mean(m_gw$SDF$Intercept)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>sp <- mean(m_sp$SDF$X.Intercept.)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>cbind(obs, sp, gw)<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>Sent from Mail for Windows 10<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>[[alternative HTML version deleted]]<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>_______________________________________________<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>R-sig-Geo mailing list<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>R-sig-Geo@r-project.org<o:p></o:p></span></p></div></div></blockquote><p class=MsoNormal style='mso-margin-top-alt:0cm;margin-right:36.0pt;margin-bottom:0cm;margin-left:6.0pt;margin-bottom:.0001pt'><span style='font-size:10.5pt;font-family:"Microsoft YaHei UI",sans-serif;color:black'>https://stat.ethz.ch/mailman/listinfo/r-sig-geo<o:p></o:p></span></p><p class=MsoNormal><span class=DefaultFontHxMailStyle><o:p> </o:p></span></p></div></body></html>