<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Walmes eu consegui criar o loop no entanto não estou conseguindo definir o critério de parada <font face="Courier New">dif</font>. <div><br></div><div>CMR:<br><div><br></div><div>dia<-c(2,5,7,10,14,19,26,31,34,38,45,52,53,60,65)</div><div><br></div><div>diag<-c(54,50,45,37,35,25,20,16,18,13,8,11,8,4,6)</div><div><br></div><div>d3 <- deriv3(~a*exp(b*x), c("a", "b"), function(x, a, b){ NULL })</div><div><br></div><div>bi<-c(55,-0.02)</div><div><br></div><div>sqresi<-crossprod(diag-c(d3(dia,bi[1],bi[2]))) </div><div><br></div><div>i <- 1</div><div><br></div><div><span style="font-size: 10pt;">while(i < 10){</span></div><div><br></div><div>est <- d3(dia,bi[1],bi[2])</div><div><br></div><div>fx <- c(est)</div><div><br></div><div>d <- attr(est, "gradient")</div><div><br></div><div>sqresf <- crossprod(diag-fx)</div><div><br></div><div>bf <- bi + (solve(t(d)%*%d)%*%(t(d)%*%(diag-fx)))</div><div><br></div><div>bi <- c(bf)</div><div><br></div><div><span style="font-size: 10pt;">sqresi <- sqresf</span></div><div><br></div><div>cat(paste(formatC(c(sqresf, bi), digits=6, format="f"), collapse="\t"), "\n")</div><div><br></div><div>i <- i + 1</div><div><br></div><div>}</div><div><br></div><div>E a proposito dei uma lida no post que você fala a respeito da anova para modelos não lineares. E fiquei com uma dúvida, seria possível testar o efeito dos betas como nos modelos</div><div><br></div><div>lineares?</div><div><br></div><div>Agradeço desde já pela ajuda.</div><div><br></div><div>Att.</div><div><br></div><div>Tiago.<span style="font-size: 10pt;"> </span></div><div><br><div><div id="SkyDrivePlaceholder"></div><hr id="stopSpelling">Date: Tue, 22 Jan 2013 16:49:16 -0200<br>From: walmeszeviani@gmail.com<br>To: r-br@listas.c3sl.ufpr.br<br>Subject: Re: [R-br] MÉTODO DE GAUSS-NEWTON<br><br><font face="trebuchet ms,sans-serif">Isso não vai ser problema não, pode haver zeros em alguns elementos da hessiana. Vários modelos se comportam como esse seu.<br><br>À disposição.<br>Walmes.<br><br clear="all"></font><div>
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br">walmes@ufpr.br</a><br>skype: walmeszeviani<br style="font-family:trebuchet ms,sans-serif"></span><span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/~walmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br></div><br>
<br>_______________________________________________
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.</div></div></div> </div></body>
</html>