<html><head></head><body><div style="color:#000; background-color:#fff; font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px"><div id="yui_3_16_0_ym19_1_1511268305376_9111"><span>Olá pessoal.</span></div><div id="yui_3_16_0_ym19_1_1511268305376_9110"><span><br></span></div><div id="yui_3_16_0_ym19_1_1511268305376_9109"><span>OI Elias.</span></div><div id="yui_3_16_0_ym19_1_1511268305376_9112"><span><br></span></div><div id="yui_3_16_0_ym19_1_1511268305376_9108" dir="ltr">Recentemente estou estudando a função spreg para tentar resolver o problema de singularidade que dá. Ocorre apenas quando coloco modelos do tipo: Y=beta1X1+beta2WX1+beta3X2+ erro onde erro pode seguir o modelo SEM. Também ocorre em modelos Sarar.</div><div id="yui_3_16_0_ym19_1_1511268305376_9108" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1511268305376_9108" dir="ltr">O importante é entender que o erro ocorre sempre que ponho na equação X1 e WX1 em conjunto, pois se não coloco os dois ao mesmo tempo não dá erro. </div><div id="yui_3_16_0_ym19_1_1511268305376_9108" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1511268305376_9108" dir="ltr">O erro que ocorre em minha base de dados é similar a esse do exemplo abaixo, pois a mensagem de erro é sempre a mesma, estranho que no outro programa do Stata roda. No R, em meus dados, dá sempre: <b id="yui_3_16_0_ym19_1_1511268305376_19547">system is exactly singular</b></div><div id="yui_3_16_0_ym19_1_1511268305376_9108" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1511268305376_9108" dir="ltr">Pode ser que tenho começar a "abrir" a função no R para entendê-la melhor....</div><div id="yui_3_16_0_ym19_1_1511268305376_9108" dir="ltr"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9481">library(spdep)</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9482">require(sphet)</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9483">data(columbus)</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9484">listw<-nb2listw(col.gal.nb)</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9485">head(columbus)</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9486">WINC<-lag.listw(listw,columbus$INC,NAOK=TRUE,zero.policy = TRUE)</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9486">columbus$WINC<-WINC<br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9487">res<- spreg(CRIME~HOVAL + INC+WINC+PLUMB, data=columbus , listw= listw, het = TRUE,</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9488"> verbose = FALSE,model = "sarar")</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9489"><br id="yui_3_16_0_ym19_1_1511268305376_9490"></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9491">Error in solve.default(HH, Hye) : </div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"> Lapack routine dgesv: system is exactly singular: U[7,7] = 0</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492">Bem, faço vários tipos de modelos, SAR, SEM, DURBIN, etc etc, mas necessito de modelo que use tenha a variável X1 e WX1 conjuntamente. </div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492">Obrigada pessoal , obrigada GRANDE Elias.</div><div dir="ltr" id="yui_3_16_0_ym19_1_1511268305376_9492"><br></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 13px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"><font size="2" face="Arial"> Em Terça-feira, 21 de Novembro de 2017 9:11, Elias T. Krainski via R-br <r-br@listas.c3sl.ufpr.br> escreveu:<br></font></div> <br><br> <div class="y_msg_container"><div id="yiv7115735874"><div>
<div>Oi Edimeire, <br clear="none">
</div>
Acho que a implementação no 'spdep' não é muito restritiva e você
pode usar pesos (adequados) que você mesmo definir. Como são esses
pesos "CLIFF-ORD"? <br clear="none">
<br clear="none">
Elias<br clear="none">
<br clear="none">
<div class="yiv7115735874yqt3929151011" id="yiv7115735874yqt93730"><div class="yiv7115735874moz-cite-prefix">On 20/11/2017 18:10, Edimeire Alexandra
Pinto via R-br wrote:<br clear="none">
</div>
<blockquote type="cite">
</blockquote></div></div><div class="yiv7115735874yqt3929151011" id="yiv7115735874yqt88986"><div><div>
<div id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550">Boa tarde,
pessoal.</div>
<div id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><br clear="none">
</div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550">Alguém
por acaso sabe com se trabalha com modelos "parecidos" com
CLIFF-ORD usando spreg do pacote sphet? Sei que o Stata faz,
mas gostaria de fazer no R. Estou trabalhando com vários
tipos de modelos espaciais (na realidade, quase todos do
pacote sphet) e esse me gerou muita dúvida, pois preciso de
que uma das variaveis explicativas com e sem peso espacial da
matriz W.</div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><br clear="none">
</div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550">Por
exemplo, acrescentando peso espacial a uma das variaveis
explicativas do modelo.</div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><br clear="none">
</div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_5827">library</span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_5828">(spdep)
data(columbus)
listw<-nb2listw(col.gal.nb)</span></div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_20356"><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_20348">res1<-
spreg(CRIME~HOVAL + INC+ INC, data=columbus , listw=
listw, het = TRUE, verbose = FALSE,model="sarar")
res2<- spreg(CRIME~HOVAL + INC+ WINC, data=columbus ,
listw= listw, het = </span></span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_5829">TRUE</span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_5830">, verbose = </span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_5831">FALSE,model="sarar"</span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_5832">) ### WINC é o
peso espacial multiplicado #pelo variavel INC, especie de
W*INC</span><br clear="none" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_13872">
</div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><span><br clear="none">
</span></div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_21938">ou mesmo</span></div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><span><br clear="none">
</span></div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_15477">res3<-
spreg(CRIME~HOVAL + INC+WINC, data=columbus , listw= listw,
het = </span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_15478">TRUE</span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_15479">, verbose = </span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_15480">FALSE,model="error"</span><span id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_15481">) </span><span><br clear="none">
</span></div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><span><br clear="none">
</span></div>
<div dir="ltr" id="yiv7115735874yui_3_16_0_ym19_1_1511207855320_2550"><br clear="none">
</div>
</div>
<br clear="none">
<fieldset class="yiv7115735874mimeAttachmentHeader"></fieldset>
<br clear="none">
<pre>_______________________________________________
R-br mailing list
<a rel="nofollow" shape="rect" class="yiv7115735874moz-txt-link-abbreviated" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a rel="nofollow" shape="rect" class="yiv7115735874moz-txt-link-freetext" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a rel="nofollow" shape="rect" class="yiv7115735874moz-txt-link-freetext" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</pre>
<br clear="none">
</div></div></div><div class="yqt3929151011" id="yqt27137">_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">Leia o guia de postagem (<a shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div><br><br></div> </div> </div> </div></div></body></html>