Julho 2012 - R-br - DINF UFPR

Cálculo entre todos os pares de distâncias em um objeto espacial
by ASANTOS 14 Ago '12

14 Ago '12

Boa noite pessoal, Gostaria de calcular todas as distâncias de 6 coordenadas em UTM, à partir de um objeto da classe espacial, para isto fiz: ### Cálculo da distancia entre todos os pontos em UTM ninhos<-c("N01","N02","N03","N04","N05","N06") lat<-c(0503252,0503319,0503716,0503682,0503688,0503707) long<-c(7651950,7651694,7652672,7652670,7652627,7652662) ##Cálculo que compara a primeira coordenada com a segunda e assim sucessivamente .... lat.diff<-diff(lat) long.diff<-diff(long) d=(long.diff^2+lat.diff^2)^.5 d ##Porém não era bem isto que eu queria ainda ##Transformando em um objeto da classe espacial require(sp) pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long)) Agora, gostaria de baseado no objeto espacial pontos, calcular todas as combinações de distâncias entre as coordenadas fornecidas. Alguém saberia me indicar uma função em algum pacote que faz isso, Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

4 4

Combinar colunas de um banco de dados...
by andrebvs 03 Ago '12

03 Ago '12

5 16

barplot + lines com dois eixos Y
by Marcelo Claro de Souza 30 Jul '12

30 Jul '12

Preciso plotar um gráfico de variáveis climáticas sendo o eixo Y esquerdo referente à precipitação e o eixo Y direito referente às temperaturas min, med e max. A precipitação deve ser apresentada em barras e as temperaturas em linhas. Como posso fazer isso? x=c(seq(1:12)) rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas Muito obrigado pela ajuda. Abraço. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Importar dados EpiData
by Pedro Emmanuel Alvarenga Americano do Brasil 30 Jul '12

30 Jul '12

Amigos de R, Há um software de entrada de dados que gosto muito chamado EpidData. ( www.epidata.dk) O Epidata utilizava até o ano passado o mesmo formato de dados que o antigo EpiInfo 6.4 em DOS. Os arquivos .rec. Então, para importar dados basta: library('foreign') dados <- read.epiinfo(...) Acontece que desde o ano passado o EpiData 3.1 esta sendo abandonado e substituído por um novo EpiData que está tentando se adequar a boas praticas clínicas com acesso hierarquizado aos dados entre outras coisas. O novo formato de arquivos são versoes XML do velho .rec... os epx e epz Alguem ja conseguiu importar dados nesse formato a partir do R? Isso está causando algumas dores de cabeça, pois não estou conseguindo migrar do velho para o novo EpiData devido a essa limitação! Cheguei mandar uma mensagem numa lista da turma do core do R mas o moderador rejeita as minhas mensagens. :-( Abraço, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 2

Encontrar caminho de um arquivo no computador
by Luis Iván Ortiz Valencia 27 Jul '12

27 Jul '12

Prezados Gostaria de saber se é possível no R encontrar o caminho de um arquivo específico no computador. Por exemplo, existe um arquivo indicadores.pdf na pasta C:\ivan\R, existe uma função no R que responda com o caminho "C:\ivan\R" dando como entrada indicadores.pdf? atte IVAN

5 9

[Dúvida] Salvando data.frame em DB
by Pedro Rafael 27 Jul '12

27 Jul '12

Pessoal estou usando o comando dbWriteTable para salvar em um banco de dados PostgreSQL. O problema que mesmo colocando o argumento row.names = FALSE ele está salvando o banco com um ID. Só que eu gostaria que ele não criasse o ID pois ele já está criado nas tabelas do BD. Alguem sabe como resolver esse problema? [ ], Pedro Rafael Diniz Marinho.

1 0

regressão polinomial
by Tito Conte 26 Jul '12

26 Jul '12

Pessoal tenho o seginte conjunto de dados: f=c(-26.0000,-25.9862,-25.9343,-25.8822,-25.8433,-25.8054,-25.7948,-25.7872,-25.7668,-25.7284,-25.7015,-25.6282,-25.4612,-25.3016,-25.2564,-25.2412,-25.2412,-25.2232,-25.0869,-25.0000,-25.0000,-24.9856,-24.9397,-24.8976,-24.8533,-24.6587,-24.6373,-24.5740,-24.5406,-24.4379,-24.3934,-24.3628,-24.3277,-24.3042,-24.2972,-24.2973,-24.3020,-24.3068,-24.3068,-24.2956,-24.2920) w=c(-45.8529,-45.8302,-45.7575,-45.6680,-45.6107,-45.5050,-45.4574,-45.4058,-45.3538,-45.2913,-45.2634,-45.2264,-45.1667,-45.0719,-45.0213,-45.0000,-45.0000,-44.9748,-44.8726,-44.8405,-44.8405,-44.8351,-44.8259,-44.7972,-44.7764,-44.6867,-44.6674,-44.5912,-44.5587,-44.4927,-44.4515,-44.4115,-44.3314,-44.2664,-44.2171,-44.1729,-44.1288,-44.0820,-44.0509,-44.0118,-44.0000) se vocês plotarem eles, obteram uma curva, gostaria de obter a função desta curva ou uma regressão polinomial de grau elevado dela já tentei usar a função lm e não deu um bom resultado, quando chego em graus elevados e calculo utilizando os coeficientes obtidos os valores dão 10 vezes maiores (não sei bem porque) mas quando ploto ela aparece bonitinha no lugar que deveria aprecer em outro caso tentei na raça mesmo por multiplicação matricial mas também deu errado pois os valores da multiplicação da inversa são muito pequenos alguém pode me dar uma sugestão ou uma luz do que eu estou fazendo de errado? Tito Conte

4 14

Tabelas em sql para verficar.
by Edson Lira 25 Jul '12

25 Jul '12

Benilton e Leonardo, só agora pude voltar aquele problema do sql. Abaixo o script usado para ler o sql. rm(list = ls()) require(RODBC) options(digits=3) options(max.print=800000) getOption("max.print") channel<-odbcConnect("saddw",uid="saddw",case="tolower") #DOAÇÃO GERAL doacao.dados<- sqlQuery(channel, "select doacao_dw.numerobolsa, doacao_dw.datahoradoacao,doacao_dw.tipodoacao, doacao_dw.tipodoacaodirigida,doacao_dw.localcoleta,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.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('27/08/1992','DD/MM/YYYY') and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') <= to_date('28/08/1992','DD/MM/YYYY')" ) No link abaixo o script que gera os dados em sql. http://www.datafilehost.com/download-054b9e31.html Abaixo erro ao rodar no R. + and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') >= to_date('28/08/1992','DD/MM/YYYY') + and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') <= to_date('28/08/1992','DD/MM/YYYY')" + ) Erro em as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 2

Utilização dos pacotes snow e rmpi em cluster
by Junior Beleti 24 Jul '12

24 Jul '12

Boa noite.. Gostaria de saber se alguém possui algum material sobre utilização dos pacotes snow e rmpi em cluster, além é claro dos próprios manuais dos pacotes. Att, Junior.

1 0

Geoprocessamento
by Edmar Caldas 24 Jul '12

24 Jul '12

Olá pessoal, gostaria de saber se alguém tem algum material sobre geo pode ser o básico, grato. estou querendo aprender. Edmar

3 2