Segue uma ideia inicial, que você pode ajustar às suas necessidades.
### code r
a - c(3.4, 4.3)
b - c(.53, .56)
aRange - c(0, 6)
bRange - c(0, .6)
a1 - t(scale(a-aRange[1], center=F, scale=diff(aRange)))
b1 - t(scale(b-bRange[1], center=F, scale=diff(bRange)))
plot1 - barplot(-a1, col=8,
Caro Alisson, boa tarde!
Embora eu ainda esteja engatinhando no Python, o código me pareceu
relativamente simples (depois de pronto, parece fácil mesmo! :D). Isso
reforça minha convicção de que deve haver uma forma fácil de fazer isso no
R também...
Vou insistir um pouco mais... e não tendo
Teodoro, bom dia!
Procure pelos pacotes httr, XML e RCurl.
Se o site não exigir login/senha, será bem mais simples e você poderá
operar com os dois primeiros.
Caso seja necessário passar o login/senha vale o tópico indicado
anteriormente (RCurl).
Prezado Elias,
Fornecendo dados junto com seu código, permitindo torná-lo reproduzível
(CMR), você obtém ajuda mais facilmente.
Suponho que o 17 seja em função da variável que você quer representar no
grid que está exportando (predição). Estou supondo que 17 seja o índice da
última coluna
Prezado Elias, boa noite!
Pra ajudar de fato seria necessário ter acesso aos dados que você está
manipulando (ou uma amostra). Mesmo assim vou arriscar um palpite.
Você removeu NA's com base somente na coluna 14 (gride[,14]) e por fim,
normaliza de 4 a 18 (grid1[,4:18]). É provável que outras
Prezado Thiago, boa tarde!
Acredito que o link a seguir possa lhe ser útil:
http://stackoverflow.com/questions/11141406/reshaping-an-array-to-data-frame
São sugeridas soluções por adply {plyr}, melt {reshape} e ftable {stats}.
Se não for exatamente isso que procura, ainda dá pra trabalhar em
Caro Tito, bom dia!
Apenas cogitando... Os arquivos podem ser HDF4 e não HDF5. Não sei se foi
corrigido, mas no passado o pacote hdf5() não abria arquivos HDF4
diretamente e era necessário convertê-los primeiro.
No endereço que segue abaixo tem inclusive o link pra um conversor:
[R] hdf files
of shorter object length
Obrigada!
Em 10 de agosto de 2013 16:41, Eder Comunello ecomu...@gmail.comescreveu:
Boa tarde,
Não entendi muito bem, mas para a saída que você apresenta seria algo
assim:
n - which(dados$spp==Dalbergia ecastophyllum)
x - dados$num[n]
dados[dados$num==x,]
Sugiro
Boa tarde,
Não entendi muito bem, mas para a saída que você apresenta seria algo
assim:
n - which(dados$spp==Dalbergia ecastophyllum)
x - dados$num[n]
dados[dados$num==x,]
Sugiro dar um dput no objeto de dados e colar na próxima mensagem. Espero
que ajude,
At.te,
Em 10 de agosto de 2013
Alexandre e colegas, bom dia!
Meu script está um pouco diferente, mas o que vai valer mesmo é o último
bloco (reprojetar RasteLayer). Espero que ajude.
Além disso, tenho um bookmark pra indicar:
Spatial data in R: using R as a GIS (old version)
Prezado Alexandre,
Quanto ao erro de spTransform(), faltou dar espaço antes do argumento
'+south'.
### INICIO
### bnds fornecido no CMS
bnds - structure(c(772259.1, 775259.1, 775259.1, 772259.1, 772259.1,
7843046, 7843046, 7840046, 7840046, 7843046), .Dim = c(5L, 2L
),
Senhores, bom dia!
Além das dicas dos colegas, pode ser interessante dar uma olhada no pacote
{aRT} em http://www.leg.ufpr.br/doku.php/software:art.
O pacote {aRT} faz a integração do R com a biblioteca TerraLib, permitindo
assim integrar também o TerraView (INPE).
Ainda estou para estudar
Hélio e colegas, bom dia!
Tenho interesse no tema, mas também tenho dúvidas em relação à forma de
implementar este grupo. Pra mim é bastante cômodo que seja aqui na lista
mesmo, pois aproveito pra aprender sobre o uso do R em outras vertentes que
também me interessam.
Taynãna, boa tarde!
Poderia enviar a saída do console no email?
Para mim a saída do último comando aparece assim:
dados[order(dados$local),] ### ordenando
dia local dado1 dado2
1 1 1 0.1015
3 2 1 0.1516
5 3 1NANA
7 4 1 0.2020
2 1 2 0.27
Olá, pessoal! Testei as dicas passadas e esta variação também parece
atender a pedida.
dia - c(1,2,4,1,4)
local - c(1,1,1,2,2)
dado1 - c(0.10,0.15,0.20,0.27,0.15)
dado2 - c(15,16,20,23,25)
dados - data.frame(dia,local,dado1,dado2); dados
tmp -
Daniel, boa noite!
Segue uma sugestão, alojando as conversões em novas colunas.
O formato de saída é configurável e o uso de as.factor() funcionaria, mas
iria cria levels para cada 'data' (n=103).
a- data.frame(Date.end = as.Date(seq(from=as.POSIXct('2012-06-26',
tz=UTC),
Olá, bom dia!
Sugiro verificar se o seu arquivo dados.r é realmente um arquivo de dados
com file.show('dados.r'). Com essa extensão era esperado ser o script de
comandos. Se for realmente o arquivo de dados com uma extensão
inapropriada, verifique os argumentos sep e header de read.table().
Assumindo o conteúdo do dput armazenado em um objeto dados...
dados2 - dados$afia3[-which(is.na(dados$afia3))]
Pra funcionar com essa solução deve-se endereçar o elemento da lista
'afia3'.
Éder Comunello
Ph.D. Student in Agricultural Systems
Carlos e colegas, bom dia!
Vou aproveitar um script anterior, fazendo algumas alterações. Como da
outra vez, creio que dá pra adaptar facilmente pro seu caso.
Vou trabalhar extraindo o vetor de interesse, substituindo este vetor no
objeto origem ao final...
# begin #
require(sp);
Edson e colegas, bom dia!
Ao invés de deixar o campo 'motivo' em branco (NA), experimente atribuir um
valor, p. ex. 'sem resposta', e rode novamente. Assim estes casos serão
considerados e aparecerão na saída. Vale pros outros campos.
Éder
Fernando e colegas, boa noite!
Para um conjunto de dados nominado 'dados', com 4 colunas:
head(dados)
X1 Y1MO ALT
1 726023.6 7534573 29.50 387.1
2 726042.6 7534567 31.60 386.7
3 726014.5 7534544 33.71 387.1
4 726033.6 7534538 34.06 386.6
5 726005.5 7534516 32.65 387.1
6
Mauro e colegas, boa tarde!
As operações com polígonos do {maptools} dependem do pacote {gpclib}
que tem licença restrita.
Ele é desativado por default, por isso usa-se o gpclibPermit().
Verifique a presença do pacote {gpclib}. Pra funcionar gpclibPermit()
deve retornar TRUE.
PS: peço enviar as
Boa tarde, Samuel e colegas!
Segue uma sugestão de código...
Testei com a delimitação dos bairros de Curitiba, do material do Prof.
Paulo Justiniano
Para testar será necessário baixar os 3 arquivos disponibilizados a seguir:
https://dl.dropbox.com/u/117618178/ctba/bairros.shp'
Olá, pessoal!
Tentei dar um 'up' no CMR, seguindo a dica do Walmes...
###
#install.packages(c(geoR, maptools, sp), dep=T)
require(geoR); require(maptools); require(sp)
### download direto! ;D
### joga na pasta/dir de trabalho ativa
Olá, Ligia!
Tem um material que encontrei hoje que está bem didático. É do Prof.
Landim, também de Rio Claro. Você consegue o texto completo no link que
consta na referência que segue.
LANDIM, Paulo M. Barbosa. Sobre Geoestatística e mapas.* Terrae didat.*,
Campinas, v. 2, n. 1, 2006 .
Lucas e colegas, bom dia!
Tem que operar com as datas no formato interno do R, só transformando pra
visualizar.
Defini as datas de corte nas vars data.ini e data.fim...
Alterei todos os valores de p1$num para '1', então se a soma estiver
correta corresponderá ao número de eventos selecionados.
Exemplo de operação com datas com seleção por subset
dados-data.frame(
data=c(01/01/2004,01/02/2004,01/03/2004,01/04/2004,01/05/2004,01/06/2004),
tmin=c(17.82966,18.05029,22.47399,28.20188,13.46510,18.82721),
tmax=c(44.23301,21.11526,22.79882,35.31366,35.84991,39.14286)); dados
Natalia e colegas, bom dia!
Ainda não estou muito familiarizado com o R, mas acho que este exemplo que
rodei irá ajudar...
Não utilizei o project(), mas sim o spTransform()...
#
#install.packages(rgdal, dep=T)
#install.packages(sp,
Kin, boa tarde!
O erro listado é por atribuir ao parâmetro pch o objeto 'membros' que não
existe ou não está acessível (veja seu arquivo de dados e declarações).
Acredito que a ideia seria diferenciar os pontos (tipos) de acordo com o
valor em membros...
Acredito que o exemplo que segue faça
29 matches
Mail list logo
