Provavelmente seu Nelder-Mead não convergiu ou convergiu para algum minimo
local. Vc tentou usar valores iniciais diferentes? Pra sair da dúvida pegue
as estimativas obtivas pelo newton-raphson e coloque como inicial do
Nelder-Mead e veja o que acontece. Se o Nelder-Mead não convergir para o
mesmo
que tal um bom chute inicial?
hist(y, freq=FALSE)
plot(function(x) dweibull(x, 1, 1e6), 0, 2e6, add=T, col=1)
plot(function(x) dweibull(x, 1.2, 1e6), 0, 2e6, add=T, col=2)
essa aplicacao do ks.test() parece que nao faz muito sentido aqui...
___
seria bom mudar a escala dos dados para evitar problemas numericos
### estima para y/100
o - optim(c(1.2, 1), function(p)
-sum(dweibull(y/1e6, p[1], p[2], log=TRUE)),
method='BFGS', hessian=TRUE)
### estima para 1000 reamostras (bootstrap)
oo - t(replicate(1000, optim(o$par,
Oi Andre,
No link abaixo tem alguma coisa de passo-a-passo para isso
http://www.leg.ufpr.br/doku.php/pessoais:elias:cursos:uem2012
mas nao esta atualizado...
Sugiro alguns links:
dados tabulados IBGE: http://www.sidra.ibge.gov.br/
microdados censo 2012:
Boa recomedacao!
On 07/05/15 20:16, Wagner Bonat wrote:
perfil de verossimilhança.
--
In God we trust, all others bring data.
– William Edwards Deming (1900-1993)
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
No IBGE, acho que o ftp realmente é a forma mais prática. E ainda tem a
vantagem de programar o R para fazer o download direto (vide www.asdfree.com
)
O Datasus também tem um ftp, mas é menos organizado. O caminho geral é
ftp://ftp.datasus.gov.br/, mas acho que as informações mais buscadas se
Olá Éder conheço o trabalho, na verdade me baseei nele para usar o mtsdi,
pretendo fazer um parecido, só que com dados de vazão
Obrigado
Em 7 de maio de 2015 12:46, Éder Comunello comunello.e...@gmail.com
escreveu:
Wagner,
Nesse trabalho, do qual faz parte o colega de lista Elias Krainski, dá
Tente padronizar o índice dos objetos, convertendo do formato Date pra
POSIX, algo como:
index(objeto) - as.POSIXct(index(objeto), America/Manaus)
Utilize o timezone adequado pro seu caso...
Éder Comunello c comunello.e...@gmail.comomunello.e...@gmail.com
Dourados, MS - [22 16.5'S, 54 49'W]
Depende muito da situação, valores iniciais é sempre uma coisa critica. O
que eu faço é iniciar o algoritmo com diferentes valores e olho se o
algoritmo acha sempre o mesmo maximo ou minimo. As vezes tento diferentes
algoritmos e pra ter certeza que eu achei o ponto certo faço um perfil de
Aline, tentou com a opção each=21?
rep (seq ( chron ( times = c(00:00:00) ), chron ( times = c(00:10:00)
),
chron ( times = c(00:00:15))), each=21 )
Éder Comunello c comunello.e...@gmail.comomunello.e...@gmail.com
Dourados, MS - [22 16.5'S, 54 49'W]
Em 7 de maio de 2015 13:28,
Oi Éder,
A opção each = 21 não funcionou, mas assa opção usando a função strptime deu
certo.
Grata.
Date: Thu, 7 May 2015 16:31:25 -0400
From: comunello.e...@gmail.com
To: r-br@listas.c3sl.ufpr.br
Subject: Re: [R-br] Sequência e repetição de série temporal (minutos)
Outra opção:
t0 -
Olá, os microdados da PNAD podes encontrar no site abaixo. O IBGE somente
disponibiliza os dados a partir de 2001, os dados anteriores podem ser
adquiridos na loja virtual.
http://www.ibge.gov.br/home/estatistica/populacao/trabalhoerendimento/pnad2011/microdados.shtm
Att.,
Luciane Pilotto
Acho que o jeito mais fácil de pegar os microdados do IBGE é pelo
ftp://ftp.ibge.gov.br/ , é bem organizado. E para as malhas digitais tem o
ftp://geoftp.ibge.gov.br/. Tem também os microdados de:
- INEP/Educação: http://portal.inep.gov.br/basica-levantamentos-acessar
/
Bom dia, André!
Tentou instalar como administrador?
[image: Imagem inline 1]
Éder Comunello c comunello.e...@gmail.comomunello.e...@gmail.com
Dourados, MS - [22 16.5'S, 54 49'W]
2015-05-06 22:17 GMT-04:00 Andre Oliveira andreolso...@yahoo.com.br:
Pessoal,
não uso Windows, mas hoje precisei
Edilson, veja se o exemplo ajuda...,
require(xts)
data(sample_matrix)
xts1 - as.xts(sample_matrix[1:30,])
xts2 - as.xts(sample_matrix[15:45,])
merge.xts(xts1, xts2, join=inner)
Éder Comunello c comunello.e...@gmail.comomunello.e...@gmail.com
Dourados, MS - [22 16.5'S, 54 49'W]
Em 6 de maio de
A questão está relacionada com a permissão de escrita na pasta program
files do windows. A partir do windows 7 apenas quando esta rodando os
programas como administrador você pode escrever nestas pastas.
A solução que acho mais confortável é criar uma pasta c:/Programas/ e
instalar R, Rstudio,
Outra opção:
t0 - strptime(00:00:00, %H:%M:%S)
t1 - strptime(00:10:00, %H:%M:%S)
times0 - format(seq(t0, t1, by=15 sec), %H:%M:%S)
rep(times0, each=21)
Éder Comunello c comunello.e...@gmail.comomunello.e...@gmail.com
Dourados, MS - [22 16.5'S, 54 49'W]
Em 7 de maio de 2015 16:30, Éder Comunello
Pessoal,gostaria de deixa aqui meus agradecimentos nas sugestões do grupo
sobre banco de dados públicos e no problema do Rstudio não citarei ninguém para
não cometer injustiça de esquecimento. No Rstudio de fato era permissão que o
aluno da escola adicionou permissão daquela pasta do Rstudio.
Wagner,
Nesse trabalho, do qual faz parte o colega de lista Elias Krainski, dá pra
ter uma ideia de métodos e pacotes a utilizar.
http://www.sbagro.org.br/anais_congresso_2011/cba2011/trabalhos/04/cba04_382_682.pdf
Éder Comunello c comunello.e...@gmail.comomunello.e...@gmail.com
Dourados, MS -
No IBGE eu costumo baixar os microdados a partir da página da pesquisa
em questão.
O Datasus oferece dados administrativos de vários sistemas de informação
desconexos. Em geral, o que você vai ter acesso é a dados agregados e/ou
com número reduzido de variáveis. Visite
Éder,
Agradeço a sua sugestão, mas não resolveu.
a função str() retorna para cada arquivo o seguinte:
str(Shanghai)
An ‘xts’ object on 2011-01-04/2014-12-31 containing:
Data: num [1:970, 1] 2853 2839 2824 2839 2792 ...
- attr(*, dimnames)=List of 2
..$ : NULL
..$ : chr Shanghai
Indexed by
Entendi, realmente não deu certo essa minha ideia, vou proceder com o
mtsdi.
Obrigado pela ajuda Éder. Abraço
Em 5 de maio de 2015 16:53, Éder Comunello comunello.e...@gmail.com
escreveu:
Wagner, boa tarde!
Eu não recomendo usar splinefun() porque a calibração do modelo feita na
execução de
Olá Wagner,
segui o seu conselho e coloquei para o caso Nelder-Mead as estimativas
iniciais como sendo as estimativas encontradas no caso de Newton-Raphson, o
valor foi realmente quase em cima! Entaõ o problema realmente era esse, a
seleção dos pontos iniciais. Fica então uma dúvida: Quando tenho
Senhores, olá novamente!
Na instalação onde gerei o printscreen abaixo (Windows 8), não tenho os
direitos de Administrador. Pedi ao setor de TI pra instalar o R e RStudio
com a opção Run as Adminstrator... e feito isso, tenho utilizado
normalmente como um usuário padrão. A diferença é que os
Éder esta também é uma possibilidade. Mas prefiro ter uma pasta com
permissão de leitura e escrita para todos os usuário fora. Fica mais
organizado e se precisar instalar algum outro programa fica MUITO mais
fácil. Hoje posso atualizar o R sem depender da TI.
Em 7 de maio de 2015 12:37, Éder
Gostei da ideia! Vou ver se consigo adotar... Valeu pela dica!
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo
E valeu por essa ftp do IBGE!
Elias
On 07/05/15 13:24, Rodrigo Coster wrote:
o jeito mais fácil de pegar os microdados do IBGE é pelo
ftp://ftp.ibge.gov.br/
--
In God we trust, all others bring data.
– William Edwards Deming (1900-1993)
___
R-br
Boa tarde prezados colegas,
Eu preciso criar um fator dena série temporal entre 0 a 10 minutos, com
intervalo de 15 segundos e repetir essa sequência 21 vezes em ordem crescente.
Para isso usei o seguinte script:
library ( chron )
sort (rep (seq ( chron ( times = c(00:00:00) ), chron ( times =
28 matches
Mail list logo