este é um modelo nao linear e vc pode usar o funcao nls() para plota a curva é só etimar os parametros, e plotar com os parametros estimados como faria o grafico de qq funcao On Thu, 16 Jun 2011, Gilson Sanchez wrote:

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.