Julho 2013 - R-br - DINF UFPR

Gráfico de barras múltiplo
by Aline Lipsky 19 Jul '13

19 Jul '13

Boa tarde, Eu tenho quatro vetores referentes a biomassa de pescado mensal de quatro municípios.Quero fazer um gráfico de barras múltiplo. Se alguém puder me ajudar, Grata, Aline

2 3

Importar no R um arquivo .txt que não possui separadores
by Adriano Borges Costa 18 Jul '13

18 Jul '13

Caros, Boa tarde! Estou tendo dificuldades para importar no R um arquivo .txt bem pesado que não possui separadores. As colunas são delimitadas pelo número de caracteres em cada linha. No dicionário estão definidas que as variáveis são: VARIAVEL TAMANHO INÍCIO TIPO Variável 1 1 1 numérico Variável 2 14 2 numérico Variável 3 4 16 numérico Variável 4 10 20 numérico Ou seja, eu possuo o tamanho de cada variável (em número de caracteres) e elas estão todas uma seguida da outra. Como faço para importar esse dataset? Desde já obrigado pela ajuda! -- Adriano Borges Costa

4 5

Problema para visualizar abrir geoTif do TOPODATA no Ubuntu
by Alexandre Santos 16 Jul '13

16 Jul '13

Boa tarde Pessoal, Estou tentando abri um GeoTif de elevação obtido no banco de dados do TOPODATA no site do INPE no R usando Ubuntu OS e não consigo, dando o seguinte erro: `21S465HN.tif' does not exist in the file system, and is not recognised as a supported dataset name. Segue meu CRM: ## Obtenção da imagem no banco do Topodata url=("http://www.webmapit.com.br/inpe/topodata/download/?ms=137398652719842&fbn=2…") download.file(url, destfile = "21S465HN.zip") ## Após descoprimir o arquivo .. # #Visualizar a imagem 21S465dem<-readGDAL(21S465HN.tif) image(21S465dem) Obrigado, Alexandre

3 2

exportar para SAGAGIS classificação via writeGDAL
by Waldir De Carvalho Junior 16 Jul '13

16 Jul '13

boa tarde a Todos Estou precisando exportar o resultado da classificação via randomforest de um SpatialPixelsDataFrame, onde a predição da classe é a 17ª coluna do objeto. Quero exportar para o SAGAGIS e estou tentando pelo "writeGDAL", mas a resposta é a seguinte: > class(grid_cls) [1] "SpatialPixelsDataFrame" attr(,"package") [1] "sp" > names(grid_cls) [1] "ID" "ascii.mde30" "Flow.Direct" "Watershed.B" "Analytical" "Slope" "Aspect" [8] "Plan.Curvat" "Profile.Cur" "Convergence" "Catchment.A" "Topographic" "LS.Factor" "Channel.Net" [15] "Vertical.Di" "Relative.Sl" "prediction" > grid_cls$prediction[1:15] [1] GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ GSJ Levels: CX GM GSJ GX LA LVA PA PVA RL RY > writeGDAL(grid_cls[,17], "Soil_Class.sdat", drivername="SAGA", mvFlag=-99999) Error in create2GDAL(dataset = dataset, drivername = drivername, type = type, : Numeric bands required O erro indica que o grid (grid_cls[,17]) é uma categoria, mas é isso mesmo, não é numérico, é classe de solo. Qualquer ajuda será bem-vinda. obrigado __________ Dr. Waldir de Carvalho Junior Pesquisador PhD - Mapeamento Digital Embrapa Solos Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Rio de Janeiro/RJ Waldir.Carvalho(a)embrapa.br Telefone: +55 (21) 2179-4614 | Fax: +55 (21) 2274-5291 www.cnps.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.

3 2

[Erro] $ operator is invalid for atomic vectors
by Pedro Rafael 16 Jul '13

16 Jul '13

Pessoal estou criando uma função em R em que em determinada linha ocorre essa mensagem: $ operator is invalid for atomic vectors O trecho do código é: myintegrate = function(...) tryCatch(integrate(...), error=function(e) NA) value_int = myintegrate(f=fdp,par=starts,lower=domain[1],upper=domain[2])$value if(isTRUE(is.na(value_int))==TRUE) warning("Make sure that fdp is a probability density function. The integral in the domain specified is not convergent.") [ ], Pedro Rafael Diniz Marinho.

2 4

Re: [R-br] gráfico de barras
by Roberto Leoni 15 Jul '13

15 Jul '13

Como faço para que os valores do vetor indices apareçam em cima da barra? Obrigado pelo auxílio. att, Leoni indices <- c(33.2179623,5.9131289, 4.6236819 , 2.3701842, 2.2887414, 2.2575243, 1.5577082, 1.3616533, 1.1463521, 0.8632573, 0.8389924, 0.8232612, 0.7109162, 0.2240344, 0.1145878, -1.5093512, -3.8103876, -4.0246571) itens <- c("Tubérculos, raízes e legumes","Pescados","Farinhas, féculas e massas","Leite e derivados","Enlatados e conservas","Sal e condimentos","Artigos de limpeza","Frutas","Hortaliças e verduras", "Carnes e peixes industrializados","Bebidas e Infusões","Aves e ovos","Cereais, leguminosas e oleaginosas Higiene pessoal","Panificados","Carnes","Açúcares e Derivados","Óleos e gorduras") barplot(indices,width=0.87,legend.text=itens,main="Índices por Item - Abril/13 (período base Mar/13)",col=c("gray10","gray15","gray20", "gray25","gray30","gray35","gray40","gray45","gray50","gray55","gray60","gray65","gray70","gray75","gray80","gray85","gray90","gray95"))

3 3

Dados quantitativos em uma análise de variança
by Alisson Lucrecio 15 Jul '13

15 Jul '13

Caro Colegas r-br, Bom dia. Eu estou com uma duvida sobre dados quantitativos numa análise de variância se eu os considero como quantitativo ou como fator? Obrigado. Alisson Lucrécio da Costa

2 2

Mapas tematicos
by Natalia Martins 15 Jul '13

15 Jul '13

Prezados, bom dia. Gostaria de saber como posso trabalhar com mapas tematicos no R, pois trabalhava com o software geoda e neste conseguia criar mapas tematicos, calcular o indice de moran, LISA, entre outros. Encontrei na internet comandos para calcular o indice de moran, mas a construção dos mapas tematicos não consegui. Segue os comando que estou utilizando: ## lendo dados tipo shapefile brasil <- readShapePoly("G:/natalia/BRASIL.shp", IDvar="UF") plot(brasil) head(brasil@data) #para ler a primeira linha como cabeçalho brasil@data$NOME <- as.character(brasil@data$UF) dados <- read.table("G:/natalia/dados.csv", sep=";", head=T, enc="latin1") head(dados) UF Estado Analfabetismo Homicidio Renda 1 AC Acre 15.19 21.55768 663.2902 2 AL Alagoas 22.52 59.18695 378.0110 3 AM Amapa 7.89 30.48163 521.1194 4 AP Amazonas 9.60 26.99392 485.6758 5 BA Bahia 15.39 36.51580 463.8209 6 CE Ceara 17.19 25.34000 427.8397 #junta o shape com os dados pela varivel comum ind ind <- match(brasil@data$UF, dados$UF) dados <- dados[ind,] row.names(dados) <- brasil$UF #linca dados com formato espacial brasil <- spCbind(brasil, dados) names(brasil) [1] "UF" "ESTADO" "REGIAO" "NOME" "UF.1" [6] "Estado" "Analfabetismo" "Homicidio" "Renda" head(brasil@data) UF ESTADO REGIAO NOME UF.1 Estado Analfabetismo Homicidio Renda AC AC Acre NO AC AC Acre 15.19 21.55768 663.2902 AL AL Alagoas NE AL AL Alagoas 22.52 59.18695 378.0110 AM AM Amazonas NO AM AM Amapa 7.89 30.48163 521.1194 AP AP Amap NO AP AP Amazonas 9.60 26.99392 485.6758 BA BA Bahia NE BA BA Bahia 15.39 36.51580 463.8209 CE CE Cear NE CE CE Ceara 17.19 25.34000 427.8397 Att, Natalia

2 1

Livro - Indicação - Medidas repetidas
by Marcelo Laia 15 Jul '13

15 Jul '13

Preciso de indicação de bom livro sobre medidas repetidas no tempo. Encontrei esse: Models for Repeated Measurements (Oxford Statistical Science Series) J. K. Lindsey 1999 2ed Eu trabalho com árvores. Você conhece o livro acima? Poderia indicar um outro? Obrigado! -- Marcelo

2 1

função matlines()
by Fernando Antonio de souza 12 Jul '13

12 Jul '13

Boa tarde Estou com uma dúvida! por que ao utilizar a função matlines() a linha adcionda ao gráfico é quebrada e não contínua (parece um conjunto de retas ligando os pontos). Como eu posso solucionar isso? grato

3 7