From 9b7f6e609b21902e817b1156d2f9d8b50c66d59d Mon Sep 17 00:00:00 2001
From: Walmes Zeviani <walmes@ufpr.br>
Date: Fri, 20 Oct 2017 18:53:50 -0200
Subject: [PATCH] usa o pacote boot.
---
scripts/ce089-06.R | 27 +++++++++++++++++++++++++++
1 file changed, 27 insertions(+)
diff --git a/scripts/ce089-06.R b/scripts/ce089-06.R
index 003ffea..39742e4 100644
--- a/scripts/ce089-06.R
+++ b/scripts/ce089-06.R
@@ -150,3 +150,30 @@ with(as.list(coef(n1)), {
})
#-----------------------------------------------------------------------
+# Usando o pacote boot.
+
+library(boot)
+
+help(boot, help_type = "html")
+
+str(turk0)
+
+# Criar uma função com dois argumentos: o data.frame original e um vetor
+# que vai representar o índice das linhas usado para reamostrar dos
+# dados.
+fitmodel <- function(dataset, index) {
+ n0 <- nls(Gain ~ Int + Ass * (1 - 2^(-A/Mei)),
+ data = dataset[index, ],
+ start = list(Int = 625, Ass = 180, Mei = 0.1))
+ c(coef(n0), s2 = deviance(n0)/df.residual(n0))
+}
+
+b0 <- boot(turk0,
+ statistic = fitmodel,
+ R = 999)
+
+pairs(b0)
+
+# TODO:
+# usar o pacote boot.
+# reproduzir os resultados do boot() e boot.ci().
--
GitLab