Skip to main content
1 de 3
M Aurélio
  • 222
  • 1
  • 4

A forma como eu gosto de fazer isso é criando uma lista de modelos no próprio objeto para fácil controle e organização com o comando nest do pacote tidyr.

library(dplyr)
library(tidyr)
library(purrr)

urlfile <- 'https://raw.githubusercontent.com/StephannyEgito/RNA/main/Amostra_Git_csv.csv'
dt <-read.csv2(urlfile)

modelo_padrao <- function(dados)
  nls(MR ~ k1*(confinante^k2)*(desvio^k3), 
      data = dados, 
      start= list(k1=2000, k2=0.4, k3=-0.2))

dt_new <- 
dt %>% 
  nest(data = -amo) %>% 
  mutate(teste = map(data, modelo_padrao),
         tidy_mod = map(teste, broom::tidy)) 

Dessa forma podemos acessar tanto os teste individualmente em cada elemento da lista

head(dt_new$teste, 3)  
#> [[1]]
#> Nonlinear regression model
#>   model: MR ~ k1 * (confinante^k2) * (desvio^k3)
#>    data: dados
#>        k1        k2        k3 
#> 1525.8538    0.5458   -0.1743 
#>  residual sum-of-squares: 32784
#> 
#> Number of iterations to convergence: 6 
#> Achieved convergence tolerance: 7.96e-06
#> 
#> [[2]]
#> Nonlinear regression model
#>   model: MR ~ k1 * (confinante^k2) * (desvio^k3)
#>    data: dados
#>        k1        k2        k3 
#> 8.115e+03 7.308e-01 2.655e-02 
#>  residual sum-of-squares: 408876
#> 
#> Number of iterations to convergence: 7 
#> Achieved convergence tolerance: 1.048e-06
#> 
#> [[3]]
#> Nonlinear regression model
#>   model: MR ~ k1 * (confinante^k2) * (desvio^k3)
#>    data: dados
#>        k1        k2        k3 
#> 8863.3252    0.6966    0.0106 
#>  residual sum-of-squares: 64422
#> 
#> Number of iterations to convergence: 7 
#> Achieved convergence tolerance: 4.794e-07

Quanto em um formato limpo para comparação direta através de broom::tidy e tidyr::unnest

unnest(dt_new, tidy_mod)
#> # A tibble: 30 × 8
#>      amo data               teste  term   estimate std.error statistic  p.value
#>    <int> <list>             <list> <chr>     <dbl>     <dbl>     <dbl>    <dbl>
#>  1     1 <tibble [18 × 18]> <nls>  k1    1526.      156.         9.76  6.84e- 8
#>  2     1 <tibble [18 × 18]> <nls>  k2       0.546     0.0605     9.02  1.90e- 7
#>  3     1 <tibble [18 × 18]> <nls>  k3      -0.174     0.0471    -3.70  2.15e- 3
#>  4     2 <tibble [18 × 18]> <nls>  k1    8115.     1431.         5.67  4.45e- 5
#>  5     2 <tibble [18 × 18]> <nls>  k2       0.731     0.102      7.15  3.31e- 6
#>  6     2 <tibble [18 × 18]> <nls>  k3       0.0265    0.0732     0.363 7.22e- 1
#>  7     3 <tibble [18 × 18]> <nls>  k1    8863.      500.        17.7   1.78e-11
#>  8     3 <tibble [18 × 18]> <nls>  k2       0.697     0.0328    21.3   1.30e-12
#>  9     3 <tibble [18 × 18]> <nls>  k3       0.0106    0.0238     0.446 6.62e- 1
#> 10     4 <tibble [18 × 18]> <nls>  k1    1313.       40.8       32.2   2.90e-15
#> # … with 20 more rows

Created on 2023-02-07 with reprex v2.0.2

M Aurélio
  • 222
  • 1
  • 4