- R-br - DINF UFPR

Ajuda - scatterplot (car)
by Núbia Batista 03 Ago '13

03 Ago '13

Caros, boa tarde. Utilizando o comando > scatterplot(alcalinidade~hora | amostra, data=dados, smooth=F, legend.plot=F, reg.line=F,pch=c(16,17,15),col=c("blue","green","red")) alcalinidade - y hora = horário de coleta - x amostra - fator Gera um gráfico de "dispersão", com boxplots invés de pontos, isso pq minha variável hora está sendo considerada fator, já mudei mas não estou conseguindo. Alguém sabe como forçar a utilização de pontos? Grata desde já, Núbia

2 1

truncar eixo y do gráfico
by Aline Lipsky 01 Ago '13

01 Ago '13

Boa tarde, estou fazendo um gráfico com os seguintes comandos: sp1 <- c(300, 69, 347, 42, 319, 81, 1, 100, 166, 279, 895, 334)sp2 <- c(0, 24, 0.6, 40, 640, 0, 0, 55, 75, 139, 967, 351)sp3 <- c(0, 0, 40, 6, 135, 238, 83, 150, 37, 107, 10)sp4 <- c(0, 0, 1, 14, 32, 412, 139, 38, 52, 31, 45)sp5 <-c(0, 0, 48, 0, 0, 0, 0, 41, 3, 102, 77, 0) captura.mensal <- rbind(sp1, sp2, sp3, sp4)colnames(captura.mensal) <- c("jan", "fev","mar", "abr", "mai", "jun", "jul", "ago", "set", "out", "nov", "dez") par(cex.lab = 0.8, cex.axis = 0.8, "usr", cex.main = 1.2, font.main=4)grafico.1.1 <- barplot(apply(captura.mensal, 2, sort(decreasing = TRUE), beside=TRUE, ylim = c(0,1000),ylab = "Captura (t)", col = c("seagreen4", "red3", "purple4", "violetred", "tomato2"), font.main = par ("font.main"), cex.axis=par("cex.axis"), cex.lab=par("cex.lab"), border=NA)) Eu preciso truncar o eixo y do gráfico.Gostaria de saber qual função que devo usar. Agradeço a atenção, Aline F. Lipsky

1 0

Escolha metodo de ajuste do p-value
by Fernando Antonio de souza 01 Ago '13

01 Ago '13

Caros amigos Gostaria de uma resposta para uma dúvida sobre a escolha do método de ajuste do taxa de erro (p-value) Qual critério devo utilizar para definir qual o melhor método? Ja li em discussões passadas que este método deve ser definido no momento que se pensa o projeto e não deve ser escolhido após a análise do dado para se encontrar o que se quer. No entanto como saber no momento de elaboração do projeto que por exemplo a false discovery rate (FDR) é mais adequada a meus dados que a family wise rate. Algúem sabe alguma literatura que possa me ajudar? att

2 1

Modelo não linear
by Tiago Souza Marçal 01 Ago '13

01 Ago '13

Boa noite pessoal, estou querendo modelar uma curva de embebição de sementes através do seguinte modelo não linear: y=A*(1-B*exp(-C*tempo)) + exp(D + E*(tempo-t0)) t0 = 6 Gostaria de saber se existe um procedimento parecido com o do pacote manipulate do Rstudio no R para ajustar este modelo. CMR tempo=1:12 y=c(8.81790512674818,8.83724295007473,8.84652971800299,8.85625726150244, 8.86597721721444,8.86647798384983,8.87625601672093,8.89503521901844,8.91478268591035,8.93378059512321,8.97371264813241,9.00355384268932) 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 #################################################################

5 5

Re: [R-br] Aumentar as dimensoes do dispositivo grafico
by Ivan Bezerra Allaman 29 Jul '13

29 Jul '13

Já tentou x11(w = x, h =y)?? (s,f,p) Allaman

1 0

Re: [R-br] Título do eixo x de gráfico em barras
by Ivan Bezerra Allaman 29 Jul '13

29 Jul '13

Tentei reproduzir o seu código com o intuito de ajudá-la e obtive a seguinte mensagem: Erro em barplot(media.capturada, beside = TRUE, main = "Média capturada por espécie ", : objeto 'media.capturada' não encontrado Por gentileza, leia o guia de postagem para que possamos ajudá-la. (s,f,p) Allaman

1 0

Substituir caracteres
by Luciano F. Sgarbi 29 Jul '13

29 Jul '13

Boa tarde a todos. Estou com uma série de nomes em um conjunto de dados e alguns destes nomes apresentam alguns "caracteres" estranhos. Como no exemplo abaixo "/". Eu preciso de um teste lógico que me diga se cada posição do meu vetor tem ou não o determinado caracter e o substitua. O exemplo ficaria assim: exemplo<-c("a/b","cd_mj","hoje","luciano/sgarbi")#dados de entrada #preciso substituir os "/" por "-" resultado<-c("a-b","cd_mj","hoje","luciano-sgarbi")#dados de saída Desde já, muito obrigado a todos! Abraços, Luciano -- *Luciano F. Sgarbi * * Currículo Lattes <http://lattes.cnpq.br/9788198135883354> Mestre em Ecologia e Evolução Dep. Ecologia, Instituto de Ciências Biológicas I* *Universidade Federal de Goiás, campus II, Goiânia-GO, Brasil*

2 2

Re: [R-br] Novo problema para mudar projeção de um raster [Resolvido]
by Éder Comunello 27 Jul '13

27 Jul '13

Alexandre, boa noite! Você pode tentar primeiro reprojetar o dem importado para utm e então recortar usando o limite utm (SP.utm). Desse modo, é bem provável que os limites coincidam e o SP.geo não será necessário. T+, -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

1 0

Re: [R-br] Novo problema para mudar projeção de um raster
by Eder Comunello 27 Jul '13

27 Jul '13

Alexandre e colegas, bom dia! Meu script está um pouco diferente, mas o que vai valer mesmo é o último bloco (reprojetar RasteLayer). Espero que ajude. Além disso, tenho um bookmark pra indicar: Spatial data in R: using R as a GIS (old version) < https://sites.google.com/site/rodriguezsanchezf/news/usingrasagis> Vi que saiu uma nova versão em < https://sites.google.com/site/rodriguezsanchezf/resources/rgis> Outro dia, alguém perguntava pela geração de mapas temáticos no R! Esses links podem ajudar bastante! At.te, ### <BEGIN> ### require(sp); require(rgdal); require(raster) ### Região de interesse pc <- cbind(X1=773759.1,Y1=7841546) dx <- 1500*c(-1, 1, 1,-1,-1); dx dy <- 1500*c( 1, 1,-1,-1, 1); dy reg <- NULL; reg <- rbind(cbind(pc[,1]+dx,pc[,2]+dy)); reg plot(reg, asp=1); points(pc, col=2) ### CRS UTM SP.utm <- SpatialPolygons(list(Polygons(list(Polygon(reg)), "1"))) proj4string(SP.utm) = CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs") SP.utm@bbox; SP.utm@proj4string plot(SP.utm, col=3, axes=T) ### CRS GEO (GeoTiff) SP.geo<- spTransform(SP.utm, CRS("+proj=longlat +datum=WGS84")) SP.geo@bbox; SP.geo@proj4string plot(SP.geo, col=5, axes=T) ### TOPODATA #getwd() #url=("http://www.dsr.inpe.br/topodata/data/geotiff/19S435SN.zip") #download.file(url, destfile = "19S435SN.zip") #unzip("19S435SN.zip") dem <- raster('19S435SN.tif') ### Abre o slope TOPODATA ### Recorte TOPODATA dem.crop <- crop(dem, extent(SP.geo), snap='out') image(dem.crop, asp=1) # Reprojetar RasterLayer (S4) projection(dem.crop) <- CRS("+proj=longlat +datum=WGS84") alt.utm <- projectRaster(dem.crop, crs="+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs") image(alt.utm, asp=1) ### <END> ### Em 25 de julho de 2013 22:32, ASANTOS <alexandresantosbr(a)yahoo.com.br>escreveu: > Boa noite Pessoal, > > Consegui resolver um problema e consegui outro, transformei as > coordenadas latlong de um geotiff do Topodata em utm usando spTransform(). > Como não consegui transformar diretamente um SpatialPixelsDataFrame do > topodata apenas colocando o novo CRS e em vários posts não vi nenhuma > solução para isso, tentei fazer da seguinte maneira, transformei apenas o > @coords do topodata de geo para utm e depois tentei inserir com replace() > os novos valores em utm dentro do SpatialPixelsDataFrame e ocorreu o > seguinte erro: > > Erro em replace(dem.sp@coords[1:10026, 1:2], dem.sp@coords, Utm@coords[1:10026, > : > valores negativos não são permitidos na subscrição de matriz > >

2 1

Re: [R-br] Mapas tematicos
by Éder Comunello 27 Jul '13

27 Jul '13

Olá, Natália! Bom dia! Se ainda estiver precisando das referências sobre geração de mapas temáticos, segue uma sugestão: Spatial data in R: using R as a GIS (old version) < https://sites.google.com/site/rodriguezsanchezf/news/usingrasagis> Nova versão em < https://sites.google.com/site/rodriguezsanchezf/resources/rgis> At.te, -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

3 2