- R-br - DINF UFPR

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