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
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.