- 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

Re: [R-br] Calculo de correlcao entre coluas que apresentam NAS de varias matrizes
by Éder Comunello 21 Mar '15

21 Mar '15

Olá, Michelle! Inicialmente gostaria de sugerir que você disponibilizasse os links dos arquivos do Dropbox diretamente no teu script, evitando as diversas mensagens de compartilhamento na lista. Para isso você pode consultar o tópico que segue: < http://r-br.2285057.n4.nabble.com/R-br-Problema-em-abrir-conexao-com-source… > Quanto ao problema em si, não entendi muito bem, mas deixo um trecho de código que pode te ajudar: ### <code r> urls1 <- c(" https://www.dropbox.com/s/kgxo3s8qrrv74p3/AAparaCalculodaMedia.txt", " https://www.dropbox.com/s/fhutciprhgt077x/AIGparaCalculodaMedia.txt", " https://www.dropbox.com/s/kvubgji1tce8sm8/AXPparaCalculodaMedia.txt") ### alterar links públicos para carregamento direto! urls2 <- sub("^.*/s/", "http://dl.dropbox.com/s/", urls1) AA <- read.table(urls2[1], header=T, as.is=T) AIG <- read.table(urls2[2], header=T, as.is=T) AXP <- read.table(urls2[3], header=T, as.is=T) sapply(c("AA", "AIG", "AXP"), function(x) length(get(x))) ### tamanhos diferentes sapply(c("AA", "AIG", "AXP"), function(x) tail(names(get(x)))) COR <- NULL for (i in 2:394) { AA_AIG <- cor(AA[,i], AIG[,i], use = "pairwise.complete.obs") AA_AXP <- cor(AA[,i], AXP[,i], use = "pairwise.complete.obs") AIG_AXP <- cor(AIG[,i], AXP[,i], use = "pairwise.complete.obs") row <- data.frame(i, AA_AIG, AA_AXP, AIG_AXP) COR <- rbind(COR, row) } head(COR) ### </code> Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] 2015-03-18 13:38 GMT-04:00 Michelle Bau Graczyk <mbgraczyk(a)gmail.com>: > Boa tarde, > > Eu tenho tres matrizes de tres empresas diferentes que estao dispostas da > seguinte maneira: > > Empresa AA > data 9:31 > data > 09:31 09:32 09:33 09:34 09:35 09:36 09:37 09:38 09:39 09:40 09:41 09:42 > 09:43 09:44 09:45 09:46 09:47 > > > > > > > > > > > > > > > > > > > > > > 2004-07-01 700 200 7000 9400 1300 1500 6300 1800 4300 1300 4800 3900 2500 > 4300 4300 > > > > > > > 2004-07-02 NA NA 68900 2800 4900 1200 12300 NA NA 6500 NA 2900 8900 11600 > 26900 > > > > > > > 2004-07-06 800 4400 13700 10700 4700 4000 18700 9000 2100 2300 700 8900 > 10500 3900 8400 > > > > > > > 2004-07-07 NA 5600 1000 7300 5600 11200 3900 5800 5700 4500 34400 4100 > 500 15300 500 > > > > > > > 2004-07-08 NA NA NA NA NA NA NA NA NA NA 154100 42400 16700 4700 5900 > > > > > > > 2004-07-09 NA 96900 1000 20400 18300 18000 22500 4100 6000 1000 1900 3700 > 600 2900 14600 > > > > > > > 2004-07-12 1400 100 4500 3000 1300 4200 6500 4500 1600 1600 5600 7600 > 1900 4600 8300 > > > > > > > 2004-07-13 NA 5600 3400 16900 9000 5700 16500 1800 21000 6300 5900 12500 > 51000 7500 13100 > > > > > > > 2004-07-14 NA NA NA NA NA NA 164400 35000 41600 56600 27800 36500 55300 > 29900 13100 > > > > > > > 2004-07-15 NA 8100 2000 6100 4900 4500 2500 1400 5200 4100 5400 4000 > 13400 3200 11300 > > > > > > > 2004-07-16 282300 18100 30400 NA 15900 40800 9000 5300 5200 27600 7600 > 12100 8400 1000 8400 > > > > > > > 2004-07-19 28700 1700 2900 9500 400 2800 1100 6200 11500 17900 11100 > 16000 6800 5200 7100 > > > > > > > > Onde as datas vao de 2004-07-01 ate 2004-12-31(ao todo 128) e os horarios > vao de 09:31 ate 16:04( ao todo 397) > > Eu quero calcular a correlacao entre essas tres empresas para cada hora > fixa,ou seja, para a empresa AA eu calculo com ela mesma e com a AIG e a > AXP para cada hora fixa, e por isso montei este programa: > > A<-read.table(file="/home/usuario/Documentos/AAparaCalculodaMedia.txt", > header=TRUE) > AA.<-as.matrix(AA) > ##write.table(AA.,"AA..txt") > > AXP<-read.table("/home/usuario/Documentos/AXPparaCalculodaMedia.txt",header=TRUE) > AXP.<-as.matrix(AXP) > ##write.table(AXP.,"AXP..txt") > > AIG<-read.table("/home/usuario/Documentos/AIGparaCalculodaMedia.txt",header=TRUE) > AIG.<-as.matrix(AIG) > > Symbols<-c("AA.","AIG.","AXP.") > matrizAA<-matrix() > j<-1 > for (j in 1:length(Symbols)){ > empresa<-paste0(Symbols)#,"[",j,"]") > > for (i in 2:397){ > empresa1<-as.numeric(paste0(empresa,"[,",i,"]")) > head(empresa1) > matrizAA[,i-1]<-cor(as.numeric(AA.[,i]), as.numeric(empresa1), > use="pairwise.complete.obs") > > } > > } > > Mas da o erro: Erro em cor(as.numeric(AA.[, i]), as.numeric(empresa1), use > = "pairwise.complete.obs") : > incompatible dimensions > Além disso: Mensagens de aviso perdidas: > NAs introduzidos por coerção > > Eu achava que usando o use = "pairwise.complete.obs" as correlacoes seriam > calculadas apenas entre pares de numeros, mas acho que isso nao esta > acontecendo. > Os tres arquivos vou deixar em dropbox. > > Muito obrigada, > > Michelle > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. >

3 5

GLM
by ana paula coelho madeira 19 Mar '15

19 Mar '15

Boa tarde! Tenho um experimento que visa analisar como transgenia e controle químico (tratamentos) afetam número de parasitóides (Chelonus). Construí o modelo: m1<-glm.nb(n_parasit_Chelonus~trat) > summary(m1) Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 2.54553 0.16046 15.864 < 2e-16 *** tratC -0.24295 0.23851 -1.019 0.3084 tratD 0.17905 0.21973 0.815 0.4152 tratF -21.84812 4713.30843 -0.005 0.9963 tratG 0.46262 0.21032 2.200 0.0278 * tratI -1.22378 0.31394 -3.898 9.69e-05 *** tratJ 0.12862 0.22165 0.580 0.5617 tratL -21.84812 4713.30843 -0.005 0.9963 tratM -0.40547 0.24759 -1.638 0.1015 tratN 0.28768 0.21586 1.333 0.1826 tratP -3.23868 0.72931 -4.441 8.96e-06 *** tratQ -0.04001 0.22868 -0.175 0.8611 tratS -21.84812 4713.30843 -0.005 0.9963 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for Negative Binomial(40.7143) family taken to be 1) Null deviance: 376.329 on 51 degrees of freedom Residual deviance: 39.605 on 39 degrees of freedom AIC: 234.09 Como posso interpretar esses resultados? Agradeço pela ajuda. Ana

2 1

Re: [R-br] Valores repetidos no eixo de um gráfico
by Paulo Pimenta 19 Mar '15

19 Mar '15

Boa saída esta apresentada, mas a série de dados por ser dia varia bastante e não é possível realizar manualmente. Existe algum comando que pega dados contendo os dias (de 1 a 31 e alguns de 1 a 29 ou 28, dependendo do ano) e plota no eixo "x" ou "y" esta mesmos valores mantendo sua devida repetição, ou seja, 0 1 2 3 4 5 6 ... 31 1 2 3 4 5 6 7 8 9 ... 28 ... 1 2 3 4 5 6 7 8 9 ... 30? Abraços. Paulo Henrique Pimenta. > ------------------------------ > > Message: 11 > Date: Thu, 19 Mar 2015 06:04:02 +0000 > From: Paulo Pimenta <paulopimenta6(a)hotmail.com> > To: "=?utf-8?Q?r-br(a)listas.c3sl.ufpr.br?=" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Valores repetidos no eixo de um gráfico > Message-ID: <BAY404-EAS22211DD7A124D389B0F770AF4010(a)phx.gbl> > Content-Type: text/plain; charset="utf-8" > > Olá Listeiros, > > > Possuo uma dúvida que consiste no seguinte: Possuo uma série de dados que estão listadas por dias ordenados da forma: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,? sendo que estes correspondem cada a um valor de uma determinada variável física. > > O problema é que ao plotar o gráfico os valores em ?x? não são mostrados de forma seguida, ou seja: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,? mas sim como de 1 a 10 e os valores de cada dia aparecem acumulados. > > Como plotar um eixo com os valores seguidos? > > Abraços. > > Paulo Henrique Pimenta. > > > > > > > Sent from Windows Mail > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150319/906996c9/atta…> > > ------------------------------ > > Message: 12 > Date: Thu, 19 Mar 2015 08:01:32 -0300 > From: "Luis G. S. e Silva" <lgsilvaesilva(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Valores repetidos no eixo de um gráfico > Message-ID: > <CALruoAF9DBBfmYw-zEjDoNZ_QMGNTb76rQnHtU7aq5aS3CECKQ(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > x <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10) > y <- rnorm(length(x)) > plot(y, xaxt = 'n') > axis(1, at = 1:length(x), labels = x) > > 2015-03-19 3:04 GMT-03:00 Paulo Pimenta <paulopimenta6(a)hotmail.com>: > > > Olá Listeiros, > > > > Possuo uma dúvida que consiste no seguinte: Possuo uma série de dados que > > estão listadas por dias ordenados da forma: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, > > 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,? sendo que estes correspondem cada a um > > valor de uma determinada variável física. > > > > O problema é que ao plotar o gráfico os valores em ?x? não são mostrados > > de forma seguida, ou seja: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, > > 7, 8, 9, 10,? mas sim como de 1 a 10 e os valores de cada dia aparecem > > acumulados. > > > > Como plotar um eixo com os valores seguidos? > > > > Abraços. > > > > Paulo Henrique Pimenta. > > > > Sent from Windows Mail

3 2

Valores repetidos no eixo de um gráfico
by Paulo Pimenta 19 Mar '15

19 Mar '15

Olá Listeiros, Possuo uma dúvida que consiste no seguinte: Possuo uma série de dados que estão listadas por dias ordenados da forma: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,… sendo que estes correspondem cada a um valor de uma determinada variável física. O problema é que ao plotar o gráfico os valores em “x” não são mostrados de forma seguida, ou seja: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,… mas sim como de 1 a 10 e os valores de cada dia aparecem acumulados. Como plotar um eixo com os valores seguidos? Abraços. Paulo Henrique Pimenta. Sent from Windows Mail

2 1

Construção de Gráficos
by Thiago Touya 18 Mar '15

18 Mar '15

Boa noite! Gostaria de saber se é possível a criação de gráficos, por exemplo, de densidades de distribuições de probabilidade, sem a necessidade de um vetor de valores X, ou seja, se é possível criar gráficos apenas a partir das restrições de parâmetros e suporte. Se for possível, gostaria de ver um exemplo. Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

3 2

Duvida modo de uso da função paste
by Michelle Bau Graczyk 18 Mar '15

18 Mar '15

Boa tarde, estou com uma dúvida de como usar a função paste ou se tem como eu fazer de outro modo. Eu tentei primeiramente assim: matrizAA<-matrix(nrow=395, ncol=31) Symbols<-c("AA","AIG","AXP") for(j in 2:3) { empresa<-Symbols[j] for (i in 2:4){ matrizAA[i-1,j]<-cor(AA[,i],Symbols[j][,i], use="pairwise.complete.obs") } } write.table(matrizAA,"mAAteste23.txt") Entâo surge o erro: Erro em show(symbols[j][, 2]) : erro na avaliação do argumento 'object' na seleção do método para a função 'show': 'Erro em symbols[j] : objeto de tipo 'closure' não possível dividir em subconjuntos Ai tentei usar a funçao paste0 desta maneira: matrizAA<-matrix(nrow=395, ncol=31) Symbols<-c("AA","AIG","AXP") for(j in 1:3) { empresa<-Symbols[j] for (i in 2:4){ empresa1<-paste0(empresa,"[,",i,"]") ou empresa1<-paste0(empresa) matrizAA[i-1,j]<-cor(AA[,i],empresa1, use="pairwise.complete.obs") ou matrizAA[i-1,j]<-cor(AA[,i],empresa1[,i], use="pairwise.complete.obs") } } write.table(matrizAA,"mAAteste23.txt") e os erros são: Erro em cor(AA[, i], empresa1, use = "pairwise.complete.obs") : 'y' deve ser numérico ou Erro em empresa1[, i] : número incorreto de dimensões O que me leva a crer que usando o paste ele passa a não mais entender como um vetor AIG[,i] por exemplo, e sim por uma string. Alguem poderia me dizer como posso resolver isso? Muuuito obrigada, Michelle

1 0

Michelle Baú Gaczyk compartilhou "AXPparaCalculodaMedia.txt" com você.
by Michelle Baú Gaczyk via Dropbox 18 Mar '15

18 Mar '15

Michelle usou o Dropbox para compartilhar um arquivo com você! Clique aqui para visualizar.: https://www.dropbox.com/l/iUzfT7j35DiOHPtoHtQPvi?text=1

1 0

Michelle Baú Gaczyk compartilhou "AIGparaCalculodaMedia.txt" com você.
by Michelle Baú Gaczyk via Dropbox 18 Mar '15

18 Mar '15

Michelle usou o Dropbox para compartilhar um arquivo com você! Clique aqui para visualizar.: https://www.dropbox.com/l/oqyecuqKWf5ydyT1wHijAg?text=1

1 0

Michelle Baú Gaczyk compartilhou "AAparaCalculodaMedia.txt" com você.
by Michelle Baú Gaczyk via Dropbox 18 Mar '15

18 Mar '15

Michelle usou o Dropbox para compartilhar um arquivo com você! Clique aqui para visualizar.: https://www.dropbox.com/l/c2Ba0YfWF90baTnAarcYzg?text=1

1 0

Calculo de correlcao entre coluas que apresentam NAS de varias matrizes
by Michelle Bau Graczyk 18 Mar '15

18 Mar '15

Boa tarde, Eu tenho tres matrizes de tres empresas diferentes que estao dispostas da seguinte maneira: Empresa AA data 9:31 data 09:31 09:32 09:33 09:34 09:35 09:36 09:37 09:38 09:39 09:40 09:41 09:42 09:43 09:44 09:45 09:46 09:47 2004-07-01 700 200 7000 9400 1300 1500 6300 1800 4300 1300 4800 3900 2500 4300 4300 2004-07-02 NA NA 68900 2800 4900 1200 12300 NA NA 6500 NA 2900 8900 11600 26900 2004-07-06 800 4400 13700 10700 4700 4000 18700 9000 2100 2300 700 8900 10500 3900 8400 2004-07-07 NA 5600 1000 7300 5600 11200 3900 5800 5700 4500 34400 4100 500 15300 500 2004-07-08 NA NA NA NA NA NA NA NA NA NA 154100 42400 16700 4700 5900 2004-07-09 NA 96900 1000 20400 18300 18000 22500 4100 6000 1000 1900 3700 600 2900 14600 2004-07-12 1400 100 4500 3000 1300 4200 6500 4500 1600 1600 5600 7600 1900 4600 8300 2004-07-13 NA 5600 3400 16900 9000 5700 16500 1800 21000 6300 5900 12500 51000 7500 13100 2004-07-14 NA NA NA NA NA NA 164400 35000 41600 56600 27800 36500 55300 29900 13100 2004-07-15 NA 8100 2000 6100 4900 4500 2500 1400 5200 4100 5400 4000 13400 3200 11300 2004-07-16 282300 18100 30400 NA 15900 40800 9000 5300 5200 27600 7600 12100 8400 1000 8400 2004-07-19 28700 1700 2900 9500 400 2800 1100 6200 11500 17900 11100 16000 6800 5200 7100 Onde as datas vao de 2004-07-01 ate 2004-12-31(ao todo 128) e os horarios vao de 09:31 ate 16:04( ao todo 397) Eu quero calcular a correlacao entre essas tres empresas para cada hora fixa,ou seja, para a empresa AA eu calculo com ela mesma e com a AIG e a AXP para cada hora fixa, e por isso montei este programa: A<-read.table(file="/home/usuario/Documentos/AAparaCalculodaMedia.txt", header=TRUE) AA.<-as.matrix(AA) ##write.table(AA.,"AA..txt") AXP<-read.table("/home/usuario/Documentos/AXPparaCalculodaMedia.txt",header=TRUE) AXP.<-as.matrix(AXP) ##write.table(AXP.,"AXP..txt") AIG<-read.table("/home/usuario/Documentos/AIGparaCalculodaMedia.txt",header=TRUE) AIG.<-as.matrix(AIG) Symbols<-c("AA.","AIG.","AXP.") matrizAA<-matrix() j<-1 for (j in 1:length(Symbols)){ empresa<-paste0(Symbols)#,"[",j,"]") for (i in 2:397){ empresa1<-as.numeric(paste0(empresa,"[,",i,"]")) head(empresa1) matrizAA[,i-1]<-cor(as.numeric(AA.[,i]), as.numeric(empresa1), use="pairwise.complete.obs") } } Mas da o erro: Erro em cor(as.numeric(AA.[, i]), as.numeric(empresa1), use = "pairwise.complete.obs") : incompatible dimensions Além disso: Mensagens de aviso perdidas: NAs introduzidos por coerção Eu achava que usando o use = "pairwise.complete.obs" as correlacoes seriam calculadas apenas entre pares de numeros, mas acho que isso nao esta acontecendo. Os tres arquivos vou deixar em dropbox. Muito obrigada, Michelle

1 0

Merge com mais de 2 bancos
by Vanessa Felix 18 Mar '15

18 Mar '15

Pessoal, eu consigo juntar 3 ou 4 bancos de uma vez só? Não queria ter que fazer vários merges... Abç, -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestranda em Saúde Coletiva - IMS/UERJ Assistente de pesquisa I - Instituto de Pesquisa Econômica Aplicada – IPEA (32) 9142-8428 (21) 97557-3956 -------------------------------------------------------------------------------------------

5 6

Extrair linhas específicas de um data.frame
by Alexandre Loures 17 Mar '15

17 Mar '15

Boa tarde! Tenho um arquivo em que uma das variáveis tem o seguinte formato: data 1958Q3 2001Q1 1989M7 1948M12 1949 2011 . . . E assim por diante. Ou seja, tenho informações quadrimestrais, mensais e anuais. Todavia, gostaria de obter um novo data.frame apenas com os dados anuais. Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures* Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE e-Mail: alexandre.loures(a)ymail.com Site: www.ccsa.ufpb.br/ppge

6 7

loop lento
by Roney Fraga Souza 17 Mar '15

17 Mar '15

Olá pessoal, Eu tenho uma lista (ver2) composta por com 14 elementos, cada elemento é um data.frame() contendo as mesmas variáveis, o que varia é apenas a quantidade de linhas. Cada elemento da lista é um ano específico do data.frame(), exemplo, ver2[[1]] é o ano de 1977, ver2[[2]] é 1978, ver2[[3]] é 1979 até o último ano ver2[[14]] que é 1990. Nos data.frame() tenho informações sobre artigos científicos: nome (código que não se repete), grupo (extraído via processo de clusterização), PY (ano que o artigo foi publicado) e one (variável apenas com valor 1 que uso em outros procedimentos). Os data.frame() são cumulativos, no ano de 1978 estão todos os artigos publicados até 1978, em 1979 estão todos os artigos publicados até 1978 mais os artigos publicados em 1979, no em de 1990 estão todos os artigos publicados até 1990. O procedimento de clusterização é realizado de modo independente para cada ano, logo o artigo 'Baker, 1976, V4, P5' pode estar em grupos diferentes para cada ano. Quero saber a origem dos artigos que estão em um determinado grupo. Exemplo, grupo 1 de 1990, é composto por 100 artigos, desses 25 faziam parte do grupo 1 em 1989, 30 do grupo 2 em 1989, 35 do grupo 3 em 1989, e os 10 restantes eu não tenho interesse, pois, foram publicados no ano de 1990 logo não existim nos anos anteriores. Segue exemplo do resultado que eu obtenho o código atual: label qtde tm1.grupo tm1.ano t.grupo t.ano 1-2 6 1 1989 2 1990 2-1 5 2 1989 1 1990 3-1 3 3 1989 1 1990 4-5 3 4 1989 5 1990 5-4 2 5 1989 4 1990 6-9 2 6 1989 9 1990 7-6 2 7 1989 6 1990 8-8 2 8 1989 8 1990 9-10 2 9 1989 10 1990 qtde = quantidade de artigos em cada grupo tm1.grupo = ano t menos 1 determinado grupo tm1.ano = ano t menos 1 t.grupo = grupo t.ano = ano Para chegar a tal resultado eu utilizo o seguinte código: # ------------------------------ # início do código # ------------------------------ # baixar o arquivo 'ver2.rds' nesse link: # https://db.tt/13dT9XQI <https://db.tt/13dT9XQI> ver2 <- readRDS("~/Downloads/ver2.rds") grupo <- list() for(k in 2:length(ver2)){ grp <- sort(unique((ver2[[k-1]]$grupo))) ### qtde de grupos no ano t-1 grp2 <- sort(unique((ver2[[k]]$grupo))) ### qtde de grupos no ano t RES <- LAB <- list() for (i in grp) { for (j in grp2) { RES <- append(RES,list(length(intersect(ver2[[k-1]]$name[ver2[[k-1]]$grupo==i], ver2[[k]]$name[ver2[[k]]$grupo==j])))) LAB <- append(LAB, paste(i,j, sep='-')) } } grupo[[k]] <- data.frame(label=sapply(LAB, "["), qtde=sapply(RES, sum)) grupo[[k]] <- subset(grupo[[k]],qtde>0) grupo[[k]]$tm1.grupo <- as.numeric(gsub('-.*','',grupo[[k]]$label)) grupo[[k]]$tm1.ano <- as.numeric(max(ver2[[k]]$PY)-1) grupo[[k]]$t.grupo <- as.numeric(gsub('^.-','',grupo[[k]]$label)) grupo[[k]]$t.ano <- as.numeric(max(ver2[[k]]$PY)) } # grupo # ------------------------------ # fim do código # ------------------------------ Meu problema é: esse código funciona mas é muito lento, muito mesmo, inviabilizando sua aplicação para o volume de dados que trabalho no dia a dia. Alguém conhece uma forma mais rápida de fazer isso. Obs.: parte desse código foi contribuição de Eder Comunello através dessa lista, mas naquela altura o problema estava estruturado de outra forma. Abraço Roney

4 7

Uso do laço for para impessão de mais de um arquivo gerado.
by Michelle Bau Graczyk 16 Mar '15

16 Mar '15

Boa noite, Com ajuda de colegas do grupo montei um programa para pegar dados do yahoo e organiza-los. Porém, eu preciso pegar para vårias empresas e por isso fiz um laço com o uso da função for para imprimir a saida. Porem, eu não estou conseguindo fazer com que ele imprima todos e sim, somente o ultimo elemento. No caso abaixo, ele só imprime e salva o arquivo referente a empresa DIS enquanto eu também quero das outras. Já tentei colocar o write.table no final da última chaves mas também não imprimiu de todos os arquivos. Muito obrigada! Michelle symbol<-c("XOM","AAPL", "DIS") for(i in seq_along(symbols)){ symbol <- symbols[i] URL <- paste0("http://chartapi.finance.yahoo.com/instrument/1.0/", symbol, "/chartdata;type=quote;range=1d/csv") tmp <- readLines(URL) tabela<-tmp[18:length(tmp)] tabela_split<-strsplit(tabela,",") # matriz para receber os dados dados<-matrix(NA,ncol=length(tabela_split[[1]]),nrow=length(tabela_split)) for(j in 1:length(tabela_split)) { dados[j,]<-as.numeric(tabela_split[[j]]) } # organizar os nomes de colunas linhas dimnames(dados)<-list(paste(symbol,1:length(tabela_split)),paste(c("timestamp","close","high","low","open","volume"),1:length(tabela_split[[1]]))) SAIDA<-paste0("Dados",symbol,"yahoo20150312.txt") write.table(dados,file= SAIDA) }

6 11

Render {rmarkdown} não funciona em emacs
by Fernando Antonio de souza 15 Mar '15

15 Mar '15

Caros amigos Preciso rodar um arquivo *.rmd e gostaria de utilizar a função Render do pacote rmarkdown. Acontece que ao utilizá-la ela retorna a seguinte mensagem "Error: pandoc version 1.12.3 or higher is required and was not found." Eu já atualizee a versão do pandoc mas a mensagem continau. Fiz uma busca na internte e encontrei esta postagem do Walmes que teve o mesmo problema http://r.789695.n4.nabble.com/Apply-rmarkdown-render-outside-the-RStudio-do… Segundo o Walmes ele solucionou criando um soft-link. Sou usuário novo em linux. Como posso fazer isso? abraços

1 0

Knitr-emacs dúvida
by Fernando Antonio de souza 13 Mar '15

13 Mar '15

Caros amigos Estou criando um relatório utilizando as funcionalidades do knitr+ html+markdown. Eu desejo após gerar o arquivo *.rmd convertê-lo em documento word/odt utilizando a função pandoc{knitr}. O meu problema é o seguinte. ao fazer a conversão para docx/odt o código chunk é impresso de forma literal no documento. No chunk abaixo mesmo a especificação "echo=FALSE" não foi atendida. O meu objetivo é que no documento doc/odt apareça apenas a tabela gerada pela função kable{knitr} e não o código utilizado para gerá-lo Utilizo o editor emacs+ess. Rodei os comandos no RStudio e nele o comando funcionou perfeitamente. Por isso acredito que seja algo do emacs que desconheço. Para trabalhar com knitr instalei o pacote polymode. ```{r,echo=FALSE} opts_chunk$set(echo=FALSE,message=FALSE,results = "asis") ``` ```{r} #install.packages("repmis") #install.packages("knitr") library(pander) library(repmis) #funçao para importar dados do dropbox library(knitr) dados <- source_DropboxData(file = "DQL.csv",key = "xtp04eaa1316kt2", sep = ";", header = TRUE) dados <- transform(dados,VACA=factor(VACA),PERIODO=factor(PERIODO),TRAT=factor(TRAT)) modelo <- aov(LEITE~VACA+PERIODO+TRAT,data=dados) #install.packages("multcomp") #library(multcomp) comp <- glht(modelo,linfct=mcp(TRAT="Tukey")) nometrat <- colnames(summary(comp)$test$pvalues) medias<- aggregate(LEITE~TRAT,data=dados,FUN=mean) comparacao <-unname( cld(summary(comp))$mcletters$Letters) tabela <- cbind(medias,comparacao) colnames(tabela) <- c("TRATAMENTO","LEITE(Kg)","Tukey (<0,05)") ``` ```{r echo=FALSE,results="asis"} kable(tabela,format="markdown",digits=3,align="c") ```

1 0

Chamada de Trabalhos - VIII Reunião da ABAVE
by Leandro Marino 11 Mar '15

11 Mar '15

*VIII Reunião da ABAVE * *Avaliação de Larga Escala no Brasil: * *ensinamentos, aprendizagens e tendências* Florianópolis/SC – 19 a 21 de agosto de 2015 A Diretoria e os membros da Diretoria Científica da ABAVE convocam os interessados a enviar os trabalhos originais, contendo resultados de suas atividades de pesquisa científica ou de execução de projetos de avaliação educacional tanto da educação básica, quanto da educação superior, para apresentação durante a reunião. O formato da reunião foi pensado para favorecer a discussão e troca de experiências. Os autores deverão submeter seus artigos em duas categorias: - *Apresentação oral:* Serão aceitos textos que apresentem pesquisas com conclusão parciais ou finais, abordando temáticas novas ou já estabelecidas na área de avaliação educacional, preferencialmente utilizando os dados disponíveis das avaliações educacionais em larga escala, que evidenciem elaboração teórica e rigor conceitual e metodológico na análise. - *Apresentação pôster: *Serão aceitos textos que apresentem projetos de pesquisa em andamento ou relatos de experiências na área de avaliação educacional. Os trabalhos deverão ser submetidos eletronicamente até 30 de março de 2015. Em breve serão divulgados os procedimentos para a submissão eletrônica dos trabalhos. No ato da submissão, os próprios autores deverão classificar o trabalho em uma das categorias: apresentação oral ou pôster. *Regras para Submissão de Trabalhos:* Os artigos que serão submetidos à análise devem conter as seguintes especificações: - Deve estar no formato PDF, na fonte Arial, corpo 12. - Os artigos podem ser escritos em português, espanhol ou inglês. - A página deverá ser configurada em papel A4, utilizando os seguintes parâmetros: margem superior 2,5 cm; inferior 2,5 cm; lateral esquerda 3,5 cm; lateral direita 2,5 cm. - O trabalho deve conter apenas o título no corpo do texto. NENHUMA informação de autoria deve ser inserida no texto do artigo ou do seu resumo. - A primeira página do artigo deverá conter APENAS as seguintes informações de identificação: 1. Título do texto 2. Nome do autor – sigla da Instituição 3. Nome(s) do co/autor(es) – sigla da(s) instituição(es) (quando houver) 4. Agência ou Instituição Financiadora (quando houver) *Para Apresentação Ora**l* O texto deve conter no máximo 20 páginas, incluindo resumo, figuras, diagramas,bibliografias e anexos. Colocar resumo em português com, no máximo, 200 palavras, terminado pelas palavras-chave, em negrito. As notas de rodapé, em caso de serem indispensáveis, devem ser colocadas no final dotrabalho, antes das referências bibliográficas. Os títulos das seções internas devem estar em negrito e posicionado à margemesquerda. As figuras e tabelas devem ter legendas e estas, como as referências, devem seguir opadrão ABNT. *Para Apresentação de Pôster*  O texto deve conter no máximo 10 páginas, incluindo resumo, figuras, diagramas, bibliografias e anexos.  Colocar resumo em português com, no máximo, 200 palavras, terminado pelas palavras-chave, em negrito.  As notas de rodapé, em caso de serem indispensáveis, devem ser colocadas no final do trabalho, antes das referências bibliográficas.  Os títulos das seções internas devem estar em negrito e posicionado à margem esquerda.  As figuras e tabelas devem ter legendas e estas, como as referências, devem seguir o padrão ABNT. *Avaliação dos Trabalhos*  Os artigos submetidos serão avaliados pelos membros da Diretoria Científica e por pareceristas ad hoc indicados pela Diretoria.  Os trabalhos serão avaliados pelo método do double blind review, que possibilita a análise inominada dos artigos, garantindo a imparcialidade da avaliação.  Serão selecionados, para a apresentação na VIII Reunião Anual da ABAVE, 50 trabalhos na categoria Apresentação Oral e 70 na categoria Pôster.  Cabe à Comissão Científica a possibilidade de reclassificar o trabalho, a partir de critérios da natureza da pesquisa desenvolvida, como também por disponibilidade de espaço para alocação das sessões apresentação oral. *Datas Relevantes*  30 de março: data final para submissão de trabalhos completos, nas categorias apresentação oral e apresentação de pôster;  01 de abril a 15 de maio: julgamento dos trabalhos pela Comissão Científica;  20 de maio: Divulgação do resultado final sobre os trabalhos aprovados. *Outras Informações*  Os artigos aprovados e apresentados serão publicados eletronicamente nos anais do evento no sítio da ABAVE.  O autor que apresentará o artigo submetido deverá estar inscrito no evento para que o artigo seja considerado para apresentação. Outras Informações devem ser solicitadas à: Secretaria da ABAVE (21) 2103-9635 abave(a)abave.org.br

1 0

Re: [R-br] Shiny App - apenas uma dúvida sobre single-file
by Rafael Tieppo 11 Mar '15

11 Mar '15

Fernando, não tenho muita experiência com o shiny, fiz apenas um trabalho com o shiny, sendo este para um trabalho final do coursera. Não criei nenhum arquivo com o nome "run.r", apenas criei ui.r e o server.r . Espero ter ajudado. Saudações. Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: http://www2.unemat.br/rafaeltieppo blog: https://sistemasagricolas.wordpress.com "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente". On Wednesday, March 11, 2015 11:00 AM, "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> wrote: Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. [R -br] Shiny App - apenas uma dúvida sobre single-file (Fernando Gama) 2. Re: [R -br] Shiny App - apenas uma dúvida sobre single-file (Fernando Gama) ---------------------------------------------------------------------- Message: 1 Date: Tue, 10 Mar 2015 16:34:09 -0300 From: Fernando Gama <f.fabiogama88(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] [R -br] Shiny App - apenas uma dúvida sobre single-file Message-ID: <CAOwqRyu3QApkACvpTn1Xp-kVMkvh4URr9Sv3nDf4buJen_QBiQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Prezados, quem trabalhou/trabalha algum tempo com o shiny notou que recentemente foi incorporado no framework a possibilidade que temos de rodar a aplicação utilizando ui.r e server.r em um unico arquivo app.r . A duvida que tenho é: não precisamos mais do run.r para rodar a aplicação? Porque rodei o app.r e notei que a aplicação foi executada. Alguém poderia me tirar essa pequena duvida? -- Att, Fernando da Mata -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150310/5f64c984/atta…> ------------------------------ Message: 2 Date: Tue, 10 Mar 2015 18:37:14 -0300 From: Fernando Gama <f.fabiogama88(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] [R -br] Shiny App - apenas uma dúvida sobre single-file Message-ID: <CAOwqRysu9kaOPXbkZXMbeS_K5S=bFssPRg2b9aj7aSmVSG5mWw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Pessoal, antes que alguém tenha o trabalho de responder. Fiz alguns testes e encontrei uma diferença bem sutil: quando eu rodo o script ''app.R'' a aplicação é executada da mesma maneira quando rodamos o script run.r . Este, carrega o diretório que guarda todos os scripts e quando atualizamos o navegador tudo o que for alterado no script é atualizado no meu browser, caso eu ignore a execução do script run.r os valores alterados no script não são alterados no browser e serei obrigado a interromper a execução e rodar novamente o app.R :) Por uma questão de lógica, penso que o script 'não carrega o diretório na qual ele foi chamado' ele simplestemente executa o que está dentro dele e por isso, o browser nao entende a mudança, diferentemente do run.R ?? Obrigado de qualquer maneira, se alguem quiser acrescentar algo fique a vontade. Boa noite. Em 10 de março de 2015 16:34, Fernando Gama <f.fabiogama88(a)gmail.com> escreveu: > Prezados, > > quem trabalhou/trabalha algum tempo com o shiny notou que recentemente foi > incorporado no framework a possibilidade que temos de rodar a aplicação > utilizando ui.r e server.r em um unico arquivo app.r . > > A duvida que tenho é: não precisamos mais do run.r para rodar a aplicação? > Porque rodei o app.r e notei que a aplicação foi executada. Alguém poderia > me tirar essa pequena duvida? > > -- > Att, > > Fernando da Mata > > > -- Att, Fernando da Mata -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150310/44825689/atta…> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------ Fim da Digest R-br, volume 51, assunto 12 *****************************************

1 0

[R -br] Shiny App - apenas uma dúvida sobre single-file
by Fernando Gama 10 Mar '15

10 Mar '15

Prezados, quem trabalhou/trabalha algum tempo com o shiny notou que recentemente foi incorporado no framework a possibilidade que temos de rodar a aplicação utilizando ui.r e server.r em um unico arquivo app.r . A duvida que tenho é: não precisamos mais do run.r para rodar a aplicação? Porque rodei o app.r e notei que a aplicação foi executada. Alguém poderia me tirar essa pequena duvida? -- Att, Fernando da Mata

1 1

Exibir categoria quando não há elemento - summarize {Hmisc}
by Leandro Marino 09 Mar '15

09 Mar '15

Caros, boa noite. Tudo bem? Vejam o CMR: dados = data.frame(valor=rnorm(100),var=sample(LETTERS[c(1,2,3,5)],100,replace=T),peso=rpois(100,2)) dados[1:10,] dados$var <- factor(dados$var,levels=LETTERS[1:5]) table(dados$var) g1 <- function(y) wtd.mean(y[,1],y[,2]) summarize(dados[,c(1,3)], llist(var=dados$var), g1,stat.name = 'med') Quando faço o table, aparece que a categoria D não teve ninguém, ou seja 0. Como posso fazer para forçar ela aparecer no summarize? Tenho depois uma serie de rotinas que utilizo as funcoes wtd... do Hmisc. No entanto preciso que todas as categorias que estão no fator sejam sempre exibidas. obrigado pela ajuda. abs Leandro

4 4

nloptr no openSUSE 13.2 (64 bits)
by Vinícius Lionel Mateus 09 Mar '15

09 Mar '15

Prezados useRs, Venho enfrentando duras batalhas para instalar o nloptr, devido a problemas na instalação do nlopt. Por fim, consegui instalar o nlopt, mas ainda tenho problemas para instalar o nloptr. Seguindo a dica do autor do nloptr, executei a linha: - R CMD INSTALL --configure-args="--with-nlopt-cflags='-I/usr/local/include' --with-nlopt-libs='-L/usr/local/lib64 -lnlopt'" nloptr No entanto, o output é: - Warning: invalid package ‘nloptr’ Error: ERROR: no packages specified Alguém tem alguma dica? Abs, Vinícius -- Best regards, Vinícius Lionel Mateus, MSc (http://lattes.cnpq.br/6501001637020665) Chemistry PhD Student Atmospheric Chemistry Laboratory - Dep. Chemistry PUC-Rio - Pontifical Catholic University of Rio de Janeiro Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brazil CEP: 22451-900 Telefone: (+45) 26 28 28 51 (+55) (21) 3527-1327 (+55) (21) 993 - 588 - 051 Skype: vinicius.lionel http://www.qui.puc-rio.br/index.html

3 5

Somar intervalo de uma coluna
by Alexandre Loures 08 Mar '15

08 Mar '15

Tenho um data frame da seguinte forma: ano mês exp 1997 JAN 431 1997 FEV 550 1997 MAR 20 . . . 2014 JAN 220 2014 FEV 100 2014 MAR 80 Gostaria de somar os valores para cada um dos anos da base. Isto é, gostaria de obter os valores por ano e não por ano/mês. Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures* Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE e-Mail: alexandre.loures(a)ymail.com Site: www.ccsa.ufpb.br/ppge

2 2

Extrair coordenadas do polígono do Brasil - maptools
by Augusto Ribas 07 Mar '15

07 Mar '15

Ola pessoal, estou com uma duvida aqui, como eu retiro as coordenadas que formam o polígono do Brasil de um spatialpolygondataframe. No pacote maptools tem o wrld_simpl que é um mapa do mundo, um dos polígonos é o mapa do Brasil, como eu consigo somente os pontos, coordenadas na ordem para fazer o polígono, não estou conseguindo entender como extrai essa informação ###### library(maptools) data(wrld_simpl) which(wrld_simpl$NAME=="Brazil") plot(wrld_simpl[21,]) wrld_simpl[21,] ####### Eu gostaria ter os pontos para escrever polygon(pontos) ou plot(pontos,type="l") e desenhar o mapa em um plot igual acima, tem como fazer isso? -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

4 4

krige.conv
by Felinto COSTA 07 Mar '15

07 Mar '15

Prezados. Tenho dois modelos teóricos (exponencial e esférico ) que foram ajustados aos dados por máxima verossimilhança. As coordenadas numéricas (UTM) dos dados foram reduzidas pela subtração da coordenada média segundo cada eixo. Em ambos considero a existência de uma tendência polinomial de 2ª ordem na média. Eles foram ajustados sem problemas, assim como a validação cruzada de ambos pode ser feita. O problema está ao querer estimar algum ponto(s) usando a função krige.conv. Apenas o modelo exponencial é aceito. Ao usar o esférico a função retorna a seguinte mensagem: krige.conv: model with mean given by a 2nd order polynomial on the coordinates krige.conv: anisotropy correction performed krige.conv: performing the Box-Cox data transformation Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 1.76896e-16 Alguém saberia como contornar esse problema? Antecipadamente agradeço a colaboração Felinto COSTA

3 3

[R] Uso da função if
by Michelle Bau Graczyk 06 Mar '15

06 Mar '15

Boa tarde, Estou com um problema de lógica para o uso do if. Meu intuito com este programa é separar por horário os valores do volume. O horário está escrito de maneira numérica e vai de 1 até 396. O valor máximo de volumes por horário é 126. Só que quando tento fazer rodar, vem a mensagem : Erro em if (AA.data$hora[k] == i) { : valor ausente onde TRUE/FALSE necessário Abaixo segue o programa. AA<-read.table("/Users/bau/Documents/Mestrado/teste/AXP_NYS_USD_1min_EST5EDT_200407s12e.dat", header=FALSE, col.names=c("data","hora","preco","volume")) hora<-as.numeric(AA$hora) AA.data<-data.frame(data=AA$data,hora=hora, preco=AA$preco, volume=AA$volume) write.table(AA.data,"testeteste.txt") matriz.dadosAA<-matrix(nrow=128,ncol=4) for(i in 1:4){ j<-1 for(k in 1:48935){ for(j in 1:126){ if(AA.data$hora[k]==i){ matriz.dadosAA[j,i]=AA.data$volume[k] j=j+1 } else { matriz.dadosAA[j,i]=0 j=j+1} k=k+1 } } } a<-matriz.dadosAA[,1] write.table(a,"a.txt") Muito obrigada, Michelle

4 7

Extrair coordenadas de um raster dada uma condição
by ASANTOS 05 Mar '15

05 Mar '15

Caros Listeiros, Fiz a seleção da coordenada dos pixels (p.cen) em um raster que continham pontos de interesse (pontos) no CRM abaixo, gostaria de saber como impor uma condição para selecionar as coordenadas dos pixels que não contém nenhum ponto. Tentei criar um novo objeto com todas as coordenadas no raster e remover as coordenadas p.cen sem sucesso, alguém poderia me ajudar? Obrigado, ### <code r> require(raster); require(sp) ## Criando um raster r <- raster(nc=10, nr=10) r <- setValues(r, round(runif(ncell(r))* 255)) ##Pontos selecionados x <- c(44.25,33.06,63.22,71.35) y <- c(22.00,81.90,15.18,71.03) pontos <- SpatialPoints(cbind(x,y)) cells <- cellFromXY(r, pontos); cells coordinates(r)[cells,] p.cen <- xyFromCell(r, cells); p.cen plot(r) points(pontos, col="red") points(p.cen) values <- extract(r,p.cen) ##Extrai as coordendas dados<-cbind(p.cen,cells,values) dados ### </code> -- ====================================================================== 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 4

Re: [R-br] [R] Uso da função if
by Michelle Bau Graczyk 05 Mar '15

05 Mar '15

Oi Fernando e Rafael , eu acho que, pelo que eu vi, se eu só mandat o dput nao vai ficar tão claro. Por isso vou mandar os dados aqui em baixo, espero que você não ache ruim. Porque assim acho que fica mais claro. Eu preciso separar por horário os valores de volume. Assim eu quero gerar uma matriz as quais as colunas vao de 09:31 até 16:03 com as linhas sendo preenchidas pelos volumes do certo dia. O horário eu deixei em numero que vai de 1 ate a 396 ao invés de 09:31 até 16:03 porque achei mais fácil, ainda tenho dificuldades com character. Bom, seguem os dados! Muito obrigada! 01.07.2004 09:31:00 51.37 27900 01.07.2004 09:32:00 51.32 100 01.07.2004 09:33:00 51.45 700 01.07.2004 09:34:00 51.46 1000 01.07.2004 09:35:00 51.50 800 01.07.2004 09:36:00 51.54 3000 01.07.2004 09:37:00 51.54 1600 . . . 31.12.2004 15:45:00 56.50 4100 31.12.2004 15:46:00 56.50 7400 31.12.2004 15:47:00 56.40 8400 31.12.2004 15:48:00 56.44 15800 31.12.2004 15:49:00 56.41 4700 31.12.2004 15:50:00 56.44 6900 31.12.2004 15:51:00 56.42 9800 31.12.2004 15:52:00 56.43 8700 31.12.2004 15:53:00 56.43 7400 31.12.2004 15:54:00 56.40 12400 31.12.2004 15:55:00 56.44 11600 31.12.2004 15:56:00 56.45 8900 31.12.2004 15:57:00 56.44 14300 31.12.2004 15:58:00 56.44 16600 31.12.2004 15:59:00 56.43 21100 31.12.2004 16:00:00 56.43 34500 31.12.2004 16:02:00 56.37 46000 31.12.2004 16:03:00 56.37 111400 Em 4 de março de 2015 13:52, Michelle Bau Graczyk <mbgraczyk(a)gmail.com> escreveu: > Oi Fernando e Rafael > , eu acho que, pelo que eu vi, se eu só mandat o dput nao vai ficar tão > claro. Por isso vou mandar os dados em anexo, espero que você não ache > ruim. Porque assim acho que fica mais claro. > Eu preciso separar por horário os valores de volume. Assim eu quero gerar > uma matriz as quais as colunas vao de 09:31 até 16:03 com as linhas sendo > preenchidas pelos volumes do certo dia. O horário eu deixei em numero que > vai de 1 ate a 396 ao invés de 09:31 até 16:03 porque achei mais fácil, > ainda tenho dificuldades com character. > > Bom, seguem os dados em anexo. > > Muito obrigada! > > Em 4 de março de 2015 13:34, Fernando Antonio de souza < > nandodesouza(a)gmail.com> escreveu: > > Michelle é que para que os usuários da lista possam lhe ajudar é >> necessário que você forneça um CMR(código mínimo reproduzível) isso permite >> que possamos acessar seus dados e refazer os seus comandos e assim lhe >> orientar melhor. >> >> Você forneceu o código, no entanto eu não tenho como acessar seus dados, >> pois o endereço que vc forneceu no read.table não é o mesmo em meu >> computador. >> >> Por isso estou lhe solicitando que faça um dput(seusbancodedados). Após >> dar o enter esta função retornara seu banco de dados em um formato que >> todos os usuários poderam ter acesso apenas por copiar e colar esta saída >> no R. Você copia e cole no email essa saída e assim todos poderão utilizar >> seus dados e lhe ajudar de maneira mais rápida >> >> Você pode fornecer todo o banco de dados ou apenas parte dele (no caso de >> banco de dados muito grande). O importante é que os dados fornecidos possam >> ser reproduzidos pela pessoa que está lhe ajudando. >> >> Em 4 de março de 2015 13:25, Michelle Bau Graczyk <mbgraczyk(a)gmail.com> >> escreveu: >> >> Oi Fernando, >>> >>> Me desculpe mas você pode ser um pouco mais claro? Sou iniciante em R e >>> por isso ainda tenho muita dificuldade. >>> >>> Em 4 de março de 2015 13:20, Fernando Antonio de souza < >>> nandodesouza(a)gmail.com> escreveu: >>> >>> OLÁ >>>> >>>> Faça um dput(seubancodedados) copie e cole a saída para que seus >>>> comandos possa ser replicados e possa lhe ajudar. Acredito que o que você >>>> deseja fazer pode ser feito com funções do pacote plyr ou lapply, o que >>>> torna muito mais eficiente. >>>> >>>> Em 4 de março de 2015 13:09, Michelle Bau Graczyk <mbgraczyk(a)gmail.com> >>>> escreveu: >>>> >>>>> Boa tarde, >>>>> >>>>> Estou com um problema de lógica para o uso do if. Meu intuito com este >>>>> programa é separar por horário os valores do volume. O horário está escrito >>>>> de maneira numérica e vai de 1 até 396. O valor máximo de volumes por >>>>> horário é 126. >>>>> Só que quando tento fazer rodar, vem a mensagem : >>>>> Erro em if (AA.data$hora[k] == i) { : >>>>> valor ausente onde TRUE/FALSE necessário >>>>> >>>>> Abaixo segue o programa. >>>>> >>>>> AA<-read.table("/Users/bau/Documents/Mestrado/teste/AXP_NYS_USD_1min_EST5EDT_200407s12e.dat", >>>>> header=FALSE, col.names=c("data","hora","preco","volume")) >>>>> hora<-as.numeric(AA$hora) >>>>> AA.data<-data.frame(data=AA$data,hora=hora, preco=AA$preco, >>>>> volume=AA$volume) >>>>> write.table(AA.data,"testeteste.txt") >>>>> >>>>> matriz.dadosAA<-matrix(nrow=128,ncol=4) >>>>> for(i in 1:4){ >>>>> >>>>> j<-1 >>>>> for(k in 1:48935){ >>>>> for(j in 1:126){ >>>>> if(AA.data$hora[k]==i){ >>>>> >>>>> >>>>> matriz.dadosAA[j,i]=AA.data$volume[k] >>>>> j=j+1 >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> } else { >>>>> >>>>> matriz.dadosAA[j,i]=0 >>>>> j=j+1} >>>>> >>>>> k=k+1 >>>>> >>>>> } >>>>> } >>>>> } >>>>> a<-matriz.dadosAA[,1] >>>>> write.table(a,"a.txt") >>>>> >>>>> >>>>> >>>>> Muito obrigada, >>>>> >>>>> Michelle >>>>> >>>>> _______________________________________________ >>>>> R-br mailing list >>>>> R-br(a)listas.c3sl.ufpr.br >>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>>> código mínimo reproduzível. >>>>> >>>> >>>> >>>> _______________________________________________ >>>> R-br mailing list >>>> R-br(a)listas.c3sl.ufpr.br >>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>> código mínimo reproduzível. >>>> >>> >>> >>> _______________________________________________ >>> R-br mailing list >>> R-br(a)listas.c3sl.ufpr.br >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>> código mínimo reproduzível. >>> >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> > >

3 7

Michelle Baú Gaczyk compartilhou "AXP_NYS_USD_1min_EST5EDT_200407s12e.dat" com você.
by Michelle Baú Gaczyk via Dropbox 04 Mar '15

04 Mar '15

Michelle usou o Dropbox para compartilhar um arquivo com você! Clique aqui para visualizar.: https://www.dropbox.com/l/Hixo6d02ckZPbWnEbH54yh?text=1

1 0

como salvar gráfico após o iso do comando identify
by Antonio Silva 03 Mar '15

03 Mar '15

Olá a todos Quero salvar um gráfico após o uso do comando identify mas não consegui. Utilizo o comando identify para adicionar etiquetas aos pontos de forma manual, para que os números não fiquem embolados ou sobre as linhas. Exemplo x<-c(1,2,3,4,5,6) y<-c(20,30,15,7,25,40) plot(x,y,type="b",ylim=c(0,45)) identify(x,y,labels=y) No entanto, quando escrevo png(filename="teste.png",width=20,height=20,units="cm",res=300) x<-c(1,2,3,4,5,6) y<-c(20,30,15,7,25,40) plot(x,y,type="b",ylim=c(0,45)) identify(x,y,labels=y) dev.off() o gráfico é gravado sem aparecer na tela e consequentemente sem a identificação dos pontos. Utilizo o R integrado ao gedit no UBUNTU Trusty Tahr 64 bits. Gostaria de fazer o procedimento pela linha de comando, sem a utilização de GUIs. Há alguma forma de gravar o gráfico com os pontos identificados? Obrigado! Antônio Olinto

3 3

Mover os vértices de um polígono
by ASANTOS 03 Mar '15

03 Mar '15

Caros Listeiros, Criei um polígono (c1) abaixo com os dados dos vértices da área em UTM e gostaria de saber se alguém conhece alguma função para mover todos os vértices 30m para dentro e reduzir a área do polígono? #Contorno da área -------------------------------------------------------------- x.coords1<-c(371299.9, 371266.4, 371205.6, 371111.8, 371047.6, 371018.2, 371014.0, 371009.3, 370983.1, 370919.7, 370853.6, 370785.6, 370748.8, 370711.8, 370687.8, 370696.4, 370785.9, 370885.5, 371035.8, 371148.1, 371205.2, 371231.7, 371236.5, 371240.3, 371285.8, 371326.5, 371397.2, 371417.1, 371432.9, 371445.0, 371455.7, 371466.4, 371476.6, 371502.6, 371536.0, 371550.0, 371546.8, 371528.3, 371470.0, 371393.3, 371299.9, 371299.9) y.coords1<-c(8246589, 8246560, 8246508, 8246428, 8246373, 8246349, 8246348, 8246352, 8246385, 8246465, 8246551, 8246638, 8246685, 8246732, 8246764, 8246771, 8246846, 8246932, 8247062, 8247160, 8247209, 8247230, 8247224, 8247221, 8247160, 8247107, 8247016, 8246991, 8246967, 8246939, 8246914, 8246892, 8246875, 8246846, 8246821, 8246809, 8246802, 8246785, 8246735, 8246669, 8246589, 8246589) # Cria um poligono t05<-unique(data.frame(x=rev(x.coords1),y=rev(y.coords1))) c1 = cbind (t05$x, t05$y) r1 = rbind (c1, c1[1, ]) # junta plot(c1) polygon(c1) # -- ====================================================================== 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

2 1

Gerar coluna a partir das informações da 1ª Coluna
by Alexandre Loures 02 Mar '15

02 Mar '15

Boa tarde! Tenho uma base cuja primeira coluna é a seguinte: ARG BRA FRA A partir dessa primeira coluna gostaria de gerar a seguinte base 1ª COLUNA 2ª COLUNA ARG BRA ARG FRA BRA ARG BRA FRA FRA ARG FRA BRA Alguém saberia me dizer como conseguir esse resultado? Desde já muito obrigado! -- *Alexandre Rodrigues Loures* Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE e-Mail: alexandre.loures(a)ymail.com Site: www.ccsa.ufpb.br/ppge

2 2

Selecionar dados específicos de um data.frame
by Aline Lipsky 02 Mar '15

02 Mar '15

Bom dia, Eu tenho um data.frame chamado desembarques com as seguintes variáveis: id Embarcação Localidade Municipio Esse data.frame tem a dimensão (280,4). E tenho um vetor embarcação com 154 elementos. Quero selecionar os dados do data.frame desembarques de acordo com as embarcações existentes no vetor embarcação. Para isso tenho que comparar a variável Embarcação do data.frame com o vetor embarcação.Gostaria de sugestões de como fazer essa seleção. Grata pela atenção, Aline

3 4

[Dúvida] Submetendo código R em um cluster
by Pedro Rafael 28 Fev '15

28 Fev '15

Caros, eu tenho uma dúvida. Atualmente tenho um código paralelo que preciso submeter a um cluster. Estou utilizando o cluster Gauss do CESUP em que os detalhes dos seus hardwares podem serem encontrados no link: http://www.cesup.ufrgs.br/recursos-e-servicos/hardware-1/sgi-altix Basicamente em meu código eu faço uso de OpenMPI via biblioteca doSNOW que faz uso da SNOW e consequentemente da Rmpi. Uso também o pacote foreach onde uso a diretiva para paralelismo de loops %dopar%. O código simplesmente trata-se de simulações de Monte Carlo (MC) em que em cada loop de MC faço um bootstrap. Dentro de cada bootstrap faço uma otimização global em uma função objetivo via Particle Swarm Optimization que é bastante intensivo. Dessa forma, me pareceu razoável paralelizar o laço refere ao bootstrap. Matematicamente são independentes e um laço não depende de nenhum outro. O cluster em questão (cluster Gauss), o gerenciamento de filas é feito pelo PBS Pro 11. Estou um pouco de dificuldades de como construir o arquivo .sh para submissão do job. Não sei se devo usar o sistema FILL_UP ou ROUND_ROBIN como regras de alocação. Obrigado, Pedro Rafael.

1 0

Manipulação de tabela obtida com o getURL
by Michelle Bau Graczyk 27 Fev '15

27 Fev '15

Boa noite, Estou tentando coletar dados intraday do yahoo finance e manipular os dados obtidos excluindo todo cabeçalho, ou seja as 18 primeiras linhas, e tentando separar em colunas os dados restantes nas outras linhas. Já fiz várias tentaivas usando read.table, write,table e suas complementares e não obtive êxito. Alguém poderia por favor me ajudar? Abaixo segue o código e o que ele me fornece. Muito obrigada! Michelle library(RCurl) tmp <- getURL(' http://chartapi.finance.yahoo.com/instrument/1.0/AAPL/chartdata;type=quote;… ') tmp2<-paste(tmp, sep=",", dec='.') write.table(tmp2, "AAPL.txt") "x" "1" "uri:/instrument/1.0/AAPL/chartdata;type=quote;range=1d/csv ticker:aapl Company-Name:Apple Inc. Exchange-Name:NMS unit:MIN timezone:EST currency:USD gmtoffset:-18000 previous_close:133.0000 Timestamp:1424788200,1424811600 labels:1424790000,1424793600,1424797200,1424800800,1424804400,1424808000,1424811600 values:Timestamp,close,high,low,open,volume close:131.3100,133.5873 high:131.5000,133.6000 low:131.1700,133.5000 open:131.2700,133.5800 volume:0,2250400 1424788200,132.7800,132.9400,132.5100,132.9400,2077100 1424788260,132.6467,132.8000,132.6350,132.7300,503100 1424788379,132.5601,132.6500,132.3600,132.6500,625100 1424788380,132.0966,132.5700,132.0900,132.5600,553200 1424788440,132.2300,132.2600,132.0500,132.0900,579800 ….

4 5

MCA gráfico
by Edson Lira 26 Fev '15

26 Fev '15

Bom dia caros amigos. Novamente retomo um problema que ainda não consegui resolver. Download MCA COMUNICACAO.docx # nesse link tem o script o gráfico e as dúvidas | | | | | | | | | | | Download MCA COMUNICACAO.docxDownload MCA COMUNICACAO.docx free. File hosted free by DataFileHost. | | | | Visualizar em www.datafileho... | Visualizado por Yahoo | | | | | Download mj1.csv # os dados | | | | | | | | | | | Download mj1.csvDownload mj1.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafileho... | Visualizado por Yahoo | | | | | No geral, necessito melhorar a aparência e a escala do gráfico e ainda não consegui. [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 1

ddply - obter ganho de peso médio
by Fernando Souza 26 Fev '15

26 Fev '15

Caros amigos Estou tentando criar um função para calcular o ganho de peso médio. Eu quero calcular o ganho de peso por animal, sendo o ganho de peso a média dos ganhos obtidas a cada duas medidas sucessivas. Eu estou utilizando a função ddply() do pacote {plyr} onde divido o dataframe em grupos por ANIMAL e aplico a seguinte função a cada grupo. No entanto eu preciso que o ganho de peso obtido entre duas medidas consecutivas para um animal seja armazenado em um vetor, para em seguida fazer a média dos ganhos de pesos. No entanto não consigo criar esse vetor.Quando utilizo o comando abaixo, a função ddply retorna uma mensagem de erro indicando que não foi possível encontrar o objeto ganho. Como posso resolver esse problema? install.packages("plyr") library(plyr) A <- function(x){ for (i in 1:length(x$Pesos)-1){ difPeso <- x$Pesos[i+1] -x$Pesos[i] dias <-as.numeric( x$DATA[i+1] - x$DATA[i]) ganho [i]<-difPeso/dias # Gostaria de guardar este resultado em um vetor e retornar somente a média } return(mean(ganho)) } ddply(GPD,.(ANIMAL),A) GPD<- structure(list(DATA = structure(c(16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16454, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463, 16425, 16437, 16451, 16463), class = "Date"), ANIMAL = structure(c(1L, 1L, 1L, 1L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 10L, 10L, 10L, 10L, 15L, 15L, 15L, 15L, 19L, 19L, 19L, 19L, 2L, 2L, 2L, 2L, 7L, 7L, 7L, 7L, 13L, 13L, 13L, 13L, 16L, 16L, 16L, 16L, 17L, 17L, 17L, 17L, 24L, 24L, 24L, 24L, 3L, 3L, 3L, 3L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 12L, 12L, 12L, 12L, 20L, 20L, 20L, 20L, 22L, 22L, 22L, 22L, 6L, 6L, 6L, 6L, 11L, 11L, 11L, 11L, 14L, 14L, 14L, 14L, 18L, 18L, 18L, 18L, 21L, 21L, 21L, 21L, 23L, 23L, 23L, 23L), .Label = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "19", "20", "21", "22", "23", "24", "25" ), class = "factor"), TRAT = structure(c(4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 4L, 4L, 4L, 4L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 4L, 4L, 4L, 4L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 4L, 4L, 4L, 4L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("C", "T1", "T2", "T3"), class = "factor"), CONSUMO = c(0.6847039559, 0.659115454, 0.8601238438, 0.9662238017, 0.653160515, 0.76204388, 0.829147492, 1.0781859662, 0.5080284272, 0.7916373952, 0.7176244537, 1.2145475842, 0.5901562027, 0.5354316571, 0.9090091044, 0.8230424453, 0.800990936, 0.8258115771, 1.0350092142, 1.1778239333, 0.4218995736, 0.748182177, 0.9008559089, 0.7793604982, 0.8503107855, 0.8743924455, 1.0149210147, 1.2117686308, 1.0008721911, 0.9224992, 1.0628987486, 1.3235606054, 0.5847562394, 0.5971266995, 0.9015356047, 1.1422017785, 0.5506020875, 0.8996851452, 1.1010903098, 1.2777703226, 0.7757312968, 0.7327742822, 0.9957085833, 1.299725333, 0.6184034276, 0.7019017179, 0.8252818991, 0.7725886442, 0.57584, 0.8274993578, 0.9568347616, 1.0673499251, 0.7403614923, 0.7388590798, 0.7659735658, 1.0623794866, 0.6109657028, 0.683632284, 0.8188398169, 0.9643152772, 0.6384021189, 1.0495757074, 1.0262077037, 1.3312530279, 0.9202033913, 0.8721430158, 1.038471733, 1.1844675851, 0.8142008896, 1.0393409306, 1.2542648171, 1.3544221328, 0.5135583976, 0.5789023337, 0.7218802192, 1.0082128242, 0.7826673895, 0.9077870689, 0.9303001507, 1.4221483678, 0.8357668107, 0.8984098957, 1.0042573313, 1.5336138147, 0.407296, 0.6941073063, 0.8535336465, 0.9517448657, 0.5418470982, 0.5011396696, 0.6771379951, 0.6793475586, 0.7004441364, 0.9387915881, 1.0943207226, 1.2500940543 ), Pesos = c(14.75, 15.35, 16.8, 18.4, 17.6, 17.95, 19.2, 21, 19.4, 20.05, 21.4, 23.4, 14.55, 14.4, 15.1, 16.5, 19.1, 19.3, 20.4, 21.9, 16.3, 17.45, 18.9, 20.6, 22.25, 22.5, 23.5, 25.2, 22.8, 21.9, 24.2, 26.6, 18.45, 18.3, 19.5, 21.2, 21.8, 23.2, 22, 24.9, 18.6, 19.05, 20, 21, 17, 17.9, 18.5, 19.45, 20.6, 21.3, 22.6, 24.5, 19.1, 19.4, 20.1, 22.25, 16.65, 16.95, 17.8, 19.4, 22.8, 23.8, 24.4, 26.3, 19.3, 20.15, 21.9, 23.6, 22.9, 24.1, 25.6, 26.45, 16.65, 16.55, 17.3, 18.85, 22.4, 23.1, 25.4, 27.3, 22.9, 22.95, 25, 26.9, 16.7, 18.15, 19.7, 20.6, 14.05, 14.9, 15.6, 16.4, 18.9, 19.85, 21.8, 23.2)), .Names = c("DATA", "ANIMAL", "TRAT", "CONSUMO", "Pesos"), row.names = c(NA, -96L), class = "data.frame") -- Fernando Antônio de Souza Zootecnista,Dsc. Nutrição e alimentação animal (UFMG) Lattes: http://lattes.cnpq.br/6519538815038307 e-mail:nandodesouza@gmail.com (82)9618-6457 (tim) (82)8113-8781 (vivo) (82)3313-8781 (net) (31)9784-5273 (vivo-MG)

3 6

aumentar a velocidade
by Ze Henrique 26 Fev '15

26 Fev '15

Boa tarde pessoal, estou trabalhando com a função disponível no link < http://morotalab.org/Mrode2005/relmat/quass.txt> e gostaria de aumentar a velocidade de execução da funcão. Qualquer opinião é bem vinda. Segue abaixo um CMR: sire <- rep(1:50,each=50) dam <- rep(1:50,50) s <- c(rep(0,50),sire) d = c(rep(0,50),dam) system.time(quass(p,m)) Desde já obrigado -- *José Henrique Soler Guilhen* Mestrando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

1 0

Duvida em estatística espacial - Kernel em duas dimensões levando em conta prevalências.
by Augusto Ribas 26 Fev '15

26 Fev '15

Ola. Gostaria de saber se alguém poderia me sugerir algo para avaliar o risco, ou aonde uma doença ocorre espacialmente, mais comumente. Meus dados são constituídos de coordenadas geográficas de cidades, na forma de latitude e longitude. Dai para cada cidade eu tenho o número de casos de uma doença para um total de vaquinhas observadas. E calculo a prevalência como os casos pelo total de vaquinhas observadas. Então eu tenho esses dados, e gostaria de fazer um mapinha, mostrando espacialmente se existe uma área que tem maior ocorrência, maior chance de ocorrência. Apenas para explorar os dados a principio. Segue um pequeno exemplo de como estão os dados. ######################################################### library(MASS) set.seed(1) dados<-data.frame(lat=sample(1:10,20,replace=T), lon=sample(1:10,20,replace=T), casos=rpois(20,5), total=rpois(20,60) ) dados$prevalencia<-dados$casos/dados$total #Eu gostaria de fazer uma figura mais ou menos como essa, mas levando em consideração as prevalências. contorno<-kde2d(dados$lat,dados$lon) contour(contorno) ############################################################# Alguém tem alguma sugestão? Eu não sei nem o que começar a ler exatamente para Agradeço desde ja a atenção. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

3 3

Re: [R-br] GGPLOT
by Rafael Tieppo 23 Fev '15

23 Fev '15

DF <- data.frame(A = rnorm(40), B=rnorm(40)) #tamanho normal ggplot()+geom_point(data= DF, aes(x=A, y=B)) +ggtitle("Teste") #aumentar fonte ggplot()+geom_point(data= DF, aes(x=A, y=B))+theme (axis.title = element_text (size =25)) +ggtitle("Teste") +theme (title = element_text (size =25)) Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: http://www2.unemat.br/rafaeltieppo blog: https://sistemasagricolas.wordpress.com On Monday, February 23, 2015 4:48 PM, "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> wrote: Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: Vídeo Aula R (Carlos Andrade) 2. Re: Vídeo Aula R (Téo Calvo) 3. Re: Vídeo Aula R (Flavia Batista Gomes - Embrapa Amazonia Ocidental - CPAA) 4. GGPLOT (Edson Lira) 5. Re: GGPLOT (Fernando Gama) 6. Re: GGPLOT (Tropidurus Torquatus) 7. Re: GGPLOT (Edson Lira) ---------------------------------------------------------------------- Message: 1 Date: Mon, 23 Feb 2015 14:04:49 -0300 From: Carlos Andrade <prf.cantonio(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Vídeo Aula R Message-ID: <CA+_U+LzVNaaa7m5-f0W+P7q-v0Aw7AqnDH0Js9O1ou7ngHhwxg(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Pedro, Obrigado pelas dicas. Em 23 de fevereiro de 2015 10:56, Pedro Emmanuel Alvarenga Americano do Brasil <emmanuel.brasil(a)gmail.com> escreveu: > Carlos, > > Nunca consegui migrar totalmente do EpiData 3.1 para o EpiData versao XML. > No windows, ele fica travando em horas inapropriadas, não consigo fazer > algumas tarefas que julgo relevantes, e nem de perto tem os mesmos recursos > que o Epidata 3.1 tem quanto aos controles. A grande desvantagem do 3.1 é > que o acesso aos dados não são hierarquizados, ou seja, qualquer um com > acesso ao banco de dados pode modifica-lo. Ha cerca de dois anos comecei > uma migração para o RedCap (http://project-redcap.org/) que é o que eu > sugiro para todos atualmente que tenham um suporte de TI e um servidor que > possa instalar. Para os alunos de pos-graduação, que geralmente conduzem > sozinhos todos os passos do caminho dos dados, eu continuo recomendando o > EpiData 3.1. Salvo em algumas situações eu passo a recomendar o RedCap, > porque tenho acesso a esse programa na instituição que trabalho. > > Recentemente dei uma olhada no sítio do CSpro, e parece que la mesmo eles > já colocaram uns tutoriais em vídeos para a versão laptop e para a versão > android. > > http://www.census.gov/population/international/training/csprovideos.html > > Enquanto a gerar um formulario em papel que não seja tosco, o único que já > experimentei e ficou de fato muito bom, foi o Teleform. Mas o teleform > possui uma dificuldade de inserção dos dados a distância, é bem complicado > configurar pra funcionar em estudos multicentricos e é também mais > complicado fazer o data management porque o programa fica instalado num > computador local. Mas nada impede de fazer um tratamento "embeleze" nos > seus formularios com o editor de texto de sua preferência, inclusive > colocando o timbre a da sua instituição e coisas que usualmente são > relvantes para o controle no papel como "data da digitação" "nome do > profissional que preencheu", um cabeçalho com o identificador em todas as > folhas e coisas assim. > > Mais uma vez, pesando os pros e os contras, atualmente eu venho > recomendando o RedCap (ahhh e funciona bem com R gerando os relatórios do > andamento da pesquisa, ou seja, o R consegue conectar no RedCap e baixar > uma cópia dos dados, dos dicionários, da lista de usúarios, da estrutura > dos dados e etc). > > Espero que tenha ajudado. Abraço forte, > > > Pedro Brasil > > Em 23 de fevereiro de 2015 10:21, Carlos Andrade <prf.cantonio(a)gmail.com> > escreveu: > > Olá, Pedro. >> >> Já experimentei a integração EpiData com R. Na versão EpiData 3.1 >> consegui fazer um formulário para digitação de dados coletados em papel e >> dar entrada com os dados no R. Ficou muito tosco. Com a versão EpiData >> -manager/entry verifiquei que sofisticou mas não consegui montar um >> formulário (máscara de entrada). >> Você poderia indicar links de materiais mais "passo a passo" de criação >> de máscara de entrada de dados seja para EpiData ou CSPro, que também achei >> interessante, mas só tem para Windows. >> >> Agradeço antecipadamente. >> >> Em 23 de fevereiro de 2015 08:37, Pedro Emmanuel Alvarenga Americano do >> Brasil <emmanuel.brasil(a)gmail.com> escreveu: >> >> Ei teo, >>> >>> Eu leciono R ha dois anos para pos graduandos da area da saude. Por >>> acaso eu estava gravando video aulas de R com a intençao de mudar o formato >>> de oficinas presenciais para uma um curso semi-presencial quando vi o seu >>> email. Assisti somente a primeira aula sua e achei bem legal e bem parecido >>> com que estou fazendo. As video aulas ainda nao estao la, mas estarao em >>> breve, por outro lado o material escrito esta la desde 2013. >>> >>> http://blog.ipec.fiocruz.br/lapclin-chagas/ensino/computacao/ >>> >>> Parabens pela iniciativa. >>> >>> Abraço forte, >>> >>> Pedro Brasil via Android (:)= >>> Em 17/02/2015 18:42, "Téo Calvo" <teocalvo2(a)gmail.com> escreveu: >>> >>>> Olá caros, tudo bem? >>>> >>>> Dei início a um novo projeto, estou gravando video aulas de R para >>>> iniciantes. >>>> >>>> Foco principal é apresentar a linguagem a alunos de 1o e 2o Ano de >>>> Estatística, assim como curisosos e pesquisadores de outras áreas. >>>> >>>> >>>> https://www.youtube.com/watch?v=qtmeG_26_Pw&list=PLOUU7axoXx0r1n-ZSgE-iYNJW… >>>> >>>> Espero que colaborar com o desenvolvimento de futuros pesquisadores e >>>> também profissionais do mercado. >>>> >>>> Conto com sugestões e críticas. >>>> >>>> Att Téo Calvo. >>>> >>>> _______________________________________________ >>>> R-br mailing list >>>> R-br(a)listas.c3sl.ufpr.br >>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>>> código mínimo reproduzível. >>>> >>> >>> _______________________________________________ >>> R-br mailing list >>> R-br(a)listas.c3sl.ufpr.br >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>> código mínimo reproduzível. >>> >> >> >> >> -- >> Atenciosamente, >> >> Prof. Carlos A. S. de Andrade >> LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de >> Produção >> Universidade Federal de Campina Grande. >> Centro de Humanidades >> Unidade Acadêmica de Economia >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150223/c3246adc/atta…> ------------------------------ Message: 2 Date: Mon, 23 Feb 2015 16:29:52 -0300 From: Téo Calvo <teocalvo2(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Vídeo Aula R Message-ID: <CAGeP_1K5HfsNSVW2wRTLD1b0-D9zTPfuhqMm6SDgC=MT6aR8-w(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Olá a todos novamente. Fico contente com os comentários dos colegas. Todo este apoio tem sido importante para a motivação em continuar desenvolvendo material para as aulas e pensar em exemplos práticos para os alunos que estão começando a ter contato com R. Continuo a postar os vídeo sempre que posso. Peço que os colegas se inscrevam no canal para que recebam a notificação de novos vídeos, assim não preciso gerar Spawn na lista e em outras redes sociais. Muito bacana Pedro! Espero que minhas aulas possam ajuda-lo de alguma forma. Muito obrigado novamente. Att Téo Calvo. Em 23 de fevereiro de 2015 08:37, Pedro Emmanuel Alvarenga Americano do Brasil <emmanuel.brasil(a)gmail.com> escreveu: > Ei teo, > > Eu leciono R ha dois anos para pos graduandos da area da saude. Por acaso > eu estava gravando video aulas de R com a intençao de mudar o formato de > oficinas presenciais para uma um curso semi-presencial quando vi o seu > email. Assisti somente a primeira aula sua e achei bem legal e bem parecido > com que estou fazendo. As video aulas ainda nao estao la, mas estarao em > breve, por outro lado o material escrito esta la desde 2013. > > http://blog.ipec.fiocruz.br/lapclin-chagas/ensino/computacao/ > > Parabens pela iniciativa. > > Abraço forte, > > Pedro Brasil via Android (:)= > Em 17/02/2015 18:42, "Téo Calvo" <teocalvo2(a)gmail.com> escreveu: > >> Olá caros, tudo bem? >> >> Dei início a um novo projeto, estou gravando video aulas de R para >> iniciantes. >> >> Foco principal é apresentar a linguagem a alunos de 1o e 2o Ano de >> Estatística, assim como curisosos e pesquisadores de outras áreas. >> >> >> https://www.youtube.com/watch?v=qtmeG_26_Pw&list=PLOUU7axoXx0r1n-ZSgE-iYNJW… >> >> Espero que colaborar com o desenvolvimento de futuros pesquisadores e >> também profissionais do mercado. >> >> Conto com sugestões e críticas. >> >> Att Téo Calvo. >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150223/e9692112/atta…> ------------------------------ Message: 3 Date: Mon, 23 Feb 2015 17:04:23 -0300 (BRT) From: Flavia Batista Gomes - Embrapa Amazonia Ocidental - CPAA <flavia.b.gomes(a)embrapa.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Vídeo Aula R Message-ID: <1621180510.1586048.1424721863692.JavaMail.root(a)embrapa.br> Content-Type: text/plain; charset="iso-8859-1" Oi Téo, Ainda não vi seus vídeos mas acessarei com certeza, já sugiro vc gravar um sobre análise de Probit e cálculo de CL50. Será muito útil. Flávia Batista ----- Mensagem original ----- De: "Téo Calvo" <teocalvo2(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Segunda-feira, 23 de fevereiro de 2015 15:29:52 Assunto: Re: [R-br] Vídeo Aula R Olá a todos novamente. Fico contente com os comentários dos colegas. Todo este apoio tem sido importante para a motivação em continuar desenvolvendo material para as aulas e pensar em exemplos práticos para os alunos que estão começando a ter contato com R. Continuo a postar os vídeo sempre que posso. Peço que os colegas se inscrevam no canal para que recebam a notificação de novos vídeos, assim não preciso gerar Spawn na lista e em outras redes sociais. Muito bacana Pedro! Espero que minhas aulas possam ajuda-lo de alguma forma. Muito obrigado novamente. Att Téo Calvo. Em 23 de fevereiro de 2015 08:37, Pedro Emmanuel Alvarenga Americano do Brasil < emmanuel.brasil(a)gmail.com > escreveu: Ei teo, Eu leciono R ha dois anos para pos graduandos da area da saude. Por acaso eu estava gravando video aulas de R com a intençao de mudar o formato de oficinas presenciais para uma um curso semi-presencial quando vi o seu email. Assisti somente a primeira aula sua e achei bem legal e bem parecido com que estou fazendo. As video aulas ainda nao estao la, mas estarao em breve, por outro lado o material escrito esta la desde 2013. http://blog.ipec.fiocruz.br/lapclin-chagas/ensino/computacao/ Parabens pela iniciativa. Abraço forte, Pedro Brasil via Android (:)= Em 17/02/2015 18:42, "Téo Calvo" < teocalvo2(a)gmail.com > escreveu: <blockquote> Olá caros, tudo bem? Dei início a um novo projeto, estou gravando video aulas de R para iniciantes. Foco principal é apresentar a linguagem a alunos de 1o e 2o Ano de Estatística, assim como curisosos e pesquisadores de outras áreas. https://www.youtube.com/watch?v=qtmeG_26_Pw&list=PLOUU7axoXx0r1n-ZSgE-iYNJW… Espero que colaborar com o desenvolvimento de futuros pesquisadores e também profissionais do mercado. Conto com sugestões e críticas. Att Téo Calvo. _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem ( http://www.leg.ufpr.br/r-br-guia ) e forneça código mínimo reproduzível. _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem ( http://www.leg.ufpr.br/r-br-guia ) e forneça código mínimo reproduzível. </blockquote> _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150223/c99fb618/atta…> ------------------------------ Message: 4 Date: Mon, 23 Feb 2015 20:48:03 +0000 (UTC) From: Edson Lira <edinhoestat(a)yahoo.com.br> To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] GGPLOT Message-ID: <1173242903.5170236.1424724483802.JavaMail.yahoo(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Meus caros amigos, estou tentando usar um exemplo de aplicação de MCA, script abaixo: #Usando MCA# load packagesrequire(FactoMineR)require(ggplot2) # load data teadata(tea)teanames(tea)# select these columnsnewtea = tea[, c("Tea", "How", "how", "sugar", "where", "always")] # take a peekhead(newtea)names(newtea)str(newtea)cats = apply(newtea, 2, function(x) nlevels(as.factor(x)))cats # apply MCAmca1 = MCA(newtea, graph = FALSE) # list of resultsmca1 # table of eigenvaluesmca1$eig #We can use the package "ggplot2()" to get a nice plot:# data frame with variable coordinatesmca1_vars_df = data.frame(mca1$var$coord, Variable = rep(names(cats), cats)) # data frame with observation coordinatesmca1_obs_df = data.frame(mca1$ind$coord) par(cex.axis=0.5)# plot of variable categoriesggplot(data=mca1_vars_df,cex.axis=0.5, cex.lab=0.5, aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df))) + geom_hline(yintercept = 0, colour = "gray70") + geom_vline(xintercept = 0, colour = "gray70") + geom_text(aes(colour=Variable)) + ggtitle("MCA plot of variables using R package FactoMineR") # MCA plot of observations and categoriesggplot(data = mca1_obs_df, aes(x = Dim.1, y = Dim.2)) + geom_hline(yintercept = 0, colour = "gray70") + geom_vline(xintercept = 0, colour = "gray70") + geom_point(colour = "gray50", alpha = 0.7) + geom_density2d(colour = "gray80") + geom_text(data = mca1_vars_df, aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df), colour = Variable)) + ggtitle("MCA plot of variables using R package FactoMineR") + scale_colour_discrete(name = "Variable") Como faço para mudar o tamanho da fonte do título, dos eixos, legenda, etc.? Tentei os comandos cex, cex.axis, cex.lab, mas não consegui. Alguém tem alguma sugestão? [ ]'s.Edson Lira Estatístico Manaus-Amazonas -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150223/d95b2c4f/atta…> ------------------------------ Message: 5 Date: Mon, 23 Feb 2015 17:55:49 -0300 From: Fernando Gama <f.fabiogama88(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Edson Lira <edinhoestat(a)yahoo.com.br> Subject: Re: [R-br] GGPLOT Message-ID: <CAOwqRysTxYL0da4xHbOo3rUU68dJ8z3VTFP9Db=LuhQeH2Ahog(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Olá Edson boa tarde, meu caro, não sei se esssa é a melhor solução mas você poderia utilizar text(). Utilizei nos meus scripts e atendeu minhas expectativas: a library da função pode ser encontrada neste link: http://stat.ethz.ch/R-manual/R-devel/library/graphics/html/text.html abraços 2015-02-23 17:48 GMT-03:00 Edson Lira <edinhoestat(a)yahoo.com.br>: > Meus caros amigos, estou tentando usar um exemplo de aplicação de MCA, > script abaixo: > > #Usando MCA > # load packages > require(FactoMineR) > require(ggplot2) > > # load data tea > data(tea) > tea > names(tea) > # select these columns > newtea = tea[, c("Tea", "How", "how", "sugar", "where", "always")] > > # take a peek > head(newtea) > names(newtea) > str(newtea) > cats = apply(newtea, 2, function(x) nlevels(as.factor(x))) > cats > > # apply MCA > mca1 = MCA(newtea, graph = FALSE) > > # list of results > mca1 > > # table of eigenvalues > mca1$eig > > #We can use the package "ggplot2()" to get a nice plot: > # data frame with variable coordinates > mca1_vars_df = data.frame(mca1$var$coord, Variable = rep(names(cats), cats > )) > > # data frame with observation coordinates > mca1_obs_df = data.frame(mca1$ind$coord) > > > par(cex.axis=0.5) > # plot of variable categories > ggplot(data=mca1_vars_df,cex.axis=0.5, cex.lab=0.5, > aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df))) + > geom_hline(yintercept = 0, colour = "gray70") + > geom_vline(xintercept = 0, colour = "gray70") + > geom_text(aes(colour=Variable)) + > ggtitle("MCA plot of variables using R package FactoMineR") > > > # MCA plot of observations and categories > ggplot(data = mca1_obs_df, aes(x = Dim.1, y = Dim.2)) + > geom_hline(yintercept = 0, colour = "gray70") + > geom_vline(xintercept = 0, colour = "gray70") + > geom_point(colour = "gray50", alpha = 0.7) + > geom_density2d(colour = "gray80") + > geom_text(data = mca1_vars_df, > aes(x = Dim.1, y = Dim.2, > label = rownames(mca1_vars_df), colour = Variable)) + > ggtitle("MCA plot of variables using R package FactoMineR") + > scale_colour_discrete(name = "Variable") > > > Como faço para mudar o tamanho da fonte do título, dos eixos, legenda, > etc.? Tentei os comandos cex, cex.axis, cex.lab, mas não consegui. > > > Alguém tem alguma sugestão? > > > > [ ]'s. > Edson Lira > Estatístico > Manaus-Amazonas > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -- Att, Fernando Gama da Mata Vale of Institute Technology -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150223/6145791f/atta…> ------------------------------ Message: 6 Date: Mon, 23 Feb 2015 18:06:54 -0300 From: Tropidurus Torquatus <t.torquatus(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Edson Lira <edinhoestat(a)yahoo.com.br> Subject: Re: [R-br] GGPLOT Message-ID: <32AED325-B4B3-4AF9-84E1-CDB3560C9F92(a)gmail.com> Content-Type: text/plain; charset="windows-1252" Edson, Para construir meus gráficos com ggplot2 normalmente eu uso e abuso desses 2 sites: Acho que esse vai atender a suas dúvidas: http://www.cookbook-r.com/Graphs/ <http://www.cookbook-r.com/Graphs/> Esse abaixo é completo? mas pode gerar mais dúvidas: http://ggplot2.org <http://ggplot2.org/> Se não me falha a memória, podem ser usados duas formas? alteração direta na inserção dos dados mas aparente não seria a forma que você está utilizando (xlab/ylab), outra forma alterando o ?theme?. > Em 23/02/2015, à(s) 17:48, Edson Lira <edinhoestat(a)yahoo.com.br> escreveu: > > Meus caros amigos, estou tentando usar um exemplo de aplicação de MCA, script abaixo: > > #Usando MCA > # load packages > require(FactoMineR) > require(ggplot2) > > # load data tea > data(tea) > tea > names(tea) > # select these columns > newtea = tea[, c("Tea", "How", "how", "sugar", "where", "always")] > > # take a peek > head(newtea) > names(newtea) > str(newtea) > cats = apply(newtea, 2, function(x) nlevels(as.factor(x))) > cats > > # apply MCA > mca1 = MCA(newtea, graph = FALSE) > > # list of results > mca1 > > # table of eigenvalues > mca1$eig > > #We can use the package "ggplot2()" to get a nice plot: > # data frame with variable coordinates > mca1_vars_df = data.frame(mca1$var$coord, Variable = rep(names(cats), cats)) > > # data frame with observation coordinates > mca1_obs_df = data.frame(mca1$ind$coord) > > > par(cex.axis=0.5) > # plot of variable categories > ggplot(data=mca1_vars_df,cex.axis=0.5, cex.lab=0.5, > aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df))) + > geom_hline(yintercept = 0, colour = "gray70") + > geom_vline(xintercept = 0, colour = "gray70") + > geom_text(aes(colour=Variable)) + > ggtitle("MCA plot of variables using R package FactoMineR") > > > # MCA plot of observations and categories > ggplot(data = mca1_obs_df, aes(x = Dim.1, y = Dim.2)) + > geom_hline(yintercept = 0, colour = "gray70") + > geom_vline(xintercept = 0, colour = "gray70") + > geom_point(colour = "gray50", alpha = 0.7) + > geom_density2d(colour = "gray80") + > geom_text(data = mca1_vars_df, > aes(x = Dim.1, y = Dim.2, > label = rownames(mca1_vars_df), colour = Variable)) + > ggtitle("MCA plot of variables using R package FactoMineR") + > scale_colour_discrete(name = "Variable") > > > Como faço para mudar o tamanho da fonte do título, dos eixos, legenda, etc.? Tentei os comandos cex, cex.axis, cex.lab, mas não consegui. > > > Alguém tem alguma sugestão? > > > > [ ]'s. > Edson Lira > Estatístico > Manaus-Amazonas > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150223/9b691c0e/atta…> ------------------------------ Message: 7 Date: Mon, 23 Feb 2015 21:09:55 +0000 (UTC) From: Edson Lira <edinhoestat(a)yahoo.com.br> To: Tropidurus Torquatus <t.torquatus(a)gmail.com>, "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] GGPLOT Message-ID: <390107215.4582006.1424725795582.JavaMail.yahoo(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Ok vou tentar. [ ]'s.Edson Lira Estatístico Manaus-Amazonas Em Segunda-feira, 23 de Fevereiro de 2015 17:07, Tropidurus Torquatus <t.torquatus(a)gmail.com> escreveu: Edson, Para construir meus gráficos com ggplot2 normalmente eu uso e abuso desses 2 sites: Acho que esse vai atender a suas dúvidas: http://www.cookbook-r.com/Graphs/ Esse abaixo é completo? mas pode gerar mais dúvidas: http://ggplot2.org Se não me falha a memória, podem ser usados duas formas? alteração direta na inserção dos dados mas aparente não seria a forma que você está utilizando (xlab/ylab), outra forma alterando o ?theme?. Em 23/02/2015, à(s) 17:48, Edson Lira <edinhoestat(a)yahoo.com.br> escreveu: Meus caros amigos, estou tentando usar um exemplo de aplicação de MCA, script abaixo: #Usando MCA# load packagesrequire(FactoMineR)require(ggplot2) # load data teadata(tea)teanames(tea)# select these columnsnewtea = tea[, c("Tea", "How", "how", "sugar", "where", "always")] # take a peekhead(newtea)names(newtea)str(newtea)cats = apply(newtea, 2, function(x) nlevels(as.factor(x)))cats # apply MCAmca1 = MCA(newtea, graph = FALSE) # list of resultsmca1 # table of eigenvaluesmca1$eig #We can use the package "ggplot2()" to get a nice plot:# data frame with variable coordinatesmca1_vars_df = data.frame(mca1$var$coord, Variable = rep(names(cats), cats)) # data frame with observation coordinatesmca1_obs_df = data.frame(mca1$ind$coord) par(cex.axis=0.5)# plot of variable categoriesggplot(data=mca1_vars_df,cex.axis=0.5, cex.lab=0.5, aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df))) + geom_hline(yintercept = 0, colour = "gray70") + geom_vline(xintercept = 0, colour = "gray70") + geom_text(aes(colour=Variable)) + ggtitle("MCA plot of variables using R package FactoMineR") # MCA plot of observations and categoriesggplot(data = mca1_obs_df, aes(x = Dim.1, y = Dim.2)) + geom_hline(yintercept = 0, colour = "gray70") + geom_vline(xintercept = 0, colour = "gray70") + geom_point(colour = "gray50", alpha = 0.7) + geom_density2d(colour = "gray80") + geom_text(data = mca1_vars_df, aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df), colour = Variable)) + ggtitle("MCA plot of variables using R package FactoMineR") + scale_colour_discrete(name = "Variable") Como faço para mudar o tamanho da fonte do título, dos eixos, legenda, etc.? Tentei os comandos cex, cex.axis, cex.lab, mas não consegui. Alguém tem alguma sugestão? [ ]'s.Edson Lira Estatístico Manaus-Amazonas_______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150223/8d9d81b3/atta…> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------ Fim da Digest R-br, volume 50, assunto 15 *****************************************

1 0

Vídeo Aula R
by Téo Calvo 23 Fev '15

23 Fev '15

Olá caros, tudo bem? Dei início a um novo projeto, estou gravando video aulas de R para iniciantes. Foco principal é apresentar a linguagem a alunos de 1o e 2o Ano de Estatística, assim como curisosos e pesquisadores de outras áreas. https://www.youtube.com/watch?v=qtmeG_26_Pw&list=PLOUU7axoXx0r1n-ZSgE-iYNJW… Espero que colaborar com o desenvolvimento de futuros pesquisadores e também profissionais do mercado. Conto com sugestões e críticas. Att Téo Calvo.

8 12