- R-br - DINF UFPR

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

(sem assunto)
by Diogo Braga Mendes 17 Jun '12

17 Jun '12

O problema está provavelmente na linha: *dataset = 0 * Troque por (use uma lista): *dataset = vector("list",vetZ)* e troque: dataset*[[i]]* = data.frame(read.csv(myDir[i],sep=";",header=TRUE)); -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Diogo Braga Mendes Bacharel em Estatística, CONRE 9558 Fone: (91) 8239-7799 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ “O âmago da ciência, que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata, nascida da infatigável curiosidade, flexibilidade e força da razão humana” Andrei Sakharov In Análise Harmônica de Processos Estocásticos

1 0

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

17 Jun '12

tente isso (ao invés de cbind): banco1 <- data.frame(id, ent, sai, grupo) banco2 <- data.frame(id, dia, teste) Você não precisará converte as variáveis (dia vai como fator para o banco3). -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Diogo Braga Mendes Estatístico do Instituto de Previdência e Assistência do Município de Belém - IPAMB Mestre em Estatística-PPGME/UFPA Bacharel em Estatística, CONRE 9558 Fone: (91) 8239-7799 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ “O âmago da ciência, que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata, nascida da infatigável curiosidade, flexibilidade e força da razão humana” Andrei Sakharov In Análise Harmônica de Processos Estocásticos

1 0

Dúvida sobre read.csv
by Caio Deustch 17 Jun '12

17 Jun '12

Boa noite, estou tentando fazer um contador para alocar bancos de dados em variáveis diferentes. Quando faço para 1 variável ( dataset ), esta variável aloca o banco corretamente. Quando ponho do contador, o meu vetor dataset[i] não aloca corretamente o banco. Alguma outra ideia de como fazer isso? Obrigado. myDirDefault = 'C:/Users/Noely/Documents/Caio/Bolsa/R/Caio/'; myFileConfig = paste(myDirDefault, 'Config.cfg', sep=""); z = read.table(myFileConfig, header=FALSE, sep="=") vetZ = length(z[,1]) - 2 myDir = 1:vetZ for (i in 1:(vetZ)) { myDir[i] = paste(z[1,2],z[i+1,2]) } myDir = sub("/ ","/", myDir, ignore.case=T) dataset = 0 for (i in 1:vetZ) { dataset[i] = data.frame(read.csv(myDir[i],sep=";",header=TRUE)); #quando tiro o i, funciona corretamente com a última interação do i. } [image: Imagem inline 1] -- Caio Deutsch Bacharelado em Matemática Aplicada a Negócios Faculdade de Economia Administração e Contabilidade de Ribeirão Preto - FEA -RP Universidade de São Paulo - USP Cel: 16 8812 4376

1 0

ggplot2 Browser
by Eder David Borges da Silva 16 Jun '12

16 Jun '12

Pessoal, A um tempo postei um tópico sobre gráficos animados para web, encontrei um gráfico do mesmo modo que quero fazer o meu: http://www.babynamewizard.com/voyager#ms=true&sw=f&exact=false Porem pesquisando na ggplot2 e na gridSVG não encontrei nenhum código de exemplo, se alguém tiver feito algo assim, poderia me dar um help? Att

1 0

Árvore de decisão
by Vinicius Brito Rocha 16 Jun '12

16 Jun '12

Senhores, estou trabalhando com um problema de árvore de decisão onde eu obrigatoriamente tenho que começar meu processo decisório por variáveis específicas. Por exemplo: necessariamente meu modelo precisa ser resolvido na ordem em que "inputo" as variáveis independentes model=tree(Resposta~ Var1+Var2+Var3,data=base) Primeira variável Var1, segunda variável Var2, e terceira Var3. Tenho como fazer algo parecido com isso? -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* 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.

2 2

intervalo de confiança bayesiana
by Bernardo Rodrigues 16 Jun '12

16 Jun '12

Boa noite, Sabem como posso obter intervalo de confiança bayesiano utilizando os metodos de monte carlo sendo por exemplo a distribuição a priori uniforme? valeu

2 2

Função - agrupamento de grupos
by Luciano F. Sgarbi 15 Jun '12

15 Jun '12

#olá pessoal, #Estou querendo montar uma função que me mostre pontos que estão á uma distância menor que uma pré-estabelecida. #o script e explicações seguem abaixo: ############################ ##Conjunto de dados para testar a<-rep(c(1,4),each=4) b<-c(1,2,4,5,1,2,4,5) cc<-c("p1","p2","p3","p4","p5","p6","p7","p8") mat<-data.frame(cc,a,b) rownames(mat)<-cc coord<-mat[,c(2,3)] plot(coord,pch="",xlim=c(0,6),ylim=c(0,6));text(coord[,"a"],coord[,"b"],rownames(coord))#esse é o meu conjunto de pontos. #note que se eu for juntar os pontos que estão a uma distância menor ou igual á dois, formarei dois grupos. #um formado pelos pontos: p1, p2, p3 e p4 e, outro grupo com os pontos p5, p6, p7 e p8 #################### ##iniciando os cálculos de distância library(vegan) (mat.dist<-vegdist(coord,"euclidean"))#calculo de distância entre os pontos resu.col1<-numeric(0) resu.col2<-numeric(0) (distancias<-as.numeric(mat.dist)) locais<-rownames(coord) n.linhas<-nrow(coord) ##Agora vou criar uma tabela para comparar as distâncias par-a-par for(i in 2:n.linhas){ linhas<-c(i:n.linhas) resu.col2<-c(resu.col2,locais[linhas]) resu.col1<-c(resu.col1,rep(locais[i-1],length(linhas))) } (mat.resu<-data.frame(resu.col1,resu.col2,distancias))#tabela com as distâncias par-a-par dist.corte<-2##argumentos usado para estabelecer distância para juntar pontos (min.dists<-mat.resu[which(mat.resu[,"distancias"]<=dist.corte),])##tabela com todos os pontos que estão a uma distância menor, entre si, que á pré estabelecida. #note que se for juntando os pontos forma-se dois grupos # resu.col1 resu.col2 distancias *#1 p1 p2 1#p1 com p2* *#8 p2 p3 2#o p2 que estava juntado com o p1 junta-se ao p3* *#14 p3 p4 1**#o p3 que estava juntado com o p2 junta-se ao p4* *#23 p5 p6 1#idem* *#26 p6 p7 2* *#28 p7 p8 1* ######################################## *#Alguém, por favor, sabe me dizer como de agora em diante eu faço para saber quais que são os pontos que formam cada um dos grupos?* Atenciosamente, Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

1 0

Alteração de pontos em Biplot
by Marcelo Claro de Souza 15 Jun '12

15 Jun '12

Pessoal, Normalmente eu gero meus biplots de pca segundo script abaixo: require(vegan) dados.pca=rda(dados,scale=T) simbolos=c(rep(19,X),rep(2,X),.....) plot(dados.pca,display="species",xlab="PC1 = Y%",ylab="PC2 = Y%") points(dados.pca,display="sites",pch=simbolos) legend("topright",c("nome 1","nome 2"....),pch=c(19,2,....)) Entretanto dessa vez eu preciso colocar o nome das espécies ao invés dos símbolos. Como são várias plantas da mesma espécie não da p importar o arquivo utilizando o comando row.names. Como eu faço p criar um vetor com "texto" que possa substiur o vetor simbolos? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1