Abril 2012 - R-br - DINF UFPR

Ajuda para evitar um loop com elegância (agora correto)
by Jose Claudio Faria 14 Abr '12

14 Abr '12

Pessoal, Gostaria de evitar um loop no procedimento abaixo da forma mais elegante possível: # R help m <- matrix(runif(24), ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(m) # média das colunas mr <- rowMeans (m) # médias das variáveis for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res Toda ajuda é bem vinda... -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 7

SciELO Books
by Benilton Carvalho 14 Abr '12

14 Abr '12

Uma sugestao de um modo legal de distribuicao de livros (agradecam tbm `a licenca Creative Commons): http://books.scielo.org/ e algo um pouco mais "especifico": http://books.scielo.org/id/p2qh6 benilton

4 3

Ajuda para evitar loop com elegância
by Jose Claudio Faria 14 Abr '12

14 Abr '12

Pessoal, Gostaria de evitar um loop no procedimento abaixo da forma mais elegante possível: m <- matrix(1:24, ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(matrixdata) # média das colunas mr <- rowMeans (matrixdata) # médias das variáveis for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res Toda ajuda é bem vinda... -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 1

Juntando os dados de dois objetos um abaixo do outro
by ASANTOS 14 Abr '12

14 Abr '12

Bom dia Pessoal, Criei um objeto com duas variáveis: v1<-1:11 v2<-rnorm(11) tab1<-cbind(v1,v2) e outro: v1<-1:11 v2<-rnorm(11) tab2<-cbind(v1,v2) e gostaria de juntar o tab 1 e tab 2 um abaixo do outro em um novo objeto. Fiz algumas pesquisas e encontrei o pacote plyr que permite especificar o local onde se deseja a fusão, mas não funcionou não sendo: require(plyr) tab12<- join(tab1,tab2,type="inner") > tab12 v1 v2 [1,] 1 2.10533290 [2,] 2 -1.85441266 [3,] 3 0.95619644 [4,] 4 1.15970056 [5,] 5 -0.60433624 [6,] 6 0.73799947 [7,] 7 -1.28049285 [8,] 8 -2.63621583 [9,] 9 1.36940938 [10,] 10 -0.02570176 [11,] 11 0.99026465 Que não me dá o resultado esperado, pois esta faltando o tab2. Alguém teria alguma sugestão, Obrigado, Alexandre -- 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

5 8

Problemas r-br com yahoo mail
by ASANTOS 14 Abr '12

14 Abr '12

Caros colegas do r-br, Alguém saberia me dizer se existe algum problema do local onde esta hospedado o r-br@listas e o yahoo mail, não sei se alguém esta passando por isso, mas estou recebendo as 22h mensagens postadas as 11h e algumas mensagens minhas estão chegando a lista duplicadas. 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 1

Converter arquivos RData em csv dentro de um loop
by Marcos Silva 14 Abr '12

14 Abr '12

Olá Pessoal, Eu tenho, em um diretório, um conjunto de mais ou menos 30 arquivos RData que armazenam data frames. O que eu gostaria de fazer é escrever um scrip que lesse os nomes dos arquivos RData no referido diretorio e os exportasse no formato .csv. A abordagem que estou tentanto desenvolver é a seguinte: arquivos <- list.files() for(nomes in arquivos){ dframe <- load(nomes) load(nomes) write.csv2(????, file=paste(dframe, ".csv", sep=""), row.names=FALSE) rm(list=ls()) } O problema para o qual ainda não consegui vislumbrar uma solução é como usar o write.csv2() nesse contexto. Os arquivos RData tem nomes no seguinte padrão: "Emp_DETRAN_2009.2012fev.RData" e o data frame contido no arquivo tem nome "empDETRAN". Basicamente o que muda de um arquivo para o outro é o nome do do órgão. Possivelmente exista outra abordagem menos macarronica e mais direta ao ponto... Abs. -- Marcos F. Silva http://sites.google.com/site/marcosfs2006

4 3

arredondamento de números no R
by Vitor Aguiar 13 Abr '12

13 Abr '12

Caros, quero fazer o seguinte cálculo: pvalue = 1 - pchisq(sum(X2, na.rm = TRUE), dfreedom(X2)) #dfreedom( ) é uma função escrita por mim mas usando os meus dados isso dá sempre 0 (zero), porque a parte do pchisq( ) é um valor muito próximo de 1, então o R arredonda para 1, ficando 1 - 1 = 0. Mesmo que o verdadeiro valor desse cálculo seja algo muitíssimo pequeno, como zero seguido de vários zeros e eventualmente algum outro número, eu preciso saber qual é. Como fazer isso? Desde já agradeço a quem puder ajudar. Vitor Aguiar

5 10

Rscript
by Leandro Marino 13 Abr '12

13 Abr '12

Caros, alguem saberia me dizer aonde fica localizada a pasta Rscript em uma instalação padrão em Linux? Como posso localizá-la? Atenciosamente, Leandro Marino http://www.leandromarino.com.br (Fotógrafo) http://est.leandromarino.com.br/Blog (Estatístico) Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

6 8

ajuda com integral
by DANIELA CUNHA DA SE 13 Abr '12

13 Abr '12

Olá a todos, Estou com o seguinte problema, quero calcular a integral de 0 até HT( valor da altura total dde cada arvore presente na base de dados em anexo). Rotininha: rm(list=ls(all=T)) setwd("C:/Documents and Settings/Daniela/Desktop/Ajustes no R") mdados<-read.table("INTEGRAL.txt", h=T, dec=",");mdados # funçao de afilamento fx<-function(hi) 0.91-(1.28*(hi/HT)^1)+(2.98*(hi/HT)^2)-(4.50*(hi/HT)^3)+(2.58*(hi/HT)^4)-(0.69*(hi/HT)^5) for (i in 1:nrow(mdados))##numero total de arvores do arquivo { integral<-integrate(fx,lower=0,upper=mdados$HT) volume<-with(mdados, ((pi*(mdados$DAP^2))/40000)*integral) mdados<-data.frame(mdados, volume) } Mas não está dando certo. A equação (polinomio do 5 grau) gerada é: 0.91-(1.28*(hi/HT)^1)+(2.98*(hi/HT)^2)-(4.50*(hi/HT)^3)+(2.58*(hi/HT)^4)-(0.69*(hi/HT)^5) Objetivo, calcular somente a integral para cada árvore (ex: 33 observações), para nao precisar faze-lo manualmente ao introduzir árvore a árvore (HT). O arquivo está no excel, a integral é de 0 até HT. O DAP servirá somente para o calculo do Volume. Alguém pode me ajudar? Meu computador é Windows. Agradeço a todos, desde já! Daniela

2 1

Problema para armazenar vértices de parcelas
by ASANTOS 13 Abr '12

13 Abr '12

Boa tarde pessoal, Tenho os centróides em UTM de 11 parcelas quadradas de 40x40m e estou tentando calcular os vértices, inclusive com a repetição do primeiro vértice para fechar o polígono. Ao final, quero armazenar os resultados. Para tanto fiz: xcc<-c(672036,672332,672674,672079,671386,671364,670753,671854,669986,670860,671366)##latitude dos centróides ycc<-c(6783650,6784105,6783470,6782354,6782416,6782418,6782199,6783073,6781853,6781791,6781911)##longitude dos centróides Fazendo o cálculo dos quatro vértices para cada par de coordenada: res=NULL for(i in 1:(length(xcc))){ for(j in 1:(length(ycc))){ xa<-xcc[i]-20 ya<-ycc[j]+20 v1<-cbind(xa,ya) xb<-xcc[i]+20 yb<-ycc[j]+20 v2<-cbind(xb,yb) xc<-xcc[i]+20 yc<-ycc[j]-20 v3<-cbind(xc,yc) xd<-xcc[i]-20 yd<-ycc[j]-20 v4<-cbind(xd,yd) lim.p<-rbind(v1,v2,v3,v4,v1) write.table(lim.p, file ="parc.mod", col.names=T, row.names=F,append=T) result=read.table("parc.mod",header=F) }} O arquivo result esta saindo todo bagunçado e com mais de 60 linhas que deveriam corresponder a cabeçalho e cinco pares de coordenadas para todos os centróides, porém quando inspeciono o objeto lim.p, os cálculos estão corretos para o último par de coordenadas do centróide. Algum poderia me dizer onde posso estar errando. Obrigado, Alexandre -- 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

1 0