
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.