<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt"><div><br></div><div 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;">Elias, Walmes e Wagner, antes de mais nada muito obrigado <span style="font-size: 12pt;">pela ajuda.</span></div><div><br></div><div>Elias muito bem observado o problema de identificabilidade <span style="font-size: 12pt;">acarretado pelo confundimento de longitude com blocos, já que </span><span style="font-size: 12pt;">os blocos "estão em linha". </span></div><div><br></div><div>Estudarei o material sugerido por você, Wagner.</div><div><br></div><div>Walmes muito obrigado pelas linhas de comando que resolveram <span style="font-size: 12pt;">o modelo 3, com
meus dados </span><span style="font-size: 12pt;">reais, já que no CMR coloquei apenas alguns dados para </span><span style="font-size: 12pt;">possibilitar a reprodução. Disponho de um grupo de experimentos </span><span style="font-size: 12pt;">em blocos ao acaso com 6 tratamentos, 5 blocos e 4 locais. A </span><span style="font-size: 12pt;">variável resposta é altura de plantas de cana-de-açúcar. Os blocos </span><span style="font-size: 12pt;">estão "em linha" o que acaba confundindo "longitude" com "blocos" </span><span style="font-size: 12pt;">mas sua sugestão para meu M3 (seu M4) resolveu muito bem. Pesquisarei esta </span><span style="font-size: 12pt;">literatura que você citou a respeito da questão de blocos ser </span><span style="font-size: 12pt;">considerado fixo ou aleatório.</span></div><div><br></div><div>Walmes, construí diversos modelos alternando blocos (fixo, aleatório <span
style="font-size: 12pt;">ou sem efeito), variâncias de locais (a mesma ou diferentes) e </span><span style="font-size: 12pt;">função espacial para a matriz R (exponencial, esférica e gaussiana), </span><span style="font-size: 12pt;">e optei pelo modelo exponencial sem efeito de blocos e com </span><span style="font-size: 12pt;">homogeneidade de variâncias como o "melhor" observando os valores de AIC.</span></div><div><br></div><div>Você acha que o critério de AIC (ou BIC - apesar de nem sempre apontarem <span style="font-size: 12pt;">o mesmo modelo) é suficiente para a escolha de um modelo? </span><span style="font-size: 12pt;">Seria necessário o uso do critério de validação cruzada neste caso?</span></div><div><br></div><div>A partir do modelo selecionado (análise conjunta), posso extrair os resíduos condicionais <span style="font-size: 12pt;">estudentizados, separá-los segundo locais, para em cada local
construir </span><span style="font-size: 12pt;">semivariogramas destes resíduos, verificando a "força" da dependência, e </span><span style="font-size: 12pt;">realizar krigagem em cada área (trabalhando com os resíduos para "limpar" </span><span style="font-size: 12pt;">os dados de seus efeitos de tratamentos e demais efeitos). </span></div><div><br></div><div>Muito grato pela força,</div><div>Cássio Dessotti.</div><div><br></div><div class="yahoo_quoted" style="display: block;"> <br> <br> <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 dir="ltr"> <font size="2" face="Arial"> Em Quinta-feira, 19 de Dezembro de 2013 18:35, Wagner Bonat <wbonat@gmail.com> escreveu:<br> </font> </div> <div class="y_msg_container"><div
id="yiv6914124313"><div><div dir="ltr"><br clear="none"><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 clear="none">
</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 clear="none">
<br clear="none"></div><div><a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/mcie">www.leg.ufpr.br/mcie</a><br clear="none"></div><div><br clear="none"></div><div><br clear="none"></div><br clear="none"></div><div class="yiv6914124313gmail_extra"><br clear="none"><br clear="none"><div class="yiv6914124313gmail_quote">Em 18 de dezembro de 2013 13:07, Luiz Roberto Martins Pinto <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:luizroberto.uesc@gmail.com" target="_blank" href="mailto:luizroberto.uesc@gmail.com">luizroberto.uesc@gmail.com</a>></span> escreveu:<br clear="none">
<blockquote class="yiv6914124313gmail_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="yiv6914124313gmail_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="yiv6914124313gmail_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:<div><div class="yiv6914124313h5"><br clear="none"><blockquote class="yiv6914124313gmail_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="yiv6914124313gmail_extra"><br clear="none"><br clear="none"><div class="yiv6914124313gmail_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="yiv6914124313gmail_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 clear="none"></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 clear="none"></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 clear="none">
</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 clear="none">
</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 clear="none"></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 clear="none"></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 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 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><font color="#888888"><br clear="none"></font></span></blockquote>
</div><span><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>
<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.<div class="yiv6914124313yqt3416196322" id="yiv6914124313yqtfd20398"><br clear="none"></div></blockquote></div></div></div><div class="yiv6914124313yqt3416196322" id="yiv6914124313yqtfd54557"><br clear="none"></div></div><div class="yiv6914124313yqt3416196322" id="yiv6914124313yqtfd20942">
<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="yiv6914124313yqt3416196322" id="yiv6914124313yqtfd41543"><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á
</div></div></div></div><br><div class="yqt3416196322" id="yqtfd11135">_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" 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 clear="none">Leia o guia de postagem (<a shape="rect" 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.</div><br><br></div> </div> </div> </div> </div></body></html>