Julho 2012 - R-br - DINF UFPR

Importar tabela com células vazias
by Luciano F. Sgarbi 13 Jul '12

13 Jul '12

Boa tarde, Estou com um problema para importar uma tabela em formato ".txt" com várias células vazias. Para tal, estou usando o comando "read.table("name.file.txt",h=F,fill=T)". O que está acontecendo é que quando eu importo, o arquivo só é importado com as duas primeiras colunas. Alguém poe me ajudar? A tabela está no formato do exemplo abaixo: 1 125 2 151 3 -59,8184 4 -14,009 5 0,00167 6 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 Abraços, Luciano F. Sgarbi -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

3 2

Acessar banco de dados
by Emerson Cotta Bodevan 12 Jul '12

12 Jul '12

Prezados membros da lista. Tenho um arquivo ".MDB" gerado pelo Epi Info. Estou tentando acessá-lo no R. Com o comando abaixo eu consegui fazer uma conexao com o arquivo. automed<-odbcConnectAccess("Estudo de automedicação em crianças.mdb") O que eu gostaria agoraera poder acessar as variáveis individualmente (como se fosse um data frame), para que eu possa analisar as variáveis separadamente. Alguém sabe como fazê-lo? Estou com o pacote RODBC instalado. Agradeço qualquer ajuda, Emerson. -- "A educação não muda o mundo. A educação muda as pessoas. As pessoas mudam o mundo." * Paulo Freire.*

3 2

Séries temporais [OFF-TOPIC]
by ASANTOS 12 Jul '12

12 Jul '12

Bom dia pessoal, Alguém tem conhecimento de algum livro que tenha em seu conteúdo algo sobre séries temporais inflacionadas com zeros usando o R. Procurei no Times series analysis do Shumway e Stoffer que tenho como referência e não encontrei nada sobre abordagens inflacionada com zeros. Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 1

Re: [R-br] Mapa Manaus por zonas e bairros:acentuação e legendas
by Benilton Carvalho 12 Jul '12

12 Jul '12

Copiando a lista novamente.... plot(mao_bz, col = mao_bz$ZONA) legend('bottomright', unique(as.character(mao_bz$ZONA)), fill=unique(mao_bz$ZONA), bty='n') b 2012/7/12 Edson Lira <edinhoestat(a)yahoo.com.br>: > Valeu Benilton!!! > E as legendas? Como faço? São 6 zonas e em torno de 60 bairros. > > Edson Lira > Estatístico > Manaus-Amazonas > ________________________________ > De: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > Para: Edson Lira <edinhoestat(a)yahoo.com.br> > Enviadas: Quarta-feira, 11 de Julho de 2012 20:54 > > Assunto: Re: [R-br] Mapa Manaus por zonas e bairros:acentuação e legendas > > o encoding original dos dados a q vc se refere e' ISO-8859-1. > > use > > toUTF8(mao_bz$NOME, 'ISO-8859-1') > > ou qq um dos compativeis (CP819, CP822... talvez ate' o 1252) > > b > > 2012/7/12 Edson Lira <edinhoestat(a)yahoo.com.br>: >> Link para o shx >> >> http://www.datafilehost.com/download-3ef3a7f9.html >> >> >> sessionInfo() >> R version 2.14.2 (2012-02-29) >> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) >> >> locale: >> [1] pt_BR.UTF-8/pt_BR.UTF-8/pt_BR.UTF-8/C/pt_BR.UTF-8/pt_BR.UTF-8 >> >> attached base packages: >> [1] stats graphics grDevices utils datasets methods base >> >> other attached packages: >> [1] classInt_0.1-17 e1071_1.6 class_7.3-3 descr_0.9.7 >> shapefiles_0.6 maptools_0.8-14 lattice_0.20-0 sp_0.9-98 >> foreign_0.8-49 >> >> loaded via a namespace (and not attached): >> [1] grid_2.14.2 tools_2.14.2 xtable_1.7-0 >> >> >> [ ]'s. >> >> Edson Lira >> Estatístico >> Manaus-Amazonas >> ________________________________ >> De: Benilton Carvalho <beniltoncarvalho(a)gmail.com> >> Para: r-br(a)listas.c3sl.ufpr.br; Edson Lira <edinhoestat(a)yahoo.com.br> >> Enviadas: Quarta-feira, 11 de Julho de 2012 20:25 >> Assunto: Re: [R-br] Mapa Manaus por zonas e bairros:acentuação e legendas >> >> e tenho a impressao de q faltou o o link pro shx? >> >> > >

1 0

Mapa Manaus por zonas e bairros:acentuação e legendas
by Edson Lira 12 Jul '12

12 Jul '12

Caros amigos, estou gerando o mapa de Manaus por zona, minha máquina é um Mac i7, 4gb ram. Estou com o problema da acentuação, a linha: mao_bz$NOME <- toUTF8(mao_bz$NOME, "IBM850") me dá a saida abixo: > mao_bz$NOME [1] PONTA NEGRA [2] LĂ˘Â\u0088Â\u009aĂ\u0083ÂŽĂ˘Â\u0088Â\u009aĂ˘Â\u0089Â¤Ă˘Â\u0088Â\u009aĂ\u0083Â˘RIO DO VALE [3] COL TERRA NOVA [4] TARUMĂ˘Â\u0088Â\u009aĂ\u0083ÂŽĂ˘Â\u0088Â\u009aĂ˘Â\u0088Â\u0082Ă\u0082ÂŹĂ\u0082ÂŁ [5] SANTA ETELVINA [6] MONTE DAS OLIVEIRAS [7] FLORES [8] CIDADE NOVA [9] NOVO ISRAEL [10] COL STO ANTĂ˘Â\u0080Â\u009aĂ\u0083ÂŽĂ\u0083ÂşĂ˘Â\u0088Â\u009aĂ˘Â\u0084Â˘NIO [11] REDENĂ˘Â\u0080Â\u009aĂ\u0083ÂŽĂ\u0083ÂşĂ˘Â\u0088Â\u009aĂ\u0082Â˘Ă˘Â\u0088Â\u009aĂ\u0083ÂŽĂ˘Â\u0088Â\u009aĂ˘Â\u0088Â\u0082Ă\u0082ÂŹĂ\u0082ÂŁO . . . O script que estou usando é o abaixo: mao_bz<-readShapePoly("/Users/edsonlira/Edinho/Edinho/Bases_MAO/Setores_MAO/Bairros_e_Zonas.shp") summary(mao_bz) plot(mao_bz) library(descr) mao_bz$NOME <- toUTF8(mao_bz$NOME, "IBM850") mao_bz$NOME # palette(c("#ccddff", "#ffddcc", "#ccffcc", "#ffffcc", "#ccffff")) palette(rainbow(6)) plot(mao_bz, col = mao_bz$ZONA) title("Mapa de Manaus por bairros e zonas geográficas") Os links abaixo tem as base de dados usadas. http://www.datafilehost.com/download-00405c50.html #bairros e zonas.shp http://www.datafilehost.com/download-c36f359a.html #bairros e zonas.dbf http://www.datafilehost.com/delete-a54c1f16-16120426.html #bairros e zonas shx Outra dúvida: Como faço para criar uma legenda com os nomes das zonas e/ou com o nome dos bairros? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 3

Como escrever o modelode Michaelis-Menten na linguagem BUGS - R2WinBUGS
by Augusto Ribas 11 Jul '12

11 Jul '12

Opa. Alguem usa o winbugs com o R? Eu estava tentando ajustar um Modelo de michelis-menten e não consegui. Se alguem puder dar uma olhadinha e ver onde estou errado. Segue um exemplo do que estou fazendo: library(R2WinBUGS) ################################################################# #Relação Riqueza~estudos ################################################################# #Modelo do que eu acho que acontece #riqueza~(c*estudos)/(d+estudos) #Fazendo um levantamento do numero de especies de parasitas que parasitam algumas especies de sapos, eu supuz que qd uma especie tinha poucos estudos #tava sempre bem subestimado a riqueza de parasitas associadas aquela especie. #dai como a unica informação que fico é com o número de estudos, pensei que deveria ter uma relação onde qt mais estudos, mais especies eu devo encontrar #até um limite, o limite é o parametro c nesse caso, e d representa como os estudos vão contribuindo pra eu chegar ao valor limite de especies. #existem alguns problemas acredito, como especies podem ser mais faceis ou mais dificies de estudar, e o limite pra cada especie depende de mais coisas e tipos de estudos. #mas ficando so nesse caso, eu fiz esse exemplo set.seed(13) #gerando dados exemplo estudos<- rpois(48,5) #parametros originais #então o limite de especies de parasitos seria 10 c<-10 #e eu preciso de mais ou menos 5 estudos pra ver legal a riqueza de especies, d<-2.5 #resposta riqueza<-round(((c*estudos)/(d+estudos))+rnorm(length(estudos),0,1),0) plot(riqueza~estudos,xlim=c(0,15),ylim=c(0,15),xlab="Riqueza de espécies",ylab="Número de estudos") curve(c*x/(d+x),0,15,add=T,col="black") #ajustando com modelo não linear modelo.estudos<-nls(riqueza~((a*estudos)/(b+estudos)),start=list(a=8,b=1.5)) summary(modelo.estudos) points(predict(modelo.estudos)[order(predict(modelo.estudos))]~estudos[order(estudos)],type="l",col="red",lty=2) #eu li que é possivel lineariza essa função e ajusta com glm da seguinte forma: #ajustando com glm modelo.menten<- glm(riqueza~c(1/estudos),family=gaussian(link = "inverse")) summary(modelo.menten) a1<-1/modelo.menten$coef[1] b1<-a1*modelo.menten$coef[2] a1 b1 curve(a1*x/(b1+x),0,15,add=T,col="blue",lty=2) ### Ajuste com winbugs # Definir modelo sink("michaelis.menten.test.txt") cat(" model { # Priors alpha ~ dnorm(0,0.001) beta ~ dnorm(0,0.001) sigma ~ dunif(0, 100) # Likelihood for (i in 1:n) { riqueza[i] ~ dnorm(mu[i], tau) mu[i] <- alpha + beta * estudos[i] } # Derived quantities tau <- 1/ (sigma * sigma) } ",fill=TRUE) sink() # Juntar os dados - Note que eu não entendi como fazer a função de ligação dentro do modelo winbugs #dai eu registrei os dados como 1/riqueza e 1/estudos aqui, acho que isso que o pepino deve ta nisso win.data <- list(riqueza=c(1/riqueza),estudos=c(1/estudos),n=length(estudos)) # Valores iniciais inits <- function(){ list(alpha=rlnorm(1), beta=rlnorm(1),sigma = rlnorm(1))} # Paramstros a estimar params <- c("alpha","beta","tau") # MCMC settings nc <- 3 ni <- 3000 nb <- 1000 nt <- 2 # Começa o winbugs out <- bugs(data=win.data, inits=inits, parameters.to.save=params, model.file="michaelis.menten.test.txt", n.thin=nt, n.chains=nc, n.burnin=nb, n.iter=ni, debug = TRUE) ###Olhar o modelo print(out, dig = 3) #retornoar os parametros como no glm acima a2<-1/out$summary[1,1] b2<-a*out$summary[2,1] a2 b2 curve(a2*x/(b2+x),0,15,add=T,col="green",lty=2) legend("topright",legend=c("Original","Ajustado com nls","Ajustado com glm","Ajustado com WinBUGS"), lty=c(1,2,2,2),col=c("black","red","blue","green")) #Minhas duvidas são: #01. como montar o modelo na linguagem bugs pra fazer esse tipo de analise? #procurando pela internet eu não achei nenhum exemplo pra esse caso especifico. #02. como eu colocaria dentro do modelo esse função de ligação inversa. #eu tentei usar a função pow(mu[i],-1) mas não rolou, #03. e o que estou fazendo no modelo winBUGS tem muito problemas? entrar com os dados daquele jeito, ao invez de especificar no modelo? #não consigo imaginar a implicação de não definir no modelo a função de ligação. Bom se alguem tiver esperiencia com winbugs e puder dar uma olhadinha no que estou fazendo errado Ou alguma dica se o que estou fazendo tem algum problema critico -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 0

Integração de C com R
by Jodavid Ferreira 11 Jul '12

11 Jul '12

Alguém conhece alguma apostila sobre integração de C com R? estou precisando estudar sobre isso, estou procurando nos sites de buscas, mas não estou encontrando resultados muito satisfatórios! -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

3 2

Re: [R-br] Como gerar dados de um experimento
by Walmes Zeviani 10 Jul '12

10 Jul '12

CAR e SAR já não são minha praia mais. Mas você pode usar a função spdep::cell2nb() para obter essa matriz de proximidade. À disposição. 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 twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

help para constrastes
by Clenilso Sehnen Mota 10 Jul '12

10 Jul '12

Prezados (as) Estou analisando os resultado de um trabalho. Minha intensão é verificar a comparação conforme a tabela 2. A tabela 1 mostra os tratamentos e as combinações dos 13 componentes. Já discuti pessoalmente com algumas pessoas, se a forma com estimei esses contrastes esta correta ou não. E não chegamos a um consenso. Gostaria de saber a resposta de você se esta correto ou não. Desde já agradeço pela ajuda. Att Clenilso Tabela 1- Composição dos tratamentos: Trat Combinações 1 A 2 B 3 B+C (3:1) 4 B+C (1:1) 5 B+C (1:3) 6 B+D+E(3:1:6) 7 B+D+C(2:1:2) 8 B+D+C (1:1:1) 9 F 10 G+C(7:3) 11 D+C(1:1) 12 D+H (1:1) 13 I 14 D+J+K(2:2:1) 15 D+K(3:1) 16 L Tabela 2 Composição dos contrastes. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Contraste 1 (A e L vs Demais) -7.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 -7.0 Contraste 2 (C vs s/C) 0.0 1.0 -1.0 -1.0 -1.0 1.0 -1.0 -1.0 1.0 -1.0 -1.0 1.0 1.0 1.0 1.0 0.0 Contraste 3 (B vs s/B) 0.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 0.0 Contraste 4 (E, F, H e K vs s/E, F, H e K) 0.0 1.0 1.0 1.0 1.0 -1.8 1.0 1.0 -1.8 1.0 1.0 -1.8 1.0 -1.8 -1.8 0.0 Contraste 5 (D vs s/D) 0.0 1.0 1.0 1.0 1.0 -1.0 -1.0 -1.0 1.0 1.0 -1.0 -1.0 1.0 -1.0 -1.0 0.0 Clenilso Sehnen Mota Eng. Agronômo - Udesc DSc Fisiologia Vegetal - UFV MSc produção Vegetal Udesc/Embrapa E-mail: <mailto:csm.sehnen@gmail.com> csm.sehnen(a)gmail.com Celular: (64) 8150 8808

2 2

[Dúvida] Linguagem R e Julia.
by Pedro Rafael 08 Jul '12

08 Jul '12

Pessoal, está existindo esforços para tornar a linguagem R mais eficiente no que diz respeito eu desempenho computacional de determinadas tarefas, como por exemplo nos laços que em R são mais demorados que algumas outras linguagens interpretadas como Lua, Python, JavaScript, etc? Atualmente vejo a linguagem Julia se promovendo com comparativos de tempo de processamento com outras linguagens como Matlab, Python, Octave e inclusive o R. Acho que daqui a uns 5 ou 10 anos poderemos utilizar essa linguagem como uma boa ferramenta para computação científica geral. Contudo, acredito que ela irá precisar de muitos outros fatores para que seja uma alternativa para estatistica como o R. Posso estar falando cedo demais, mas acredito que Julia não só vislumbra a estatística. Vejo que será uma linguagem alternativa compilada para Scilab, Matlab, Octave e demais outros softwares de computação numérica. Dessa forma, se esse cenário permanecer vejo que o R ainda será o ícone de software livre para estatística. Sem falar na existência de um local (CRAN) que reúne boa parte das contribuições de colaboradores de toda parte do mundo que é de extrema importância para evolução de um software livre. [ ], Pedro Rafael Diniz Marinho.

1 0