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

 

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

---
Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente.