- R-br - DINF UFPR

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 13 Abr '12

13 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 13 Abr '12

13 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

Pergunta boba
by Sérgio Henrique almeida da silva ju 13 Abr '12

13 Abr '12

No meu conjunto de dados têm muitos dados em branco e o R não reconhece como NA, como posso fazer para que ele transforme o "em branco" para NA? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

3 12

Regressão não-linear
by Marcelo Claro de Souza 13 Abr '12

13 Abr '12

Existe algum pacote onde o R teste automaticamente vários modelos não-lineares sugerindo o mais adequado ao meiu conjunto de dados ou preciso testar manualmente diversos modelos como Gaussiano, Bolztman, Exponencial, etc... até encontral o ideal? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 2

cut
by Vinicius Brito Rocha 13 Abr '12

13 Abr '12

Pessoal, estou usando o comando cut, para criar a legenda de uma mapa do brasil que estou plotando. O problema é que minha variável tem uma escala muito alta e por isso ao executar o comando cut estou ficando com algo parecido com isso: cut(mapa$dados),10) [1] (397,1.32e+05] (397,1.32e+05] (397,1.32e+05] [4] (397,1.32e+05] (3.95e+05,5.27e+05] (3.95e+05,5.27e+05] [7] (1.32e+05,2.64e+05] (397,1.32e+05] (1.32e+05,2.64e+05] [10] (1.32e+05,2.64e+05] (397,1.32e+05] (397,1.32e+05] O que preciso é que de uma notação mais amigável para o cut. Alguém pode me ajudar? -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 2