Abril 2014 - R-br - DINF UFPR

Ler .MAP no R - Shapefile
by Sérgio Henrique almeida da silva ju 12 Abr '14

12 Abr '14

Olá Alguém sabe como faço para ler arquivos .MAP no R? Esse formato é usado pelo tabwin para ler shapefiles. Sei que existe o script abaixo http://dev-pg2.datasus.gov.br/DATASUS/tabwin/rx/Importa_mapa.htm Porém quando faço > mapa = read.map("BR.MAP") Erro em readChar(zz, 25) (from #33) : entrada UTF-8 inválida em readChar() Não roda. 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) 968463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

2 5

Erro na leitura de base em SQL
by Edson Lira 11 Abr '14

11 Abr '14

Caros amigos, estou lendo uma base de dados em SQL com a rotina abaixo #DOAÇÃO GERAL channel<-odbcConnect("saddw",uid="saddw",case="tolower") doacao.dados<- sqlQuery(channel, "select doacao_dw.numerobolsa, doacao_dw.datahoradoacao,doacao_dw.tipodoacao, doacao_dw.tipodoacaodirigida,doacao_dw.localcoleta,doacao_dw.localespecificocoleta, doacao_dw.apto, doacao_dw.codigoinaptidao,doador_dw.tipagemsanguinea,doador_dw.rh, doador_dw.idade, doador_dw.sexo, doador_dw.data_primeira_doacao, doador_dw.data_ultima_doacao,bairroresidencia, bairro_nao_tabelado,doador_dw.qtde_doacoes,doador_dw.municipioorigem, doador_dw.naturalidade, doador_dw.numerogente from saddw.doador_dw, saddw.doacao_dw where doador_dw.numerogente = doacao_dw.numerogente and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') >= to_date('01/01/1992','DD/MM/YYYY') and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') <= to_date('31/12/1999','DD/MM/YYYY')" ) Que está me gerando a mensagem abaixo: Erro em as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format Os anos com problema são: 1992, 1995 e 1996 - rodei cada ano separadamente. Como consequência, não gera o data frame doacao.dados Alguém teria alguma sugestão para corrigir esse erro? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

1 0

R+SIG
by Emerson Cotta Bodevan 11 Abr '14

11 Abr '14

Pessoal, alguém poderia me indicar um bom material sobre uso do R com SIG. Como ler e trabalhar com arquivos shape. Como converter de um sistema de coordenadas geográfica para outro. (p. ex. mudar do South America 69 para o SIRGAS 2000, ou ainda do South America para o SIRGAS UTM 23, etc). Agradeço qualquer ajuda. Att., *Emerson*

3 3

Problema "Desconhecido no R"
by Rodrigo Muniz 11 Abr '14

11 Abr '14

Caros membros da lista, gostaria mais uma vez de solicitar a ajuda de vocês. Estou tendo um problema, que não sei identificar o que é. O problema é o seguinte: Tenho um data frame 20 colunas chamado "Carbono", onde um das colunas é um fator identificando os tratamentos do experimento, e tenho um outra coluna com as classes texturais do solo. Eu fiz uma análise por classe textural, então extraí um outro data frame chamado "Arenosa", utilizando a função "subset" somente para a classe arenosa, desta forma neste segundo data frame não há todos os tratamentos que estão presentes no data frame "Carbono". Porém qualquer gráfico que eu faça como por exemplo, plot(tapply(arenosa, arenosa$Tratamento, mean), aparece no gráfico todos os tratamentos do data frame "Carbono". Esse é um exemplo, mas tenho tido problemas recorrentes referente a isto. Alguém tem alguma ideia do que possa ser? A título de informação, meu sistema operacional é Windows 7 e eu utilizo o Rstudio. Desde já agradeço... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

5 7

Repetição dentro da parcela
by ROSIANA RODRIGUES ALVES 10 Abr '14

10 Abr '14

Olá Pessoal, Estou com um problema em que minha parcela é um aquario com 15 peixes, com 3 tratamentos e 3 aquários por tratamento. Como faço para identificar no meu modelo essa repetição dentro da parcela? Att, Rosiana Rodrig ues Alves Estat ística e Experi mentação Agropecuária Em brapa Pesca e Aquicu ltura Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Palmas/TO +55 (63) 3229-7800/7850 | Skype: cnpasa.rosiana http://cnpasa.sede.embrapa.br | twitter.com/embrapa Confira também: www.facebook.com/agrosustentavel ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure.

2 1

Somar intervalo de valores nas colunas
by Alexandre Loures 10 Abr '14

10 Abr '14

Olá pessoal! Tenho uma base de dados com 216 observações de 10 variáveis. Sendo que essas observações são mensais, mas tenho que trabalhar com o total anual. Ou seja, eu tenho dados para os meses de 1997, 1998, 1999, 2000, ..., 2014. Mas quero os totais desses anos em cada coluna. Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures*

6 11

Caracteres especiais como string
by Rodrigo Muniz 10 Abr '14

10 Abr '14

Caros membros da lista, gostaria mais uma vez de solicitar a ajuda de vocês. Estou fazendo um gráfico boxplot onde eu gostaria de nomear os rótulos com o argumento ylab="δ13C", porém o resultado no plot é "d13C". Outra coisa é que o número "13" precisa estar em sobrescrito. Como eu poderia fazer esses ajustes? Desde já agradeço... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

4 4

Re: [R-br] Somar intervalo de valores nas colunas
by Robert Iquiapaza 09 Abr '14

09 Abr '14

Alexandre, O objeto que enviou não parece ter a estrutura que vc falou. Mas se seus dados estiverem ordenados por meses, > sp[c(1,5,89,200,216),] Sp1 Sp2 Sp3 Sp4 1 101.20315 96.83025 95.74570 104.38100 5 99.13003 98.05711 98.44344 97.72637 89 98.88583 101.97051 99.31058 102.20602 200 95.72225 100.38504 101.06742 102.75147 216 100.06324 100.50797 100.69345 101.10212 O código a seguir pode funcionar> mes=rep(1:12,(2014-1996)) > ano=rep(1997:2014,each = 12) > (Tot <-aggregate(sp, by=list(ano), FUN=sum, na.rm=TRUE)) Group.1 Sp1 Sp2 Sp3 Sp4 1 1997 1203.948 1197.491 1203.537 1207.432 2 1998 1201.671 1225.573 1191.417 1192.366 3 1999 1206.173 1206.108 1191.918 1218.957 4 2000 1199.445 1214.414 1197.911 1190.456 Att From: Alexandre Loures Sent: Tuesday, April 08, 2014 8:37 PM To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Somar intervalo de valores nas colunas Olá Rodrigo Coster! É isto mesmo que você solicitou? > dput(head(sp)) structure(c(38795229L, 63924605L, 0L, 0L, 0L, 0L, 33680798L, 51477786L, 0L, 0L, 0L, 0L, 32400472L, 39645654L, 0L, 0L, 0L, 0L, 68022854L, 56066200L, 0L, 0L, 0L, 0L, 31881633L, 67045744L, 0L, 0L, 0L, 0L, 452118806L, 500158705L, 0L, 0L, 0L, 0L, 21795038L, 13562364L, 0L, 0L, 0L, 0L, 2730838L, 4259594L, 0L, 0L, 0L, 0L, 103132447L, 85622631L, 0L, 0L, 0L, 0L, 60353896L, 61121454L, 0L, 0L, 0L, 0L, 32403301L, 40519152L, 0L, 0L, 0L, 0L, 116570336L, 110079499L, 0L, 0L, 0L, 0L, 166832925L, 149309546L, 0L, 0L, 0L, 0L, 78060254L, 101650482L, 0L, 0L, 0L, 0L, 25985390L, 36699369L, 0L, 0L, 0L, 0L, 69222831L, 80306324L, 0L, 0L, 0L, 0L, 4392027L, 4567085L, 0L, 0L, 0L, 0L, 33752865L, 41537988L, 0L, 0L, 0L, 0L, 80216802L, 88982746L, 0L, 0L, 0L, 0L, 1489612L, 1465175L, 0L, 0L, 0L, 0L, 60973115L, 54078115L, 0L, 0L, 0L, 0L, 28441802L, 50033001L, 0L, 0L, 0L, 0L, 60420736L, 17987780L, 0L, 0L, 0L, 0L, 28824716L, 48507914L, 0L, 0L, 0L, 0L, 444712L, 328451L, 0L, 0L, 0L, 0L, 786654L, 838397L, 0L, 0L, 0L, 0L, 12933000L, 10051098L, 0L, 0L, 0L, 0L, 39827569L, 45256447L, 0L, 0L, 0L, 0L, 58803226L, 76290276L, 0L, 0L, 0L, 0L, 46468076L, 38537192L, 0L, 0L, 0L, 0L, 53314093L, 33200142L, 0L, 0L, 0L, 0L, 34642085L, 23556561L, 0L, 0L, 0L, 0L, 119213241L, 121005345L, 0L, 0L, 0L, 0L, 274982262L, 244713051L, 0L, 0L, 0L, 0L, 74376824L, 35888744L, 0L, 0L, 0L, 0L, 2144991L, 2440378L, 0L, 0L, 0L, 0L, 114384520L, 74615518L, 0L, 0L, 0L, 0L, 54130383L, 62022477L, 0L, 0L, 0L, 0L, 3340357L, 4785296L, 0L, 0L, 0L, 0L, 3396350L, 3538952L, 0L, 0L, 0L, 0L, 7356480L, 6493192L, 0L, 0L, 0L, 0L, 30067652L, 34811828L, 0L, 0L, 0L, 0L, 10236714L, 14744233L, 0L, 0L, 0L, 0L, 21064904L, 28239710L, 0L, 0L, 0L, 0L, 5987416L, 9771117L, 0L, 0L, 0L, 0L, 16779486L, 11925052L, 0L, 0L, 0L, 0L, 9060866L, 13985017L, 0L, 0L, 0L, 0L, 15448503L, 24449350L, 0L, 0L, 0L, 0L, 8407209L, 20023742L, 0L, 0L, 0L, 0L, 38288341L, 43777071L, 0L, 0L, 0L, 0L, 610282312L, 508090709L, 0L, 0L, 0L, 0L), .Dim = c(6L, 51L), .Dimnames = list( NULL, c("sp1", "sp2", "sp3", "sp4", "sp5", "sp6", "sp7", "sp8", "sp9", "sp10", "sp11", "sp12", "sp13", "sp14", "sp15", "sp16", "sp17", "sp18", "sp19", "sp20", "sp21", "sp22", "sp23", "sp24", "sp25", "sp26", "sp27", "sp28", "sp29", "sp30", "sp31", "sp32", "sp33", "sp34", "sp35", "sp36", "sp37", "sp38", "sp39", "sp40", "sp41", "sp42", "sp43", "sp44", "sp45", "sp46", "sp47", "sp48", "sp49", "sp50", "sp51"))) Em 08-04-2014 18:37, Rodrigo Coster escreveu: Sem saber como são os teus dados fica complicado a gente de dar um comando que certamente funcionará... da um dput(head(dados)) e nos envia o que aparecer na tela 2014-04-08 18:34 GMT-03:00 Daniel C Bezerra <danielcbezerra(a)gmail.com>: apply(nomedobanco,2,sum) 2014-04-08 18:07 GMT-03:00 Alexandre Loures <alexandre.loures(a)ymail.com>: Olá pessoal! Tenho uma base de dados com 216 observações de 10 variáveis. Sendo que essas observações são mensais, mas tenho que trabalhar com o total anual. Ou seja, eu tenho dados para os meses de 1997, 1998, 1999, 2000, ..., 2014. Mas quero os totais desses anos em cada coluna. Alguém saberia como fazer isso? Desde já muito obrigado! -- Alexandre Rodrigues Loures _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. -- Alexandre Rodrigues Loures -------------------------------------------------------------------------------- _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Trocar caracteres específicos repetidos
by Waldelene Moura 07 Abr '14

07 Abr '14

Olá Sou iniciante no R, preciso retirar de uma string um caractere específico que esteja repetido, trocando estes por apenas um, Ex: a = "23345556333388" gostaria de retirar o 3 especificamente que estão repetidos, logo a string ficaria a = "2345556388" - só altero onde o 3 aparece, os outros continuam como estão. outro exemplo texto = "BOLLA CARRO LLATA CASSINO " aplicando somente para a letra L texto = "BOLA CARRO LATA CASSINO" Alguém poderia me dar uma dica de qual função usar? Waldelene Moura Analista de Sistemas Universidade Federal do Maranhão - PROEN

2 2

table with >= 2^31 elements
by Roney Fraga Souza 07 Abr '14

07 Abr '14

Pessoal, Utilizando o pacote survey resolvi estratificar a amostra de uma base de dados de com as seguintes dimensões: > dim(mydata) [1] 20274412 4 com o seguinte comando: monom00.svy <- svydesign(id=~V0300, strata=~AREAP, weights=~P001p, data=mydata) obtive a seguinte mensagem de erro: Error in table(ids[, 1], strata[, 1]) : attempt to make a table with >= 2^31 elements Com o comando traceback() tenho: 6: stop("attempt to make a table with >= 2^31 elements") 5: table(ids[, 1], strata[, 1]) 4: is.data.frame(x) 3: rowSums(table(ids[, 1], strata[, 1]) > 0) 2: svydesign.default(id = ~V0300, strata = ~AREAP, weights = ~P001p, data = mydata) 1: svydesign(id = ~V0300, strata = ~AREAP, weights = ~P001p, data = mydata) Alguém tem alguma sugestão de como lidar com isso? Informações úteis: > .Machine$integer.max [1] 2147483647 > sessionInfo() R version 3.1.0 beta (2014-03-28 r65330) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=pt_BR.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=en_US.UTF-8 LC_PAPER=pt_BR.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] splines grid tools stats graphics grDevices utils datasets methods base other attached packages: [1] plotrix_3.5-5 ggplot2_0.9.3.1 doBy_4.5-10 MASS_7.3-31 pastecs_1.3-18 [6] boot_1.3-11 psych_1.4.3 weights_0.80 gdata_2.13.3 Hmisc_3.14-3 [11] Formula_1.1-1 survival_2.37-7 lattice_0.20-27 survey_3.29-5 reshape_0.8.4 [16] plyr_1.8.1 quantreg_5.05 SparseM_1.03 ETLUtils_1.2 ffbase_0.11.3 [21] ff_2.2-12 bit_1.1-11 sqldf_0.4-7.1 RSQLite.extfuns_0.0.1 RSQLite_0.11.4 [26] DBI_0.2-7 gsubfn_0.6-5 proto_0.3-10 vimcom.plus_0.9-93 setwidth_1.0-3 [31] colorout_1.0-2 loaded via a namespace (and not attached): [1] chron_2.3-45 cluster_1.15.2 colorspace_1.2-4 dichromat_2.0-0 digest_0.6.4 [6] fastmatch_1.0-4 gtable_0.1.2 gtools_3.3.1 labeling_0.2 latticeExtra_0.6-26 [11] lme4_1.1-5 Matrix_1.1-3 minqa_1.2.3 munsell_0.4.2 nlme_3.1-117 [16] RColorBrewer_1.0-5 Rcpp_0.11.1 RcppEigen_0.3.2.1.1 reshape2_1.2.2 scales_0.2.3 [21] stringr_0.6.2 tcltk_3.1.0 Atenciosamente Roney

1 0