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 

Desde já muito obrigado.

Bach. Agr.Gilson Sánchez Chia, MSc.

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