<div dir="ltr"><div>O que seria essa matrix de projecao? Como eu mudaria na funcao log.verossimilhancao, abaixo:  <br></div><div><br></div><div>montaSigma <- function(s2, t2, phi, Umat){<br>    Sigma <- as.matrix(s2 * exp(-Umat/phi))<br>    Sigma <- ifelse(Umat==0, diag(Sigma)+t2, Sigma)<br>    return(Sigma)<br>}<br><br></div><div><br></div><div>ll.geo <- function(s2, t2, phi, modelo, Umat, dados, logpars = F) {<br>    if (logpars) {<br>        s2 <- exp(s2)<br>        t2 <- exp(t2)<br>        phi <- exp(phi)<br>    }<br>    mf <- model.frame(modelo, dados)<br>    y <- model.response(mf)<br>    D <- model.matrix(modelo, mf)<br>    Sigma <- montaSigma(s2 = s2, t2 = t2, phi = phi, Umat = Umat)<br>    R <- chol(Sigma)<br>    invRD <- backsolve(R, D, transpose = TRUE)<br>    invRy <- backsolve(R, y, transpose = TRUE)<br>    bhat <- solve(crossprod(invRD), crossprod(invRD, invRy))<br>    invRe <- invRy - invRD %*% bhat<br>    nll <- drop(length(y) * log(2 * pi)/2 + sum(log(diag(R))) + crossprod(invRe)/2)<br>    return(nll)<br></div><div><br></div><div>Abraco<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 31 de julho de 2018 16:22, Elias T Krainski via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div text="#000000" bgcolor="#FFFFFF">
    <p>Escreva <br>
    </p>
    <p> y = X\beta + As + erro</p>
    <p>em que A e' uma matriz de projecao do efeito espacial s. Disso
      temos que <br>
    </p>
    <p>Cov(Y) = Cov(Ax + erro) = A'Cov(s)A + Cov(Erro)<br>
    </p>
    Elias. <br><div><div class="h5">
    <br>
    <div class="m_1270534479705122277moz-cite-prefix">On 31/07/2018 11:12, Wagner Wolff via
      R-br wrote:<br>
    </div>
    </div></div><blockquote type="cite"><div><div class="h5">
      
      <div dir="ltr">
        <div>Olá pessoal da lista!</div>
        <div><br>
        </div>
        <div>Estou tentando estimar os parâmetros de um modelo
          geoestatistico no qual a matriz de covariancia é obtida por
          meio de uma matriz de distância entre linhas. Ou seja, linhas
          interceptadas tem distância zero, assim não somente a diagonal
          da matriz de distância é zero. Isso acaba acarretando em
          problema na hora de decompor a matriz de covariância usando
          Cholesky. Alguém tem alguma solucão ou dica? Já tentei algumas
          coisas como aproximacões para ser positiva definida, entre
          outras, mas na hora de estimar os parâmetros nada converge.</div>
        <div><br>
        </div>
        <div>Abraco <br>
        </div>
      </div>
      <br>
      <fieldset class="m_1270534479705122277mimeAttachmentHeader"></fieldset>
      <br>
      </div></div><pre>______________________________<wbr>_________________
R-br mailing list
<a class="m_1270534479705122277moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a class="m_1270534479705122277moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="m_1270534479705122277moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-<wbr>guia</a>) e forne�a c�digo m�nimo reproduz�vel.</pre>
    </blockquote>
    <br>
  </div>

<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-<wbr>guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="font-size:medium"><b><i>Wagner Wolff, </i></b></span><i><b>PhD</b></i><br>"<b>Luiz de Queiroz</b><b><span>"</span> College of Agriculture,</b><br>University of São Paulo<br>Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil<br>Phone:  <a href="tel:+55%2019%2098238-5582" value="+5519982385582" target="_blank">+55 19 982385582</a>  <br><span><span><a href="http://orcid.org/0000-0003-3426-308X" target="_blank">http://orcid.org/0000-0003-3426-308X</a><br><a href="https://github.com/wwolff7" target="_blank">https://github.com/wwolff7</a><br><a href="http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1" target="_blank">http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1</a></span></span></div></div>
</div>