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@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@listas.c3sl.ufpr.br <mailto:R-br@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@listas.c3sl.ufpr.br <mailto:R-br@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.