- R-br - DINF UFPR

Re: [R-br] Problema com decimais
by Éder Comunello 11 Ago '13

11 Ago '13

Caro Sérgio, boa noite! Também havia me deparado com isso há algum tempo atrás e pelo que pude ver é devido ao default de dígitos significantes do comando summary(). No meu caso é 4 (7-3). ## Default S3 method: summary(object, ..., digits = max(3, getOption("digits")-3)) Você pode mudar o valor, mas me parece que não tem efeito além de 7 digitos. summary(dados, digits=7) Tem um tópico na lista internacional do R que relaciono abaixo... [R] inconsistent behavior of summary function https://stat.ethz.ch/pipermail/r-help/2011-October/291602.html At.te, -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

2 1

Função str()
by Tiago Souza Marçal 11 Ago '13

11 Ago '13

Prezados, andei dando uma olhada na saída do str() do exemplo da função aov() e fiquei com uma dúvida. Como posso definir a saída de dados para que eles apareçam de forma resumida no formato Named num como mostrado abaixo: List of 13 $ coefficients : Named num [1:13] 54.88 1.71 1.68 -1.82 -1.01 ... ..- attr(*, "names")= chr [1:13] "(Intercept)" "block1" "block2" "block3" ... $ residuals : Named num [1:24] -1.39 4.47 -5.03 1.94 -5.3 ... ..- attr(*, "names")= chr [1:24] "1" "2" "3" "4" ... $ effects : Named num [1:24] -268.83 4.84 8.23 -12.63 -9.07 ... 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 1

Problema com decimais
by Sérgio Henrique almeida da silva ju 11 Ago '13

11 Ago '13

Prezados Estou trabalhando com um banco que tem uma variável numérica com duas decimais. Quando leio o banco ele mostra essa variável com as decimais, porém toda análise que faço os resultados não mostram as decimais. O que pode estar dando errado? Ex. dados <- read.csv("BANCO.csv",dec=",", sep=";", stringsAsFactors=FALSE) dput(dados) structure(list(VAR = c("PJ", "PJ", "PF", "PF", "PJ", "PJ", "PJ", "PF", "PJ", "PJ", "PF", "PJ", "PJ", "PJ"), Saldo = c(99514,04, 49781,73, 49731,17, 49723,11, 398253,84, 49805,66, 49764,19, 100363,2, 49739,54, 99580,98, 49781,73, 198958,16, 49739,54, 49830,54)), .Names = c("VAR", "Saldo"), class = "data.frame", row.names = c(NA, -14L)) > summary(dados) VAR Saldo Length:14 Min. : 49723 Class :character 1st Qu.: 49746 Mode :character Median : 49794 Mean : 96040 3rd Qu.: 99564 Max. :398254 > tapply(dados$Saldo, list(VAR=dados$VAR), sum, na.rm=TRUE) VAR PF PJ 6111262 82400812 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) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

2 1

extrapolacao de numero de espécies
by rodrigo koblitz 11 Ago '13

11 Ago '13

Buenos, Ecologia. estou querendo fazer a previsão de espécies de uma comunidade com a 'unconditional variance'. No 'vegan' os prórpios autores advertem:´The functions are based on assumption that there is a species pool: The community is closed so that there is a fixed pool size *S_P*. Such cases may exist, although I have not seen them yet. All indices are biased for open communities.´ Procurei no BiodiversityR, mas nada. Se alguém tiver trabalhando ou souber de alguém que tenha algum script que gere essa variancia e puder diponibilizar.... bom estou interessado. abs, Koblitz

1 0

R em pesquisa operacional
by Zhu Sha Zang 09 Ago '13

09 Ago '13

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA256 Bom dia, alguém tem material específico do R para trabalhar com otimização e "tomada de decisão"? De fato, algum material onde se compara o uso do R no lugar do Solver®? Att - -- - --- Rodolfo T. da Silva -----BEGIN PGP SIGNATURE----- Version: GnuPG v2.0.20 (GNU/Linux) Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/ iF4EAREIAAYFAlIEwiYACgkQSYE5fAFtI6skegD+I9ed8J3HyPSWkPwsjxK3fzlt FzHbk7wDhTuYlrHwVPoA/3NkvC3mr49CAasl+ZjzSFOjj0+cXGMUoE6OsNcn7kX7 =BxU5 -----END PGP SIGNATURE-----

4 3

Popular um data.frame em dois loops
by Thiago V. dos Santos 09 Ago '13

09 Ago '13

Pessoal, Estou fazendo um loop para variáveis (tair, tsoil, biomass, lai) e dentro desse loop estou fazendo outro com anos (2005, 2006, 2007) etc. Ao final dos dois loops, gostaria de ter um data.frame assim: > d var date value 1 tair 2005-01-01 10 2 tair 2006-01-01 20 3 tair 2007-01-01 30 4 tair 2008-01-01 40 5 tsoil 2005-01-01 10 6 tsoil 2006-01-01 20 7 tsoil 2007-01-01 30 8 tsoil 2008-01-01 40 etc etc Qual é o meio mais eficiente para fazer isso? Talvez até exista uma forma mais rápida de fazer isso do que usar loops, então qualquer sugestão será muito bem-vinda. Segue o CMR (que não está eficiente, é só para fins de compreensão do problema): ------------------------------------------------------------------------------------------------------------- library(raster) # Create random rasterstack and set names r1 <- r2 <- r3 <- r4 <- raster(nrow = 50, ncol = 100) r1[] <- 10 r2[] <- 20 r3[] <- 30 r4[] <- 40 s <- stack(r1, r2, r3, r4) # Set years names(s) <- c('X2005.01.01', 'X2006.01.01', 'X2007.01.01', 'X2008.01.01') # Create random variable names vars <- c('tempsoi', 'tempair', 'biomass', 'lai') for( i in 1:length (vars) ) { # Loop over years (names), extract raster average values and put in the data.frame d <- NULL for ( j in 1:length (names(s)) ) { var <- vars[i] date <- format (as.Date (gsub (pattern='\\.', replacement='-', substring (names(s)[j],2)))) value <- mean (getValues (s[[j]]), na.rm = TRUE) d <- rbind (d, data.frame(var, date, value)) } # end loop over names } # end loop vars ------------------------------------------------------------------------------------------------------------- Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

3 3

Correlação não linear
by Alisson Lucrecio 09 Ago '13

09 Ago '13

Caro Colegar r-br, Boa tarde. Qual teste uso para fazer uma correlação não linear? Obrigado Alisson Lucrécio da Costa

2 1

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