Junho 2013 - R-br - DINF UFPR

[Dúvida] Estimação por Máxima Verossimilhança.
by Pedro Rafael 07 Jun '13

07 Jun '13

Senhores tenho uma dúvida. Na verdade não é dúvida, apenas quero sugestões. Tenho uma função densidade de probabilidade "complicada". Trate-se de uma distribuição chamada Kwmarashwamy Weibull Poisson. Tenho alguns bancos de dados e gostaria de verificar o ajustamento dessa distribuição à estes bancos de dados. Estou estimando os parâmetros pelo método de máxima verossimilhança. Essa distribuição tem suporte nos reais positivos (x>0) e todos os seus parâmetros são positivos. Optei em utilizar o método L-BFGS-G para restringir a busca nos reais positivos. Segue abaixo o comando. Nesse exemplo não houve convergência. Percebi que os chutes iniciais influenciam muito as estimativas dos parâmetros no caso em que há convergência. Usando métodos de maximização diferentes em muitos casos há grandes diferenças nas estimativas. Existe alguma forma mais tranquila e direta para encontrar as estimativas pelo método de máxima verossimilhança em R? O código segue abaixo: vero <- function(par,x){ a = par[1] b = par[2] c = par[3] lambda = par[4] beta = par[5] -sum(log((a*b*c*lambda*(beta^c)*(x^(c-1))*((1-exp(-(x*beta)^c))^(a-1)) * ((1-(1-exp(-(beta*x)^c))^a)^(b-1)) * exp(-lambda*(1-(1-(1-exp(-(beta*x)^c))^a)^b) - (beta*x)^c))/(1-exp(-lambda)))) } dados = c(17.23, 28.92, 33.00, 41.52, 42.12, 45.60, 48.80, 51.84, 51.96, 54.12, 55.56, 67.80, 68.64, 68.64,68.88, 84.12, 93.12, 98.64, 105.12, 105.84, 127.92, 128.04, 173.40) optim(par=c(1,1,1,1,1),fn=vero, method="L-BFGS-B",x=dados/1000, lower=c(0.001,0.001,0.001,0.001,0.001), upper=c(Inf,Inf,Inf,Inf,Inf)) O engraçado é que essa distribuição é encaixada com a distribuição Weibull. Quando tento ajustar a Weibull à esses dados há convergência. O que vocês acham que devo fazer? [ ], Pedro Rafael Diniz Marinho.

6 12

Inserir NA em objeto do tipo lista
by Aline Lipsky 07 Jun '13

07 Jun '13

Bom dia, Eu tenho um objeto do tipo lista com mais de 6 mil vetores com tamanhos diferentes.Como faço para acrescentar NA nos vetores de tamanho menor e transformar esse objeto em um data frame? Grata, Aline

2 1

encoding de páginas web readLines XML
by walmes . 07 Jun '13

07 Jun '13

Saudações, Estou fazendo um estudo com vídeos do youtube. A partir dos links eu retiro informações deles, como o título. No entanto, não estou conseguindo acertar o encoding, seja para ler corretamente a página ou para tranformar as strings após a leitura. Com o CMR abaixo eu leio o código de um vídeo do youtube e pego seu título, o problema é que os acentos não são corretamente representados, e eu preciso que sejam. Tentei atribuir um encoding na readLines() mas nenhum dos usados (latin1 e utf-8) funcionaram. Tentei usar a iconv() mas sem sucesso. Uso ubuntu 12.04. Sessão R está no CMR. library(XML) # lê a página do vídeo com esse endereço lines <- readLines("http://www.youtube.com/watch?v=P3m6rJR2yDI") # lines <- readLines("http://www.youtube.com/watch?v=P3m6rJR2yDI", encoding=?) # organiza h <- htmlTreeParse(lines, asText=TRUE, useInternalNodes=TRUE) summary(h) # pega o título do vídeo n <- unlist(getNodeSet(doc=h, path="//title", fun=xmlValue)) n Encoding(n) # retira os caracteres estranhos substituindo por vazio iconv(nap, from="UTF-8", to="ASCII", sub="") # iconv(n, from=?, to=?) *# Como sai # "MatemÃ¡tica - EstatÃstica: Conceito e frequÃªncias - YouTube" # Como eu gostaria # "Matemática - Estatística: Conceito e frequências - YouTube"* ## > sessionInfo() ## R version 3.0.1 (2013-05-16) ## Platform: i686-pc-linux-gnu (32-bit) ## locale: ## [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C ## [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 ## [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8 ## [7] LC_PAPER=C LC_NAME=C ## [9] LC_ADDRESS=C LC_TELEPHONE=C ## [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C ## attached base packages: ## [1] stats graphics grDevices utils datasets methods base ## other attached packages: ## [1] XML_3.96-1.1 ## loaded via a namespace (and not attached): ## [1] compiler_3.0.1 tools_3.0.1 ## > Alguém saberia solucionar esse problema? Desde já grato. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

4 4

Adicionar Na em uma lista
by Aline Lipsky 06 Jun '13

06 Jun '13

Boa tarde, Eu tenho um objeto do tipo lista com mais de 6 mil vetores com tamanhos diferentes. Eu consigo acrescentar o NA em um elemento da lista: especies[[1]][[7]] <- NA Se alguém puder me ajudar, como faço para repetir esse passo em todos os elementos da minha lista? Tentei com a função sapply, mas não dá certo. Grata, Aline

5 5

MANTEL
by Cézar Borges 05 Jun '13

05 Jun '13

Olá caros,Realizei um teste de mantel para 19 fragmentos para verificar a dependência espacial da comunidade arbórea, o resultado (r=-0.01; p=0.55) me mostrou que os fragmentos são independentes espacialmente. No entanto, estou com a seguinte dúvida: posso interpretar este resultado como se todos os fragmentos fossem diferentes entre si??? ou se não, como posso estar verificando quais os fragmentos que são iguais e diferentes entre si???Abraços! Att, Cézar Fonseca Borges Mestrando em Ciências Ambientais UFPA/MPEG/EMBRAPA - Amazônia Oriental Geógrafo Email: cezar-fonseca(a)hotmail.com/ caborges(a)museu-goeldi.br Lattes: http://lattes.cnpq.br/9805721240137485

1 0

teste de Jonckheere
by Cristina 04 Jun '13

04 Jun '13

Olá pessoal! Gostaria de saber se alguém já usou, no R, o teste de Jonckheere. tentei usar o teste mas consigo encontrar qual pacote está instalado. Obrigada Teresa Cristina

3 3

Coordenadas UTM para Lat Lon
by Tito Conte 04 Jun '13

04 Jun '13

Alguém poderia me explicar como faço a converção de utm para lat lon, estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão Tito Conte

4 4

erro na validação cruzada
by Hélio Gallo Rocha 04 Jun '13

04 Jun '13

Caros da Lista, Está dando este erro ao tentar rodar a validação cruzada na direção de 45 graus: xvalid: performing cross-validation at location ... 1, Erro em solve.default(Vcov, trend.d) : sistema é computacionalmente singular: condição recíproca número = 8.6266e-21 Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Remover compilação
by Sérgio Henrique almeida da silva ju 04 Jun '13

04 Jun '13

Amigos Estou usando a versão mais nova do R (3.0.1) no ubuntu. Tive que compilar a versão 2.14, para fazer um trabalho, agora quero remover essa compilação, como faço? Sempre quando chamo o R no terminal aparece a versão mais antiga, como faço pra voltar a chamar o R 3.0.1? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

2 2

Separação de duas distribuições de dados
by ASANTOS 03 Jun '13

03 Jun '13

Boa tarde pessoal, Simulei duas distribuições normais e representei em um histograma que obviamente assumiu a forma bimodal, sendo: ### set.seed(10) dist1<-rnorm(n=400,mean=10,sd=5) dist2<-rnorm(n=400,mean=35,sd=8) dist.c<-rbind(dist1,dist2) # hist(dist.c) # Gostaria de saber se existe algum procedimento que me permita separar dist1 e dist2 com seus respectivos intervalos de confiança, imaginando que eu não conheça os parâmetros media e desvio padrão que os criaram? Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 1