Junho 2012 - R-br - DINF UFPR

Contar atraso de um número numa matrix...
by andrebvs 18 Jun '12

18 Jun '12

3 5

Buscar por palavras no site da ABL
by Lucas Barbosa 18 Jun '12

18 Jun '12

Boa tarde, pessoal! No site da ABL ( http://www.academia.org.br/abl/cgi/cgilua.exe/sys/start.htm?sid=23) é possível buscar palavras no Vocabulário Ortográfico da Língua Portuguesa. Mas a busca não pode ser definida a partir de um endereço html, é necessário preencher o campo e clicar em "Pesquisar". Como faço para automatizar esse tipo de busca usando R? É possível? Encontrei, na stackoverflow <http://stackoverflow.com/>, problemas semelhantes resolvidos com as funções postForm ou getForm (pacote RCurl), mas não estou entendendo muito bem como elas funcionam. O html da caixa de texto e do botão "Pesquisar": <input id="palavra" size="30" name="palavra" maxlenght="30"> <input onclick="ajaxport(palavra.value,'')" value="Pesquisar" name="pesquisar" type="button"> Obrigado! Lucas

2 3

Divulgação curso: "Business Intelligence - A partir de uma perspectiva Bayesiana"
by nani.design＠bol.com.br 18 Jun '12

18 Jun '12

1 0

regressão linear
by Bernardo Rodrigues 18 Jun '12

18 Jun '12

Alguem me pode ajudar como é que se determina distribuições posteriori da regressão linear a partir de distribuição a priori nao informativa?

4 3

Atualização do R no Linux Ubuntu 10.04
by Andre Oliveira 18 Jun '12

18 Jun '12

Boa noite, estou com a dúvida uso ubuntu 10.04 depois de alguns anos considero que em minha máquina está tudo configurado como preciso criei ate uma versão persolizada em um DVD via garfio. Em abril de 2013 o 10.04 deixará de ter suporte. Minha dúvida é a seguinte eu receberei atualização do R caso eu decida ficar com o 10.04 após abril de 2013? obrigado André Oliveira Souza

3 2

Diferença entre duas abordagens de comparação usando GLM
by ASANTOS 18 Jun '12

18 Jun '12

Boa noite pessoal, Gostaria de saber se do ponto de vista da interpretação estatística do método empregado, existe diferença entre comparar 3 tratamento hipotéticos usando: 1) require(multcomp) y<- c(mapply(rpois, lambda=c(5,20,5), MoreArgs=list(n=100)))##Criação da variável resposta trat <- gl(3,100) ##Criação dos tratamentos m.1 <- glm(y~trat,family="poisson") anova.m.1 <- anova(m.1, test="Chi") anova.m.1 summary(glht(m.1, linfct=mcp(trat="Tukey"))) #Os tratamentos 1 e 3 são semelhantes e diferentes do tratamentos 2 ou 2) m.1 <-glm(y~trat,family="poisson") m.0 <- glm(y~1,family="poisson") anova(m.1,m.0,test="Chi")####Existe diferença sort(tapply(y,trat,mean))## Ordenando as médias levels(trat) ##Modelo T1 e T3 t1t3<-trat levels(t1t3) levels(t1t3)[1]<-"1+3" levels(t1t3)[3]<-"1+3" levels(t1t3) m.2<- glm(y~t1t3,family="poisson") anova(m.1,m.2,test="Chi")####Semelhantes ## Mas se faço qualquer uma das outras combinações "1+2" ou "2+3" os valores do teste Chi são significativos, podendo-se concluir que os tratamentos 1 e 3 são semelhantes e diferentes do tratamentos 2. 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 Fone: +55 (35) 9223-0304

2 2

Trabalhando com equações
by Diego Leal 18 Jun '12

18 Jun '12

Oi pessoal, gostaria de saber como faço pra trabalhar com equações no R project (como represento uma equação, maximar função, derivar e integrar) Existe algum tutorial, ou coisa do tipo onde eu possa me informar melhor sobre o assunto, não sei nem por onde começar?

4 7

[OFF-TOPIC]Usar "for" via terminal linux?
by Ivan Bezerra Allaman 17 Jun '12

17 Jun '12

Boa noite senhores! Alguém sabe se é possível usar algum comando semelhante ao "for" via terminal no linux? Eu tenho uma sequência de imagens para recortar e estou utilizando o camando: convert crn_01.png -crop 970x474+25+60 +repage crn_01.png Algo como: for(i in 1:200){ convert crn_i.png -crop 970x474+25+60 +repage crn_i.png } já ajudaria! Alguém sabe se é possível fazer algo semelhante no terminal? Desde já muito obrigado! (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}

4 5

como juntar dois bancos em função de um intervalo
by Paula Mendes Luz 17 Jun '12

17 Jun '12

oi pessoal eu normalmente utilizo o merge() para juntar bancos em função de 1 ou mais variáveis que são iguais nos dois bancos. também normalmente tenho bancos com medidas repetidas para cada indivíduo de modo que novas linhas são criadas no ato do merge() em função das informações dos bancos originais. agora preciso juntar dois bancos em função de um intervalo de números (na verdade são datas mas isso não importa): um exemplo do que preciso fazer p o banco original com > 3000 ids seria banco1 tem o id, a entrada, a saída, e o grupo id <- c(1, 1, 1, 2, 2) ent <- c(1, 10, 20, 1, 30) sai <- c(9, 19, 40, 29, 50) grupo <- c(1, 2, 3, 1, 2) banco1 <- cbind(id, ent, sai, grupo) banco2 tem o id, o dia, tipo do teste id <- c(1, 1, 1, 1, 1, 2, 2, 2, 2) dia <- c(5, 8, 12, 18, 25, 22, 25, 29, 41) teste <- c("A", "A", "B", "B", "A", "B", "B", "C", "D") banco2 <- cbind(id, dia, teste) e o banco 3 seria o merge() dos dois em função da avaliação se o DIA do banco2 está contido no intervalo (entre ENT e SAI) do banco1 o resultado seria o banco 3 com 5 linhas para o id1 com as 5 vars da seguinte forma id ent sai grupo dia teste 1 1 9 1 5 A 1 1 9 1 8 A 1 10 19 2 12 B 1 10 19 2 18 B 1 20 40 3 25 A e 4 linhas para o id 2 alguém pode ajudar???? eu li a respeito da função findInterval() mas não consegui fazer funcionar.... agradeço desde já! abraços, paula.

2 4

simulação da normal bivariada
by Natalia Martins 17 Jun '12

17 Jun '12

Boa noite pessoal, estou tentando simular dados diários de temperatura e radiação solar usando uma normal bivariada e após a simulação eu tentei fazer uma reta de regressão para verificar se minha simulação esta de acordo, ou seja, eu esperaria uma reta pois os valores deveriam ser iguais. No entanto quando a faço o meu R^2 ajustado é ruim (sei que tambem devo fazer testes e não apenas olhar o R), alguem pode me ajudar, pois acredito que estou fazendo algo errado? a programação que estou utilizando segue a baixo: # gerando os dados normal bivariada set.seed(1234) #semente sigma <- matrix(c(19260.71398,28.92129,28.92129,2.173789), ncol=2) sim_bnorm <- rmvnorm(n=1000, mean=c(mean(RST$RS),mean(RST$Temp)), sigma=sigma) bnorm <- data.frame(sim_bnorm) simulado <- sample(bnorm$X2,31) modelo <- lm(simulado~RST$Temp) summary(modelo) plot( RST$Temp , simulado,xlab="Valores reais",ylab="Valores simulados") abline(modelo,col="red") Muito obrigada --

4 8