Re: [R-br] sobre conjuntos

[salah via R-br]

segue sugestão


SP1 = c("a","b","c","d","g","h","i")
SP2 = c("a","b","c","d","e","f","j")
SP3 = c("a","b","g","h","e","f","k")

A = setdiff(SP1, SP2)
B = setdiff(SP2, SP3)
C = setdiff(SP3, SP1)

setdiff(A, SP3)
setdiff(B, SP1)
setdiff(C, SP2)

intersect(SP1, B)
intersect(SP2, C)
intersect(SP3, A)

intersect(intersect(SP1, SP2), SP3)


Em 31/08/2018 16:49, Antonio Silva via R-br escreveu:

Obrigado Jasmine, ajudou bastante.

Meu desafio agora é obter o que é exclusivo para uma combinação de 
dois conjuntos. Por exemplo:


> inner_join(SP1,SP2)
Joining, by = "value"
# A tibble: 4 x 1
  value
  
1 a
2 b
3 c
4 d

No entanto eu queria como resposta apenas o "c" e o "d" (exclusivos de 
SP1 e SP2) pois o "a" e o "b" também está em SP3.


Mais uma vez obrigado pela atenção.

A.O.

Em sex, 31 de ago de 2018 às 15:37, Jasmine Moreira 
> escreveu:


Olá Antônio,

Veja se ajuda:


library(dplyr)

SP1 <- c("a","b","c","d","g","h","i")
SP2 <- c("a","b","c","d","e","f","j")
SP3 <- c("a","b","g","h","e","f","k")
SP1 <- as_tibble(SP1)
SP2 <- as_tibble(SP2)
SP3 <- as_tibble(SP3)

#Para SP1
setdiff(setdiff(SP1,SP2),SP3)

#Para SP2
setdiff(setdiff(SP2,SP1),SP3)


#Para SP3
setdiff(setdiff(SP3,SP1),SP2)


O resto vc pode fazer com inter_join()

inner_join(SP1,SP2)
inner_join(SP1,SP3)
inner_join(inner_join(SP1,SP2),SP3)


Abraços,
Jasmine


> Em 31 de ago de 2018, à(s) 15:01, Antonio Silva via R-br
mailto:r-br@listas.c3sl.ufpr.br>> escreveu:
>
> Car@s
>
> Gostaria de saber se há alguma função no R, em algum pacote, que
retorne as interseções e  valores únicos de diversos conjuntos ou
grupos de subconjuntos. Por exemplo para:
> SP1 <- c("a","b","c","d","g","h","i")
> SP2 <- c("a","b","c","d","e","f","j")
> SP3 <- c("a","b","g","h","e","f","k")
>
> teríamos a saída:
> SP1: i
> SP2: j
> SP3: k
> SP1-SP2: c,d
> SP1-SP3: g,h
> SP2-SP3: e,f
> SP1-SP2-SP3: a,b
>
> Encontrei funções como intersect, setdiff e Reduce, Nem mesmo os
valores exclusivos de cada conjunto eu consegui obter.
>
> Agradeço qualquer indicação.
>
> Obrigado,
>
> --
> Antônio Olinto Ávila da Silva
> Instituto de Pesca
> São Paulo, Brasil
> ___
> R-br mailing list
> R-br@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@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@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.

Re: [R-br] Eixo x

[salah via R-br]

segue sugestões

qtidade=c(61,53,44,35,24,17,15,16,18,12,8,3,3,4,2,2,1,2,1,1,1,1,1)
Anos=c(2018,2017,2016,2015,2014,2013,2012,2011,2010,2009,2008,2007,2006,2005,2003,2002,2001,1999,1998,1997,1992,1991,1981)

plot(Anos, qtidade, xaxt = "n")
axis(1, at=Anos, labels=Anos, las=2)

## 1980-2018 = 19
plot(Anos, qtidade, xaxt = "n")
axis(1, xaxp=c(1980, 2018, 19), las=2)


Em 16/07/2018 13:19, Fabiana Goncalves Barbosa via R-br escreveu:

Olá Pessoal

Fiz um gráfico (plot) com os dados abaixo, porém no eixo x só 
apareceram alguns anos. Quero que apareçam todos os anos no eixo x. 
Como faço?


qtidade=c(61,53,44,35,24,17,15,16,18,12,8,3,3,4,2,2,1,2,1,1,1,1,1)
Anos=c(2018,2017,2016,2015,2014,2013,2012,2011,2010,2009,2008,2007,2006,2005,2003,2002,2001,1999,1998,1997,1992,1991,1981)

plot(Anos,qtidade)

Grata
Fabiana





___
R-br mailing list
R-br@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@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.

Re: [R-br] Extrair Período de datas

[salah via R-br]

olá

Bem correto seus comentários Cesar, não compreendi bem o problema.

segue sugestão:

data = seq(as.Date('2016-01-01'), as.Date('2016-12-31'), by = 1)
DF =  data.frame(data, N = 1:366)

DF[DF$data > "2016-02-05" & DF$data < "2016-02-10",]

ou

subset(DF, data > "2016-02-05" & DF$data < "2016-02-10")


saudações

___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Extrair Período de datas

[salah via R-br]

Olá

Uma das alternativas é usar a função  da biblioteca {data.table}

exemplo adaptado do site:

http://stackoverflow.com/questions/17244077/select-subset-by-date-in-r#17246437

library(data.table)

## dados
n.child = as.numeric(c(1,2,3,4,5,6))
sex = as.factor(c("f","f","f","m","m","f"))
date = as.Date(c("01/01/2002", "01/12/2002", "13/02/2003", "17/02/2003", 
"03/01/2004", "09/09/2004"), format="%d/%m/%Y")

DF = data.frame(n.child, sex, date)

## - dataframe -
fev = DF[month(DF$date) == 2,]
fev

## - data.table -
## converte para "data.table"
DT = data.table(DF)

## seleciona mes fevereiro, mes 2
fev = DT[month(date) == 2]
fev

## converte para dataframe
setDF(fev)
fev

o data.table é muito mais eficiente para manipulação de uma base de dados

saudações

___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] trabalhando com números pequenos

[salah via R-br]

segue sugestão

library(Rmpfr)

mpfr(500^(-10), precBits = 255) - mpfr(1, precBits = 255)
# 1 'mpfr' number of precision  255   bits
# [1] 
-0.998975914663505297831951000646311908486498


## caracteristica numerica da maquina
.Machine
___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Unindo colunas de um data-frame

[salah via R-br]

Caro

Segue sugestões:

## usando data.frame
DD = data.frame(A = c(3L, 1L, 3L), B = c(4L, 2L, 2L), C = c(6L, 9L, 5L))

## sugestão 1 - usando nome das colunas
DD$ABC = as.numeric(do.call(paste0, c(DD[c("A", "B", "C")])))
DD
## sugestão 2 - usando a posição das colunas
DD$ABC = as.numeric(do.call(paste0, c(DD[c(1:3)])))
DD

## usando o poderoso data.table
library(data.table)

setDT(DD) ## converte para data.table

## sugestão 3 - separado da base
ABC = DD[ ,list(ABC = as.numeric(paste0(A,B,C)))]
ABC

## sugestão 4 - inserindo na base
DD[ ,ABC := as.numeric(paste0(A,B,C))]
DD

setDF(DD) ## converte para data.frame

saudações

Em Sex, Set 2, 2016 em 11:59 , Emerson Cotta Bodevan via R-br 
 escreveu:

Pessoal, boa noite.

Tenho um conjunto de dados armazenados em um data.frame com três 
colunas, como segue:


A B C
3 4 6
1 2 9
3 2 5
. . .

Gostaria de gerar um vetor com os seguintes números:
346 129 325 ...

Ou seja, os números do vetor possuem três algarismos. Eles são 
formados pela união dos trẽs números de cada linha do data frame.


Obrigado por qualquer ajuda,

Emerson
___
R-br mailing list
R-br@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.

Re: [R-br] Sub-amostragem

[salah via R-br]

Caro

Não consegui reproduzir seu erro. Seu código com loop  funcionou 
perfeitamente

Fiz uma amostragem com 2 observações.
Encapsulei numa função simples onde se tem a opção do n amostral, 
com e sem reposição


segue sugestões:

Id = 1:2
X = 
rep(c("A","B","A","C","A","A","B","A","C","A","A","B","A","C","A","A","B","A","C","A"), 
1000)

Dados = data.frame(Id,X) ## Dados

## sugestão 1
amostragem1 = function(Dados, n=5, replace=TRUE)
{
IdDados = sample(Dados[,1], n, replace = replace)
DadosAm = NULL
for (i in 1:n) DadosAm = rbind(DadosAm, Dados[Id == IdDados[i],])
return(DadosAm)
}##end amostragem1

## amostragem com n = 5 e sem repetição
amostragem1(Dados, n=5)

## sugestão 2
amostragem2 = function(Dados, n=5, replace=TRUE)
{
IdDados = sample(Dados[,1], n, replace = replace)
return(Dados[Dados[,1] %in% IdDados,])
}##end amostragem2

## amostragem com n = 5 e sem repetição
amostragem2(Dados, n=5)

saudações

Em Seg, Ago 29, 2016 em 4:56 , Julimar Pinto via R-br 
 escreveu:
Dependendo do tipo de amostragem, costumo utilizar a package 
"sampling". Que, no caso de amostras complexas (com ou sem 
reposição), trabalha de modo bastante conjugado com a package 
"survey".


Saudações,

Julimar

Muito bom isso: Nunca votar no PSDB/DEM, muito bom mesmo!

Em 28/08/2016 20:29, "Marcos Bissoli via R-br" 
 escreveu:

Prezados,

Peço desculpas de antemão se meu problema é de extrema 
trivialidade ou se o tema já fora aqui debatido. No entanto, 
pesquisei os arquivos do fórum e não encontrei, dessa vez, uma 
solução.


Gerei um código simplificado de meu problema. Meu intuito é 
extrair uma sub-amostra de uma amostra de dados já coletados para 
fins de análise com melhor equilíbrio entre os grupos observados. 
Assim, elaborei algo semelhante com o abaixo.


Id <- 1:20
X <- 
c("A","B","A","C","A","A","B","A","C","A","A","B","A","C","A","A","B","A","C","A")

Dados <- data.frame(Id,X)
IdDados <- sample(Dados$Id,5,replace = FALSE)
DadosAm <- Dados[id==IdDados[1],]
for (i in 2:5) DadosAm <- rbind(DadosAm,Dados[id==IdDados[i],])
DadosAm

Este código funciona, e consigo extrair uma sub-amostra aleatória 
de n1=5 a partir de uma amostra inicial de n=20 em data.frame.


No entanto, meu real problema é gerar uma sub-amostra de n1=88 em 
uma amostra inicial de n=1668. Mas, quando tento fazer com tais 
dimensões a sub-amostra gera uma série de NA's, que não existem 
na amostra original.


Chequei o funcionamento na amostra real e maior, e percebi que o 
primeiro Id de IdDados não corresponde ao Id de Dados adicionado 
já no primeiro comando de criação de DadosAm.


Creio que seja algum erro meu de implementação, mas cheguei a 
fazer testes com simulações de n e n1 maiores e deram certo. Mas 
quando vou para o meu banco real permanece o problema.


Desde já, agradeço qualquer ajuda, e reitero minhas desculpas pela 
possível trivialidade da dúvida.


Saudações acadêmicas,


--
MARCOS BISSOLI

Faculdade de Nutrição
Universidade Federal de Alfenas

Blog: bocademiamaldita.blogspot.com/
E-mail: mbiss...@gmail.com
Twitter: #mbissoli

Alfenas, Minas Gerais, Brasil


*Pense na Natureza antes de Imprimir*
Divulgue ON-LINE

Eu apoio a ENEN "na luta por um Brasil sem fome"

"por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la 
esperantan"
(para cada povo sua própria língua, para todos os povos o 
Esperanto)


E nunca votarei no PSDB/DEM!

___
R-br mailing list
R-br@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@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.

Re: [R-br] Sub-amostragem

[salah via R-br]

Caro,

segue sugestão:

Id <- 1:20
X <- 
c("A","B","A","C","A","A","B","A","C","A","A","B","A","C","A","A","B","A","C","A")

Dados <- data.frame(Id,X)
IdDados <- sample(Dados$Id,5,replace = FALSE)

DadosAm = subset(Dados, Id %in% IdDados)
DadosAm

saudações

Em Dom, Ago 28, 2016 em 8:29 , Marcos Bissoli via R-br 
 escreveu:

Prezados,

Peço desculpas de antemão se meu problema é de extrema 
trivialidade ou se o tema já fora aqui debatido. No entanto, 
pesquisei os arquivos do fórum e não encontrei, dessa vez, uma 
solução.


Gerei um código simplificado de meu problema. Meu intuito é extrair 
uma sub-amostra de uma amostra de dados já coletados para fins de 
análise com melhor equilíbrio entre os grupos observados. Assim, 
elaborei algo semelhante com o abaixo.


Id <- 1:20
X <- 
c("A","B","A","C","A","A","B","A","C","A","A","B","A","C","A","A","B","A","C","A")

Dados <- data.frame(Id,X)
IdDados <- sample(Dados$Id,5,replace = FALSE)
DadosAm <- Dados[id==IdDados[1],]
for (i in 2:5) DadosAm <- rbind(DadosAm,Dados[id==IdDados[i],])
DadosAm

Este código funciona, e consigo extrair uma sub-amostra aleatória 
de n1=5 a partir de uma amostra inicial de n=20 em data.frame.


No entanto, meu real problema é gerar uma sub-amostra de n1=88 em 
uma amostra inicial de n=1668. Mas, quando tento fazer com tais 
dimensões a sub-amostra gera uma série de NA's, que não existem na 
amostra original.


Chequei o funcionamento na amostra real e maior, e percebi que o 
primeiro Id de IdDados não corresponde ao Id de Dados adicionado já 
no primeiro comando de criação de DadosAm.


Creio que seja algum erro meu de implementação, mas cheguei a fazer 
testes com simulações de n e n1 maiores e deram certo. Mas quando 
vou para o meu banco real permanece o problema.


Desde já, agradeço qualquer ajuda, e reitero minhas desculpas pela 
possível trivialidade da dúvida.


Saudações acadêmicas,


--
MARCOS BISSOLI

Faculdade de Nutrição
Universidade Federal de Alfenas

Blog: bocademiamaldita.blogspot.com/
E-mail: mbiss...@gmail.com
Twitter: #mbissoli

Alfenas, Minas Gerais, Brasil


*Pense na Natureza antes de Imprimir*
Divulgue ON-LINE

Eu apoio a ENEN "na luta por um Brasil sem fome"

"por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la 
esperantan"

(para cada povo sua própria língua, para todos os povos o Esperanto)

E nunca votarei no PSDB/DEM!
___
R-br mailing list
R-br@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.

Re: [R-br] Contador por ordem de entrada

[salah via R-br]

sugestão:

usando data.table

library(data.table)

DD = data.frame(
id = c(101L, 101L, 101L, 101L, 102L, 102L, 103L, 104L, 104L, 
104L),
nome = c("jose", "jose", "jose", "jose", "ana", "ana", "mara", 
"ze", "ze", "ze"),

data_inicio = c(
"01/11/2013", "01/05/2014", "01/09/2014", "01/04/2015", 
"01/09/2015",
"01/09/2016", "01/09/2014", "01/04/2015", "01/09/2015", 
"01/09/2016"),

data_fim = c(
"01/03/2014", "01/08/2014", "01/10/2014", "01/06/2015", 
"01/03/2016",
"01/03/2017", "01/10/2014", "01/06/2015", "01/03/2016", 
"01/03/2017"),

contador = c(1L, 2L, 3L, 4L, 1L, 2L, 1L, 1L, 2L, 3L))

setDT(DD)##converte para data.table

## o argumento (.N) conta o numero de elementos (by) pela coluna id
DD[, .N, by = id]

## o atributo (:=) insere uma nova coluna contando de 1 a .N por id
DD[, Cont := 1:.N, by = id]
DD

setDF(DD)##converte para data.frame
DD

saudações

Em 13/07/2016 13:42, Wagner Tassinari via R-br escreveu:

Olá pessoal,

Gostaria de criar um contador para saber quantas vezes o registro está 
entrando no banco, por exemplo


id   nome  data_inicio   data_fim  contador
101  jose  01/11/2013   01/03/2014  1
101  jose  01/05/2014   01/08/2014  2
101  jose  01/09/2014   01/10/2014  3
101  jose  01/04/2015   01/06/2015  4
102  ana  01/09/2015   01/03/2016   1
102  ana   01/09/2016   01/03/2017  2
103  mara 01/09/2014 01/10/2014  1
104  ze 01/04/2015 01/06/2015  1
104  ze 01/09/2015 01/03/2016  2
104  ze 01/09/2016   01/03/2017  3

Obrigado
-
-
Wagner S. Tassinari
Departamento de Matemática
Universidade Federal Rural do Rio de Janeiro.
BR-465, Km 7 - Seropedica, RJ - Brasil
CEP: 23890-000
Cel: (21) 96488-5982 (WhatsApp)
Skype: wagner.tassinari
wtassin...@gmail.com 
--



___
R-br mailing list
R-br@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@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.

Re: [R-br] Simbolo

[salah via R-br]

sugestão:

x = 1:10
y = 25:34
plot(x ~ y, ylab = expression(theta~degree~C))
ou
plot(x ~ y, ylab = expression(Temperatura~degree~C))

saudações

Em 13/07/2016 10:50, Marcus Nunes via R-br escreveu:

Fabiana,

Tente rodar o código abaixo e veja se ele funciona. Se funcionar, 
adapte-o para o seu contexto.


set.seed(123)
x <- rnorm(25)
y <- x+rnorm(25)
plot(y ~ x, xlab="Temperatura (ºC)", ylab="Riqueza")

2016-07-13 10:26 GMT-03:00 Fabiana Goncalves Barbosa via R-br 
>:


Olá Pessoal

Fiz um diagrama de dispersão (Riqueza /versus/ temperatura) e
coloquei os nomes dos eixos y = Riqueza e x = temperatura.
A unidade da temperatura é °C, assim, como faço para inserir o
símbolo °?

__


___
R-br mailing list
R-br@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.




--
Marcus Nunes
http://marcusnunes.me/


___
R-br mailing list
R-br@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@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.

Re: [R-br] Instalação de pacote RCurl no R [Linux 16.04]

[salah via R-br]

olá

segundo a linha de erro: checking for curl-config... no
você precisa instalar o |curl-config |
no terminal linux digite:

sudo apt-get install libcurl4-gnutls-dev

No R instale o Rcurl

saudações|
|||

Em 11/07/2016 12:39, Paulo Henrique Pimenta via R-br escreveu:

Olá Listeiros,

Estou tendo dificuldades para instalar o pacote Rcurl no meu R [usando 
a IDE RStudio e o Linux Ubuntu 16.04].


Ao fazer a instalação no RStudio recebo a seguinte mensagem:

2016-07-11 12:19:33 (494 KB/s) - 
“/tmp/RtmpeRnbpy/downloaded_packages/RCurl_1.95-4.8.tar.gz” salvo 
[916934/916934] * installing *source* package ‘RCurl’ ... ** package 
‘RCurl’ successfully unpacked and MD5 sums checked checking for 
curl-config... no Cannot find curl-config ERROR: configuration failed 
for package ‘RCurl’ * removing 
‘/home/paulo/R/x86_64-pc-linux-gnu-library/3.2/RCurl’ Warning in 
install.packages : installation of package ‘RCurl’ had non-zero exit 
status The downloaded source packages are in 
‘/tmp/RtmpeRnbpy/downloaded_packages’


>


E ao tentar fazer a instação do pacote usando o terminal do Linux 
[baseado no tutorial: 
http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase36.html 
] 
recebo a seguinte mensagem:



root@tolstoi:/tmp/RtmpeRnbpy/downloaded_packages# R install 
RCurl_1.95-4.8.tar.gz

ARGUMENTO 'install' __ignorado__

ARGUMENTO 'RCurl_1.95-4.8.tar.gz' __ignorado__

Espero que vcs possam me ajudar!

Um grande abraço a todos.

Paulo Henrique de A. S. Pimenta.

Graduando em Meteorologia (Bacharelado) - IAG/USP.
Fone: +5511981318435.


___
R-br mailing list
R-br@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@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.

Re: [R-br] subset de observações específicas de uma variável

[salah via R-br]

sugestão:

DADOS[nchar(DADOS$HScode)  == 4,]

ou

subset(DADOS, nchar(HScode) == 4)

saudações

Em 07/07/2016 23:36, Alexandre Loures via R-br escreveu:

Boa noite!

Tenho uma base de dados em que uma das variáveis (HScode) possui três 
níveis. Por exemplo,


01
0101
010122
0102
010233
02
0201
020179

Contudo preciso de um subset apenas com as observações com quatro 
dígitos. Ou seja,


0101
0102
0201

structure(list(reporter = c("Albania", "Albania", "Albania",
"Albania", "Albania", "Albania", "Albania", "Albania", "Albania",
"Albania", "Albania", "Albania", "Albania", "Albania", "Albania",
"Albania", "Albania", "Albania", "Albania", "Albania", "Albania",
"Albania", "Albania", "Albania", "Albania", "Albania", "Albania",
"Albania", "Albania", "Albania", "Albania", "Albania", "Albania",
"Albania", "Albania"), HScode = c("01", "0101", "010110", "010190",
"0102", "010210", "010290", "0103", "010310", "010391", "010392",
"0104", "010410", "010420", "0105", "010511", "010512", "010519",
"010594", "010599", "0106", "010611", "010612", "010619", "010620",
"010631", "010632", "010639", "010690", "02", "0201", "020110",
"020120", "020130", "0202"), boundtariff = c("6.571969696969697",
"8.75", "7.5", "10", "7.2916667", "5", "9.584",
"5", "5", "5", "5", "5", "5", "5", "9.5", "5", "5", "5", "12.5",
"20", "5", "5", "5", "5", "5", "5", "5", "5", "5", "10.27914614121511",
"10", "10", "10", "10", "10")), datalabel = "", time.stamp = " 7 Jul 
2016 22:53", .Names = c("reporter",
"HScode", "boundtariff"), formats = c("%54s", "%9s", "%20s"), types = 
c(54L,

7L, 20L), val.labels = c("", "", ""), var.labels = c("", "",
""), version = 12L, row.names = c("1", "2", "3", "4", "5", "6",
"7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17",
"18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28",
"29", "30", "31", "32", "33", "34", "35"), class = "data.frame")


Alguém saberia como fazer isso?


Desde já muito obrigado!



--
*Alexandre Rodrigues Loures*
Doutorando em Economia Aplicada
Universidade Federal da Paraíba - UFPB
Centro de Ciências Sociais Aplicadas - CCSA
Programa de Pós-Graduação em Economia - PPGE
Site: www.ccsa.ufpb.br/ppge
orcid www.orcid.org/-0002-1288-0135


___
R-br mailing list
R-br@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@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.

Re: [R-br] Inputando registros

[salah via R-br]

olá

segue sugestões:

## usando data.frame
DF = data.frame(
var1 = c(20L, 40L, 40L, 40L, 50L, 50L, 60L),
var2 = c(25L, NA, 45L, NA, 55L, NA, 60L))

DF[is.na(DF$var2),] = DF[is.na(DF$var2),][1]
DF

## usando data.table
library(data.table)

##  data.table
DT = data.table(
var1 = c(20L, 40L, 40L, 40L, 50L, 50L, 60L),
var2 = c(25L, NA, 45L, NA, 55L, NA, 60L))

DT[is.na(var2), var2 := var1]
DT

saudações

Em 30/06/2016 13:54, Wagner Tassinari via R-br escreveu:

Olá pessoal, tenho o seguinte banco de dados:

var1 var2
20 25
40
40 45
40
50 55
50
60 60


Quero criar uma var3 que seja igual a var2 mas onde está em branco 
impute os valores de var1, ex:


var3
25
40
45
40
55
50
60

Estou fazendo da seguinte forma:

banco$var1=as.character(banco$var1)
banco$var2=as.character(banco$var2)

banco$var3= banco$var2
banco$var3[banco$var2 == " "] = banco$var1

Mas não esta dando certo !!

Obrigado pessoal

-
Wagner S. Tassinari
Departamento de Matemática
Universidade Federal Rural do Rio de Janeiro.
BR-465, Km 7 - Seropedica, RJ - Brasil
CEP: 23890-000
Cel: (21) 96488-5982 (WhatsApp)
Skype: wagner.tassinari
wtassin...@gmail.com 
--



___
R-br mailing list
R-br@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@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.

Re: [R-br] RES: Operação entre datas

[salah via R-br]

Boa ideia Éder!

Considerando que se deseja apenas o cálculo dos anos

segue outra ideia

## data.frame
DADOS = data.frame(
NASCIMENTO = c("26/11/97", NA, "14/05/99","06/05/98"),
DATA_FATO  = c("23/02/15", "28/03/15", "08/04/15", "08/04/15"))

within(DADOS,
{
## converter para Date
NASCIMENTO = as.Date(NASCIMENTO, "%d/%m/%y")
DATA_FATO = as.Date(DATA_FATO, "%d/%m/%y")
IDADE = as.numeric(format(DATA_FATO, "%Y")) - 
as.numeric(format(NASCIMENTO, "%Y"))

## opcional
## modificando o formato de saida da data
NASCIMENTO = format(NASCIMENTO, "%d/%m/%Y")
DATA_FATO = format(DATA_FATO, "%d/%m/%Y")
})##end within

saudações

Em 24/06/2016 09:41, Éder Comunello escreveu:

Uma ideia...

### 
if(!require(eeptools)){install.packages("eeptools"); require(eeptools)}

DADOS <- data.frame(
  NASCIMENTO = c("26/11/97", "20/03/99", "14/05/99","06/05/98", 
"03/01/00", "19/05/97", "01/02/01", "28/11/97", "10/02/00"),
  DATA_FATO =  c("26/11/07", "20/04/09", "14/06/09","06/10/08", 
"13/01/10", "29/05/07", "01/02/11", "28/11/07", "10/02/10"))


DADOS$NASCIMENTO <- as.Date(DADOS$NASCIMENTO, "%d/%m/%y")
DADOS$DATA_FATO  <- as.Date(DADOS$DATA_FATO,  "%d/%m/%y")

DADOS$IDADE <-  round(age_calc(DADOS$NASCIMENTO, DADOS$DATA_FATO, 
units = "years"), 2)

DADOS

# Os warnings são devidos à forma de implementação da função, que só 
permite verificar a condição do primeiro elemento (enddate < dob)

if (DADOS$DATA_FATO < DADOS$NASCIMENTO) {print("ok")}

# A forma abaixo suprime os warnings por realizar o cálculo 
"independente" de cada elemento
sapply(1:9, function(x) round(age_calc(DADOS$NASCIMENTO[x], 
DADOS$DATA_FATO[x], units = "years"), 2))


### simulando o erro...
DADOS$NASCIMENTO[c(3,6,7)] <- NA
DADOS$IDADE <- NA
DADOS$IDADE <- round(age_calc(DADOS$NASCIMENTO, DADOS$DATA_FATO, units 
= "years"), 2) # ERRO!

# Error in seq.int <http://seq.int>(r1$year, to$year, by) :
#   'from' cannot be NA, NaN or infinite
# In addition: Warning message:
#   In if (enddate < dob) { :
#   the condition has length > 1 and only the first element will 
be used

DADOS

### solução
sel <- complete.cases(DADOS[,1:2])
DADOS$IDADE[sel] = round(age_calc(DADOS$NASCIMENTO[sel], 
DADOS$DATA_FATO[sel], units = "years"), 2)

DADOS
# NASCIMENTO  DATA_FATO IDADE
# 1 1997-11-26 2007-11-26 10.00
# 2 1999-03-20 2009-04-20 10.08
# 32009-06-14NA
# 4 1998-05-06 2008-10-06 10.42
# 5 2000-01-03 2010-01-13 10.03
# 62007-05-29NA
# 72011-02-01NA
# 8 1997-11-28 2007-11-28 10.00
# 9 2000-02-10 2010-02-10 10.00
### 




Éder Comunello
Researcher at Brazilian Agricultural Research Corporation (Embrapa)
DSc in Agricultural Systems Engineering (USP/Esalq)
MSc in Environ. Sciences (UEM), Agronomist (UEM)
---
Embrapa Agropecuária Oeste, Dourados, MS, Brazil ||

GEO, -22.2752, -54.8182, 408m
UTC-04:00 / DST: UTC-03:00




Em 23 de junho de 2016 22:42, Amikobh via R-br 
<r-br@listas.c3sl.ufpr.br <mailto:r-br@listas.c3sl.ufpr.br>> escreveu:


Ei! Eis os scripts e a saída de erro:

>if(!require(eeptools)){install.packages("eeptools");
require(eeptools)}

>banco$NASCIMENTO = as.Date(banco$NASCIMENTO, "%d/%m/%y")
>banco$DATA_ATO = as.Date(banco$DATA_ATO, "%d/%m/%y")
## calcula o intervalo em anos
>banco$IDADE2 = floor(age_calc(banco$NASCIMENTO, banco$DATA_ATO,
units = "years"))

    Error in seq.int <http://seq.int>(r1$year, to$year, by) :
  'from' cannot be NA, NaN or infinite
Além disso: Warning message:
In if (enddate < dob) { :
  a condição tem comprimento > 1 e somente o primeiro elemento
será usado

Existem linhas do banco$NASCIMENTO  vazias.
Quando uso os scripts sugeridos pelo Eder dá certo sendo que
aparece NA em banco$idade2 na correspondente que contém NA em
banco$NASCIMENTO...

Muito obrigado.


Em 22-06-2016 13:55, salah via R-br escreveu:

Olá Amiko

É necessário que você disponibilize um trecho do seu script ou
pelo menos  a saída de erro
Executei o meu script e do Eder e ambos obtiveram exito

Na conversão de datas usando a função <as.Date()> o padrão de
conversão tem que ser exatamente com o formato da data se não
retorna NA

exemplo:

as.Date(c("26/11/97"), "%d/%m/%y"

Re: [R-br] RES: Operação entre datas

[salah via R-br]

Olá Amiko

É necessário que você disponibilize um trecho do seu script ou pelo 
menos  a saída de erro

Executei o meu script e do Eder e ambos obtiveram exito

Na conversão de datas usando a função <as.Date()> o padrão de conversão 
tem que ser exatamente com o formato da data se não retorna NA


exemplo:

as.Date(c("26/11/97"), "%d/%m/%y")## padrão da data como "dia/mes/ano"
## [1] "1997-11-26"

as.Date(c("26/11/97"), "%d-%m-%y")## padrão da data como "dia-mes-ano"
## NA

usando o lubridate() como sugeriu o Eder

library(lubridate)
dmy("26/11/97")## não é necessário um padrão, mas tem que estar como dia 
mes e ano

## [1] "1997-11-26"

dmy("26-11-97")
##[1] "1997-11-26

dmy("11/26/97")
## NA

saudações

Em 22/06/2016 09:28, Amiko Bh via R-br escreveu:


Prezado,

Quando executo o script para calcular o período retorna erro por conta 
dos NA’S.


Obrigado pela ajuda.

Enviado do meu telefone Windows 10

*De: *salah via R-br <mailto:r-br@listas.c3sl.ufpr.br>
*Enviado:*quinta-feira, 16 de junho de 2016 23:23
*Para: *r-br@listas.c3sl.ufpr.br <mailto:r-br@listas.c3sl.ufpr.br>
*Assunto: *Re: [R-br] Operação entre datas

segue sugestão

library(eeptools)

DADOS = data.frame(
NASCIMENTO = c("26/11/97", "20/03/99", "14/05/99","06/05/98", 
"03/01/00", "19/05/97", "01/02/01", "28/11/97", "10/02/00"),
DATA_FATO = c("23/02/15", "28/03/15", "08/04/15", "08/04/15", 
"08/04/15", "08/05/15", "17/05/15", "03/06/15", "03/06/15"))


## converter para Date
DADOS$NASCIMENTO = as.Date(DADOS$NASCIMENTO, "%d/%m/%y")
DADOS$DATA_FATO = as.Date(DADOS$DATA_FATO, "%d/%m/%y")

## calcula o intervalo em anos
DADOS$IDADE = floor(age_calc(DADOS$NASCIMENTO, DADOS$DATA_FATO, units 
= "years"))


DADOS

saudações

Em 16/06/2016 14:45, Amikobh via R-br escreveu:

Prezados,

Desde já, obrigado pela ajuda.

Tendo duas variáveis com datas referentes a determinada pessoa:

NASCIMENTO



DATA_FATO

26/11/97



23/02/15

20/03/99



28/03/15

14/05/99



08/04/15

06/05/98



08/04/15

03/01/00



08/04/15

19/05/97



08/05/15

01/02/01



17/05/15

28/11/97



03/06/15

10/02/00



03/06/15

Como procedo para criar mais uma coluna que calcula a idade no dia
da ocorrência de determinado fato? De tal forma que obtenha:

NASCIMENTO



DATA_FATO



IDADE

26/11/97



23/02/15



17

20/03/99



28/03/15



16

14/05/99



08/04/15



15

06/05/98



08/04/15



16

03/01/00



08/04/15



15

19/05/97



08/05/15



17

01/02/01



17/05/15



14

28/11/97



03/06/15



17

10/02/00



03/06/15



15

Valeu, gente!




___

R-br mailing list

R-br@listas.c3sl.ufpr.br <mailto:R-br@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@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@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.

Re: [R-br] Operação entre datas

[salah via R-br]

segue sugestão

library(eeptools)

DADOS = data.frame(
NASCIMENTO = c("26/11/97", "20/03/99", "14/05/99","06/05/98", 
"03/01/00", "19/05/97", "01/02/01", "28/11/97", "10/02/00"),
DATA_FATO = c("23/02/15", "28/03/15", "08/04/15", "08/04/15", 
"08/04/15", "08/05/15", "17/05/15", "03/06/15", "03/06/15"))


## converter para Date
DADOS$NASCIMENTO = as.Date(DADOS$NASCIMENTO, "%d/%m/%y")
DADOS$DATA_FATO = as.Date(DADOS$DATA_FATO, "%d/%m/%y")

## calcula o intervalo em anos
DADOS$IDADE = floor(age_calc(DADOS$NASCIMENTO, DADOS$DATA_FATO, units = 
"years"))


DADOS

saudações

Em 16/06/2016 14:45, Amikobh via R-br escreveu:


Prezados,

Desde já, obrigado pela ajuda.

Tendo duas variáveis com datas referentes a determinada pessoa:

NASCIMENTO
DATA_FATO
26/11/9723/02/15
20/03/9928/03/15
14/05/9908/04/15
06/05/9808/04/15
03/01/0008/04/15
19/05/9708/05/15
01/02/0117/05/15
28/11/9703/06/15
10/02/0003/06/15

Como procedo para criar mais uma coluna que calcula a idade no dia da 
ocorrência de determinado fato? De tal forma que obtenha:


NASCIMENTO  DATA_FATO   IDADE
26/11/9723/02/1517
20/03/9928/03/1516
14/05/9908/04/1515
06/05/9808/04/1516
03/01/0008/04/1515
19/05/9708/05/1517
01/02/0117/05/1514
28/11/9703/06/1517
10/02/0003/06/1515


Valeu, gente!








___
R-br mailing list
R-br@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@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.

Re: [R-br] Windows Server 2008 (64-bit)

[salah via R-br]

Olá Daniel

Sugestão:

Instale o Windows Server 2008 (64-bit) numa máquina virtual

tenho apreço por essa:
https://www.virtualbox.org/

windows:
https://www.microsoft.com/en-us/download/details.aspx?id=23163

saudações

Em 09/06/2016 15:17, Daniel Marcelino via R-br escreveu:

Windows Server 2008 (64-bit)


___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Leitura de Arquivo

[salah]

sugestão:

 dataset = 
read.csv("C:\\Users\\Anselmo\\Desktop\\trab\\BC-pOct89.TL-1", sep="", 
header = FALSE)


saudações

Em 02/06/2016 00:15, Anselmo Alves de Sousa escreveu:

Boa noite listeiros

Estou com a missão de auxiliar um colega no trabalho escolar. A 
primeira missão é ler o dado. Como begginer no R não estou conseguindo 
ler o arquivo ftp://ita.ee.lbl.gov/traces/BC-pOct89.TL.Z. Sabemos que 
tem duas colunas, uma com o timestamp e outra com o tamanho em bytes.


Usando
 dataset<-scan("C:\\Users\\Anselmo\\Desktop\\trab\\BC-pOct89.TL\\BC-pOct89.TL")

no diretório apropriado

setwd("C://Users//Anselmo//Desktop//trab_vinicius//BC-pOct89.TL")
e str(dataset)
 num [1:200] 0.0177 87 0.0368 142 0.0368 ...

Ou seja, obtenho apenas as duas colunas numa só. Como posso ler 
adequadamente esse arquivo?


[]'s
Anselmo


___
R-br mailing list
R-br@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@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.

Re: [R-br] Erro ao fazer download de PDFs: failed, no valid url links detected

[salah]

Caro Elias

Você pode baixar o pacote fonte metagear (a beleza do opensource) aqui:

https://cran.r-project.org/src/contrib/Archive/metagear/metagear_0.1.tar.gz

Examinando a função PDF_download() notamos que ela usa três arquivos:
PDF_download.R, PDF_utils.R, PDF_collect.R e isPDF.R

para o doi "10.1371/journal.pone.0123511" a função está extraindo como alvo:

"content=\"http://journals.plos.org/plosone/article/asset?id=10.1371/journal.pone.0123511.PDF;

e o correto seria o link com atributo href:

"href=\"/plosone/article/asset?id=10.1371%2Fjournal.pone.0123511.PDF"

resolvendo em: 
"http://journals.plos.org/plosone/article/asset?id=10.1371%2Fjournal.pone.0123511.PDF;


Segue um pequeníssimo ROBOT, veja que ele não prevê https e necessita 
das urls dos jornais alvo.


# Author: salah
###
require(RCurl)
require(XML)

rm(list=ls())

## seleciona o link alvo e converte em caracter
baixaURL = function(doi)
{
## link doi
urlDOI = paste0("http://dx.doi.org/;, doi)

## download url
web = getURLContent(urlDOI)
tc = textConnection(web)
web = readLines(tc)
close(tc)

## captura os links
lnk = getHTMLLinks(web)

## converte para chr
doc = htmlParse(lnk)
w = capture.output(doc)

return(w)
}##end baixaURL

## prepara a url para download
urlPDF = function(w, journal)
{
a = unlist(strsplit(w, split = "\\s"))
b = unique(grep("pdf|PDF", a, value = TRUE))
cc = grep("href+", b, value = TRUE)
d = gsub("href\\=|\"|>", "", cc)
h = gsub(".*pdf|.*PDF", "", d)
outPDF = gsub(h[1], "", d)

return(paste0("http://;, journal, outPDF))
}##end naturePDF

## faz o download do pdf
baixaPDF = function(urls, dest = "~")
{
j = unlist(strsplit(urls, "/"))
namePDF = grep(".pdf|.PDF", j, value = T)

## se houver pdf extra
for(i in 1:length(urls))
{
dir_name = paste0(dest, namePDF[i])
download.file(urls[i], dir_name)
}##end for
}##end baixaPDF

doiNature = c("10.1038/nutd.2016.11", "10.1038/srep17841", 
"10.1038/srep25762")

doiPLOS = c("10.1371/journal.pone.0123511", "10.1371/journal.pbio.1002461")

j1 = "www.nature.com/"
j2 = "journals.plos.org"

DATA = data.frame(DOI = c(doiNature, doiPLOS), Journal = c(j1, j1, j1, 
j2, j2), stringsAsFactors = FALSE)


destDir = "~/"

for(n in 1:nrow(DATA))
{
w = baixaURL(DATA$DOI[n])
urls = urlPDF(w, DATA$Journal[n])
print(urls)
baixaPDF(urls, destDir) ## faz o download
}


saudações

Em 22/05/2016 23:23, Elias Carvalho escreveu:
Desculpe Leonardo, segue o código mínimo, tem dois DOIs para poder 
testar:


*library(metagear)
*
*
*
*DOI <- "10.1371/journal.pone.0123511" # Disponível em: 
http://goo.gl/rhtvFx*
*# DOI <- "10.1038/srep17841"# Disponível em: 
http://goo.gl/rNwOpk*

*filename <- "ArtigoTeste.pdf"
*
*folder <- "~"
*
*PDF_download(DOI, directory = folder, theFileName = filename, 
validatePDF = TRUE, quiet = FALSE)*



Em 22 de maio de 2016 09:41, Elias Carvalho <ecaca...@gmail.com 
<mailto:ecaca...@gmail.com>> escreveu:


Bom dia Pessoal

Estou baixando tentando baixar alguns pdf de artigos científicos
pelo DOI usando o pacote metagear, porém a maioria apresenta a
mensagem "failed, no valid url links detected", ou
"cannot open: HTTP status was '404 Not Found' PDF download...
skipped", mesmo ligando o proxy da universidade.

De 1800 artigos que preciso consegui por volta de uns 200 download
apenas.

Alguém poderia me dar uma ajuda em como resolver este problema ?

-- 
Best regards... 8^)


“The mind that is open to new ideas never come back
to its original size” /Albert Einstein/


-- 
Obrigado

Elias




--
Best regards... 8^)

“The mind that is open to new ideas never come back
to its original size” /Albert Einstein/


--
Obrigado
Elias


___
R-br mailing list
R-br@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@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.

Re: [R-br] Contar casos em cada variável e em todas as variáveis

[salah]

Caro

Segue sugestões

library(reshape2)
ATO$id = 1:7 ##acrescenta uma coluna de controle para o melt
##remodelando sua tabela
ATOMelt = melt(ATO, id="id", variable.name="ATO_n", value.name="ART")

##usando dplyr
library(dplyr)
ATOMelt %>% count(ART)

##usando o poderoso data.table
library(data.table)
setDT(ATOMelt) ##converte para data.table
##contando
ATOMelt[, .N, by = .(ART)]

saudações

Em 26/05/2016 16:28, Amiko Bh escreveu:

Pessoal,

Dado o banco:

*ATO_1*



*ATO_2*



*ATO_3*

ART 121



ART 155



ART 180

ART 33



ART 28



ART 10

ART 121



ART 15



ART 155

ART 155



ART 33



ART 121

ART 28



ART 33



ART 33

ART 121



ART 121



ART 121

ART 15



ART 121



ART 180



a) Preciso contar cada elemento em cada variável, para que vislumbre o 
seguinte resultado:


*ATO_1*

ART 121 = 3
ART 33 = 1
ART 155= 1
ART 28= 1
ART 15 = 1

e assim com *ATO_2 e **ATO_3.*

b) preciso saber o total de cada elemento em todos as colunas, assim:

ART 121 = 7
ART 33 = 4
ART 155= 3
ART 28= 2
ART 15 = 2
ART 180=2
ART 10=1

Muito obrigado.



___
R-br mailing list
R-br@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@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.

Re: [R-br] Multiplots usando lattice

[salah]

Olá Marcos

Talvez este exemplo lhe ajude

library(lattice)

# Data
w <- as.matrix(dist(Loblolly))
x <- as.matrix(dist(HairEyeColor))
y <- as.matrix(dist(rock))
z <- as.matrix(dist(women))

# Plot assignments
pw <- levelplot(w, scales = list(draw = FALSE))  # "scales..." removes axes
px <- levelplot(x, scales = list(draw = FALSE))
py <- levelplot(y, scales = list(draw = FALSE))
pz <- levelplot(z, scales = list(draw = FALSE))

# Plot prints
print(pw, split = c(1, 1, 2, 2), more = TRUE)
print(px, split = c(2, 1, 2, 2), more = TRUE)
print(py, split = c(1, 2, 2, 2), more = TRUE)
print(pz, split = c(2, 2, 2, 2), more = FALSE)  # more = FALSE is redundant

retirado de 
http://stackoverflow.com/questions/2540129/lattice-multiple-plots-in-one-window


saudações

Em 13/05/2016 00:14, marcos paulo escreveu:

Pessoal estou querendo plotar 4 gráficos na mesma janela gráfica,
usando o pacote lattice pra confecção desses plots. Porém não estou 
conseguindo colocar mais que dois na mesma janela, usando a função print.


Segue exemplo:

g1<- xyplot(NDVI~DAE|Cultivar, data=dados, col=1, layout=c(5,1),
   type=c("p","r"),
   xlab="DAE", ylab="NDVI",
   main='Vitrine 2014/2016',
   strip=strip.custom(bg="gray90"))

g2<-xyplot(NDVI~IAF|Cultivar, data=dados, col=1, layout=c(5,1),
   type=c("p","r"),
   xlab="IAF", ylab="NDVI",
   strip=strip.custom(bg="gray90"))


g3<-xyplot(Chl~DAE|Cultivar, data=dados, col=1, layout=c(5,1),
   type=c("p","r"),
   xlab="DAE", ylab="Clorofila",
 strip=strip.custom(bg="gray90"))

g4<-xyplot(Chl~IAF|Cultivar, data=dados, col=1, layout=c(5,1),
   type=c("p","r"),
   xlab="IAF", ylab="Clorofila",
   strip=strip.custom(bg="gray90"))
print(g1, split=c(1,1,1,4), more=TRUE)
print(g2, split=c(1,2,1,4))
print(g3, split=c(1,3,1,4))
print(g4, split=c(1,4,1,4))

Cordialmente, Marcos Paulo.

*Técnico Agrícola em Zootecnia
Bacharel em Agronomia *
*
*
*Mestrando em Produção Vegetal - Universidade Federal de Goiás*
*Estagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo 
Antônio de Goiás-GO*

*
*
*Contato: (62) 85075783*

http://lattes.cnpq.br/4322347592884852

 
	Livre de vírus. www.avast.com 
. 





___
R-br mailing list
R-br@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@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.

Re: [R-br] Pacote XTable - Erro de saida

[salah]

Olá Elias

No RStudio abra um novo arquivo R Markdown, marque a opção HTML. Abrirá 
um arquivo com exemplos, é só acrescentar:


```{r results='asis', echo=FALSE}
  library(xtable)
  data(tli)
  print(xtable(tli[1:20, ]), type='html')
```
é só rodar clicando no ícone "Knit HTML"

Não há nenhum erro, o xtable lhe retorna o código para usar em Latex ou HTML

veja ajuda: http://rmarkdown.rstudio.com/

saudações

Em 08/05/2016 19:21, Cesar Rabak escreveu:

Elias,

Tecnicamente essa questão deixa de ser sobre o "R" e passa a ser sobre 
o "RStudio". . .


Você tem essa expectativa de ver uma tabela em LaTeX sem compilá-la? O 
RStudio tem esse tipo de funcionalidade?



2016-05-08 13:28 GMT-03:00 Elias Carvalho >:


Obrigado pela resposta Felipe, mas o RStudio não deveria abrir
automaticamente o leitor html ao apresentar os dados ?

Em 8 de maio de 2016 12:55, Elias Carvalho > escreveu:

Boa Tarde Pessoal

Estou testando o pacote "xtable, e o meu ambiente é ubuntu
14.04, RStudio 0.99.896, porém estou tendo o seguinte problema:

Quando executo o código abaixo:

*
# Load library
library(xtable)

## Load example dataset
data(tli)

*
## Demonstrate data.frame
*
tli.table <- xtable(tli[1:20, ])
xtable(mtcars)
*

O meu resultado não é uma tabela apresentável:

% latex table generated in R 3.3.0 by xtable 1.8-2 package
% Sun May  8 12:53:18 2016
\begin{table}[ht]
\centering
\begin{tabular}{}
   \hline
  & mpg & cyl & disp & hp & drat & wt & qsec & vs & am & gear & carb \\
   \hline
Mazda RX4 & 21.00 & 6.00 & 160.00 & 110.00 & 3.90 & 2.62 & 16.46 & 0.00 & 1.00 
& 4.00 & 4.00 \\


Já busquei algumas soluções no google mas nada muito explicativo.

Alguém pode me ajudar a solucionar o problema ?

-- 
Obrigado

Elias




-- 
Best regards... 8^)


“The mind that is open to new ideas never come back
to its original size” /Albert Einstein/


-- 
Obrigado

Elias

___
R-br mailing list
R-br@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@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@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.

Re: [R-br] Completar linhas de uma variável

[salah]

Olá

segue:

## vetor tipo character
PROCESSO = 
rep(c("5.155.855-8","5.153.696-8","5.154.220-6",NA,NA,NA,"5.155.401-1","5.155.960-6",NA,"5.157.157-7","5.157.658-4"), 
1000)


um pequeno ajuste

## evita que transforme em factor
dados = data.frame(PROCESSO, stringsAsFactors = FALSE)
indice = which(is.na(dados$PROCESSO))
##Agora é so usar o indice para preencher a coluna
dados[indice,"PROCESSO"] <- 1:length(indice)
dados

Você também pode usar o poderoso data.table

library(data.table)
## converte para data.table
dados = data.table(PROCESSO)
dados[is.na(PROCESSO), PROCESSO := as.character(1: .N)]
dados

saudações

Em 27/04/2016 02:00, Amikobh escreveu:

indice<-which(is.na (dados$PROCESSO))

##Agora é so usar o indice para preencher a coluna
dados[indice,"PROCESSO"]<-1:length(indice)


___
R-br mailing list
R-br@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.

Re: [R-br] Casas Decimais

[salah]

Olá

Como os colegas já colocaram a conversão do número binário para um 
número real não é trivial.
Se você pretende trabalhar com alta precisão numérica  você precisa 
estudar o deus da informática - ponto flutuante - (/floating point)./


sugiro a leitura:

https://en.wikipedia.org/wiki/Floating_point

https://pt.wikipedia.org/wiki/Ponto_flutuante

http://www.burns-stat.com/pages/Tutor/R_inferno.pdf

O R sempre trabalha com o 22 dígitos, a função options(digits = ) 
controla o número de dígitos que será visualizado.

A precisão está em 16 dígitos

Verifique:

options(digits = 22)
#options(digits = 16)

(x1 = 10^-324)
(x2 = 10^-323)
(y1 = 10^+308)
(y2 = 10^+309)

(n17 = 0.1)
(n16 = 1.0001)
(n15 = 1.001)
(n14 = 1.01)
(n13 = 1.1)
(n12 = 1.0001)
(n11 = 1.001)
(n10 = 1.01)
(n09 = 1.1)
(n08 = 1.0001)
(n07 = 1.001)
(n06 = 1.01)
(n05 = 1.1)
(n04 = 1.0001)
(n03 = 1.001)
(n02 = 1.01)
(n01 = 1.1)
(n00 = 1)

n16 + n17
n15 + n16
n14 + n15

(n14 + n15) == 2.01110223
(n14 + n15) == 2.0111
(n14 + n15) == 2.011
(n14 + n15) == 2.01

0 == n17
1 == n16
1 == n15

.1
.2
.1 + .2
.3
.1 + .2 == .3

bons estudos!

Em 22/04/2016 08:19, Rodrigo Coster escreveu:
Contribuindo com a resposta do Jobenil, isso acontece porque a 
representação do teu número em binário é uma dizima periódica (agora 
não tenho certeza se é só dizima ou dizima periódica) - isso é: tem 
infinitas casas decimais. Por não ter como armazenar as infinitas 
casas, ele é arredondado e por isso esse 'resíduo'. Embora parece algo 
que possa ser ignorado, devemos nos preocupar com esses casos, como 
mostra o exemplo:


> .1 + .2
[1] 0.3
> .1 + .2 == .3
[1] FALSE



2016-04-22 8:11 GMT-03:00 Jobenil Luiz Magalhães Júnior 
>:


Bom dia Júlio,
Acho que é devido ao processo de conversão interna de decimal para
binário e vice versa. Como você exige que a representação dos
números tenham 22 dígitos na reconversão de P para decimal ele
retorna essa precisão.
Teste outro número por exemplo 665.5, penso que não encontrará
este problema pois sua representava binária será exata.

Atenciosamente
Jobenil Junior


Enviado do meu iPad

> Em 21 de abr de 2016, às 19:35, julio cesar oliveira
> escreveu:
>
> Colegas,
>
> Sou novato no R, e por isso gostaria da ajuda para esclarecer
algumas dúvidas.
>
> Preciso que os cálculos sejam realizados com máxima precisão em
relação ao número de casas decimais. Por isso tenho a seguinte dúvida:
>
> 1) Como configurar as casas decimais no R ?
>
> estou usando o comandooptions(digits=22) para apresentar os
dados com 22 casas. Além dele posso trabalhar com o pacote (Rmpfr)
>
> exemplo
>
> > mpfr(665.598,128)
> 1 'mpfr' number of precision  128   bits
> [1] 665.59799563442543148994445800781
>
> 2) por que aparece um resíduo? por exemplo
>
>  p<- 665.598
> > p
> [1] 665.59799563443
>
> Grato
>
> julio
> ___
> R-br mailing list
> R-br@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@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
fornea cdigo mnimo reproduzvel.




___
R-br mailing list
R-br@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@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.

Re: [R-br] Separação de dados em classes usando if e ifelse

[salah]

Olá Alexandre

Tente:

dados$classes <- ifelse(dados$V2 %in% GROUP1,"Group1" ,"Group2")

? match

Sugiro você usar o pacote data.table ou dplyr

V1 = c(1,1,1,1,1,12,12,12,12,12,31,31,31,31,31,31,31,142,142,142)
V2 = paste('P',  V1, sep = "")
GROUP1 = c("P12B","P12C","P12D","P12E","P12F", 
"P12G","P12H","P12I","P13A","P13B","p142", 
"P142","p148","P148","P15B","P15C","P15D", 
"P15E","P15F","P15G","P15H","P15I","P15J", 
"P16A","P16B","P16C","P16D","P16E","P16F", 
"P16G","P16H","P16I","P16J","P16K","P16", "P1A",
"P1J","P1K","P1L","P21A","P27A", "P27B","P27C",
"P27D","P27E","P2A","P2B", "P2D","P2E","P2F",
"P2G","P2","P31A", "P31B", "P34","P37","P3A",
"P3B","P3C", "P3D","P3E", "P3F","P492","P6A","P6B", 
"P6C","P7A","P7B", "P7C","P7D","P7E", "P7F",
"P7G","P7H","P12A")


library(data.table)

dados = data.table(V1, V2)
dados[V2 %in% GROUP1, classes := "Group1"]
dados[ ! V2 %in% GROUP1, classes := "Group2"]
dados

saudações

Em 22/03/2016 16:40, ASANTOS escreveu:

Caros Listeiros,

   Estou tentando empregar as dicas do Rafael para achar uma 
maneira de classificar o nome de variáveis, mas não individualmente 
com os níveis e sim com um vetor que indica quais são os níveis de 
interesse que quero classificar. É possível fazer isso com um vetor, 
no meu caso  GROUP1, ao invés de utilizar o nome de cada nível 
individualmente seguido de | na função ifelse?


dados<-NULL
dados$V1<-c(1,1,1,1,1,12,12,12,12,12,31,31,31,31,31,31,31,142,142,142)
dados$V2 <- paste('P', dados$V1,sep="")
GROUP1<-c("P12B","P12C","P12D","P12E","P12F", 
"P12G","P12H","P12I","P13A","P13B","p142", 
"P142","p148","P148","P15B","P15C","P15D", 
"P15E","P15F","P15G","P15H","P15I","P15J", 
"P16A","P16B","P16C","P16D","P16E","P16F", 
"P16G","P16H","P16I","P16J","P16K","P16", 
"P1A","P1J","P1K","P1L","P21A","P27A", "P27B",
"P27C","P27D","P27E","P2A","P2B", "P2D","P2E",
"P2F","P2G","P2","P31A", "P31B","P34","P37",
"P3A","P3B","P3C", "P3D","P3E","P3F","P492",
"P6A","P6B", "P6C","P7A","P7B","P7C","P7D",
"P7E", "P7F", "P7G","P7H","P12A")

dados$classes<- ifelse((dados$V2)==GROUP1,"Group1" ,"Group2")

Obrigado,



___
R-br mailing list
R-br@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.

Re: [R-br] Aumentar o tamanho de fontes em gráfico ggplot2

[salah]

Olá Marcelo

A função opts foi substituída por theme
O título pode ser acrescentado com ggtitle("Título do gráfico")
A data da publicação do r-bloger é de 2012, de lá prá cá o ggplot2 
sofreu muitas alterações(melhorias)

A função geo_text() é para inserir texto no gráfico
 veja aqui:
http://www.cookbook-r.com/Graphs/Fonts/

Você pode encontrar mais informações aqui:

http://docs.ggplot2.org/dev/index.html

http://www.cookbook-r.com/Graphs/Fonts/#themes-and-elementtext

segue código atualizado:

p1 = ggplot(data=Data,
aes(x=factor(year), y=proportion.tasty,
group=group,
shape=group,
color=group)) +
geom_line() +
ggtitle("Proportion Tasty by Year, Quality, and Group") +
theme(axis.text.x = element_text(angle = 90, vjust = .5)) + ## 
altera o angulo do texto - eixo x

geom_point() +
scale_x_discrete("Year") +
scale_y_continuous("Proportion Tasty") +
facet_grid(.~quality )
p1
## mudando o tema globalmente
theme_set(theme_gray(base_size = 18))
p1
## ou
theme_set(theme_bw(base_size = 18))
p1

saudações

Em 05/03/2016 17:11, Marcelo Laia escreveu:

Olá,

Rodei o exemplo de gráfico disponível em
http://www.r-bloggers.com/how-to-plot-three-categorical-variables-and-one-continuous-variable-using-ggplot2/

Tive dois problemas:

1. não consegui incluir o Título do gráfico. É informado que a função opts não
existe.

2. as fontes ficaram muito pequenas. Se eu insiro a função geom_text(size = 12,
aes()) dá erro. Alguém saberia em que local do código eu posso alterar
as fontes globalmente (de todo o gráfico)?

Obrigado.



___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Função para remover colunas

[salah]

Olá

O colega já lhe deu a solução, mas se é pra facilitar a vida do usuário 
leigo segue sugestão


(df = data.frame( a=letters[10], b=LETTERS[10], c=seq(10), 
d=letters[1:10], e=letters[1:10] ))


##  exemplo1 
remove.var = function(dataname, vars)
{
dataname[vars] <- list(NULL)
return(dataname)
}##end function

vars = c("a", "b")
remove.var(df, vars)

##  exemplo2 
remove.var = function(dataname, vars)
{
return(dataname[, ! names(dataname) %in% vars, drop = FALSE])
}##end function

vars = c("a", "b")
remove.var(df, vars)

saudações

Em 28/02/2016 08:55, Elias Carvalho escreveu:

Ola Pessoal

Para facilitar a vida de um usuário leigo, eu criei uma função para 
excluir as colunas de uma tabela, conforme código abaixo:


 1. remove.variables <- function(dataset.name ,
vars.to.remove)
 2. {
 3.   # scan the variables to be removed
 4.   for (x in 1 : length(vars.to.remove))
 5.   {
 6. # transform command.Remove in "data$variable <- NULL"
 7. command.Remove <- paste(dataset.name
,"$",vars.to.remove[x], " <- NULL",sep="")
 8. # execute the command
 9. command.Remove <- eval(parse(text=command.Remove))
10.   } #  for (x in 1 : length(vars.to.remove))
11. } # remove.variables <- function(dataset.name
, vars.to.remove)


Fazendo um teste local, ou seja, executando o procedimento abaixo com 
as variáveisdataset.name , vars.to.remove 
carregadas respectivamente com o nome do meu data frame "data" e com o 
nome de duas variáveis "var1" e "var5" a serem removidas e executando 
da linha 2 a 11 tudo tunciona bem e as colunas são removidas.


dataset.name  <- "data"
vars.to.remove<-c("var1", "var5")

No entanto se carrego a função e executo os procedimentos abaixo em 
outro script, a função é executada sem erro, mas não remove nenhuma 
coluna do data frame:


vars.to.remove = c("var1", "var5")
dataset.name  = "data"

remove.variables(dataset.name , vars.to.remove)

Eu acredito que tem a ver com o ambiente, acesso global algo assim mas 
não estou achando a solução.


Agradeço se aguem puder me ajudar.

--
Best regards... 8^)

“The mind that is open to new ideas never come back
to its original size” /Albert Einstein/


_
Prof. Elias César Araújo de Carvalho
CV: http://lattes.cnpq.br/4248328961021251


___
R-br mailing list
R-br@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@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.

Re: [R-br] Cores em mapas

[salah]

Olá Mateus

segue sugestão:

## paleta de cores
cor1 = heat.colors(15)
cor2 = terrain.colors(15)
cor3 = topo.colors(15, alpha = .5)
cor4 = palette(rainbow(9))
library(RColorBrewer)
cor5 = brewer.pal(7, "BrBG")

contour(
sort(lon),
lat,
levels = intervalos,
nlevels = 15,
media_ColumnAmountO3[ order(lon), ],
add = T,
lwd = 2,
col = cor1,
labcex = 1.3)

## usando o pacote raster
library(raster)
oz = raster('teste.nc')
xlim = c(-100, -10)
ylim = c(-60, 10)
plot(oz, xlim = xlim, ylim = ylim)
map(xlim = xlim, ylim = ylim, add = TRUE, col = "black")
title(main = "Campo médio de Ozônio Novembros" )
contour(oz, add = TRUE, xlim = xlim, ylim = ylim)

saudações

Em 19/02/2016 17:18, Mateus Dias Nunes escreveu:
Olá eu gostaria de colocar cores (com a barra de cores ao lado da 
figura) no meu mapa ao invés de usar os contornos com a função "contour";

abaixo o link da figura que consegui gerar através dessa função.

https://www.dropbox.com/s/xi622mbpqg7tgch/campo_medio_O3.png?dl=0

abaixo o script que gerou o grafico com os contornos.



# Carregando biblioteca para manipular arquivos netCDF

library(maps)
library(ncdf4)
#==


#BIBLIOTECA "ncdf4".
# PARA ESTE EXEMPLO SÃO USADOS RECURSOS PARA ABRIR A BIBLIOTECA )
#POIS RNetCDF, ncdf e ncdf4 APRESENTAM COMANDOS DIFERENTES PARA 
ABRIRMOS AS VARIÁVEIS



 dados <- nc_open('teste.nc ')
# lendo coordenadas espaço-temporal
lat <- ncvar_get( dados, 'lat' )
lon <- ncvar_get( dados, 'lon' )
time <- ncvar_get( dados, 'time' )

#===

# lendo dados coluna total de Ozônio
ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' )

# dimensoes da variavel ColumnAmountO3
dims_ColumnAmountO3 <- dim(ColumnAmountO3)

# tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos 
de grade X tempo


dim(ColumnAmountO3) <- c( 
dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2], dims_ColumnAmountO3[3] )


# calculando a média e retornado-a em 2D

media_ColumnAmountO3 <- rowMeans( ColumnAmountO3)
dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1], 
dims_ColumnAmountO3[2] )


#==
# longitude varia de 0 a 360, convertendo para -180 a 180, essa 
conversão é feita para plotagem sobre o mapa

for (i in 1:dim(lon)) { if (lon[i]>180) { lon[i] <- lon[i]-360 } }


# criando arquivo PNG que receberá o campo com o mapa

#png( filename="campo_medio_O3_jan2005.png",width=600,height=800 )

# plotando mapa da America do Sul
map( xlim=c(-100,-10), ylim=c(-60,10) )
map.axes()  # plotando eixos
title( main="Campo médio de ozonio janeiros" )   # título do gráfico

# definindo intervalo de 5 Dobson Units (DU)
intervalos = seq( trunc(min(ColumnAmountO3)), 
trunc(max(ColumnAmountO3)), 5 )


# adicionando campo de coluna de ozonio


contour( sort(lon), lat, media_ColumnAmountO3[ order(lon), ], add=T, 
levels=intervalos, lwd=2, labcex=1.3, col="black" )



# fechando arquivo PNG
#dev.off()


obrigado


MATEUS DIAS NUNES
MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET
UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL
TELEFONE: +55 (53) 81125154



___
R-br mailing list
R-br@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@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.

Re: [R-br] Gráfico Pareto

[salah]

Caro,

segue sugestão:

## dados
dd = data.frame(valor = 2*abs(f1), grupo = names(f1))

## equação
lb1 = paste("S[ef]%.%t[crit]==", s)

library(ggplot2)

## gráfico
p = ggplot(dd, aes(grupo, valor, fill=grupo)) +
theme_bw() +
geom_bar(stat = "identity") +
xlab(NULL) +
ylab(NULL) +
ggtitle("Gráfico de pareto dos efeitos") +
geom_text(aes(label = round(valor, 1)),
position = position_dodge(width = 0.9),
angle = -90,
vjust = -0.25) +
annotate("text", x = 2, y = 10, label = lb1, parse=TRUE) +
coord_flip()

p

ggsave("grafico_pareto.jpg", p, width = 6, height = 4)

saudações

Em 18/02/2016 13:20, Ari Clecius escreveu:
Prezados colegas, este assunto me interessa, gostaria de saber como 
automatizar este gráfico, fiz um script, mas é todo manual, os 
gráficos de pareto dos pacotes não formatam por cor, seria possível 
visualizar os efeitos negativos e positivos em cores diferentes?


Att,
Ari Clecius Alves de Lima
Engenheiro Químico
Me. Engenharia Civil
Dr. Engenharia Civil
(085)84174333
(085)33669042

*/#3/*
*/#A=temperatura/*
*/#B=vibração/*
*/#A<-c(20,20,60,60,20,20,60,60,20,20,60,60)/*
*/#B<-c(5,15,5,15,5,15,5,15,5,15,5,15)/*
*/#y<-c(74,50,46,34,68,48,44,32,65,51,40,26)/*
*/
/*
*/A<-c(-1,-1,1,1,-1,-1,1,1,-1,-1,1,1)/*
*/B<-c(-1,1,-1,1,-1,1,-1,1,-1,1,-1,1)/*
*/y<-c(74,50,46,34,68,48,44,32,65,51,40,26)/*
*/group<-c("A","B","C","D","A","B","C","D","A","B","C","D")/*
*/falha<-data.frame(A,B,y,group)/*
*/aov.2=aov(y~A+B+A*B, data=falha)/*
*/summary(aov.2)/*
*/
/*
*/
/*
*/efeitos<-2*(aov.2$coef[2:4])/*
*/
/*
*/#Estimativa da variância do erro experimental (s^2)/*
*/grau_L1<-2/*
*/grau_L2<-2/*
*/grau_L3<-2/*
*/grau_L4<-2/*
*/ag <- data.frame(aggregate(. ~ group, falha, function(x) c(mean = 
mean(x), sd = sd(x/*

*/s2<-(sum(2*(ag$y[,2])^2))/(grau_L1+grau_L2+grau_L3+grau_L4)/*
*/
/*
*/#syi^2=s2/n onde n é o número de replicatas(3)/*
*/#sy1^2=sy2^2=sy3^2=sy4^2=syq/*
*/#sy=raiz(syq)/*
*/n<-3/*
*/sy=sqrt(s2/n)/*
*/
/*
*/#t(ef)=ef/s(ef)/*
*/tefA=efeitos[1]/sy/*
*/tefB=efeitos[2]/sy/*
*/tefAB=efeitos[3]/sy/*
*/
/*
*/
/*
*/#Gráfico de pareto/*
*/
/*
*/f1<-c(tefA,tefB,tefAB);f1/*
*/z=data.frame(2*abs(f1))/*
*/z=t(z)/*
*/v1=sy*(qt(0.975,8))/*
*/coef3=sort(z)/*
*/coef3/*
*/
/*
*/#png(file = "pareto_ceto.png",res=300, width = 4, height = 4, units 
= 'in')/*
*/g=barplot(sort(z),horiz=TRUE,names.arg=c( 
rev(names(f1))),col=ifelse(sort(abs(f1))>v1,"blue","gray"),xlim=c(0,25))##xlim, 
ylim/*

*/box("plot", col="blue") /*
*/abline(v=v1,col="red",lty=2)/*
*/s=format(v1,digits=2)/*
*/text(6,1,label=expression("S"[ef]*" x t"[crit]*"="),cex=1)/*
*/text(9,1,label=s,cex=1)/*
*/mtext("Gráfico de pareto dos efeitos", side=3, line=2, cex=1.0, 
col="blue", outer=FALSE)/*

*/p=format(f1,scientific = TRUE)/*
*/legend(x=16,y=2,cex=0.6,bg="yellow",box.col="blue", legend=c( 
"A=-22.3","B=-16.0","AB=3.33" ), col="blue", lwd=1:2, lty=c(NA,NA,NA))/*

*/#dev.off()/*


Em Quinta-feira, 18 de Fevereiro de 2016 10:18, Paulo Dick 
 escreveu:



Bom dia Ana Paula,

Seu gráfico não aparece para mim.

Sugiro também que envie a sintaxe que utilizou.

Abraço

*
*
*Paulo Dick*
Estatístico
Mestre em Epidemiologia em Saúde Pública
Tel.: (55 21) 99591-2716

Em 17 de fevereiro de 2016 21:59, ana paula coelho madeira 
> escreveu:


Prezados,

boa noite!

Construí um gráfico de Pareto utilizando a função paretoPlot e a
saída é da forma:

Gostaria de fazer uma das seguintes alterações:

1) Apresentar as barras na horizontal;
ou
2) Colocar as estimativas na vertical ( para não ficar tão difícil
de visualizar).

Desde já agradeço.

Att.,

Ana Paula

___
R-br mailing list
R-br@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@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� 
c�igo m�imo reproduz�el.





___
R-br mailing list
R-br@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@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.

Re: [R-br] Importação de planilha de contingência

[salah]

Olá Marie

A exigência de valores únicos para row.names é inerente a formação do 
data.frame.

veja em:
?row.names
?data.frame

Provavelmente há algum nome duplicado.
Aplique e o script abaixo para ter uma ideia

## arquivo no formato CSV
## read.csv2 vem como default sep=";" e dec=","
## carregue os ddos sem indicar row.names
dados = read.csv2("C:\\endereço\\nomeArquivo.csv")

## verifica duplicatas
dados[duplicated(dados),]

## modifica os nomes para valores únicos
(rn = make.names(dados$X, unique = TRUE))

## visualiza row.names - linhas numeradas
rownames(dados)
dados

## nova row.names
rownames(dados) = rn
dados

## remove a coluna X
dados = dados[-1]
dados

## verifique
?duplicated
?make.names

saudações

Em 16/02/2016 10:43, Marie-Christine . escreveu:

Olá Salah,


Muitíssimo obrigada pela ajuda. Dessa vez deu certo e as análises 
rodaram tranquilamente ;).


Saberias, contudo, me dizer o motivo pelo qual o R só consegue 
interpretar a planilha da maneira desejada quando adicionamos
uma numeração junto aos nomes das estações? Porque no caso do meu 
exemplo, eu não tinha nome de estação duplicada...


Att.,

Marie



To: r-br@listas.c3sl.ufpr.br
From: salah3.1...@gmail.com
Date: Mon, 15 Feb 2016 20:20:02 -0200
Subject: Re: [R-br] Importação de planilha de contingência

Olá Marie

O erro: " double 'row.names' are not allowed" , significa que não é 
permitido duplicatas no row.names.
Sugiro você editar o nomes das estações amostrais numerando ou 
adicionando a data de coleta
No seu exemplo de tabela adicionei AS-B-D e o read.csv2 me deu o mesmo 
erro, quando acrescentei os números( AS-B-D6, AS-B-D7) carregou 
novamente.

Observe o row.names em aravo$spe , AR01, AR02,...,AR59



Abu sp  Abu vai Aca bi
AS-A-C1 0   0   0
AS-A-D2 1   0   0
AS-B-A3 1   0   1
AS-B-B4 0   0   0
AS-B-C5 0   0   0
AS-B-D6 0   0   0
AS-B-D7 0   0   0



O erro para read.xls: perl executable not found
ele não encontrou o executável perl, supondo que você esteja usando 
windows


## arquivo no formato XLSX
dados = read.xls("C://endereço//nomeArquivo.xlsx", sheet = 1, 
row.names = 1, perl = "C://endereço//perl.exe")



saudações

Em 15/02/2016 19:14, Marie-Christine . escreveu:

 double 'row.names' are not allowed



___ R-br mailing list 
R-br@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.


Diese E-Mail wurde von einem virenfreien Computer gesendet, der von 
Avast geschützt wird.

www.avast.com <https://www.avast.com/sig-email>



___
R-br mailing list
R-br@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@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.

Re: [R-br] Importação de planilha de contingência

[salah]

Olá Marie

O erro: " double 'row.names' are not allowed" , significa que não é 
permitido duplicatas no row.names.
Sugiro você editar o nomes das estações amostrais numerando ou 
adicionando a data de coleta
No seu exemplo de tabela adicionei AS-B-D e o read.csv2 me deu o mesmo 
erro, quando acrescentei os números( AS-B-D6, AS-B-D7) carregou novamente.

Observe o row.names em aravo$spe , AR01, AR02,...,AR59



Abu sp  Abu vai Aca bi
AS-A-C1 0   0   0
AS-A-D2 1   0   0
AS-B-A3 1   0   1
AS-B-B4 0   0   0
AS-B-C5 0   0   0
AS-B-D6 0   0   0
AS-B-D7 0   0   0



O erro para read.xls: perl executable not found
ele não encontrou o executável perl, supondo que você esteja usando windows

## arquivo no formato XLSX
dados = read.xls("C://endereço//nomeArquivo.xlsx", sheet = 1, row.names 
= 1, perl = "C://endereço//perl.exe")



saudações

Em 15/02/2016 19:14, Marie-Christine . escreveu:

 double 'row.names' are not allowed


___
R-br mailing list
R-br@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.

Re: [R-br] Importação de planilha de contingência

[salah]

Olá Marie

indique que sua primeira coluna é o row.names, caso contrário ele irá 
carregar como uma variável (coluna) e no R colunas sem nome será 
acrescentado o nome default "X" ( se tiveres várias colunas sem nome, 
será: X, X.1, X.2, ...)

Lembre-se que row.names devem ter nomes únicos(não repetidos).

segue sugestão:

## arquivo no formato CSV
## read.csv2 vem como default sep=";" e dec=","
dados = read.csv2("endereço/nomeArquivo.csv", row.names = 1)

library(gdata)
## arquivo no formato XLSX
dados = read.xls("endereço/nomeArquivo.xlsx", sheet = 1, row.names = 1)

saudações

Em 15/02/2016 16:15, Marie-Christine . escreveu:

Olá Thiago,

Desculpa. Sou nova aqui neste fórum e não sabia a maneira mais correta 
de expor a minha dúvida.
Talvez eu deva contextualizar, antes de tudo, melhor o caso de uma das 
minhas planilhas.

Abaixo exponho-lhe parte da minha planilha:


Abu sp  Abu vai Aca bi
AS-A-C  0   0   0
AS-A-D  1   0   0
AS-B-A  1   0   1
AS-B-B  0   0   0
AS-B-C  0   0   0
AS-B-D  0   0   0







As linhas da primeira coluna correspondem às minhas estações de 
amostragens (AS-A-C, etc). Já a primeira linha da 2 a 3 coluna 
corresponde ao nome das minhas espécies estudadas (Abu sp, etc), com 
suas respectivas abundâncias para cada estação de coleta.



Para importar esta planilha eu usei tanto a função read.csv() quanto 
read.table() e read.delim().


No caso da função read.csv usei os argumentos básicos (sep=";" e 
header=T). A mesma coisa fiz báscimanete com o read.table.
Para função read.delim() eu descobri que, se eu acrescentar o 
argumento check.names=FALSE, eu consigo importar a planilha da forma 
desejada.
Isto é, obtenho um data.frame na forma de tabela de contingência sem o 
''X'' no lugar da primeira célula.


O problema, contudo, é que dessa forma o R interpreta minha primeira 
coluna como sendo um fator com 23 níveis. Eu preciso que  as estações 
também sejam interpretadas como um cabecalho. Se, por exemplo, você 
der um View(aravo$spe - exemplo do pacote ade4), observe que as 
estações de coleta constam como cabecalho (as linhas não estão 
enumeradas).
No meu caso, como o R lê a minha primeira linha como um fator, quando 
dou um View() da minha planilha, minhas linhas estão enumeradas e só 
depois delas é que começa a planilha.


No caso quando eu uso a função ftable() do pacote stats, aparece a 
seguinte mensagem de erro:



Error in table(x, exclude = exclude) :
  attempt to make a table with >= 2^31 elements


E quando eu uso a função read.xls() do pacote gdata, o seguinte erro 
aparece:


Error in findPerl(verbose = verbose) :
  perl executable not found. Use perl= argument to specify the correct 
path.

Error in file.exists(tfn) : invalid 'file' argument


Eu já tentei modificar esse argumento "perl", mas de nada adiantou.

Bom, não sei se consegui ser clara o suficiente desta vez. Caso não, 
posso tentar detalhar mais e inclusive te enviar a minha planilha para 
dar uma olhada.


Desde já muito obrigada pela sua atenção




Date: Mon, 15 Feb 2016 12:03:25 -0200
From: thiago.sera...@gmail.com
To: r-br@listas.c3sl.ufpr.br
Subject: Re: [R-br] Importação de planilha de contingência

Olá Marie.
É difícil te ajudar se você não fornecer algumas informações básicas.

"Sempre que importo, o R atribui um ''X'' no lugar da primeira célula 
da minha planilha, aonde, contudo, não deveria constar nada."
Que função você usou para "importar" o arquivo? Com quais argumentos? 
Tentou usar outros argumentos? Quais?


"Já tentei usar a função ftable() e read.xls(), mas não consegui 
aplicá-las."
Não conseguiu aplicar porque? O que acontece? Também aparece X na 
primeira célula? Ocorre um erro? Se sim, qual a mensagem?


Quais pacotes você está usando? Quando executo read.xls(), recebo a 
mensagem "Error: could not find function "read.xls""


Acho que esse é um bom começo.

Atenciosamente,
Thiago


2016-02-15 11:31 GMT-02:00 Marie-Christine . >:



Olá,

Estou querendo fazer uma análise RLQ pelo pacote ade4.
Para isso, preciso fazer uso de diversas funções do pacote, tais
como dudi.coa(), dudi.hillsmith() e dudi.pca().
Para fazer esta análise necessito usar três planilhas diferentes,
cada qual está no formato de contingência (vide exemplo do pacote
com o banco de dados ''aravo'').
Entretanto estou com problemas para importar as minhas planilhas
do excel (que já estão no formato de tabela de contingência) e
aplicar essas funções em cima delas.
Sempre que importo, o R atribui um ''X'' no lugar da primeira
célula da minha planilha, aonde, contudo, não deveria constar nada.
Já tentei usar a função ftable() e read.xls(), mas não consegui
aplicá-las.

Alguém teria alguma solução para resolver o meu problema?

Att.,

Marie


  

Re: [R-br] Abrir Arquivo de Extensão .nc

[salah]

Olá

sugestão no site

http://stackoverflow.com/questions/16443211/error-when-trying-to-import-netcdf-to-r

|library(ncdf4) mycdf <- nc_open('endereço/nome_arquivo.nc') str(mycdf) |

saudações


Em 10/02/2016 14:56, Yury Duarte escreveu:

Boa tarde colegas programadores!

Estou com alguns problemas em trabalhar com arquivos de extensão .nc.
Acredito que esses são arquivos do tipo raster.
Baixei os pacotes ncdf4; ncdf.tools e o RNetCDF.
Com nenhum desses pacotes consegui fazer a leitura do arquivo, até o 
momento.
Estou selecionando o diretório de onde irei buscar os arquivos .nc e 
em seguida uso o comando library() para instalar um dos três pacotes 
que mencionei acima.
Feito isso, o próximo passo sugerido pelo R é usar um comando 
open.nc/open.ncdf  (dependendo do pacote 
escolhido), entretanto o R não reconhece os comandos open.
Alguém tem familiaridade com esse tipo de arquivo e pode me dar uma 
ajuda nessa questão?


Desde já, agradeço pela disponibilidade de todos!

Abs

Yury Duarte
Engenheiro Agrônomo - ESALQ/USP


___
R-br mailing list
R-br@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@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.

Re: [R-br] problema biblioteca netcdf

[salah]

Olá

Package ‘ncdf’ was removed from the CRAN repository.

Formerly available versions can be obtained from the archive 
.


Archived on 2016-01-11: use 'RNetCDF' or 'ncdf4' instead.

https://cran.r-project.org/web/packages/ncdf/index.html

saudações

Em 01/02/2016 18:25, Mateus Dias Nunes escreveu:

Olá

Quando tento chamar a biblioteca "library(ncdf)" resulta no erro:

> library(ncdf)
Error in library(ncdf) : there is no package called ‘ncdf’

ai quando vou instalar o pacote, apresenta-se esta informação:

package ‘ncdf’ is not available (for R version 3.2.3)

como eu poderia resolver este problema.

att



MATEUS DIAS NUNES
MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET
UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL
TELEFONE: +55 (53) 81125154



___
R-br mailing list
R-br@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@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.

Re: [R-br] Large Merge

[salah]

Tente:

## determine a chave primaria usando data.table
setkey(bancolink, NUMERO_CPF)
CadfAmerge = bancolink[fichaAlink]

saudações

Em 14/01/2016 14:13, Wagner Tassinari escreveu:
CadfAmerge = bancolink[fichaAlink] 


___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Loop para leitura de picos

[salah]

Olá

segue uma sugestão:

vp = function(onda){

pico = NULL
vale = NULL

for(i in 2:(length(onda)-1))
{

## seleciona valores maximos
if(onda[i+1] < onda[i] & onda[i-1] < onda[i])
{
p = onda[i]
pico = rbind(pico, p)
}## end if pico

## seleciona valores minimos
if(onda[i+1] > onda[i] & onda[i-1] > onda[i])
{
v = onda[i]
vale = rbind(vale, v)
}## end if vale
}##end for

return(rbind(pico, vale))
}## end function

## seus dados
amostra = c(49.936157, 48.133488, 45.775238, 43.229397, 40.021607,
35.927558, 31.365912, 26.849848, 23.233177, 20.770278, 18.572067,
15.590341, 11.936978, 8.478520, 5.572115, 3.272254, 1.642041, 0.97,
1.473886, 2.616919, 4.169691, 6.464019, 9.790810, 13.820101, 17.978423,
21.674281, 24.065947, 24.875069, 24.590811, 23.620608, 22.240339,
20.676521, 19.088229, 17.689105, 16.720912, 16.358739, 16.792281,
18.321545, 21.053896, 24.56, 29.121638, 33.551339, 37.628649,
41.272700, 44.405581, 46.899820, 48.602420, 49.253714, 48.816692,
47.380737, 44.920718)
wavelenght = seq(400, 450, by=1)

## monta um data.frame
dados = data.frame(wavelenght, amostra)

## seleciona vales e picos
a = vp(amostra)

## separa a posição dos vales e picos nos dados
pp = dados[c(match(a, dados$amostra)),]

## grafico
plot(dados, type='l')
points(pp, pch = 19, col = 'red')



Em Sáb, 2015-12-19 às 19:04 -0200, José Lucas Safanelli escreveu:
> Boa tarde, estou com o seguinte caso:
> 
> 
> DADO:
> amostra = c(49.936157, 48.133488, 45.775238, 43.229397, 40.021607,
> 35.927558, 31.365912, 26.849848, 23.233177, 20.770278, 18.572067,
> 15.590341, 11.936978, 8.478520, 5.572115, 3.272254, 1.642041,
> 0.97, 1.473886, 2.616919, 4.169691, 6.464019, 9.790810, 13.820101,
> 17.978423, 21.674281, 24.065947, 24.875069, 24.590811, 23.620608,
> 22.240339, 20.676521, 19.088229, 17.689105, 16.720912, 16.358739,
> 16.792281, 18.321545, 21.053896, 24.56, 29.121638, 33.551339,
> 37.628649, 41.272700, 44.405581, 46.899820, 48.602420, 49.253714,
> 48.816692, 47.380737, 44.920718)
> wavelength = seq(400, 450, by=1)
> 
> plot(wavelenght, amostra, xlab  = "Comp. Onda (nm)", ylab = 'Valor',
> type = "l")
> axis(side=1, at=seq(400, 450, by=5))
> 
> 
> PROBLEMA:
> Estou tentando ler os valores mínimos e máximos de cada pico para ter
> o valor da amplitude.
> Consigo fazer isso delimitando o intervalo de comp. de ondas da
> posição do pico, mas é muito laborioso
> EX: max(amostra[,1:20]) e min(amostraamostra[,1:20]). Amplitude seria
> a diferença deles.
> 
> 
> Para várias amostras tenho utilizado a função apply para dentro da
> linha, PORÉM decorrente da heterogeneidade das amostras, as posições
> variam causando erro de leitura e modelagem.
> 
> 
> Outro fato é que existem inúmeros picos dentro de uma faixa de
> interesse, e delimitando visualmente torna-se inviável. Para isso,
> estou tentando criar uma rotina que leia todos os picos dentro da
> faixa de interesse, para que possa utiliza-los diferentemente nas
> modelagens.
> 
> 
> Tentei o código abaixo, mas me retorna o erro que o argumento é de
> comprimento zero! 
> 
> 
> for(i in 1:length(amostra)){
>   if((amostra[i+1] > amostra[i])&(amostra[i-1]>amostra[i]))
>print(amostra[i])
>   }
> 
> 
> (Error in if ((amostra[i + 1] > amostra[i]) & (amostra[i - 1] >
> amostra[i])) print(amostra[i]) : argument is of length zero)
> 
> 
> Um resultado que observei foi de do seguinte código, somente com uma
> condição colocada nele:
> 
> 
> for(i in 1:length(amostra)){
>   if(amostra[i+1] > amostra[i])
>print(amostra[i])
>   }
> 
> 
> Nesse código, consegue-se ler os valores crescentes dos picos,
> desconsiderando os valores do objeto quando os valores são
> decrescentes. A partir dele que tentei colocar duas condições depois
> do IF para que ele lê-se somente os valores mínimos, mas sem sucesso. 
> 
> 
> Portanto, peço ajuda para confecção de um código que faça a leitura
> dos min e max, ou amplitude, de todos os picos da amostra, retornando
> em um objeto/dataframe para modelagem.
> 
> 
> Grato pela atenção! José Lucas.
> ___
> R-br mailing list
> R-br@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 fornea cdigo 
> mnimo reproduzvel.


___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] precisão da operação

[salah]

Esse é o grande Deus da informática: "floating point"

Faça

options(digits = 22)
k = 123456789123456789
sprintf('%1.2f', k)

Sugiro a leitura:

http://www.burns-stat.com/pages/Tutor/R_inferno.pdf

saudações

Em Ter, 2015-09-22 às 07:54 -0300, Cleber N.Borges escreveu:
>  > 12345678912345678+1
> [1] 12345678912345680
> 
>  > 1.2345678912345678+0
> [1] 1.2345678912345679
> 
> Prezados,
> pq exatamente o fato acima apresentado acontece?
> 
> se a capacidade de representação dos números estão
> na faixa de 1.xxe-308 a 1.xxe+308
> 
> 
> agradeço antecipadamente pela atenção
> 
> cleber
> 
> 
> 
> ---
> Este email foi escaneado pelo Avast antivírus.
> https://www.avast.com/antivirus
> 
> ___
> R-br mailing list
> R-br@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 fornea cdigo 
> mnimo reproduzvel.


___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Escala de um eixo

[salah]

Olá

 Veja se lhe agrada

library(reshape2)
library(ggplot2)

taxa = read.table("clipboard",header=TRUE,dec=",")

taxa.m = melt(taxa, id=1, variable.name='Local', value.name='valor')

ggplot(taxa.m, aes(factor(Ano), valor, group=Local, color=Local)) +
theme_bw() +
ylab('legenda y') +
xlab('legenda x') +
geom_line() +
geom_point(aes(shape=Local), size=2)

saudações

___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Erro instalar pacote smatr

[salah]

Olá Marcelo

  Na mensagem de erro diz: there is no package called ‘Rcpp’
  Instale o pacote Rcpp
saudações

___
R-br mailing list
R-br@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 fornea cdigo mnimo 
reproduzvel.

Re: [R-br] Encoding

[salah]

Experimente usar o pacote rgdal

library(rgdal)

mapa1 - readOGR(dsn = BRA_adm1.shp, layer = BRA_adm1, encoding =
latin1)

saudações

___
R-br mailing list
R-br@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.

Re: [R-br] Histograma em escala logarítmica

[salah]

Segue sugestões

r =
c(1017285.7 ,2989151.8 ,4059346.0 ,4256299.0 ,8376572.0 ,9560400.0 ,13007977.0 
,25303118.0 ,33621704.0 ,55951560.0 ,101155984.0 ,144322192.0 ,376711232.0 
,731957760.0 ,9444513800.0 ,9912163300.0 ,9918688300.0 ,9921105900.0 )

## sugestão 1
h = hist(r, plot = FALSE, breaks = 100)
plot(h$counts, log = xy, col = blue, type = 'h', lwd = 2, las = 1) 

## sugestão 2
logr = log(r)
par(las = 2)
h2 = hist(logr, axes = FALSE, breaks = 10, col = rainbow(8))
Axis(side = 2)
Axis(side = 1, at = h2$mids, labels = format(exp(h2$mids), digits = 1),
lty = NULL) 

## sugestão 3
library(ggplot)
qplot(r, geom = histogram, binwidth = 1/5) + scale_x_continuous(trans
= log10)

#sugestão 4
library(fdth)
lr = log(r)
f = fdt(lr, start = 10, end = 24, h = 2)
plot(f)

Em Sex, 2015-06-12 às 12:32 -0300, Romero Luiz M. Sales Filho escreveu:
 Olá pessoal, muito obrigado pela ajuda, entretanto ainda não consegui
 desenvolver esse histograma com as dicas que aparecem no link...
 Tentei de tudo e não deu certo, não sei o que pode está acontecendo...
 
 
 Os dados que eu preciso fazer o histograma são os seguintes:
 
 
 mydata-c(1017285.7 ,2989151.8 ,4059346.0 ,4256299.0 ,8376572.0 ,9560400.0 
 ,13007977.0 ,25303118.0 ,33621704.0 ,55951560.0 ,101155984.0 ,144322192.0 
 ,376711232.0 ,731957760.0 ,9444513800.0 ,9912163300.0 ,9918688300.0 
 ,9921105900.0 )
 
 
 
 
 Se alguém puder me ajudar na construção, ou entendimento do link que
 dá as dicas já me ajudará bastante (como eu disse testei tudo que tem
 lá, mas nada deu certo para esses dados)... Preciso fazer isso o mais
 rápido possível, então toda ajuda é válida!! :)
 
 
 Grande abraço,
 
 
 Romero.
 
 
 -- 
 Profº Romero Luiz Mendonça Sales Filho
 Estatístico / Mestre em Engª de Produção
 UFRPE - UAG
 Cel.: (81) 9667-4477
 romero.sfi...@gmail.com 
 romero_sfi...@yahoo.com.br
 ___
 R-br mailing list
 R-br@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@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.

Re: [R-br] Regressão espacial com variavel transformada log normal

[salah]

tente: log(x + 1)


Em Qua, 2015-06-10 às 10:34 -0300, Samuel luna de almeida escreveu:
 Muito obrigado Luis, Paulo e Luis!
 
 
 São zero mesmo (locais onde não há casos para resultarem em taxas)...
 
 
 Consegui gerar lm para a regressão espacial com os dados originais,
 com muitos zeros, porém com a variável normalizada acho q o -inf
 impede...
 
 
 Será q eu consigo considerar os -inf como zero de algum modo? Ou `e
 o caso de encontrar um modelo/distribuição mais adequado mesmo?
 
 
 Agradecido,
 
 Samuel
 
 
 Em 10 de junho de 2015 10:07, Paulo Justiniano paulo...@leg.ufpr.br
 escreveu:
 
 log(0) = -Inf
 
 Veja que voce tem muitos zeros, mais que 25% dos teus dados,
 visto que o 1o quartil ainda é zero
 
 Não tenhosugestão específica mas acredito que voce poderia:
 - ver se seus zeros são zeros mesmo ou valores censurados
 abaixo de um cento limite de detecção
 
 - sendo zeros mesmo eu procurar por algum modelo/distribuição
 que permitisse modelar esta proporção bem razoável de zeros
 que voce tem nos teus dados
 
 
 
 
 
 On Wed, 10 Jun 2015, Samuel luna de almeida wrote:
 
 Bom dia pessoal,
 
 Trabalhando com uma variável de taxas que que
 apresenta muitos valores baixos, fui aconselhado a
 transforma-la com aplicação de log normal antes de
 buscar
 correlações...
 A variável foi transformada e inserida no dataframe
 passando a apresentar comportamento normal, no entanto
 apresenta alguns valores -inf como pode ser visto no
 summary copiado abaixo:
 
 TX_J04_10 LOGNTX_J04 
 Min.   :   0.0   Min.   :-Inf 
 1st Qu.:   0.0 1st Qu.:-Inf 
 Median : 138.4Median :  -4 
 Mean   : 240.7Mean   :-Inf 
 3rd Qu.: 295.23rd Qu.:  -4 
 Max.   :9626.1Max.   :   0
 
 Acredito que seja por isso que não estou conseguindo
 gerar lm para rodar regressão espacial, veja o erro
 que aparece:
 
  baseGN04_lm - lm(LOGNTX_J04 ~ DV_TOTAL +
 idhm_N_gri, data=baseG04)
 Error in lm.fit(x, y, offset = offset, singular.ok =
 singular.ok, ...) :
   NA/NaN/Inf in 'y'
 
 Será que é por isso mesmo, alguém poderia me indicar
 uma resolução?
 
 
 ___
 R-br mailing list
 R-br@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@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@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.

Re: [R-br] Mapa / dengue

[salah]

Olá, Segue uma sugestão

library(sp)
library(maptools)
library(rgdal)
library(maps)
library(ggplot2)

##setwd('endereço de trabalho que receberá os arquivos')

## site do shapefile - neste site você pode baixar de todos os estados
url0  = 'ftp://geoftp.ibge.gov.br/malhas_digitais/municipio_2010/ba.zip'

## faz o download do shapefile
if (!file.exists(basename(url0))) download.file(url0,
dest=basename(url0), mode = wb)
unzip(basename(url0), list = T) ## mostra o conteúdo do zip
unzip(basename(url0)) ## extrai os arquivos

##
##LENDO O SHAPEFILE
##
## carrega o estado do nordeste com todos os municípios
nordeste = readOGR(dsn = 'endereço do shape/29MUE250GC_SIR.shp', layer =
'29MUE250GC_SIR', encoding = latin1)

## os municípios estão na coluna 3
str(nordeste@data)
##Nome dos municípios
nordeste@data$NM_MUNICIP

## separando Itabuna
itabuna = nordeste[nordeste@data[, 3] %in% 'ITABUNA', ]

##--
##usando o gráfico básico do R
##--
plot(itabuna, axes=TRUE)
points(-39.3, -14.9, pch=19, cex=2, col='red')
points(-39.4, -15, pch=19, cex=5, col='red')
points(-39.35, -14.8, pch=19, cex=1, col='blue')
## acrescenta a escala
map.scale(-39.27, -15, ratio = FALSE, metric = TRUE, cex = 0.8) 

##--
##usando o ggplot2
##--

## dataframe com as coordenadas e quantidade de infectados
dengue = data.frame(long = c(-39.3, -39.4, -39.35), lat=c(-14.9, -15,
-14.8), Infectados=c(30, 12, 2))
dengue

## gráfico
ggplot(itabuna, aes(long, lat, group = group)) +
geom_polygon(fill = white) + ## cor do mapa
geom_path(col = #7f7f7f, size = 0.25) + ## cor e espessura da 
borda
do mapa
## plota os pontos do dataframe dengue
geom_point(data = dengue, aes(x = long, y = lat, 
group=Infectados,
size=Infectados), col='red', alpha=0.7) +
ggtitle('Dengue - Itabuna')

Em Qua, 2015-06-03 às 21:44 -0300, David Feitosa escreveu:
 Oi Luiz,
 
 
 Também iniciei um projeto com mapas e R há pouco tempo.
 Estou usando ggmap como biblioteca de interface para o Google Maps.
 Ele trabalha com o ggplot2, dai você precisa ver que tipo de mapa você
 quer gerar.
 No meu caso é um de densidade. 
 Eu gero o mapa em cinza e ploto a densidade de acordo com as
 ocorrências de uma determinada variável
 em uma localidade.
 
 
 Se o seu caso for similar, podes tentar ver neste link:
 http://rstudio-pubs-static.s3.amazonaws.com/7025_33556c639122442d9eb74e56457a17dc.html
 
 
 
 E, como eu gero uma animação a partir de vários gráficos, o pacote
 animation .
 Eu gero mapas no tempo e uso este pacote para gerar uma animação que
 dá ideia da evolução
 quando o tempo passa.
 
 
 
 
 
 Atenciosamente,
 
 
 David F.
 
 Em 1 de junho de 2015 22:35, Luiz Roberto Martins Pinto
 luizroberto.u...@gmail.com escreveu:
 Caros,
 
 
 Itabuna-BA é um dos municípios brasileiros com maior
 incidência de dengue. Eu desejo elaborar, periodicamente,
 mapas de ocorrência de dengue, para avaliar a sua evolução.
 
 
 Alguém da lista pode me auxiliar?
 
 
 Luiz Roberto.
 
 
 
 Luiz Roberto Martins Pinto
 Prof. Pleno/DCET/UESC
 Laboratório de Estatística Computacional
 Universidade Estadual de Santa Cruz
 Ilhéus-Bahia-Brasil
 
 luizroberto.u...@gmail.com
 skype: lrmpinto
 http://lattes.cnpq.br/2732314327604831 
 
 
 
 
 ___
 R-br mailing list
 R-br@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@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@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.

Re: [R-br] RES: Impressão do objeto Date

[salah]

Olá

Está faltando um parenteses no final de:

write.csv(dfi, paste(diretorio, saida.csv, sep = /, collapse = ))


___
R-br mailing list
R-br@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.

Re: [R-br] Knitr iniciante

[salah]

Caro,

 O objeto luglikp1 deve ser definido dentro do chunk

segue exemplo:

---
title: teste
output: html_document
---


```{r}
summary(cars)
```

You can also embed plots, for example:

```{r, echo=FALSE}
plot(cars)
```
objeto luglipk1 definido dentro do chunk

```{r, echo=FALSE}
luglikp1 = c(1:10)
plot(luglikp1)
```

saudações
salah

___
R-br mailing list
R-br@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.

Re: [R-br] Leitura de arquivo .shp

[Salah]

Olá Hélio

Não entendi muito bem sua pergunta. Veja se isso lhe ajudalibrary(maptools)
library(ggplot2)

PP = data.frame(
long = c(
-45.75753,
-43.78958,
-44.86678,
-43.20956),
lat = c(
-17.45523,
-18.02578,
-19.23028,
-20.10195),
Temp  = c(
23.2,
28.5,
33.7,
39.6),
Local = c(
'A',
'B',
'C',
'D'))

## site do shapefile
##
http://www.igam.mg.gov.br/images/stories/arquivos/ArcView/Minas%20Gerais.zip
minas = readShapeSpatial(/endereço do shapefile/Mg_region.shp)

## acessando as coordenadas como matrix class: SpatialPolygonsDataFrame
minas@polygons[[1]]@Polygons[[1]]@coords

## convertendo para dataframe fica mais fácil
minas - fortify(minas)

minas$long
minas$lat

## plotando o mapa
ggplot(minas, aes(x = long, y = lat)) +
geom_polygon(fill = lightsteelblue2) +
geom_path(colour = grey40) +
geom_point(data = PP, aes(x  = long, y = lat, colour = Local, size
= Temp), alpha = 0.5)

saudações
salah

Em 20 de abril de 2015 19:18, Hélio Gallo Rocha heliogalloro...@gmail.com
escreveu:

 Paulo, obrigado pela indicação do texto, mas ele ão mostra como se extrai
 as coordenadas do polígono.
 Após plotar o polígono o comando text insere uma informação no centro do
 polígono, nada mais.

 Grato

 Em 20 de abril de 2015 16:08, abreups [via R-br] 
 ml-node+s2285057n4664284...@n4.nabble.com escreveu:

 Hélio,

 Você conhece esse artigo:
 http://www.uesc.br/editora/livrosdigitais_20140513/r_cientistas.pdf ?
 Na página 134 tem um capítulo que fala sobre mapas. Mais especificamente,
 nas páginas 138 e 139 há um código sobre como extrair coordenadas
 geográficas.
 Será que ajuda?

 Paulo Abreu



 Date: Mon, 20 Apr 2015 10:09:56 -0300
 From: Hélio Gallo Rocha [hidden email]
 http:///user/SendEmail.jtp?type=nodenode=4664284i=0
 To: r-br [hidden email]
 http:///user/SendEmail.jtp?type=nodenode=4664284i=1
 Subject: Re: [R-br] Leitura de arquivo .shp
 Message-ID:
 [hidden email] http:///user/SendEmail.jtp?type=nodenode=4664284i=2
 Content-Type: text/plain; charset=utf-8

 Caros,

 Andei procurando em alguns CRM's e na lista, mas ão achei e já que estão
 tratando do assunto de arquivos shape, pergunto:

 Como extrair as coordenadas de um polígono.

 No meu caso tenho os dados climatológicos de Minas e o shape do IBGE do
 perímetro de Minas.

 Grato

 ___
 R-br mailing list
 [hidden email] http:///user/SendEmail.jtp?type=nodenode=4664284i=3
 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.

 --
  If you reply to this email, your message will be added to the
 discussion below:

 http://r-br.2285057.n4.nabble.com/R-br-Leitura-de-arquivo-shp-tp4664267p4664284.html
  To unsubscribe from R-br, click here
 http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_codenode=3357982code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4
 .
 NAML
 http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewerid=instant_html%21nabble%3Aemail.namlbase=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespacebreadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml




 --
 Hélio Gallo Rocha
 IFSULDEMINAS - Câmpus Muzambinho

 ___
 R-br mailing list
 R-br@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.




-- 
Izi


Open Source! Use R! Use Linux! Software Free!
___
R-br mailing list
R-br@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.