- R-br - DINF UFPR

Tinn-R 2.3.7.4 pre-release
by Jose Claudio Faria 01 Nov '12

01 Nov '12

Prezados usuários do Tinn-R, Aproveitando que ando me recuperando de uma cirurgia de hérnia inguinal dupla convencional em casa, com calma e bastante tempo (não estou indo trabalhar na Universidade), andei fazendo umas alterações no Tinn-R e creio que ficaram bem boas. Contudo, gostaria de escutar as opiniões do usuários da casa antes de liberar para todos os usuários no SourceForge. Estou usando a nova versão já a alguns dias e já testei em dois sabores do WIndows 7, tudo parece estável. Mas seria bom testes em outras plataformas e máquinas. Os uauários que quiserem testar podem entrar em contato direto comigo: joseclaudio.faria(a)gmail.com que mando o link de onde pode ser baixado. As principais alterações estão abaixo (em Inglês): 2.3.7.4 (Nov/01/2012) - The graphical interface was updated with some improvements, mainly the Application options. - A new toolbar "Format" was added to the main toolbar. - A new resource allowing reformat R code (selection or whole file) using formatR package was added. The icon resource was placed in the Format toolbar, so that from this version on, the formatR package will be necessary together with the already traditional TinnR and svSocket packages. - Due to new resource related to reformat code the variable .trPaths was changed. As a result, it will be necessary to run again the R/Configure/Permanent (Rprofile.site). In this case, do not forget to remove any prior script generated by Tinn-R in the Rprofile.site file. - The R/Hotkeys interface was deeply reworked and it now has two tabs: Default and Custom. - Default: Contains the already traditional instructions of Tinn-R; - Custom: Allows the user to customize any instructions to be send to R interpreter (thanks to Philemon Lenherr for the suggestion). The instructions must be as follows: - Simple: search(). The R interpreter will receive > search(); - Replace word or small selection: View(%s, title='Test'). If the editor cursor is over the word iris or it is selected, the R interpreter will receive > View(iris, title='Test') - Replace whole file: source(%f, echo=TRUE, verbose=TRUE). The R interpreter will receive > source(.trPaths[4], echo=TRUE, verbose=TRUE). All rules related to send file will be preserved. - Sorry, due to deeply changes in R/Hotkeys, all hotkeys configured prior to this version will be lost. It will be necessary to reconfigure all. - The chapter Some secrets for an efficient use is being revised by Ricardo Pietrobon and soon will be completed. Many thanks for his hard work! Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Problemas na integração R e Tinn-R!
by Danilo Cuzzuol Pedrini 01 Nov '12

01 Nov '12

Prezados, Instalei as últimas versões do R e do Tinn-R, ao enviar uma linha de comando do Tinn-R ao R, aparece a seguinte mensagem de erro: "source(.trPaths[4]) Erro em source(.trPaths[4]) : objeto '.trPaths' não encontrado" Atenciosamente, Danilo Cuzzuol Pedrini www.danilocp.blogspot.com

2 1

[Dúvida] Função Gama Incompleta.
by Pedro Rafael 01 Nov '12

01 Nov '12

Pessoal qual o comando certo ou pacote que trabalha com a função gamma incompleta? Estou obtendo valores muito diferentes da gamma incompleta no R e no mathematica. Por exemplo: No R: pgamma(.10,1.5) = 0.0224107 No Mathematica: N[Gamma[.10, 1.5], 10] = 0.107921 Qual seria o problema que não estou enxergando? [ ], Pedro Rafael Diniz Marinho.

2 1

Problema ao ler arquivo UTF-16 com nulos
by Roberto de Pinho 01 Nov '12

01 Nov '12

Caros, Estou tentando ler o arquivo de servidores disponível no portal da transparência. ( http://arquivos.portaldatransparencia.gov.br/PortalTransparenciaEscolheTipo… - mais de 200MB) É um CSV codificado em UTF-16. No entanto, ele reconhece prematuramente fins de linhas em algumas linhas quando encontra campos com nulo [00 00]. Há valores não nulos não sendo lidos após campos nulos que provocam o reconhecimento de fim de linha. Alguém saberia como poderia fazer com que o R possa ler este arquivo? Estou lendo o arquivo com: servidores = read.table("teste.csv", header = T, sep="\t",quote="",stringsAsFactors=T,fileEncoding="UTF-16", fill=T) o fill=T estou sando para evitar que dê erro de diferente numero de campos. Envio, em anexo, um extrato das primeiras linhas do arquivo. Obrigado pela atenção. ==== Sequencia de bytes que provocam o reconhecimento de nova linha: 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 -- Roberto de Pinho robertodepinho(a)gmail.com http://www.ascoisas.com http://lattes.cnpq.br/4816166073408660

7 20

LER CSV N0 MAC OS
by Edson Lira 31 Out '12

31 Out '12

Caros amigos da lista, meu pc é MacBookPro, OS X vefrsão 10.6.8, processador I7, com 2,7 HZ, ram de 4gb. Ao ler arquivos no formato csv, tenho tido dificuldades, no último que estou tentando ler, inclusive já postei em um email anterior uma amostra do arquivo. Quais os cuidados mínimos que devo tomar para não ter problemas? Com a rotina abaixo, h<-read.csv("/Users/edsonlira/Edinho/Edinho/Consulta/Ciesa/BSC_2307.csv",row.names = 1) Está me dando o erro abaixo: Erro em read.table(file = file, header = header, sep = sep, quote = quote, : more columns than column names Edson Lira Estatístico Manaus-Amazonas

3 2

Re: [R-br] em plot.xy(xy.coords(x, y), type = type, ...) : objeto 'membros' não encontrado
by Eder Comunello 31 Out '12

31 Out '12

Kin, boa tarde! O erro listado é por atribuir ao parâmetro pch o objeto 'membros' que não existe ou não está acessível (veja seu arquivo de dados e declarações). Acredito que a ideia seria diferenciar os pontos (tipos) de acordo com o valor em membros... Acredito que o exemplo que segue faça algo parecido com o que você deseja, embora usando somente a biblioteca sp(). Inclui os pontos que você forneceu. ########### require(sp) br.df <- read.table("http://dl.dropbox.com/u/117618178/brlite.dat", sep=";", head=T) class(br.df); head(br.df); str(br.df) br.lin <- Lines(list(Line(br.df)), ID = "Brasil") br.sl <- SpatialLines(list(br.lin)) proj4string(br.sl) <- CRS('+proj=longlat +datum=WGS84') summary(br.sl); class(br.sl) est.x <- c(-53.35,-55.60,-51.80,-51.05,-50.13,-52.41,-52.10,-52.51) est.y <- c(-33.51,-30.83,-29.45,-29.68,-29.05,-31.78,-32.03,-30.53) est.id <- c(1:8) est.df <- data.frame(est.x,est.y,est.id) est.sp <- est.df coordinates(est.sp) <- c('est.x','est.y') proj4string(br.sl) <- CRS('+proj=longlat +datum=WGS84') est.sp #png(filename="c:/temp/mapa.png",width=600,height=800) plot(br.sl, axes=T) plot(est.sp, pch=21, bg=3, add=T) #dev.off() ########### -- ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00 > Date: Tue, 30 Oct 2012 21:09:26 -0200 > From: kin Huts <kim-met(a)hotmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kim-met@hotmail.com> > > > To: <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: [R-br] em plot.xy(xy.coords(x, y), type = type, ...) : > objeto 'membros' não encontrado > Message-ID: <SNT141-W57AFF1FD679C8C33FDDE18F1620(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > > > Olá, > Gostaria de uma ajuda para entender o erro encontrado e como corrigir esse > script,Estou tentando gerar um mapa para o RS com a localização das > estações que estou utilizando. > > > require(MASS)require(maps) > > # GERAÇÃO DO MAPA COM AS ESTAÇÕES > # carregando divisão política brasileiramapaBrasil <- > read.table("BASE_BR.txt",sep="") > # carregando informações de localização das estações > meteorológicaslatlonRS <- read.csv2("coordenadas.txt",h=F,sep="") > # gerando mapa em PNGpng(filename="mapa_grupos.png",width=600,height=800) > # criando grafico em branco p/ o estado do > RSmap(xlim=c(-80,-35),ylim=c(-45,10), type="n")map.axes() > # plotando mapalines(mapaBrasil[,2], mapaBrasil[,1]) > # plotando estaçõespoints( latlonRS[,2], latlonRS[,1], pch=membros) > # fechando PNGdev.off() > > dados que estou usando são: V1 V21 -33.51 -53.352 -30.83 > -55.63 -29.45 -51.84 -29.68 -51.055 -29.05 -50.136 -31.78 -52.417 -32.03 > -52.18 -30.53 -52.51 > e o erro que aparece é: > > Carregando pacotes exigidos: MASSCarregando pacotes exigidos: mapsErro em > plot.xy(xy.coords(x, y), type = type, ...) : objeto 'membros' não > encontrado > >

1 0

Re: [R-br] coordenadas geograficas
by Eder Comunello 31 Out '12

31 Out '12

Natalia e colegas, bom dia! Ainda não estou muito familiarizado com o R, mas acho que este exemplo que rodei irá ajudar... Não utilizei o project(), mas sim o spTransform()... ############################################################# #install.packages("rgdal", dep=T) #install.packages("sp", dep=T) require(sp); require(rgdal) ### long coord.x = c(-54.607845, -54.599831, -54.430241, -54.258459, -54.344408, -54.333258, -53.607359, -53.589245, -53.546762, -53.105123, -52.972362, -52.107579, -52.089572, -51.940453, -51.718151, -51.263728, -50.741531, -50.734130, -50.667779, -50.662230, -49.986102, -49.970018, -49.901474, -49.911543, -49.678260, -49.654858, -49.610676, -49.599658, -49.511570, -49.200709, -49.222508, -49.317621, -49.311839, -49.275882, -48.581791, -48.578971, -48.582604, -48.568999, -48.528527, -48.411012, -48.251622, -48.027824, -48.029886, -48.023501, -48.590796, -48.594522, -48.643601, -49.175230, -49.555033, -49.567207, -49.689784, -49.941462, -49.942167, -50.571750, -50.591636, -50.553035, -50.718453, -51.208345, -51.216170, -51.292329, -51.239319, -51.258597, -51.283372, -51.411480, -53.551638, -53.552753, -53.643791, -53.835953, -53.842011, -53.858677, -53.955544, -54.099576, -54.108952, -54.429627, -54.438295, -54.593660, -54.593683, -54.619090, -54.617907, -54.607845) ### lat coord.y = c(-25.446263, -25.434436, -25.159459, -24.357597, -24.140630, -24.116248, -22.951327, -22.936230, -22.890090, -22.621923, -22.570479, -22.516685, -22.521237, -22.567335, -22.669007, -22.669047, -22.961387, -22.962502, -22.897375, -22.895686, -22.897501, -22.910601, -22.999147, -23.051110, -23.165054, -23.210043, -23.851684, -23.857080, -23.925858, -24.344464, -24.362727, -24.566630, -24.663615, -24.693608, -24.671609, -24.685526, -25.051554, -25.051240, -25.100784, -24.980119, -24.978739, -25.221284, -25.226189, -25.230964, -25.976982, -25.977087, -25.955943, -26.001844, -26.237326, -26.233284, -26.188456, -26.008879, -26.008876, -26.003790, -26.022134, -26.055436, -26.244779, -26.297107, -26.313987, -26.404261, -26.628383, -26.625825, -26.653990, -26.717283, -26.292744, -26.292051, -26.253072, -25.971930, -25.687925, -25.680270, -25.646984, -25.495677, -25.495742, -25.695430, -25.691293, -25.592352, -25.544262, -25.457555, -25.452072, -25.446263) pr.df <- data.frame(coord.x,coord.y); head(pr.df) ### criar um objeto "SpatialPoints" pr.spgeo <- pr.df; head (pr.spgeo) coordinates(pr.spgeo) <- c('coord.x','coord.y') ### definir projeção e datum de entrada proj4string(pr.spgeo) <- CRS('+proj=longlat +datum=WGS84') class(pr.spgeo) head(pr.spgeo) plot(pr.spgeo) ### lista parâmetros disponíveis projInfo("proj") projInfo("datum") projInfo("ellps") ### transforma para utm ### atenção para o uso do parâmetro '+south' ### ter claro que o Paraná utiliza mais de uma Zona UTM e maior parte do estado em 22S pr.sputm <- spTransform(pr.spgeo, CRS('+proj=utm +zone=22 +south')); head(pr.sputm) plot(pr.sputm) ############################################################# Esse exemplo assumiu o WGS84 como datum padrão de entrada e saída. Para datums (data) diferentes tem que dar uma estudada mais a fundo, porque tem poucas opções pré-configuradas no pacote que usei (como pode ser visto usando <projInfo("datum")>). -- ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00 > ---------------------------------------------------------------------- > > > 2012/10/29 Natalia Martins <nsmbarreto(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=nsmbarreto@gmail.com> > > > > > Walmes, bom dia. > > Estou tentando transformar minhas coordenadas geograficas (das estações > > metereologicas do PR) > > em utm, procurando no help do R encontrei: > > > > data(state) > > s <- project(state.center, "+proj=merc") > > > > mas eu não estou conseguindo entender qual proj eu devo usar, no link > > > > http://www.remotesensing.org/geotiff/proj_list/transverse_mercator.html > > > > até encontri mas quando tento fazer aparece o seguinte erro > > > > s <- project(state.center, "+proj=utm") > > Erro em project(state.center, "+proj=utm") : elliptical usage required > > > > > > se vc puder me ajudar... > > > > Muito obrigada!! > > > > Em 28 de outubro de 2012 11:11, Walmes Zeviani <walmes(a)ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br>> > escreveu: > > > > Sim, para calculo de distâncias você tem que ter coordenadas decimais e > >> não exagesimais (º ' ''). > >> > >> > ========================================================================== > >> Walmes Marques Zeviani > >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 > W) > >> Departamento de Estatística - Universidade Federal do Paraná > >> fone: (+55) 41 3361 3573 > >> VoIP: (3361 3600) 1053 1173 > >> e-mail: walmes(a)ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br> > >> skype: walmeszeviani > >> twitter: @walmeszeviani > >> homepage: http://www.leg.ufpr.br/~walmes > >> linux user number: 531218 > >> > ========================================================================== > >> > > > > > > > > -- > > > > Natália da Silva Martins > > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > > Contato: (19) 8306-4743 > > > > > > Mensagem Verificada Pelo AntiSpam UFPR > > > > >

2 1

Carregar dados pelo Rcmdr_ ubuntu 12.04
by Andre Oliveira 31 Out '12

31 Out '12

Boa noite, eu sempre uso R via terminal mesmo. Mas pensando em dar uma aula para o ADM pensei no Rcmdr. Porém estou com um problema não consigo carregar arquivos .txt, .xls e .csv. Vou em carregar arquivo seleciono o arquivo e tem o aviso. 1] NOTA: Versão do R Commander 1.8-4: Fri Oct 26 20:09:58 2012 [2] ERRO: número mágico do arquivo de restauração inválido (arquivo pode estar corrompido) -- nenhum dado carregado Alguma dica? André Oliveira Souza

4 10

Modelos econométricos no R
by Gilbert Queiroz 31 Out '12

31 Out '12

Olá, amigos e amigas, Alguém poderia me indicar alguma leitura sobre Modelos Econométricos no R??? Saudações. Gilbert.

4 4

Ajuda para plotar dados (usando o ggplot2)
by Thiago Veloso 30 Out '12

30 Out '12

Pessoal, Alguém domina bem a produção de gráficos usando o pacote ggplot2? Meus dados contém a frequência de dados em 13 classes diferentes ao longo de vários anos (simulação climática). Segue o link para baixar o arquivo contendo os dados: https://www.dropbox.com/s/50loxmtd2z2eqco/rcp26.csv (a primeira coluna pode ser ignorada). Gostaria de fazer algo parecido com esse gráfico: https://www.dropbox.com/s/8p3g42c54w3dnof/Untitled.png. Meus dados contém todas as informações para montar o gráfico: as classes (eixo x) e suas frequências (ao invés de hectares - eixo y) e os anos (três barras). Por enquanto, o código mínimo reprodutível que eu consigo enviar é apenas a abertura dos arquivos: a <- read.csv ('~/Dropbox/web/rcp26.csv') b <- read.csv ('~/Dropbox/web/rcp85.csv') Preciso de ajuda para reproduzir a imagem acima usando meus dados, com os anos 1700, 2005 e 2100 no lugar de 1985, 1993 e 2005. Agradeço desde já, Thiago.

3 5