- R-br - DINF UFPR

Base de dados
by Edmar Caldas 13 Ago '13

13 Ago '13

Olá pessoal Queria saber se alguém tem alguma arquivo / Base de dados de exemplo de credit scoring? ou aonde poderia conseguir? Obrigado. Edmar

2 1

Re: [R-br] Tukey e plot
by Marcos Paulo Felizardo 13 Ago '13

13 Ago '13

Então, estas funções elas mudariam o gráfico, mas pelo que estudei nos documentos do R, não me trazem a possibilidade da separação das sete propriedades como citei acima. Preciso especificar melhor o meu problema. Tenho dados de área e diâmetro máximo, nas linhas que colocarei abaixo eu faço a comparação do teste e da replicata. Area<-data.frame(Teste=gl(2,465),Area=c(t2$Area,t3$Area)) factor(Area$Teste) Area.aov<-aov(Area$Area~Area$Teste) TukeyArea<-TukeyHSD(Area.aov,ordered=TRUE)plot(TukeyArea) ##O gráfico gerado é <http://www1.datafilehost.com/d/d4f9cd50> Dmax<-data.frame(Teste=gl(2,465),Dmax=c(t2$Dmax,t3$Dmax)) is.factor(Dmax$Teste) Dmax.aov<-aov(Dmax$Dmax~Dmax$Teste) TukeyDmax<-TukeyHSD(Dmax.aov,ordered=TRUE)plot(TukeyDmax) ##O gráfico gerado é <http://www1.datafilehost.com/d/98e6c918> Mas neste caso eu teria dois gráficos um para Área e outro para o Diâmetro máximo (Dmax). Para deixar mais concisa a minha análise eu gostaria de colocar estes dois dados no mesmo gráfico, e aí que fica o problema, pois ao tentar fazer isso consegui colocar os dois um do lado do outro particionando o dispositivo gráfico. Mas eu queria plotar essas duas análises como se truncando o eixo y e numa linha do eixo y apresentar os dados do TukeyArea e numa linha superior os dados de TukeyDmax utilizando a mesma escala do eixo x. Alguém sabe como me ajudar sobre funções que me possibilitam isso, ou que seria possível realizar algo semelhante. Por favor me ajudem, aqui no meu departamento não encontro ninguém que trabalha com R. Desde já agradeço. Marcos Paulo Felizardo RE: A funcao cld e plot.cld do pacote multcomp faz isso para vc. Veja documentacao Olá, peço ajuda, estou realizando uma análise inferencial e estou usando gráficos para mostrar o quanto a diferença entre as médias estão longe de 0. Com as seguintes linhas: Area<-data.frame(Teste=gl(2,465),Area=c(t2$Area,t3$Area)) factor(Area$Teste) Area.aov<-aov(Area$Area~Area$Teste) TukeyArea<-TukeyHSD(Area.aov,ordered=TRUE) plot(TukeyArea) Nestas linhas peguei a propriedade Área, mas são 7 propriedades do material analisado. Mas não sei como realizar a plotagem destas sete propriedades no mesmo gráfico. Já encontrei diferentes modos de colocar lado a lado, particionando o dispositivo gráfico, mas, a minha necessidade seria ter o gráfico com as sete propriedades comparadas nos dois testes (t2 e t3) utilizando o mesmo eixo x. Alguém pode me ajudar, por favor. Muito obrigado Marcos Paulo Felizardo

1 0

Boxplot ordem
by Polliana Zocche 13 Ago '13

13 Ago '13

Olá, novamente, Procurei nos tópicos sobre como ordenar boxplot e não encontrei a solução. Tenho um boxplot feito com a interação de duas variáveis. box2009 = boxplot(x~setor*trat, ano2009, ylab="Riqueza/m²") A ordem das caixas que vem são: "Fr.Adj" "Fu.Adj" "Fr.Dalb" "Fu.Dalb" E eu quero trocar para: "Fr.Adj", "Fr.Dalb", "Fu.Dalb", "Fu.Adj" Vi que mexendo em factor() e colocando: factor = factor(box2009, levels=c("Fr.Adj", "Fr.Dalb", "Fu.Dalb", "Fu.Adj")) Era pra dar certo, mas na verdade vem este erro: Error in model.frame.default(formula = x ~ factor, data = ano2009) : variable lengths differ (found for 'factor') Alguém tem uma sugestão de como poderia mudar essa ordem? Muito obrigada e boa noite! Polliana -- Polliana Zocche de Souza Bióloga/Mestre em Ecologia Doutoranda em Ecologia Departamento de Biologia Vegetal, IB, UNICAMP

4 4

Re: [R-br] Vista do OpenStreetMap
by Éder Comunello 13 Ago '13

13 Ago '13

Gilbert e colegas, boa tarde! Tenho esse trecho de código que pode servir pra testar. #install.packages("OpenStreetMap", dep=T) #Pacotes exigidos: c("OpenStreetMap", "rJava", "raster", "sp", "rgdal") require(OpenStreetMap) #Dourados-MS map1 <- openmap(c(lat=-21.5, lon=-56), c(lat=-23, lon=-54.5), minNumTiles=4, type="esri") map2 <- openmap(c(lat=-21.5, lon=-56), c(lat=-23, lon=-54.5), minNumTiles=4, type="esri-topo") plot(map1) plot(map2) #opcoes <- c("osm", "maptoolkit-topo", "waze", "mapquest", "mapquest-aerial", # "bing", "stamen-toner", "stamen-terrain", "stamen-watercolor", # "osm-german", "osm-wanderreitkarte", "mapbox", "esri", "esri-topo", # "nps", "apple-iphoto", "skobbler", "opencyclemap", "osm-transport", # "osm-public-transport", "osm-bbike", "osm-bbike-german") At.te, -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

1 0

ggplot2 geom_bar scale_y_continuous
by Alisson Lucrecio 13 Ago '13

13 Ago '13

Caro Colegas r-br, Boa tarde. Eu estou tentando usar esse script para plotar um gráfico de barra no ggplot2, mas quando coloco limits = c(2, 3) as barras não aparecem. Vocês sabem como resolver esse problema? ggplot(deployment, aes(x = Grupo, y = Means)) + facet_grid(. ~ Metodo) + geom_bar(stat = 2, position = "dodge") + geom_errorbar(data = deployment, aes(y=Means, ymin = min, ymax=max), width=0.1, position="identity") + geom_text(data = deployment, aes(x = Grupo, y = max, label = SK), parse=T, hjust=0.5, vjust=-0.5, size=5) + scale_y_continuous(limits = c(2, 3), breaks=seq(2,3,0.25)) Obrigado. Alisson Lucrécio da Costa

1 0

Tukey e plot
by Marcos Paulo Felizardo 13 Ago '13

13 Ago '13

Olá, peço ajuda, estou realizando uma análise inferencial e estou usando gráficos para mostrar o quanto a diferença entre as médias estão longe de 0. Com as seguintes linhas: Area<-data.frame(Teste=gl(2,465),Area=c(t2$Area,t3$Area))factor(Area$Teste)Area.aov<-aov(Area$Area~Area$Teste)TukeyArea<-TukeyHSD(Area.aov,ordered=TRUE)plot(TukeyArea) Nestas linhas peguei a propriedade Área, mas são 7 propriedades do material analisado. Mas não sei como realizar a plotagem destas sete propriedades no mesmo gráfico. Já encontrei diferentes modos de colocar lado a lado, particionando o dispositivo gráfico, mas, a minha necessidade seria ter o gráfico com as sete propriedades comparadas nos dois testes (t2 e t3) utilizando o mesmo eixo x. Alguém pode me ajudar, por favor. Muito obrigado Marcos Paulo Felizardo

2 1

Vista do OpenStreetMap
by Gilbert Queiroz 13 Ago '13

13 Ago '13

Olá pessoal!!! Estou tentando abrir uma vista do OpenStreetMap no R mas estou tendo problemas. Baixei o pacote, mas existe um erro de conexão, mas não é da internet. O que pode ser???? Abs.

2 7

Função de busca em um arquivo de dados
by Tiago Souza Marçal 13 Ago '13

13 Ago '13

Prezados, estou com uma dificuldade de criar uma função de busca em um data.frame(). A situação é a seguinte eu tenho um arquivo de dados com rótulos identificados em ordem crescente de 1 a n. carac id 1 A1 1 2 A2 2 3 A3 3 4 A4 4 5 A5 5 6 B2 6 7 B3 7 8 C4 8 9 C8 9 10 C9 10 Agora eu sorteei 3 números aleatórios das identificações (9,3,5). Gostaria que esta função capturasse os rotulos da coluna carac na sequencia dos números sorteados aleatoriamente (9,3,5) CMR: carac=c("A1","A2","A3","A4","A5","B2","B3","C4","C8","C9") id=1:10 id2=c(9,3,5) dados=data.frame(carac,id) dados Att. Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

2 2

Seleção dados em 2 categorias
by Polliana Zocche 12 Ago '13

12 Ago '13

Olá, Pessoal, Preciso de um subset de uns dados de acordo com 2 condições, mas não estou conseguindo. data09e13 ano num setor trat spp pmed 3450 2013 42 Fr Adj Oxipetalum tomentosum 2.5 3451 2013 42 Fr Adj Alternanthera maritima 10.0 3452 2013 42 Fr Adj Achyrocline satureioides 20.0 3458 2013 43 Fr Adj Dalbergia ecastophyllum 20.0 3459 2013 43 Fr Adj Ipomoea imperati 2.5 3460 2013 43 Fr Adj Oxipetalum tomentosum 2.5 3461 2013 43 Fr Adj Noticastrum malmei 20.0 3462 2013 43 Fr Adj Remiria maritima 37.5 3463 2013 43 Fr Adj Spartina ciliata 20.0 3464 2013 43 Fr Adj Alternanthera maritima 10.0 3472 2013 44 Fr Adj Achyrocline satureioides 10.0 3473 2013 44 Fr Adj Remiria maritima 10.0 3474 2013 44 Fr Adj Conyza sp. 2.5 3475 2013 44 Fr Adj Dichondra repens 2.5 3476 2013 44 Fr Adj Alternanthera maritima 20.0 3477 2013 44 Fr Adj Oxipetalum tomentosum 2.5 3478 2013 44 Fr Adj Noticastrum malmei 10.0 3479 2013 44 Fr Adj Senecio crassiflorus 2.5 3480 2013 44 Fr Adj Spartina ciliata 20.0 3481 2013 44 Fr Adj Hidrocotile bonariensis 2.5 Quero selecionar as espécies com o mesmo data09e13$num onde Dalbergia esteja presente. Ou seja, nesse conjunto de exemplo, eu quero somente os dados de data09e13$num 43, único num onde ela ocorre. 3458 2013 43 Fr Adj Dalbergia ecastophyllum 20.0 3459 2013 43 Fr Adj Ipomoea imperati 2.5 3460 2013 43 Fr Adj Oxipetalum tomentosum 2.5 3461 2013 43 Fr Adj Noticastrum malmei 20.0 3462 2013 43 Fr Adj Remiria maritima 37.5 3463 2013 43 Fr Adj Spartina ciliata 20.0 3464 2013 43 Fr Adj Alternanthera maritima 10.0 Como faço para acessar os dados dessa forma? Obrigada! Polliana. -- Polliana Zocche de Souza Bióloga/Mestre em Ecologia Doutoranda em Ecologia Departamento de Biologia Vegetal, IB, UNICAMP

4 7

Questão de concurso...
by andrebvs＠bol.com.br 12 Ago '13

12 Ago '13

4 5