<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt"><div id="yiv5585344374"><div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 12pt;"><div id="yiv5585344374yui_3_13_0_ym1_6_1387415804491_12">Luiz Roberto e Wagner</div><div id="yiv5585344374yui_3_13_0_ym1_6_1387415804491_12"><br clear="none"></div><div id="yiv5585344374yui_3_13_0_ym1_6_1387415804491_12" style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; background-color: transparent; font-style: normal;">Acredito que a utilização de blocos como sendo de efeito aleatório geralmente é mais adequada, porém também estou fazendo o modelo com blocos fixo a nível de comparação com os
demais (e neste não tive problemas no R - função gls), o problema está quando quero fazer o modelo (M3) com efeito de blocos aleatório, e função exponencial para a matriz de covariâncias dos resíduos. </div><div id="yiv5585344374yui_3_13_0_ym1_6_1387415804491_12" style="color:rgb(0, 0, 0);font-size:16px;font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;background-color:transparent;font-style:normal;" class="yui_3_13_0_ym1_1_1387418671154_7021">Wagner a estrutura de efeito espacial é a mesma dentro de cada local.</div><div id="yiv5585344374yui_3_13_0_ym1_6_1387415804491_12" style="color:rgb(0, 0, 0);font-size:16px;font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;background-color:transparent;font-style:normal;" class="yui_3_13_0_ym1_1_1387418671154_7021"><br></div><div id="yiv5585344374yui_3_13_0_ym1_6_1387415804491_12" style="color:rgb(0, 0,
0);font-size:16px;font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;background-color:transparent;font-style:normal;" class="yui_3_13_0_ym1_1_1387418671154_7021">Abraços.</div><div id="yiv5585344374yui_3_13_0_ym1_6_1387415804491_12" style="color:rgb(0, 0, 0);font-size:16px;font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;background-color:transparent;font-style:normal;" class="yui_3_13_0_ym1_1_1387418671154_7021">Cássio Dessotti.</div><div class="yiv5585344374yqt4850901660" id="yiv5585344374yqt03809"><div class="yiv5585344374yahoo_quoted" id="yiv5585344374yui_3_13_0_ym1_6_1387415804491_20" style="display: block;"> <br clear="none"> <br clear="none"> <div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5439" style="font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 12pt;"> <div
class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5440" style="font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 12pt;"> <div dir="ltr"> <font size="2" face="Arial"> Em Quarta-feira, 18 de Dezembro de 2013 13:07, Luiz Roberto Martins Pinto <luizroberto.uesc@gmail.com> escreveu:<br clear="none"> </font> </div> <div class="yiv5585344374y_msg_container"><div id="yiv5585344374"><div><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="yiv5585344374gmail_extra"><br clear="all"><div><div>Luiz Roberto Martins Pinto<br clear="none">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 clear="none"><br clear="none"><a rel="nofollow" shape="rect" ymailto="mailto:luizroberto.uesc@gmail.com" target="_blank" href="mailto:luizroberto.uesc@gmail.com">luizroberto.uesc@gmail.com</a><br clear="none">skype: lrmpinto</div>
<div><a rel="nofollow" shape="rect" target="_blank" href="http://lattes.cnpq.br/2732314327604831">http://lattes.cnpq.br/2732314327604831</a> <br clear="none"><br clear="none"><br clear="none"></div></div>
<br clear="none"><br clear="none"><div class="yiv5585344374gmail_quote">Em 18 de dezembro de 2013 10:37, Wagner Bonat <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:wbonat@gmail.com" target="_blank" href="mailto:wbonat@gmail.com">wbonat@gmail.com</a>></span> escreveu:<br clear="none"><blockquote class="yiv5585344374gmail_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 clear="none"><br clear="none"></div>Você vai assumir que a estrutura do efeito espacial é sempre a mesma dentro dos blocos ?<br clear="none"></div><div class="yiv5585344374gmail_extra"><br clear="none"><br clear="none"><div class="yiv5585344374gmail_quote">
Em 17 de dezembro de 2013 16:00, Cássio Dessotti <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:cassiodessotti@yahoo.com.br" target="_blank" href="mailto:cassiodessotti@yahoo.com.br">cassiodessotti@yahoo.com.br</a>></span> escreveu:<br clear="none"><blockquote class="yiv5585344374gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div class="yiv5585344374h5">
<div><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5444" 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 clear="none"></div>
<div style="display:block;">
<div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5447" style="font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 12pt;"><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5448" style="font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 12pt;">
<div><div><div><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5449" style="font-size: 12pt; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;"><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5450" style="font-style: normal; font-size: 16px; background-color: transparent; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;">
<br clear="none"></div><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5451" 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 class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5454" 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 clear="none">
</span></div><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5456" 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 class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5458" 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 clear="none">
</span></div><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5460" 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 clear="none"></div><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5462" 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 class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5463" style="font-style: normal; font-size: 16px; background-color: transparent; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;"><br clear="none"></div><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5464" 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 class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5465" 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 class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5466" 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 class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5467" style="font-style: normal; font-size: 16px; background-color: transparent; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;">
<br clear="none"></div><div>### Segue o código (com dados fictícios) para exemplificar o meu problema:<span style="background-color:transparent;"> </span></div><div class="yiv5585344374yui_3_13_0_ym1_1_1387415804491_5469" 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 clear="none"></div><div><br clear="none"></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 clear="none"></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 clear="none"></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 clear="none"></div> </div> </div> </div> </div></div><br clear="none"></div></div>_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<span class="yiv5585344374HOEnZb"><font color="#888888"><br clear="none"></font></span></blockquote>
</div><span class="yiv5585344374HOEnZb"><font color="#888888"><br clear="none"><br clear="all"><br clear="none">-- <br clear="none">Wagner Hugo Bonat<br clear="none">
LEG - Laboratório de Estatística e Geoinformação<br clear="none">UFPR - Universidade Federal do Paraná
</font></span></div>
<div class="yiv5585344374yqt0661849653" id="yiv5585344374yqtfd70566"><br clear="none">_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br clear="none"></div></blockquote></div><div class="yiv5585344374yqt0661849653" id="yiv5585344374yqtfd63176"><br clear="none"></div></div></div></div><br clear="none"><br clear="none"></div> </div> </div> </div></div> </div></div></div></div></body></html>