[R-br] adicionar curva de ajuste no grafico

[Gilson Sanchez]

Pessoal mais uma vez com minhas coisas.

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.

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?

Como eu posso adcioar esse ajuste no meu grafico... da pra fazer isso?

Desde já agradeço

dados
http://www.datafilehost.com/download-061180ff.html

CMR:

require(XLConnect)
especie<-loadWorkbook("especie1.xls", create= FALSE)
especie1<-readWorksheet(especie, sheet="especie.1")
especie1
attach(especie1)
#-------------------------------------------------------------------------
#Média
especie.me <- tapply(especie1$A.Det, INDEX = especie1$PAR, FUN = mean)
especie.me
plot(especie.me)
#Luz
luz<-c(0,25,50,75,100,500,1000,1500,2000)
#Desvio padrão
especie.sd <- tapply(especie1$A.Det, INDEX = especie1$PAR, FUN = sd)
especie.sd
#Amplitude total
especie.le <- tapply(especie1$A.Det,INDEX = especie1$PAR, FUN = length)
especie.le
especie.se <- especie.sd / sqrt(especie.le)
especie.se
plot(especie.me~luz)
stripchart(especie.me~luz, ylim=c(-1,5),
           vert = TRUE, pch = 1, method = "jitter",
           jit = 0.05, xlab = "PAR", ylab = "Fotossintese", data=especie1)
points(1:9, especie.me, pch = 16, col=2,cex = 1.0)
arrows(1:9, especie.me, 1:9, especie.me + especie.se, lwd = 1.5,
       angle = 90, length = 0.1, col=2, add=TRUE)
arrows(1:9, especie.me,1:9, especie.me - especie.se, lwd = 1.5,
       angle = 90, length = 0.1, col=2)
abline(h=0, lty=4)



MSc. Gilson Sánchez Chia

Laboratório de Fisiologia Vegetal

Embrapa Amazônia Ocidental

Fone: (92) 3303-7841
---
Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente.
Please consider the environment before printing this email.
 Renew, Reduce and Recycle. The Planet Knows.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110616/d0d518bb/attachment.html>