[R-br] Peso espacial em X usando spreg- regressão espacial

[Elias T. Krainski]

como vc tem um modelo do tipo

             y=lambda W y + X beta + u
             u=rho W u + e

automaticamente ja tem X lagged porque E(Y) = Xbeta. Por isso nao 
precisa criar XW.

Elias


On 21/11/2017 11:12, Edimeire Alexandra Pinto wrote:
> Olá pessoal.
>
> OI Elias.
>
> 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.
>
> 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.
>
> 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: *system 
> is exactly singular*
>
> Pode ser que tenho começar a "abrir" a função no R para entendê-la 
> melhor....
>
> library(spdep)
> require(sphet)
> data(columbus)
> listw<-nb2listw(col.gal.nb)
> head(columbus)
> WINC<-lag.listw(listw,columbus$INC,NAOK=TRUE,zero.policy = TRUE)
> columbus$WINC<-WINC
> res<- spreg(CRIME~HOVAL + INC+WINC+PLUMB, data=columbus , listw= 
> listw,   het = TRUE,
>     verbose = FALSE,model = "sarar")
>
> Error in solve.default(HH, Hye) :
> Lapack routine dgesv: system is exactly singular: U[7,7] = 0
>
>
>
>
>
> 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.
>
>
> Obrigada pessoal , obrigada GRANDE Elias.
>
>
>
> Em Terça-feira, 21 de Novembro de 2017 9:11, Elias T. Krainski via 
> R-br <r-br em listas.c3sl.ufpr.br> escreveu:
>
>
> Oi Edimeire,
> 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"?
>
> Elias
>
> On 20/11/2017 18:10, Edimeire Alexandra Pinto via R-br wrote:
> Boa tarde, pessoal.
>
> 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.
>
> Por exemplo, acrescentando peso espacial a uma das variaveis 
> explicativas do modelo.
>
> library(spdep) data(columbus) listw<-nb2listw(col.gal.nb)
> 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 = TRUE, verbose = 
> FALSE,model="sarar") ### WINC é o peso espacial multiplicado #pelo 
> variavel INC, especie de W*INC
>
> ou mesmo
>
> res3<- spreg(CRIME~HOVAL + INC+WINC, data=columbus , listw= listw, het 
> = TRUE, verbose = FALSE,model="error")
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br <mailto:R-br em listas.c3sl.ufpr.br>
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br <mailto:R-br em listas.c3sl.ufpr.br>
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a 
> c�digo m�nimo reproduz�vel.
>
>

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20171121/514b2ece/attachment.html>