[R-br] adicionar curva de ajuste no grafico

Paulo Justiniano paulojus em leg.ufpr.br
Quinta Junho 16 11:17:43 BRT 2011


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.
> 
> 
>


Mais detalhes sobre a lista de discussão R-br