Junho 2012 - R-br - DINF UFPR

(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 17 Jun '12

17 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 17 Jun '12

17 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 17 Jun '12

17 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

Ajuda em uma função
by Marcelo Zeuli 15 Jun '12

15 Jun '12

Bom dia, Algume conheçe como implementar função impulso, no R? Do tipo delta de Dirac: O sinal delta de Dirac δ(t) foi criado pelo físico inglês Paul A. Maurice Dirac (1902 - 1984) e pode ser visto como o equivalente contínuo do delta de Kronecker (matemático alemão,1823 - 1991) e como o limite de uma distribuição normal Grato, Marcelo.

1 0

Ajuda em uma função
by Sérgio Henrique almeida da silva ju 15 Jun '12

15 Jun '12

Olá Gostaria de fazer uma função e como não entendo nada, gostaria de pedir ajuda de vocês Eu uso o Alfa de Cronbach para analisar a consistência interna de uma escala, normalmente se calculo o alfa total e o alfa se uma variável for retirada. O alfa total é fácil de fazer, porém quando eu vou fazer o calculo de quando retiro uma variável tenho que fazer tudo manualmente. Ex.: a = c(1,3,5,4,3,2,1,4,5,3,4,5,2,1,1,3,4,5,5,3,2) b = c(4,5,5,4,4,3,2,1,4,3,4,5,3,2,3,4,5,5,4,4,1) c = c(3,3,4,4,5,4,3,2,4,5,5,3,4,4,5,4,3,2,3,4,5) d = c(2,2,3,3,4,3,2,3,4,5,5,5,4,3,2,3,4,5,4,3,2) dados = data.frame(a,b,c,d) library(psy) # Alfa Total cronbach(dados[,1:4]) alfa total = 0.5324005 Agora preciso retirar a variável "a" e rodar novamente o calculo, depois colocar ela denovo no banco e retirar a variável "b" e assim sucessivamente, dessa forma terei o valor do alfa caso eu retire uma variável. # Alfa retirando a variável a = 0.220339 b = 0.3239247 c = 0.7669716 d = 0.2198046 Qual a melhor forma de fazer isso através de uma função? 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

2 5