- R-br - DINF UFPR

leitura dados PNAD 2003 - variáveis V1701 a V1377
by Luciane Maria Pilotto 25 Mar '15

25 Mar '15

Bom dia, estou tentando ler as variáveis da PNAD de 2003 - V1701 a V1377- relacionadas ao suplemento Saúde- e as mesmas aparecem como NA. O banco está disponível no site do IBGE http://www.ibge.gov.br/home/estatistica/populacao/trabalhoerendimento/pnad2…. Também podem ser acessados no dropbox: https://www.dropbox.com/sh/d0z62bg4qqyg2x3/AAAVDYWjVlb4JpUTYlQfvq-Ka?dl=0. Estou usando os comandos abaixo: #----------------Ler microdados da PNAD 2003 (variáveis referente ao suplemento Saúde) utilizando a função read.fwf.r ----------------------# source(file.choose()) #carregar a função selectvar da pasta funcoes source(file.choose()) #carregar a função getdic da pasta funcoes # Gerar os dicionários de dados a partir dos arquivos: 'input PES2003.txt' dicPNAD2003Pes <- getdic(file.choose())# abrir o arquivo input PES2003.txt, dentro da Pasta Input head(dicPNAD2003Pes) str(dicPNAD2003Pes) nrow(dicPNAD2003Pes) dim(dicPNAD2003Pes) # Para EXCLUSÃO do "UF" (que apresenta erro na importação por começar na mesma posição que número de controle): dicPNAD2003Pes <- subset(dicPNAD2003Pes, cod !="UF") head(dicPNAD2003Pes) dim(dicPNAD2003Pes) ####Usar estes comandos para importar *apenas algumas variáveis* da base de dados das PESSOAS variaveis <- selectvar(dicPNAD2003Pes) # Com a tecla CTRL pressionada, selecione as variáveis de interesse, neste caso são as que aparecem como NA: Erro em substring(x, first, last) : invalid substring arguments ### PROBLEMA: Ao abrir a janela das variáveis utilizando o comando acima, as variáveis de interesse referentes a saúde aparecem como NA, diferentemente das demais variáveis que aparecem no banco. No dicionário as mesmas aparecerem corretamente. Atenciosamente, __________________________________________________ Luciane Maria Pilotto

5 6

metodologia para validação de tratamentos
by Fernando Antonio de souza 25 Mar '15

25 Mar '15

Caros amigos, Estou analisando dados que comparam diferentes metodologias. Os dados referem a medidas de áreas de olho de lombo e estou comparando 4 metodologias utilizadas para medila (controle, A, B, C). Eu realizei a anova e utilizei o teste de dunnet para fazer as comparações dos tratamento com o grupo controle. Gostaria de saber quais avaliações estatísticas a mais eu posso fazer para melhorar esta análise? Quais análises de resíduos eu posso utilizar? abçs -- ======================================================================= Fernando Souza Zootecnista, DSc. Produção Animal celular: (+55) 82 8113-8781 e-mail:nandodesouza@gmail.com https://producaoanimalcomr.wordpress.com/ ========================================================================

3 4

Reduzir a distância entre código e output no knitr
by Marcus Nunes 24 Mar '15

24 Mar '15

Olá Fiz esta pergunta há alguns minutos no StackExchange, mas ainda não obtive resposta. (aliás, se alguém quiser ler ela no contexto original e responder por lá ou, quem sabe, vê-la identada corretamente etc, o link é http://stackoverflow.com/questions/29240036/adjust-the-distance-between-cod… ) Estou escrevendo uma série de slides sobre as funcionalidades gráficas do R. Um dos meus slides é assim: --- \begin{frame}[fragile] \frametitle{A Região Gráfica} <<RegiaoGrafica01>>= plot(0:10, 0:10) @ \end{frame} --- O problema é que meu slide compilado fica da seguinte forma: http://i.stack.imgur.com/3Ozdd.jpg Meu gráfico está muito distante do código. Eu gostaria de saber se existe alguma opção no knitr que diminua esta distância. Já procurei no manual e não consegui encontrar nada do gênero. Segue meu CMR (i.e., meu arquivo .rnw): --- \documentclass[10pt, compress, usetitleprogressbar]{beamer} \usetheme{m} \usepackage{booktabs} \usepackage[scale=2]{ccicons} \usepackage{minted} \usepackage{amssymb} \usepackage{amsfonts} \usepackage{amsmath} \usepackage{amssymb} \usepackage{amsthm} \usepackage[portuges]{babel} \usepackage{bm} \usepackage{icomma} \usepackage[retainorgcmds]{IEEEtrantools} \usepackage{mathtools} \usepackage{multirow} \usepackage{natbib} \usepackage{shadow} \usepackage{subfigure} \usepackage{tikz} \usetikzlibrary{shapes.arrows} \usepackage{verbatim} \usepackage{xcolor} \usepgfplotslibrary{dateplot} \usemintedstyle{trac} %%% Identification \title{O Subsistema Gráfico} \subtitle{EST0091 - Introdução ao R} \date{26 de março de 2015} \author{Marcus Nunes} \institute{Universidade Federal do Rio Grande do Norte} %%% New colors \definecolor{mDarkBrown}{HTML}{604c38} \definecolor{mDarkTeal}{HTML}{23373b} \definecolor{mLightBrown}{HTML}{EB811B} \definecolor{mMediumBrown}{HTML}{C87A2F} \definecolor{mLightRed}{HTML}{EB331B} \definecolor{mLightGreen}{HTML}{14B03D} \definecolor{mLightTeal}{HTML}{137D91} \begin{document} % new commands \newcommand{\R}{\mathbb{R}} \newcommand{\Xb}{\bm{X}} \newcommand{\Yb}{\bm{Y}} %%%%%%%%%%%%%%%%%%%%%% % knitr options <<setup, include=FALSE>>= library(knitr) opts_chunk$set(fig.path='figure/beamer-',fig.align='center',fig.show='hold',size='footnotesize') @ <<options, cache=FALSE, echo=FALSE, eval=TRUE, tidy=TRUE, dpi=144>>= options(width=60) #par(cex=2) par(mar=c(0, 0, 0, 0)) opts_knit$set(out.format = "latex") opts_chunk$set(prompt=TRUE, tidy.opts=list(width.cutoff=40), tidy=FALSE, fig.width=6, fig.height=4.5) knit_theme$set("print2") @ %%%%%%%%%%%%%%%%%%%%%% \maketitle %%%%%%%%%%%%% %%% SLIDE %%% \plain{A Região Gráfica} %%%%%%%%%%%%% %%%%%%%%%%%%% %%% SLIDE %%% \begin{frame}[fragile] \frametitle{A Região Gráfica} <<RegiaoGrafica01>>= plot(0:10, 0:10) @ \end{frame} %%%%%%%%%%%%% \maketitle \end{document} --- Obrigado, -- Marcus Nunes http://marcusnunes.me/

2 2

Coincidências entre vetores
by Ze Henrique 24 Mar '15

24 Mar '15

Bom dia pessoal, estou precisando de uma ajuda de vcs. Tenho 2 vetores numéricos de tamanho diferente e gostaria de verificar as coincidências entre ambos os vetores. a=c(1,3,5,7,8) b=c(2,3,6,7) por exemplo sei que neste caso ocorreu 2 coincidências entre os vetores a e b. Desde já agradeço a ajuda de todos. -- *José Henrique Soler Guilhen* Mestrando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

4 6

friedman - agricolae - between blocks
by Daniela Recchia 24 Mar '15

24 Mar '15

Boa tarde pessoal, meu experimento a ser analisado é com medidas repetidas, por exemplo: Bloco1 Trat1 Valor1Bloco1 Trat1 Valor2Bloco1 Trat1 Valor3Bloco1 Trat2 Valor1Bloco1 Trat2 Valor2 Bloco1 Trat2 Valor3 Bloco2 Trat1 Valor1 Bloco2 Trat1 Valor2.... usando a função friedman do agricolae consigo realizar apenas a comparação dos tratamentos dentro de cada bloco, mas gostaria também de poder realizar uma comparação entre os blocos, isso seria possível? Obrigada! Daniela -- Daniela Rodrigues Recchia

1 0

Web Scraping Lattes com R
by L N 24 Mar '15

24 Mar '15

Olá Pessoal, Alguém que saiba extrair os CV Lattes (XML) via R? Tenho listagem de todos os endereços e quero montar um banco de dado com isso. Há o ScriptLattes para Python mas eu já tentei e não consegui. Abs, *Leonardo F. Nascimento*

3 4

O tamanho das bolas ajuda a descobrir quais são as bolas pretas?
by Ricardo Dobrovolski 22 Mar '15

22 Mar '15

Oi Pessoal, Tudo bem? Suponham que quero descobrir quais são as bolas pretas ("col" = 1) em um conjunto de milhares de bolas. Sei porém que dois fatores podem afetar minha capacidade de adivinhar a cor da bola. Um é conhecido, a rugosidade, bolas rugosas ("rug" = 1) tendem a ser pretas. Outro é o efeito que quero descobrir, o tamanho da bola ("siz"). Como fazer um teste de aleatorização para avaliar separadamente esses efeitos? Em anexo, está o arquivo com os dados e uma sugestão de solução do problema, baseado em um teste de aleatorização em que seleciono dentre o conjunto de bolas, aquelas que são rugosas e dentre estas vou selecionando bolas com tamanhos cada vez menores a fim de descobrir se existe efeito do tamanho e qual o tamanho máximo que a bola tem de ter para influir na minha decisão. Porém, eu tenho dois problemas, o poder do teste de aleatorização, pois quando vou diminuindo o tamanho da bola, meu universo de possibilidades é reduzido. E o outro é o fato de que faço teste para diferentes categorias de tamanho de bola, como extrair uma conclusão geral. Não sei se alguém vai ter saco de me ajudar, mas isso acontecer, agradeço. Se alguém souber de literatura que trate desses temas, também será bem-vinda. Abraço, Ricardo Dobrovolski

2 1

P-valor
by ana paula coelho madeira 22 Mar '15

22 Mar '15

Bom dia! Aplicando o teste de Kolmogorov à uma amostra "proveniente " de uma distribuição gama, obtem-se: ks.test(d, "pgamma", ahat, bhat) #p-value = 0.9779 One-sample Kolmogorov-Smirnov test data: d D = 0.0503, p-value = 0.9779 A dúvida é: como faço para encontrar o p-valor associado ao teste de uma distribuição não implementada n R? Sei que devo calcular a estatística D, mas e o p valor? Desde já agradeço Ana Paula

2 2

Empilhar variáveis em colunas para linhas
by ASANTOS 22 Mar '15

22 Mar '15

Caros Listeiros, Tenho um conjunto de dados artificiais abaixo, constituído de dados de contagem de três espécies (Esp1, Esp2 e Esp3) e duas variáveis trat e tempo. As três espécies e as variáveis estão dispostas em colunas, sendo: #Dados artificiais ----------------------------------- da <- expand.grid(trat=gl(4,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) da$Esp1<-rpois(da[,1],5) da$Esp2<-rpois(da[,1],10) da$Esp3<-rpois(da[,1],15) head(da) trat tempo Esp1 Esp2 Esp3 1 1 1 4 13 8 2 2 1 7 7 17 3 3 1 3 12 16 4 4 1 8 8 13 5 1 2 4 13 22 6 2 2 4 8 10 Porém, eu gostaria de reorganizar o banco de dados e substituir as três espécies (Esp1, Esp2 e Esp3), por duas colunas, sendo uma coluna com o nome das espécies e uma nova coluna para os dados de contagem, representada por N, de forma que mantenha a correspondência entre as variáveis trat e tempo, ficando da seguinte maneira: head(new.da) trat tempo Esp N 1 1 1 Esp1 4 2 1 1 Esp2 13 3 1 1 Esp3 8 4 2 1 Esp1 7 5 2 1 Esp2 7 6 2 1 Esp3 17 Alguém poderia me indicar alguma função ou procedimento para fazer isso? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção 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 Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

3 5

R versões 3.2 e 3.3 para windows
by Cleber N.Borges 21 Mar '15

21 Mar '15

olá a todos, Vi no site do R, que as próximas versões para Windows terá modificações, como por exemplo, na versão 3.2 prestes a ser lançada: WINDOWS BUILD CHANGES A number of changes to the Windows build system are in development. The following are currently in place. * Installation using external binary distributions of zlib, bzip, liblzma, pcre, libpng, jpeglib and libtiff is now required, and the build instructions have been revised. (mas que, até onde percebi, ainda terá essas libs (dll) empacotadas no R-3.2 já no R-3.3 (só ano que vem), isso será removido de vez... Alguém sabe falar o porquê do R-Core Team ter adotado essa estratégia? Algum problema de patente ou coisa assim? ... ou se isso implica em alguma melhoria? Agradeço antecipadamente por qualquer comentário. Grato Cleber Borges --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

1 0