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

[Edimeire Alexandra Pinto]

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
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
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/83374ef7/attachment.html>