<div dir="ltr"><br><div>O modelo de acordo com o que entendi é apenas um intercepto aleatório pros blocos e o efeito espacial, se é possível estimar na função gls() eu não sei, pelo que o Walmes fala é possível estruturar os resíduos sem problemas mas colocar mais um efeito aleatório de bloco acho que não. O Walmes pode confirmar ... <br>
</div><div>Agora sem olhar seus dados o modelo teoricamente é estimável como é tudo gaussiano, você pode escrever a verossimilhança e maximizar numericamente, deve ser um problema bem estável. No nosso livro do SINAPE tem um exemplo de geoestatístico pode te ajudar a começar. Tudo que precisa fazer é somar mais uma matriz na matriz de covariâncias marginal.<br>
<br></div><div><a href="http://www.leg.ufpr.br/mcie">www.leg.ufpr.br/mcie</a><br></div><div><br></div><div><br></div><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 18 de dezembro de 2013 13:07, Luiz Roberto Martins Pinto <span dir="ltr"><<a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Penso que 'Blocos' não pode ser considerado como efeito fixo.<div>Blocos são agrupamentos particulares das repetições. Portanto são repetições.</div>
<div>Para que os testes estatísticos sejam válidos as repetições tem que ser aleatórias, caso contrário o experimento não tem valor, posto que as repetições são tendenciosas e as conclusões não podem ser inferidas para a população.</div>

<div>Além disto o resíduo (aleatório) contém a interação blocos*tratamentos, que, portanto, deve ser aleatória. Note que no delineamento em que os tratamentos são aleatorizados nos blocos o gl do resíduo é (t-1)*(b-1). O resíduo é sempre aleatório porque o bloco é aleatório.</div>

<div> </div></div><div class="gmail_extra"><br clear="all"><div><div>Luiz Roberto Martins Pinto<br>Prof. Pleno/DCET/UESC</div><div>Laboratório de Estatística Computacional</div><div>Universidade Estadual de Santa Cruz</div>

<div>Ilhéus-Bahia<br><br><a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a><br>skype: lrmpinto</div>
<div><a href="http://lattes.cnpq.br/2732314327604831" target="_blank">http://lattes.cnpq.br/2732314327604831</a> <br><br><br></div></div>
<br><br><div class="gmail_quote">Em 18 de dezembro de 2013 10:37, Wagner Bonat <span dir="ltr"><<a href="mailto:wbonat@gmail.com" target="_blank">wbonat@gmail.com</a>></span> escreveu:<div><div class="h5"><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div dir="ltr"><div>Por que não considera os blocos como efeito fixo ?<br><br></div>Você vai assumir que a estrutura do efeito espacial é sempre a mesma dentro dos blocos ?<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">


Em 17 de dezembro de 2013 16:00, Cássio Dessotti <span dir="ltr"><<a href="mailto:cassiodessotti@yahoo.com.br" target="_blank">cassiodessotti@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div><div>
<div><div style="font-size:12pt;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif"><div><span style="background-color:transparent;font-size:12pt">Olá pessoal, boa tarde!! </span><br></div>

<div style="display:block">
<div style="font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif;font-size:12pt"><div style="font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif;font-size:12pt">


<div><div><div><div style="font-size:12pt;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif"><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif">


<br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif">Estou trabalhando com um grupo de experimentos em blocos ao acaso, com dados georreferenciados (x,y), buscando comparar modelos mistos (usuais) com modelos mistos que levam em conta a<span style="font-size:12pt"> dependência espacial, utilizando de</span><span style="background-color:transparent"> funções geoestatísticas na matriz de resíduos (R). </span></div>


<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif"><span style="background-color:transparent"><br>


</span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif"><span style="background-color:transparent">Consigo construir o modelo M1 (BA) com efeito de blocos aleatório a partir da função "lme", além do modelo M2 (Exp-H) que desconsidera o efeito de blocos, e considera uma função exponencial na matriz de covariâncias para os resíduos, além de considerar variâncias diferentes para os 2 locais.</span></div>


<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif"><span style="background-color:transparent"><br>


</span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif"><span style="background-color:transparent">1) A primeira questão está no modelo M3 (BA-Exp-H) que deve "unir" as características dos dois anteriores, considerando efeito de blocos aleatório, função exponencial para a matriz de covariâncias, além de heterogeneidade de variâncias, porém não consegui tal realização nem por meio da função "lme", nem por "gls". Seria necessário o uso de uma outra função? Como eu poderia trabalhar com este modelo no R?  </span></div>


<div><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif">2) A segunda questão é o seguinte: a partir do modelo que eu selecionar (provavelmente por AIC), desejo extrair os resíduos condicionais estudentizados, e separá-los segundo seus locais, para que em cada local eu possa construir semivariogramas destes resíduos, verificando a "força" da dependência, e realizar a krigagem de cada área (trabalhando com os resíduos para "limpar" os dados de seus efeitos de tratamentos e demais efeitos). Posso seguir esta linha? Acredito que é um ganho de informação trabalhar com todos os locais (análise conjunta) neste sentido, ao invés de fazer todo esse processo em cada local separadamente.</div>


<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif"><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif">


Muito obrigado desde já pela atenção.</div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif">

Abraço a todos.</div>
<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif">Cássio Dessotti.</div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif">


<br></div><div>### Segue o código (com dados fictícios) para exemplificar o meu problema:<span style="background-color:transparent"> </span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:HelveticaNeue,'Helvetica Neue',Helvetica,Arial,'Lucida Grande',sans-serif">


local <- as.factor(c(rep(1,6),rep(2,6))); blo <- as.factor(rep(c(rep(1,3),rep(2,3)),2))</div><div>lat <- as.factor(rep(1:3,4)) ; long <- as.factor(rep(c(rep(1,3),rep(2,3)),2))</div><div>blo_local <- as.factor(c(rep("1_1",3),rep("2_1",3),rep("1_2",3),rep("2_2",3)))</div>


<div>trat <- as.factor(c(2,3,1,1,3,2,1,3,2,3,2,1))</div><div>resp <- c(2,1.9,1.9,1.8,2,1.9,2,2.1,2.3,2,1.9,1.8)</div><div><span style="font-size:12pt">dados <-
 data.frame(local,trat,blo,lat,long,blo_local,resp)</span><br></div><div><br></div><div>require(nlme)</div><div>### M1 - BA (efeito de blocos aleatório)</div><div>M1 <- lme(resp ~ 1 + trat + local + local:trat,</div><div>


random = list(blo_local = pdIdent(~1)),</div><div>method = "REML", na.action = na.omit, data=dados, keep.data=FALSE)</div><div><br></div><div>### M2 - Exp H (sem efeito de blocos - função exponencial para R - heterogeneidade de variâncias nos locais)</div>


<div>M2 <- gls(resp ~ 1 + trat + local + local:trat, weight=varComb(varIdent(form = ~ 1|local)),</div><div>correlation=corExp(form=~as.numeric(as.character(lat))+as.numeric(as.character(long))|local,</div><div>metric="euclidean", nugget=FALSE), method="REML", na.action=na.omit, data=dados)</div>


<div><br></div><div>### M3 - BA-Exp-H (PROBLEMA - blocos aleatórios - função exponencial - heterogeneidade de variâncias) </div><div>M3 <- gls(resp ~ 1 + trat + local +
 local:trat, random = list(blo_local = pdIdent(~1)),</div><div>weight=varComb(varIdent(form = ~ 1|local)),</div><div>correlation=corExp(form = ~ as.numeric(as.character(lat))+as.numeric(as.character(long))|local,</div><div>


metric="euclidean", nugget=FALSE), method="REML", na.action=na.omit, data=dados)</div></div></div></div><br></div>  </div> </div>  </div> </div></div><br></div></div>_______________________________________________<br>



R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a 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>
Leia o guia de postagem (<a 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.<span><font color="#888888"><br></font></span></blockquote>
</div><span><font color="#888888"><br><br clear="all"><br>-- <br>Wagner Hugo Bonat<br>
LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná
</font></span></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a 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>
Leia o guia de postagem (<a 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.<br></blockquote></div></div></div><br></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a 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>
Leia o guia de postagem (<a 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.<br></blockquote></div><br><br clear="all"><br>-- <br>Wagner Hugo Bonat<br>
LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná
</div>