Maio 2012 - R-br - DINF UFPR

Bayesiana [OFF-TOPIC]
by ASANTOS 23 Mai '12

23 Mai '12

Bom dia Pessoal, Gostaria de saber se alguém quem trabalha com Bayesiana poderia recomendar algumas boas literaruras para estudo para que esta començando a usar esta abordagem, preferencialmente gostaria de material aplicado ao R e não ao winBugs. Obrigado, *Alexandre dos Santos *Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG

4 3

duvidas sobre teste de causalidade de granger
by Marcelo Zeuli 23 Mai '12

23 Mai '12

Bom dia, qual é o pacote mais adequado para teste de causalidade de granger? o FIAR? ********************************************* Type Package Title Functional Integration Analysis in R Version 0.4 Date 2011-04-01 Author Bjorn Roelstraete Maintainer Bjorn Roelstraete <Bjorn.Roelstraete(a)ugent.be> Description Contains DCM, (autoregressive) SEM, and multivariate Granger causality tests for analyzing fMRI connectivity data. Depends R (>= 2.11.1),lavaan(>= 0.4-7),Matrix,boot,accuracy,np ********************************************************** Grato, Marcelo Zeuli

2 2

Re: [R-br] Alteração em função de pacote
by Junior Beleti 22 Mai '12

22 Mai '12

Olá Benilton. No manual encontrei instruções de como se criar uma função em um novo pacote. No meu caso, onde apenas modifiquei a função eigen.f para eigem.f, gostaria de saber como fazer para carregá-la no momento da execução. Quando chamo a função eigen: eigen(cbind(c(1,-1),c(-1,1))) $values [1] 2 0 $vectors [,1] [,2] [1,] -0.7071068 -0.7071068 [2,] 0.7071068 -0.7071068 Já na chamada a eigem: eigem(cbind(c(1,-1),c(-1,1))) Erro: não foi possível encontrar a função "eigem" Att, Carlos.

3 2

Visualização interativa
by Eder David Borges da Silva 21 Mai '12

21 Mai '12

Pessoal gostaria de trocar uma ideia sobre visualização de dados. Ultimamente venho dando uma olhada do pacote googleVis, em especial na função gvisMotionChart, que possibilita criar gráficos onde podemos interagir com eles selecionando time, dando destaque a um nivel etc. O que estou procurando e algo como o gvisMotionChart, porem com uma especie de filtro a mais, por exemplo imaginem que tenho um conjunto de dados, emails da R-br, nos três últimos anos, onde uma coluna tenho o ano, uma segunda o participante, e uma terceira o numero de emails. Eu gostaria de fazer um gráfico selecionando por exemplo, no ano 0, apenas o user1 versus user2, onde essa seleção seria feia em uma caixa de seleção, selecionando ano 0, depois em um segunda caixa selecionando user1 e user2, e automaticamente o gráfico de barras por exemplo iria ser atualizado. para quem já ouviu falar seria um qlikview http://www.inteligencia-negocios.com/produtos/qlikdemos porem feito via R com a saida final de interface html. Oque quero compartilhar com vocês é o seguinte: alguém já viu algo assim feito em R?

2 2

Plots com dimensões diferentes
by Marcelo Claro de Souza 21 Mai '12

21 Mai '12

Preciso plotar 4 figuras no mesmo frame (2 x2) mas preciso que as figuras tenham dimensões diferentes. No caso preciso que as 2 figuras superiores ocupem apenas 25% do frame e as figuras inferiores 75%. Como posso fazer isso? Abraço. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Re: [R-br] modelo linear com efeito aleatório
by Ivan Bezerra Allaman 21 Mai '12

21 Mai '12

library(lme4) mod <- lmer(tempo ~ desfecho + (1|pac),data=data) summary(mod) ?lmer (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 2

modelo linear com efeito aleatório
by Pedro Emmanuel Alvarenga Americano do Brasil 21 Mai '12

21 Mai '12

Amigos de R, Pacientes foram expostos a experimentos para verificar o tempo de resposta a determinadas exposições. Para cada exposição (que são 3, mas prefiro aborda-las individualmente) verificou-se se o sujeito acertou a resposta desejada e quanto tempo levou para responder. A hipótese do estudo é que quando o sujeito acerta a resposta o tempo de resposta é menor do que quando erra. O problema é que cada paciente é exposto diversas vezes. Então eu tenho um banco com 10 pacientes, mas com mais de 300 linhas, semelhante a simulação abaixo. A minha questão é que eu imagino que haja dependencia dos pacientes, ou seja, que as respostas sejam mais parecidas dentro de cada sujeito do que entre cada sujeito. Provavelmente a variação teria dois componentes, o de dentro do sujeito e entre os sujeitos. POr isso, eu imaginei que utilizar um modelo aleatório com o paciente como o elemento aleatório fosse mais interpretável que um modelo linear simples. Mas infelizmente nao estou conseguindo fazer a função funcionar. Olhei os exemplos , li algumas vezes a documentção... mas ainda não foi. Qualquer sugestão de como fazer funcionar é bem vinda! z <- rnorm(200, mean = 120 , sd = 10) x <- rnorm(100, mean = 90, sd = 20) data <- as.data.frame(rbind(cbind(x,1,1:10),cbind(z,0,1:10))) names(data) <- c('tempo','desfecho','pac') summary(lm(tempo ~ desfecho + pac,data)) library('nlm') > lme(tempo ~ desfecho + pac,data) Erro em getGroups.data.frame(dataMix, groups) : Invalid formula for groups > lme(tempo ~ desfecho,data,random= pac) Erro em reStruct(random, REML = REML, data = NULL) : objeto 'pac' não encontrado > lme(tempo ~ desfecho,data,random= ~ pac) Erro em getGroups.data.frame(dataMix, groups) : Invalid formula for groups > lme(tempo ~ desfecho,data,random= ~desfecho|pac) Erro em lme.formula(tempo ~ desfecho, data, random = ~desfecho | pac) : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence (10) Abraço forte, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

1 0

séries em escalas diferentes com lattice
by Lucas Barbosa 20 Mai '12

20 Mai '12

Boa tarde pessoal, Como posso fazer um gráfico de duas séries em diferentes escalas (com dois eixos verticais) usando o lattice? Sei fazer com o plot normal (exemplo abaixo), mas seria muito vantajoso para mim criar uma função panel -- quero fazer vários gráficos, um para cada valor de uma outra variável ("condicionando" com aquela barra vertical). x=1:100 y1 = rpois(100, 20) y2 = rnorm(100, 5) par(mar=c(5,5,4,5)) plot(x, y1, type="l", axes=FALSE, frame=TRUE, ann=FALSE, col=4) axis(4) mtext(text="y1", 4, line=3) par(new=T) plot(x, y2, type='l', lwd=1.5)

2 2

Dúvida tapply
by Leandro Marino 19 Mai '12

19 Mai '12

Caros, estou executando o seguinte comando utilizando os dados da ProvaBrasil 2009 no RJ. > unlist(tapply(alunorj$NU_THETAT_L,list1,length,simplify=F)) [1] 94410 72130 17151 93198 55995 6091 2863 > unlist(tapply(alunorj$NU_THETAT_L,list1,length,simplify=T)) 4ª série / 5º ano EF 8ª série / 9º ano EF Matutino 94410 93198 Vespertino 72130 55995 Noturno NA 6091 Intermediário 17151 2863 Note que o NA que tem no noturno 4ª serie nao aparece quando o simplify=F, existe alguma forma de forçá-lo? Preciso que tenha uma saída em vetor destes dados. Atenciosamente, Leandro Marino Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

2 1

aplicando cálculo sobre linhas de um dataframe
by Vitor Aguiar 18 Mai '12

18 Mai '12

Olá. Postei essa pergunta ontem mas não obtive resposta. Desculpe por estar postando novamente. Tenho um dataframe com milhares de linhas com valores NAs e não-NAs. Suponhamos a seguinte linha: 1. NA NA NA NA NA NA 6 6 NA NA 8 11 NA NA 17 18 NA NA 10 10 11 13 8 9 11 11 11 11 Vamos considerar os dados em pares. Assim, meu primeiro par é NA NA, o quarto par é 6 6. Há alguma forma de eu contar quantos pares iguais eu tenho em cada linha do dataframe? Assim, nessa linha por exemplo eu teria o resultado 4, pois tenho 4 pares de valores iguais (6 6, 10 10, 11 11, 11 11). Então preciso aplicar o cálculo para todas as linhas. O caminho seria um loop onde eu contaria dados[i, 1] == dados[i, 2]? Obrigado, Vitor

3 3