[R-br] ajuste com nls

[Gilson Sanchez]

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** *
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120130/df130714/attachment.html>