<div>Pessoal mais uma vez com minhas coisas.</div><div> </div><div>O seguinte estou tentando fazer um grafico de fotossintese com pontos que representam a média e desvio padrão em diferente intesidades de luz que variam de 0 até 2000.</div>
<div> </div><div>Mas eu quero adicionar uma curva de ajuste [por exemplo dada pela equação A.Det~a+b/(1+exp(-(PAR-c)/d))]. Como estimo os valores de a,b,c e d; de esta equação?</div>
<div> </div><div>Como eu posso adcioar esse ajuste no meu grafico... da pra fazer isso?</div><div> </div><div>Desde já agradeço</div><div> </div><div>dados</div><div><a href="http://www.datafilehost.com/download-061180ff.html" target="_blank">http://www.datafilehost.com/download-061180ff.html</a></div>

<div> </div><div>CMR:</div><div><br>require(XLConnect)<br>especie<-loadWorkbook("especie1.xls", create= FALSE)<br>especie1<-readWorksheet(especie, sheet="especie.1")<br>especie1<br>attach(especie1)</div>

<div>#-------------------------------------------------------------------------<br>#Média<br><a href="http://especie.me" target="_blank">especie.me</a> <- tapply(especie1$A.Det, INDEX = especie1$PAR, FUN = mean)<br><a href="http://especie.me" target="_blank">especie.me</a><br>

plot(<a href="http://especie.me" target="_blank">especie.me</a>)<br>#Luz<br>luz<-c(0,25,50,75,100,500,1000,1500,2000)</div><div>#Desvio padrão<br><a href="http://especie.sd" target="_blank">especie.sd</a> <- tapply(especie1$A.Det, INDEX = especie1$PAR, FUN = sd)<br>

<a href="http://especie.sd" target="_blank">especie.sd</a></div><div>#Amplitude total<br>especie.le <- tapply(especie1$A.Det,INDEX = especie1$PAR, FUN = length)<br>especie.le</div><div><a href="http://especie.se" target="_blank">especie.se</a> <- <a href="http://especie.sd" target="_blank">especie.sd</a> / sqrt(especie.le)<br>

<a href="http://especie.se" target="_blank">especie.se</a></div><div>plot(<a href="http://especie.me" target="_blank">especie.me</a>~luz)</div><div>stripchart(<a href="http://especie.me" target="_blank">especie.me</a>~luz, ylim=c(-1,5),<br>
           vert = TRUE, pch = 1, method = "jitter",<br>
           jit = 0.05, xlab = "PAR", ylab = "Fotossintese", data=especie1)<br>points(1:9, <a href="http://especie.me" target="_blank">especie.me</a>, pch = 16, col=2,cex = 1.0)<br>arrows(1:9, <a href="http://especie.me" target="_blank">especie.me</a>, 1:9, <a href="http://especie.me" target="_blank">especie.me</a> + <a href="http://especie.se" target="_blank">especie.se</a>, lwd = 1.5,<br>

       angle = 90, length = 0.1, col=2, add=TRUE)<br>arrows(1:9, <a href="http://especie.me" target="_blank">especie.me</a>,1:9, <a href="http://especie.me" target="_blank">especie.me</a> - <a href="http://especie.se" target="_blank">especie.se</a>, lwd = 1.5,<br>
       angle = 90, length = 0.1, col=2)<br>
abline(h=0, lty=4)</div><div> </div><div> </div><div> </div><div>MSc. Gilson Sánchez Chia<br></div>
<p style="margin: 0cm 0cm 0pt;"><font size="3" face="Calibri">Laboratório de Fisiologia Vegetal</font></p>
<p style="margin: 0cm 0cm 0pt;"><font size="3" face="Calibri">Embrapa Amazônia Ocidental</font></p>
<p style="margin: 0cm 0cm 0pt;"><font size="3" face="Calibri">Fone: (92) 3303-7841</font></p>
<div>---<br><font color="#000099">Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente.</font></div>
<div><span lang="en-us"><font color="#000099" face="arial,helvetica,sans-serif">Please consider the environment before printing this email.</font></span></div>
<div><span lang="en-us">
<div><span lang="en-us"><font color="#000099"><font face="arial,helvetica,sans-serif">Renew, Reduce and Recycle. The Planet Knows.</font></font></span></div></span></div><br>