Pessoal mais uma duvida, eu fiz um ajuste com o nls, mas eu quero saber como utilizar a e b ou seja se tem algum comando para que eu possa utilizar os valores que foram estimadas pela equação Segue o CMR que estou trabalhando # Planta 1 planta1=foto[1:90,-3:-6] names(planta1) [1] "especie" "esp" "PAR" "Tleaf" "W" "Pn" "E" [8] "C" "int_CO2" colnames(planta1)<-c("especiea", "espa", "PARa","Tleafa", "Wa", "Pna", "Ea", + "Ca", "int_CO2a") # Obtendo PAR planta1.luz =tapply(planta1$PARa, INDEX=planta1$PARa, FUN=mean) planta1.luz 0 25 50 75 100 500 1000 1500 2000 0 25 50 75 100 500 1000 1500 2000 planta1.me.Pna=tapply(planta1$Pna, INDEX=planta1$PARa, FUN=mean) planta1.me.Pna 0 25 50 75 100 500 1000 1500 2000 -0.428 0.144 0.411 1.275 2.052 4.349 5.082 4.078 3.794 planta1.ajuste=nls((planta1.me.Pna~(a+0.428)*(1-exp(-b*planta1.luz/(a+0.428)))-0.428), + start=list(a=200, b=0.5)) planta1.ajuste Nonlinear regression model model: planta1.me.Pna ~ (a + 0.428) * (1 - exp(-b * planta1.luz/(a + 0.428))) - 0.428 data: parent.frame() a b 4.37564 0.02943 residual sum-of-squares: 1.258 Number of iterations to convergence: 5 Achieved convergence tolerance: 3.175e-06 max=c(planta1.ajuste$a) = ????? alfa=c(plnat1.ajuste$b) = ??????? Desde já muito obrigado. -- Bach. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 http://lattes.cnpq.br/3684553005529634 --- 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.*** Salva un Árbol.* * *** * PIENSA**** EN EL MUNDO PARA TUS HIJOS** *