Olá Listeiros, Tenho uma dúvida bem simples da qual gostaria de pedir ajuda a vocês. Criei um código que gera as evapotranspirações pelo método de Penman, Thronthwait e Blaney-Criddle. A partir desses resultados gostaria de gerar gráficos em dois modos: O primeiro seria pegando o resultado de Penman mostrando sua evolução com os meses do ano O segundo é pegando os métodos de Thornthwait e Blaney-Criddle juntos mostrando, também suas evoluções no ano. Sempre que eu tento plotar algo aparece no meu IDE (RStudio) uma série de caixas. Com este conjunto de dados, os .csv que possuo não conisgo criar os dados desejados. Os dados estão em: http://s000.tinyupload.com/?file_id=49893669914373744636 O código está em: http://s000.tinyupload.com/?file_id=96211791048541413747 Abraços a todos e obrigado desde já. Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435.
Paulo, Você consegue gerar um exemplo reproduzível com uma serie de 10 valores que você pode colocar na própria mensagem de email? # Gerando uma sequência diária de tamanho 100 a partir de hoje. x <- Sys.Date() x <- seq(x, x + 100, by = "1 day") # Gerando um vetor de números aleatórios normais (ruído branco). y <- rnorm(length(x)) library(lattice) class(x) class(y) # Gráfico da série. xyplot(y ~ x, type = "l") À disposição. Walmes.
Olá Walmes, Eu pensei em enviar os códigos e dados por achar melhor a visualização. O que eu fiz foi gerar os resultados das etp0, etp_1 e etp_2. Pegando uma linha que representa um determinado ano tento plotar com um vetor de strings dos meses do ano, porém o resultado não é o desejado, assim como o meu desconhecimento de como por duas curvas num mesmo gráfico. Espero que possa me ajudar! Grande abraço. Paulo Henrique Pimenta Graduando em Meteorologia - IAG USP Celular/WhatsApp/Telegram: +5511981418435 Este email foi enviado de um telemóvel Samsung Galaxy S5 mini Obter o Outlook para Android On Fri, Oct 28, 2016 at 8:20 PM -0200, "Walmes Zeviani" <walmeszeviani@gmail.com<mailto:walmeszeviani@gmail.com>> wrote: Paulo, Você consegue gerar um exemplo reproduzível com uma serie de 10 valores que você pode colocar na própria mensagem de email? # Gerando uma sequência diária de tamanho 100 a partir de hoje. x <- Sys.Date() x <- seq(x, x + 100, by = "1 day") # Gerando um vetor de números aleatórios normais (ruído branco). y <- rnorm(length(x)) library(lattice) class(x) class(y) # Gráfico da série. xyplot(y ~ x, type = "l") À disposição. Walmes. ?
## dados dados <- structure(c(2.71776792884716, 2.80971971474984, 2.65490262437138, 2.81712246434492, 2.96314794044754, 2.48998766335135, 3.90971067836738, 3.56740293299039, 3.12240458083219, 2.44655356979548, 2.60740924526302, 3.04740091582981, 2.68533220439005, 3.11200959394312, 2.79089422500219, 3.31952328998178, 2.70689315323196, 2.9800351726207, 2.78334397202765, 2.701492681034, 2.81262037000008, 2.06557636006153, 3.10254526952989, 2.6060370591154, 2.82776256991182, 2.55475281126045, 2.75669239356911, 2.17642127258464, 2.27914200068542, 2.34369002346959, 2.38569538401917, 2.20725309018804, 2.30034217251805, 2.31125265578062, 2.25004526681248, 2.82017666857997, 2.15477493636968, 2.13869253661224, 2.15582724426898, 2.04481764402476, 2.08895705916863, 1.98549711505164, 1.95236035664275, 1.85956001381077, 1.84849520585033, 1.73612037173621, 1.84606327731254, 2.06050010899242, 2.04340322596147, 1.48059542419918, 1.5897647349459, 1.83493113709698, 1.87723712290912, 1.86135270464316, 2.50520554686019, 1.49145374049223, 2.03706434893404, 1.98181583080439, 2.16558578443014, 1.92228278753351, 1.94135389144282, 1.85041568811821, 2.40033084397389, 2.22902645167497, 2.38070069752142, 2.46428303483512, 2.32698445546991, 2.97655298983297, 2.6730853845772, 2.58224425238168, 2.7792170454293, 2.35224398095467, 1.06029132663049, 1.94432371641046, 2.30470221885002, 2.95965145506887, 2.80586624373959, 2.58419197793719, 2.75901886105652, 2.38893524012871, 0, 2.52015355419062, 2.32376593912811, 2.74572912831904, 2.82729209155284, 3.15276574678139, 2.66311281766592, 3.44658279995033, 2.83634142101896, 0, 2.78936876560345, 3.10384386680356, 2.94733611308899, 3.23818940525282, 2.9025264542456, 2.92016810377702, 3.01739975150593, 2.51447627714607, 0, 3.14272836942838, 2.77071446571446, 2.82080993740916, 3.0813578745847, 3.09390782916134, 3.57884974438802, 3.4796451520399, 3.05046593961632, 0), .Dim = c(9L, 12L), .Dimnames = list(NULL, c("JAN", "FEV", "MAR", "ABR", "MAI", "JUN", "JUL", "AGO", "SET", "OUT", "NOV", "DEZ"))) ## O que você quer é isso? xyplot(ts(dados)) ## Ou isso? xyplot(ts(dados), superpose = TRUE,main='Equação de Penman', auto.key = list(space = "left")) -- Atenciosamente, ================================== Felipe E. Barletta Mendes Estatístico(UFPR) - Conre3 9766-A Mestrando em Bioestatística(UEM) +55 (41)-92077191 +55 (41)-33287216 ===================================
participantes (3)
-
Felipe -
Paulo Henrique Pimenta -
Walmes Zeviani