Março 2015 - R-br

Função identify
by ana paula coelho madeira 27 Mar '15

27 Mar '15

Caros colegas, boa noite! Ao fazer o diagnóstico de um GLM tento identificar os pontos de alavanca e pontos influentes utilizando a função identify( ), porém o R fia executando e nunca libera o resultado, daí tenho que parar a computação. Alguém tem ideia do pode estar acontecendo? Desde já agradeço. Ana Paula

4 5

Diferença entre dados coletados do yahoo usando diferentes funções.
by Michelle Bau Graczyk 26 Mar '15

26 Mar '15

Boa tarde, Eu estava verificando os diferentes modos de coletar dados intraday do yahoo e percebi que para cada função que eu uso, os dados obtidos são diferentes. Por exemplo, se eu uso: dados <- readLines(' http://chartapi.finance.yahoo.com/instrument/1.0/AAPL/chartdata;type=quote;…',sep = ",", dec = ".") eu obtenho: "18" "1427203800,127.1800,127.2600,127.0000,127.2300,631900" "19" "1427203868,126.9200,127.1900,126.8750,127.1700,202000" "20" "1427203920,126.9600,127.0600,126.8700,126.9100,189500" "21" "1427203981,126.9500,126.9900,126.8800,126.8900,198200" "22" "1427204099,126.7799,126.9900,126.7700,126.9500,217500" "23" "1427204100,126.6700,126.9000,126.6401,126.7799,286300" "24" "1427204219,126.7100,126.8000,126.6500,126.6800,222800" "25" "1427204221,126.8300,126.8700,126.5600,126.7000,238500" "26" "1427204280,126.7100,126.9200,126.6700,126.8200,167500" "27" "1427204340,126.8300,126.9300,126.7000,126.7166,114900" "28" "1427204457,126.6400,126.8300,126.6000,126.8100,156100" "29" "1427204460,126.6500,126.7900,126.6385,126.6900,102800" se eu uso: dados <- getURL(' http://chartapi.finance.yahoo.com/instrument/1.0/AAPL/chartdata;type=quote;…',sep = ",", dec = ".") obtenho: 1427203805,127.1900,127.2600,127.0796,127.2000,626100 1427203862,126.9200,127.1900,126.8750,127.1800,207800 1427203979,126.9000,127.0600,126.8700,126.9100,198000 1427203980,126.9500,126.9900,126.8800,126.9000,172200 1427204099,126.7799,126.9900,126.7700,126.9500,235000 1427204100,126.6700,126.9000,126.6401,126.7799,286300 1427204219,126.7100,126.8000,126.6500,126.6800,222800 1427204221,126.8300,126.8700,126.5600,126.7000,238500 1427204280,126.7100,126.9200,126.6700,126.8200,167500 1427204340,126.8300,126.9300,126.7000,126.7166,114900 1427204459,126.6915,126.8300,126.6100,126.8100,160900 1427204460,126.6500,126.7900,126.6385,126.6900,98000 Usando read.csv também dá dados diferentes de começo. Alguem saberia me dizer o por que? E em qual devo confiar? Porque nenhuma é igualzinha ao que aparece quando você coloca esse endereço http://chartapi.finance.yahoo.com/instrument/1.0/AAPL/chartdata;type=quote;… no navegador. O resultado deste é: 1427203800,127.1200,127.2300,127.0000,127.2300,631900 1427203900,127.0300,127.1900,127.0000,127.1700,87400 1427203978,126.9200,127.0400,126.8700,126.9050,309600 1427204038,126.9600,126.9900,126.8800,126.9399,174900 1427204099,126.7799,126.9900,126.7700,126.9600,235300 1427204100,126.6700,126.9000,126.6401,126.7799,286300 1427204219,126.7100,126.8000,126.6500,126.6800,222800 1427204221,126.8300,126.8700,126.5600,126.7000,238500 1427204280,126.7100,126.9200,126.6700,126.8200,167500 1427204340,126.8300,126.9300,126.7000,126.7166,114900 1427204459,126.6915,126.8300,126.6000,126.8100,160900 1427204460,126.6500,126.7900,126.6385,126.6900,98000 1427204520,126.7100,126.7200,126.6000,126.6500,108700 1427204581,126.7100,126.8400,126.6800,126.7100,129600 1427204641,126.7400,126.7700,126.6600,126.7300,98700 Obrigada, Michelle

3 3

Pacote QCC
by Edson Lira 26 Mar '15

26 Mar '15

Caros amigos, estou tentando replicar um exemplo do Livro do Montgomery, Introduction to Stastistical Quality Control - 7 edição, pg 240. Estou usando o pacote QCC do R. Abaixo o banco e o script Download iqcpg240.csv | | | | | | | | | | | Download iqcpg240.csvDownload iqcpg240.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Download exemplo iqc wiley pg240.R | | | | | | | | | | | Download exemplo iqc wiley pg240.RDownload exemplo iqc wiley pg240.R free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Tenho algumas dúvidas. Para os limites x de barra usa-seXbarra+-A2*Rbarra, aonde o A2 é tabelado de acordo com o tamanho da amostra. No pacote QCC, fiz algumas tentativas e não vi como ou aonde ele usa o A2, visto que quando plota os lintel ele só soma ou diminui o desvio padrão. Alguém poderia me ajudar?. No livro do Montgomery os limites para o barra são (25 amostra de tamanho 5)1,5056+-(0,577)*(0,32521) Não consegui replicar esses resultados no gráfico. Usei as 25 primeiras amostras do banco, que é o que o livro usa no exemplo. Alguém poderia me ajudar? [ ]'s.Edson Lira Estatístico Manaus-Amazonas

1 0

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 25 Mar '15

25 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 23 Mar '15

23 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