Maio 2015 - R-br - DINF UFPR

Calibração de dados preenchidos
by Wagner Wolff 07 Mai '15

07 Mai '15

Olá pessoal da lista, utilizei um pacote chamado mtsdi para preenchimento de dados em séries temporais. Gostaria de saber como posso proceder para fazer a calibração destes preenchimentos, por exemplo, eu tiro da minha série um período amostrado e com dados vizinhos correlacionados eu faço o preenchimento deste período assim eu comparo o observado com o preenchido, até aí tudo bem, mas eu gostaria de fazer uma calibração neste preenchimento e assim utilizá-la para os dados que realmente serão preenchidos. Segue um exemplo hipotético do uso do pacote. library(mtsdi) library(hydroGOF) x0 <- rnorm(500) x1 <- rnorm(500) x2 <- rnorm(500) x3 <- rnorm(500) x4 <- c(x3[1:299],rep(NA,tim=400-299),x3[401:500]) x <- data.frame(x0,x1,x2,x4) f <- ~x0+x1+x2+x4 i <- mnimput(f,x,eps=1e-3,ts=T,log=F, method="spline",sp.control=list(df=c(7,7,7,7))) imput <- data.frame(predict(i)) (pred <- imput$x4[300:400]) (obs <- x3[300:400]) ggof(pred,obs)

2 8

Banco de dados publicos
by Andre Oliveira 07 Mai '15

07 Mai '15

Prezados boa tarde, estou a procurar de dados em geral do IBGE, DATASUS entre outros órgãos. Estou com dificuldade de encontrar estes materiais mesmo sabendo que o IBGE e DATASUS disponibilizam, e me perco dentro do site. Alguém teria algum material que explica como obter dados destas instituições? Pode ser outras instituições! Onde posso obter malhas digitais, pois me compliquei também dentro do site do IBGE! Caso alguém tenha algum material resumindo sobre como obter banco de dados públicos e puder compartilhar ficarei agradecido. obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

6 7

Estimativas de parâmetros de uma Weibull utilizando duas metodologias: Newton-Raphson "vs" Nelder Mead (estimativas muito diferentes) o que pode estar acontecendo??
by Romero Luiz M. Sales Filho 07 Mai '15

07 Mai '15

Pessoal, boa noite! Tenho um conjunto de dados que provavelmente segue uma distribuição Weibull. Os dados são os seguintes: y<-c(295000, 869000 , 869900 , 1573335 , 151400 , 152000 , 183700 , 218000 , 30200 , 45100 , 46900 , 47300 ) Efetuei a estimação dos parâmetros de duas formas, uma utilizando as expressões analíticas onde tive que usar o método de Newton-Raphson... E outra utilizando a função optim em que selecionei o método de Nelder Mead... Os dois algoritmos seguem abaixo: ###### Newton Paphson Method ########## n=length(y) x=vector() ##### Função que quero obter a raiz (estimador de lambda) f <- quote((sum((y^delta)*log(y))/sum(y^delta))-(1/delta)-(1/n)*sum(log(y))) fx0 <- function(delta){ eval(f) } ##### Primeira derivada da função (estimador de lambda) f1 <-quote((1/delta^2)-((sum((log(y))*(y^delta)))^2/(sum(y^delta))^2)+(sum((log(y))^2*(y^delta))/sum(y^delta))) fx1 <- function(delta){ eval(f1) } ### Obtendo a estimativa para delta i <- 1 # valor inicial para o passo dif <- 1 # valor inical para a diferença entre duas avaliações tol <- 10^-6 # diferência mínima entre duas avaliações (tolerância) passo <- 1000 # número máximo de passos x <- 0.3 # valor inicial para a raiz while(i<passo & dif>tol){ x[i+1] <- x[i]-(fx0(x[i])/fx1(x[i])) dif <- abs(x[i+1]-x[i]) d<-x[i] i<-i+1 } ### Obtendo a estimativa para alpha alpha=((sum(y^d))/n)^(1/d) ### d é o parâmetro de forma e alpha é o parâmetro de escala d alpha ########## Teste de aderência ypop<-rweibull(1000,d,alpha) ks.test(y,ypop) ##################################################### ###### Nelder-Mead através da função Optim ########## n=length(y) vero <- function(par,x){ k = par[1] #### Shape parameter lambda = par[2] #### Scale Parameter saida<-(n*log(k/lambda) + (k-1)*sum(log(x/lambda)) - sum((x/lambda)^k)) return(-saida) } saiday<-optim(par=c(1.8,4),fn=vero, method= "Nelder-Mead",x=y ) k<-saiday[1]$par[1] lambda<-saiday[1]$par[2] ### k é o parâmetro de forma e lambda e o parâmetro de escala k lambda ########## Teste de aderência ypop<-rweibull(1000,k,lambda) ks.test(y,ypop) ############################################### O que acontece é que as estimativas através dos dois métodos estão dando muito diferentes e através de um KS.test concluo que as estimativas através de Newton-Raphson são muito boas, pois o p-valor fica realmente bem alto; já as estimativas utilizando Nelder-Mead ficam péssimas e o ajuste, claro, fica péssimo com p-valor muito pequeno.... Gostaria de saber se existe alguma explicação para que isso esteja ocorrendo... O método de Nelder-Mead não é totalmente indicado?? O que estaria provocando essa disparidade nas estimativas... Sei que são dois métodos de otimização diferentes, mas acredito que os dois eram pra fornecer vaores bem próximos, concordam??? Observe que se eu gerar uma weibull, por exemplo: Y<-rweibull(12,5,6) As estimativas serão bem próximas e o ks.test me dará um p-valor bem condizente com a realidade (isso para os dois métodos de otimização)... O que pode estar acontecendo?? Vcs tem alguma idéia??? O problema está nos dados que estou analizando?? que tipo de problema seria esse?? Estou confuso pois quando gero dados no R os dois métodos apresentam boas estimtivas e p-valor do ks.test é bem condizente... Entretanto se uso os dados que tenho, os parâmetros só são bem estimados pelo método de Newton-Raphson o outro método fica péssimo... Agradeço qualquer dica! Abraço, Romero.

3 3

Cálulo numérico
by Andre Oliveira 06 Mai '15

06 Mai '15

Prezados,meu interesse neste momento é usar o R para uma apresentação com enfoque cálculo numérico em uma escola. Fiz alguns testes aqui como pacote library(animation), com a função optim e seus argumentos e contruí um algoritmo de Newton Raphso em caso univariado. Gostaria de saber se alguém do grupo tem experiência com R em cálculo numérico e se tem algum material para optimização e algoritmos que possa compartilhar? No caso pode ser uni e multivariado! André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 1

Arquitetura do Ubuntu 14.04
by Andre Oliveira 06 Mai '15

06 Mai '15

Prezados,sempre usei Ubuntu LTS 32 bits. Estou pretendendo instalar 64 bits. Posso ter problemas por este motivo? É uma boa ideia? Quais as experiencias? obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 1

Como definir um modelo para esse experimento?
by Simone D. Sartorio 06 Mai '15

06 Mai '15

Boa tarde! Pessoal, tenho um experimento em DIC, com 6 repetições, em esquema de parcelas subdivididas, onde tenho nas parcelas um fatorial (2 variedades X 2 vinhaças) e nas subparcelas aplicou-se aluminio (4 níveis: 60, 45, 35 e 20%). Minha dúvida: como defino no R um modelo inicial para analisar isso? Seria algo talvez parecido com isso: aov(y ~ variedade*vinhaça*aluminio + Error(parcelas/subparcelas)) ? Desde já, muito obrigadaSimone D. Sartorio

2 1

Dendograma
by Mauricio Silva. 06 Mai '15

06 Mai '15

Bom dia.estou com dados quantitativo de 117 linhagens, e gostaria plotá-los em dendograma, sendo que o circular é melhor opção de visualização.Assim quais rotinas a ser feita para plotar este dendograma. Obrigado

3 2

Interagir R com wxMaxima
by aalmeidaalessandro＠gmail.com 06 Mai '15

06 Mai '15

Alguém poderia me ajudar em como posso interagir(“conversar”) o R com wxMaxima?

4 3

Erro no Debian 7
by Savano Pereira 05 Mai '15

05 Mai '15

Prezados, Estou fazendo umas simulações utilizando o R no Debian 7. Porém, está dando o seguinte erro "address 0x1001c2c82b, cause memory not mapped". O estanho é que rodei o mesmo script no Windows 8 e no Ubuntu 14 e não apresentaram este problema, a simulação gerou os resultados esperados. Alguém poderia ajudar a resolver este problema no debian 7? Att, --- Savano Pereira TEL: 21 9-8771-9063

4 3

georeferenciamento CEP Curitiba
by Elias T. Krainski 05 Mai '15

05 Mai '15

Ola, A alguns anos ja o google mudou a API de georeferenciamento. Antigamente tinhamos o codigo de acuracia. Na nova API nao ha mais. Alem disso nao ha mais o formato csv como output (meu preferido por ser mais simples). A nova API oferece output em formato xml e json. Usei a funcao pronta para essa nova API ggmap:::georef para georeferenciar uma lista de 9183 CEPs em Curitiba. Estou usando isso para analisar dados de internacoes hospitalares. Dados individuais de internacoes hospitalares via SUS estao disponiveis no DATASUS (desde 1994). Gostaria de compartilhar essa lista em http://leg.ufpr.br/~elias/tmp/cep_curitiba_cwb_georef_longlat.txt com dois objetivos: 1) para quem tiver interesse em usar 2) convidar quem trabalha com isso ao nivel profissional a checar se ha' erros Alguem ja teve experiencia com a lista do CNEFE? Se alguem tiver interesse em georeferenciar aprox. meio milhao de enderecos do CNEFE para Curitiba, poderiamos fazer isso em equipe, pois para usuario comum e' permitido apenas 2500 por dia/IP. (sozinho eu demoraria aprox. um mes, usando meia duzia de maquinas...) Elias. -- In God we trust, all others bring data. – William Edwards Deming (1900-1993)

3 5