Agosto 2013 - R-br

Gráfico
by Emerson Cotta Bodevan 14 Ago '13

14 Ago '13

Prezados membros da lista. Gostaria de saber se o R faz o gráfico que se encontra na fig 1 do artigo em anexo. International Journal of Food Microbiology Evaluation of lactic acid bacteria and yeast diversity in traditional white pickled and fresh soft cheeses from the mountain regions of Serbia and lowland regions of Croatia doi: 10.1016/j.ijfoodmicro.2013.05.032 Não sei o nome do gráfico. Agradeço qualquer ajuda. Att., * Emerson*

3 3

Experimento em parcela sub-subdividida no tempo
by Raphael Antonio Prado Dias 14 Ago '13

14 Ago '13

Olá, peço ajuda, preciso realizar uma ANOVA de um experimento no delineamento em blocos casualizados, em que as variedades foram aplicadas ao nível de parcelas, doses de hormônios foram aplicados ao nível de subparcelas e a variável responta foi medida ao longo do tempo. Alguém tem um material e/ou script para me fornecer que possa me ajudar? Obrigado e bom dia! ============================================================= Raphael Antonio Prado Dias e-mail: raphael.dias(a)muz.ifsuldeminas.edu.br <http://carpa.ciagri.usp.br/sm/src/compose.php?send_to=rapdias%40esalq.usp.br> Professor do IFSULDEMINAS - Muzambinho-MG Doutorando em Estatística e Experimentação Agronômica Universidade de São Paulo - USP Escola Superior de Agricultura “Luiz de Queiroz” – ESALQ Piracicaba - SP - Brasil =============================================================

3 2

Material sobre Regressão Logística Ordinal
by Monica Magnanini 14 Ago '13

14 Ago '13

Bom dia, Estou começando a estudar Regressão Logística Ordinal. Alguém teria um bom material para indicar? Abs Monica Magnanini IESC/UFRJ

2 2

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