Agosto 2013 - R-br

Clusters de mesmo tamanho
by Raphael Saldanha 08 Ago '13

08 Ago '13

Olá pessoal, Estou tentando procurando um modo de criar alguns clusters, mas não estou tendo sucesso. A ideia é a seguinte: tenho um dataframe chamado samp, com coordenadas geográficas (lat lon). Preciso criar n clusters com o mesmo tamanho k. Tentei algumas abordagens a partir do k-means e derivados, mas não obtive sucesso. Abaixo, coloco um exemplo para testes. library(RgoogleMaps) data(columbus) km1 <- kmeans(subset(columbus, select=c("X","Y")), 10) print(km1$size) Agradeço as ideias e sugestões. Obrigado! --- Atenciosamente,Raphael Saldanha rfsaldanha(a)outlook.com

4 4

Plot Médias diárias com ggplot2
by Jefferson Ferreira Ferreira 07 Ago '13

07 Ago '13

Pessoal, Tenho um dataframe com duas colunas mes_dia e cota, representando as médias diárias de uma varável para os 365 dias do ano (no formato mes/dia). Ao plotar ggplot(data=md, aes(x=mes_dia, y=cota, group=1)) + geom_line() os labels, obviamente ficam ilegíveis. Eu gostaria que os labels Y (cotas) fossem discretos, em intervalos, e não contendo todos os valores possíveis. Da mesma maneira que os labels X (mes_dia) contesse apenas os meses e não 365 valores. Alguma sugestão? Obrigado! -- *Jefferson Ferreira Ferreira* Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI Jefferson.ferreira(a)mamiraua.org.br *Instituto de Desenvolvimento Sustentável Mamirauá* Ministério da Ciência, Tecnologia e Inovação Telefone: +55 97 3343-9710 *Google Maps* - Mapas deste e-mail: Exibir mapa ampliado<https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> * Contatos particulares:* *(55) 9615-0100 MSN: jecogeo(a)yahoo.com.br*

3 5

Média diária de uma série temporal
by Jefferson Ferreira Ferreira 07 Ago '13

07 Ago '13

Pessoal; Estou com um probleminha aparentemente simples. Espero que alguém possa me ajudar. Tenho uma série temporal de dados, duas colunas (data e cota), de 20 anos. Estou precisando calcular a cota média de cada dia. (média dos dia 1/jan, 2/jan,....). Minha coluna de data já está as.Date() library("utils") setwd("D:/JEFFERSON/Mapeamento/Dados Hidrológicos") cotas <- read.table("cotas2.csv",sep=",",header=T,allowEscapes=T,) as.Date(cotas$data,format="%d/%m/%Y") cotas$dia <- as.Date(cotas$data,format="%d/%m/%Y") Se alguém puder me ajudar, agradeço!! Abração! -- *Jefferson Ferreira Ferreira* Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI Jefferson.ferreira(a)mamiraua.org.br *Instituto de Desenvolvimento Sustentável Mamirauá* Ministério da Ciência, Tecnologia e Inovação Telefone: +55 97 3343-9710 *Google Maps* - Mapas deste e-mail: Exibir mapa ampliado<https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> * Contatos particulares:* *(55) 9615-0100 MSN: jecogeo(a)yahoo.com.br*

3 6

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