- R-br - DINF UFPR

Re: [R-br] extrair bordas de arquivos .shp
by Eder Comunello 04 Nov '12

04 Nov '12

Boa tarde, Samuel e colegas! Segue uma sugestão de código... Testei com a delimitação dos bairros de Curitiba, do material do Prof. Paulo Justiniano Para testar será necessário baixar os 3 arquivos disponibilizados a seguir: https://dl.dropbox.com/u/117618178/ctba/bairros.shp' https://dl.dropbox.com/u/117618178/ctba/bairros.dbf' https://dl.dropbox.com/u/117618178/ctba/bairros.shx' Ou então baixar e descompactar: https://dl.dropbox.com/u/117618178/ctba/bairros.zip' ######################## #install.packages(c("geoR", maptools", "sp"), dep=T) require(geoR); require(maptools); require(sp) ctba <- readShapePoly("bairros.shp", IDvar="CODE") slotNames(ctba) npol <- as.integer(length(ctba@polygons)) ### numero de poligonos for (i in 1:npol){ nsubpol <- length(ctba@polygons[[i]]@Polygons) ### numero de poligonos internos for (ii in 1:nsubpol){ pol <- ctba@polygons[[i]]@Polygons[[ii]]@coords #assign(paste("P",i,ii,sep="-"),pol) assign(paste("P",formatC(i, flag=0, width=3),ii,sep="-"),pol) } } ######################## Neste caso aí serão originados 75 objetos correspondentes ao 75 polígonos. O código prevê sub-polígonos, pois é bem normal que eles existam. Não sei o que você pretende com esta operação, mas te confesso que fico pensando que ela não a melhor solução para "bordas". Sendo assim, acrescento outro trecho de código para unir todos os polígonos e feito isso, exportar uma borda única. ### une todos os poligonos e então exporta uma borda unica gpclibPermit() ## para permitir unionSpatialPolygons borda.spol <- (unionSpatialPolygons(ctba, ID=rep(1, length(ctba)), avoidGEOS=T)) borda.df <- data.frame(borda.spol@polygons[[1]]@Polygons[[1]]@coords) ######################## Estou usando R version 2.15.2 (2012-10-26) - Platform: x86_64-w64-mingw32/x64 (64-bit) Espero que ajude, -- ================================================ É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

2 1

extrair bordas de arquivos .shp
by Samuel Carvalho 04 Nov '12

04 Nov '12

Bom dia pessoal, Alguem poderia ajudar com a seguinta dúvida? Estou usando o pacote sp pra visualização e manipulação de dados espaciais no R. Um dos arquivos que estou trabalhando é do tipo .shp Qual procedimento para extrair as bordas dos polígonos deste arquivo que não seja tão lusitânica quanto o seguinte código: dados <- readShapePoly('dados.shp') length(dados@polygons) #supondo que o resultado deste comando seja == 2 length(dados@polygons[[1]]@Polygons) #e supondo que o resultado aqui seja == 10 length(dados@polygons[[2]]@Polygons) #e supondo que o resultado aqui seja == 10 Entao meu codigo para extrair as coordenadas ditas como bordas ou limites seria dados@polygons[[1]]@Polygons[[1]])@coords dados@polygons[[1]]@Polygons[[2]])@coords . . . dados@polygons[[1]]@Polygons[[10]])@coords e estes mesmos comandos para o poligono 2 o que totaliza 20 linhas de comando Talvez existisse algo como dados@borders mas não existe pois é possível verificar o que se pode extrair do objeto dados que é um objeto da classe "SpatialPolygonsDataFrame" pelo comando getSlots('SpatialPolygonsDataFrame') Obrigado Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

4 3

Fwd: Eclipse Juno
by Leandro Marino 03 Nov '12

03 Nov '12

*Izidro, não estou conseguindo acessar de jeito nenhum.* * * *Tentei também no juno e o erro se repete. Será que o servidor está fora do ar?* 2012/11/3 <izidro3.1416(a)gmail.com> > Olá Leandro > > Eu uso o eclipse Juno com o repositório > http://download.walware.de/eclipse-3.8 > > Izidro >

1 0

Eclipse Juno
by Leandro Marino 03 Nov '12

03 Nov '12

*Caros,* * * *formatei meu notebook e estou instalando alguns programas que considero adequados para o uso no dia a dia. Um deles é o StatET com o Eclipse.* * * *Baixei o eclipse e notei que é a versão 4.2 e ainda não há nada disponível no StatET. Baixei a versão indigo e tentei configurar e não ob**tive sucesso está fornecendo o seguinte erro.* * * * Unable to read repository at http://download.walware.de/eclipse-3.8/content.xml. Unable to read repository at http://download.walware.de/eclipse-3.8/content.xml. Cannot assign requested address: connect alguem sabe como posso resolver isto? *

1 0

Especificação de Pesos em glm(family=poisson)
by Diego Pujoni 02 Nov '12

02 Nov '12

Olá Pessoal, possuo uma variável preditora que é o número de espécies exóticas de peixes encontrados no lago i. Entretanto, para cada lago, eu possuo um esforço amostral diferente, isto é, existem lagos que foram amostrados apenas uma vez e lagos que foram amostrados 10 vezes. De forma que eu possuo uma incerteza maior para o lago que só foi amostrado uma única vez. Gostaria de incluir estas incertezas no ajuste e tentei utilizar o argumento "weights". Entretanto a escala dos meus pesos influencia muito o desvio residual. Existe uma forma correta para especificar estes pesos? Obrigado Sp = c(2,2,3,2,2,4,0,0,3,3,0,1,5,2,3,2,0,2,3) Peso = c(8,8,10,8,8,10,8,8,10,10,1,1,2,1,1,1,2,1,1) Var1 = c(1.48,1.65,2.50,2.20,1.48,2.29,2.83,3.59,2.29,1.65,2.70,4.94, 3.96,15.30,1.01,1.01,0.51,1.79,1.79) Var2 = c(1.48,0.00,0.00,0.00,0.00,0.00,2.83,3.59,0.00,0.00,14.99,4.94, 0.00,15.30,0.00,0.00,0.51,0.00,0.00) GLM = glm(Sp~Var1+Var2-1,family=poisson) summary(GLM) GLM1 = glm(Sp~Var1+Var2-1,weights = Peso, family=poisson) summary(GLM1) Diego PJ

2 1

Re: [R-br] selecionar intervalo de datas
by Eder Comunello 02 Nov '12

02 Nov '12

Exemplo de operação com datas com seleção por subset ############ dados<-data.frame( data=c("01/01/2004","01/02/2004","01/03/2004","01/04/2004","01/05/2004","01/06/2004"), tmin=c(17.82966,18.05029,22.47399,28.20188,13.46510,18.82721), tmax=c(44.23301,21.11526,22.79882,35.31366,35.84991,39.14286)); dados dados$data2 <- as.Date(dados$data,"%d/%m/%Y"); dados length(dados$tmax); mean(dados$tmax) data.cut <- as.Date("01/02/2004","%d/%m/%Y") novo <- subset(dados,data2>data.cut); novo length(novo$tmax); mean(novo$tmax) ############## Até mais, -- ================================================ É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

2 2

Re: [R-br] somar valores em intervalo pré-definido
by Eder Comunello 02 Nov '12

02 Nov '12

Lucas e colegas, bom dia! Tem que operar com as datas no formato interno do R, só transformando pra visualizar. Defini as datas de corte nas vars data.ini e data.fim... Alterei todos os valores de p1$num para '1', então se a soma estiver correta corresponderá ao número de eventos selecionados. Aqui deu certo! ;D #################### #Criando intervalo de datas Data<-seq(from=as.Date('01/01/1980',format='%d/%m/%Y'),to=as.Date('31/12/1980',format='%d/%m/%Y'),by='day') #transformando em data frame p1 <- data.frame(Data) p1$num <- 1 str(p1) data.ini <- as.Date('11/11/1980',format='%d/%m/%Y') data.fim <- as.Date('11/12/1980',format='%d/%m/%Y') data.dif <- data.fim-data.ini+1; data.dif #somando valores em um intervalo definido sum(p1$num[p1$Data >= data.ini & p1$Data <= data.fim]) #################### -- ================================================ É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

2 1

somar valores em intervalo pré-definido
by luc.souza＠usp.br 02 Nov '12

02 Nov '12

Bom dia. Estou precisando somar os valores de uma coluna em um intervalo pre-definido por mim. Por exemplo: tenho alguns valores referentes a algumas datas. Preciso somar alguns valores no período compreendio entre 11/11/1980 e 11/12/1980, conforme exemplo abaixo. Tentei fazer algumas coisas (abaixo) mas aparece o erro (tb abaixo). Só para lembrar, utilizo a data no formato americano e quando a faço a transformação, ela fica como character. Obriagdo pela ajuda Lucas #Criando intervalo de datas Data<-seq(from =as.Date('01/01/1980',format='%d/%m/%Y'),to=as.Date('31/12/1980',format='%d/%m/%Y'),by='day') #transformando em data frame p1 <- data.frame(Data) #colocando a data no formato americano que eu utilizo p1$Data <- as.Date(p1$Data, "%d/%m/%Y") p1$Data <- strftime(p1$Data, "%m/%d/%Y") #criando tabela de valores p1$num <- 1:366 #somando valores em um intervalo definido soma <- with(p1, sum(p1[p1$Data >= 11/11/1980 & p1$Data <= 11/12/1980, p1$num])) Error in `[.data.frame`(p1, p1$Data >= 11/11/1980 & p1$Data <= 11/12/1980, : undefined columns selected sum(which(p1$num[p1$Data >= '11/11/1980'] & p1$num[p1$Data <= '11/12/1980'])) erro: longer object length is not a multiple of shorter object length Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

1 0

selecionar intervalo de datas
by luc.souza＠usp.br 01 Nov '12

01 Nov '12

Boa tarde Possuo uma coluna com datas (m/d/aaaa), cujo formato é character. Preciso selecionar um intervalo entre duas datas. Já tentei subset, mas não funcionou. Alguma solução? obrigado Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

3 4

Problemas com ajuste de modelo ARIMAX
by Fatima do Nascimento Silva 01 Nov '12

01 Nov '12

Prezados, estou tentando ajustar um modelo ARIMAX no R, ao enviar as seguintes linhas de comando, aparece a seguinte mensagem de erro: "Erro em optim(init[mask], armaCSS, method = "BFGS", hessian = FALSE,: valor não finito fornecido por optim #linhas de comando (dados=read.table("H:vazao.txt",header=T,dec=","))#dados para ajuste do modelo attach(dados) (dados1=read.table("H:A.txt",header=F,dec=","))#data frame com as covariáveis #dimensão: [132,3] attach(dados1) library(forecast)#carregando pacote #Ajustando o modelo (mdl=arimax(dados, order = c(0, 2, 1), include.mean = T, xreg = dados1)) Atenciosamente, Fátima

2 4