Agosto 2013 - R-br

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

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 10 Ago '13

10 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