[R-br] RES: Dúvidas

Quarta Agosto 23 21:11:37 -03 2017

dados<-read.table(header = T,text = "
   tempo       X5      X15      X25      X35      X45      X55      X65      X75      X85      X95
1      1 46.58864 49.50402 46.24913 49.81371 49.20052 49.22530 45.96857 50.10608 48.47951 48.07773
2      2 46.58864 47.46737 50.06054 46.99793 48.66818 48.21422 49.24994 50.10608 50.11063 49.95951
3      3 47.40711 45.52298 46.23095 46.20678 45.45665 46.15861 49.75706 48.58758 48.56062 48.10503
4      4 47.40711 44.25048 46.23095 46.20678 45.45665 45.09717       NA 46.49488 46.86848 46.06307
5      5 48.68825 43.44687 45.16022 45.47870 44.19588 44.17779       NA 44.73585 45.36915 44.69637
6      7 50.07952 43.96944 44.83563 45.13914 43.65230 43.38405       NA 44.14179 44.71608 44.10601
7      9 50.07952 43.70181 44.83563 44.81555 43.65230 43.14397       NA 43.60323 43.84888 43.57078
8     13 48.77485 42.97211 44.23223 44.21402 43.39968 42.09673       NA 43.11429 43.09952 43.32199
9     19 48.77485 42.75091 43.43109 43.41536 42.50047 41.73839       NA 42.07306 42.44876 42.23848
10    25 48.77485 42.14447 42.95775 42.72354 42.10862 41.40869       NA 42.26307 42.25076 41.86852
11    31 45.90420 41.95949 41.94801 42.12054 41.58065 41.25354       NA 41.38864 41.70581 41.86852
12    37       NA 42.14447 42.52654 42.12054 41.41875 40.96084       NA 41.38864 41.53891 41.52854
13    43       NA       NA 42.73721 42.12054 41.41875 41.10437       NA 41.38864 41.53891 41.69501
14    49 42.55277       NA 42.73721 41.76025 41.26321 40.82263       NA 41.38864 41.22468 42.43592
15    55 42.55277       NA 42.13249 41.93656 40.83125 40.68944       NA 41.08584 41.07661 41.52854
16    61 42.15649       NA 42.32516 41.76025 40.69776 40.56101       NA 40.80544 40.79688 41.06740
17    67 42.15649       NA 41.94801 41.27291 40.56903 40.31740       NA 40.67288 40.66462 40.78834
18    73 42.15649       NA 41.94801 41.27291 41.11367 40.68944       NA 40.94300 41.07661 41.36870
19    79       NA       NA 41.28262 40.83990 40.32488 40.08997       NA 40.42165 40.41396 40.65639
20    91 41.45925 40.99678 41.28262 40.97877 40.20899 39.98180 40.42935 40.30250 40.29507 40.52912
21   103 41.79326 40.46814 40.71444 40.45257 39.68355 39.58209 40.19456 40.07603 39.96160 40.78834
22   116 41.15110 40.46814 40.58512 40.45257 39.68355 39.58209 40.30994 39.96832 39.85753 40.06216
23   128 41.30213       NA 40.22349 40.33237 39.58806 39.48963 39.87056 39.66513 39.65902 39.75043
24   141 40.73120       NA 40.22349 39.99534 39.49542 39.31272 39.87056 39.66513 39.75672 39.65292
25   153 40.60128       NA 40.11098 39.89021 39.49542 39.31272 39.67126 39.57020 39.56427 39.65292
26   165       NA       NA 40.00213 39.78838 39.23334 39.14571 39.48385 39.47809 39.47233 39.55835
27   177       NA       NA 39.79475 39.59403 39.07067 38.98775 39.30727 39.21745 39.21217 39.37751
28   189       NA       NA 39.79475 39.50122 39.23334 38.98775 39.39427 39.21745 39.29640 39.37751
29   213       NA       NA 39.41677 39.15604 38.84273 38.69607 39.14056 38.97801 38.97315 39.12517
30   237       NA       NA 39.41677 39.23866 38.84273 38.62773 39.06064 38.82886 38.89772 38.89299
31   261       NA       NA 39.16121 38.92145 38.43655 38.31006 38.83347 38.61917 38.61490 38.74819
32   285       NA       NA 38.92622 38.63631 38.37448 38.31006 38.76166 38.48789 38.61490 38.74819
33   309       NA       NA 38.85201 38.70487 38.31387 38.13664 38.55689 38.42461 38.54856 38.61064
34   333       NA       NA 38.70928 38.56944 38.19677 38.08139 38.55689 38.42461 38.35895 38.54440
35   406       NA       NA 38.26209 38.03414 37.68243 37.63373 38.07792 37.91621 38.01718 38.06755
")

library(minpack.lm)
regs<-lapply(dados[,2:ncol(dados)],function(x) {(return(summary(nlsLM(x~a*tempo^b,
          data = dados,start=list(a=-0.6,b=0.5),algorithm = "LM",control = list(maxiter = 500)))))})
regs

coefs<-sapply(1:(ncol(dados)-1),function(x) coef(regs[[x]])[c(1,2)])
coefs

Veja se é algo assim que precisa.

Att., Delcio Bortolanza.

________________________________
De: R-br <r-br-bounces em listas.c3sl.ufpr.br> em nome de João Pedro Domingues via R-br <r-br em listas.c3sl.ufpr.br>
Enviado: quarta-feira, 23 de agosto de 2017 12:31
Para: a lista Brasileira oficial de discussão do programa R.
Assunto: [R-br] RES: Dúvidas

Prezados,

Estou pretendendo realizar uma análise dos parâmetros alpha, betha e gama gerados a partir do HoltWinters.

Porém tenho mais de 300 séries para serem comparadas, alguém saberia um meio de gerar esses parâmetros automaticamente de todas as séries para que não tenha que fazer um por um?

Obrigado pela ajuda.

João Pedro Araujo Domingues

C +55 27 99232-9582<tel:%2B55%2027%2099232-9582>

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170824/501d6fe5/attachment.html>