Na parte de leitura do ASCII, você não queria um stack a partir do r_ascii? Me
parece que as duas funções extract estão fazendo a mesma coisa - extraindo os
valores do stack s (criado do tiff) usando as coordenadas especificadas.
Não rodei o script, mas tenho uns palpites baseados em
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).
Em geral eu evito usar esses pacotes de baixo nível para abrir netcdf (ncdf4,
RNetCDF) porque eles normalmente exigem muito código para acessar os dados que
realmente interessam.
Se você usar o pacote raster, abrir e escrever um arquivo netcdf pode ser tão
fácil quanto:
Obrigado pela resposta Walmes,
Estou tentando fazer testes two-way post-hoc para dados
multivariados, isso baseado em testes one-way post-hoc que tenho
encontrados na literatura, que usam TukeyHSD() que foi desenvolvidos
para dados univariados, mas veem sendo usados para dados
Boa noite Éder,
Montei uma rotina abaixo para iniciarmos um novo tópico e é
possível observar o que você havia comentado, sendo:
###
require(raster)
require(sp)
require(rgdal)
# RasterLayer inventado
r <- raster(nrows=1, ncols=1)
r <- setValues(r, 1:ncell(r))
r_SGDF<- as(r,
Olá,
eu gerei alguns resultados de campo médio de Ozonio, a partir de um arquivo
.nc (netCDF)
a partir deste meu resultado eu gostaria de gerar um novo arquivo no
formato ascii ou .nc porém não consegui gerar esse arquivo de saída.
abaixo a estrutura do Script
# Carregando biblioteca para
Pedro, boa tarde!
No site do Google Maps você deve utilizar latitude, longitude. Na saída da
geocode(), por padrão ,longitude vem primeiro.
require(ggmap)
adr1 <- "Dourados, MS, Brazil, 79823-460"
res1 <- geocode(adr1, output="more")
res1[c(1:4,10)]
# lon lattype
Ei Eder e Elias,
O material indicado é bem o que eu queria pra começar mesmo. Muito
obrigado. Já percebi que no pacote ggmap não há uma funçao georef, mas há
uma geocode que consegue capturar do google maps e alternativamente do
datascience toolkit o georeferenciamento a partir do endereço. Achei
Alexandre, boa tarde!
Sobre o formato dos dados, fiz uns testes aqui e dá pra tirar umas
conclusões, mas acho melhor criar um tópico novo pra não embolar.
Resumidamente, ao acessar usando o pacote raster (antes que rgdal) você
mantém a informação no disco e economiza memória. Você poderia usar um
Charlenny,
Com as escusas de praxe de eu não conhecer esse pacote especificamente,
veja se ele ao executar o comando de decomposição ele retorna um objeto
específico (deve ser um dataframe "customizado" ou coisa que o valha) e se
no pacote há um método para plotar esse objeto (ressaltando, em
Funcionou perfeitamente!
Obrigado, Thiago e Éder!
Att
Yury Duarte
Engenheiro Agrônomo - ESALQ/USP
2016-02-17 15:16 GMT-03:00 Éder Comunello :
> Thiago, boa tarde!
>
> testei aqui...
>
> as.Date(substr(names(s), 2,11), "%Y.%m.%d")
>
>
>
Thiago, boa tarde!
testei aqui...
as.Date(substr(names(s), 2,11), "%Y.%m.%d")
Éder Comunello
Agronomist (UEM), MSc in Environ. Sciences (UEM)
DSc in Agricultural Systems Engineering (USP/Esalq)
Brazilian Agricultural Research Corporation
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=","
##
Alexandre,
Objetos da classe adonis não retornam vcov e a LSmatrix não foi escrita
para receber objetos dessa classe também. Daria para escrever a matriz na
mão e lá na frente passar ela para glht estimar os contrastes, etc. Porém,
para isso, precisa-se de vcov. Além do mais, você está
Prezados,
boa tarde!
Construí uma 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
Uma solução é usar substr
as. Date(substr(X.2007.01.01, 3, 12))
Não estou perto do computador agora para testar esse código, mas com pouca ou
nenhuma alteração ele deve funcionar.
(veja tb o primeiro código q te enviei, abaixo do comentário sobre consertar o
dataframe)
Sent from Yahoo Mail
Minhas desculpas pela falta de atenção. Esse repositório não é mais público
porque alguns dos scripts têm endereços de IP e porta. Assim que eu
escondê-los, eu deixo público. O script que você precisa, no entanto, eu
coloquei em um snippet: https://gitlab.c3sl.ufpr.br/snippets/33.
À disposição.
Bom dia colegas!
Gostaria de agradecer a todos pelas ajudas com o código e também pela
orientação na utilização desse tipo de arquivo, todas as dicas estão sendo
muito bem utilizadas!
Minha última duvida acerca desse tema seria quanto a extração das datas de
leitura do arquivo.
Sempre que extraio
Prezados bom dia!
Estou utilizando o comando decompose(logdados,type="additive") e ele me
retorna o gráfico de decomposição de uma série, com as instruções no
gráfico tais como: tendência, sazonalidade, aleatoriedade, série observada
e título do gráfico em inglês, e preciso passar para português.
19 matches
Mail list logo