logica condicional e criar arquivo netCDF

Olá, Eu calculei a média e o percentil 95% (P95), da área que desejo trabalhar, porém preciso fazer duas coisas: 1 - fazer a média somente dos casos acima do P95, ou seja, não consegui achar a lógica para que mostre os valores numéricos onde casos <- UV[UV>=resultado] e fazer a média somente destes casos acima do P95 para cada ponto de grade. 2- criar um arquivo .nc que contenha os valores de "casos". Abaixo o link do meu script https://www.dropbox.com/s/ph03e1b9l6ziylg/TESTE.R?dl=0 E aqui abaixo o link do dado .nc que estou utilizando https://www.dropbox.com/s/81k73nffpt1zwq7/310times_fillmissJAN.nc?dl=0 Desde já agradeço, ____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET GRUPO DE ESTUDOS EM POLUIÇÃO ATMOSFÉRICA - GEPA UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 98112-5154 SKYPE: mateus.diasnunes ____________________________________________________________________________ Sent with Mailtrack <https://mailtrack.io/install?source=signature&lang=en&referral=nunes.mateusdias@gmail.com&idSignature=22>

Mateus, A lógica que você propõe no seu post é a correta e funcionará se em " casos <- UV[UV>=resultado]" resultado for um escalar. No seu script a linha que parece ser similar ao que você quer é "acima_95 = as.matrix(Con(UV>=Percentile95,UV,Con(UV>=Percentile95, UV, 0))) " onde 'Percentile95' não está definido, e o mais próximo dele seria 'percentile95' (lembre-se que o R é sensível a maiúsculas versus minúsculas) que é resultado de uma chamada a levelplot. . . SE seu problema for com "a lógica" para obter valores acima de um certo limiar, refaça seu CMR com menos dependências nos pacotes que você chama no seu script, por outro lado se sua dúvida for em alguma função de algum dos pacotes, por favor seja mais centrado no CMR para a gente ver qual é sua dúvida exatamente. HTH -- Cesar Rabak 2016-11-25 2:00 GMT-02:00 Mateus Dias Nunes via R-br < r-br@listas.c3sl.ufpr.br>:
Olá,
Eu calculei a média e o percentil 95% (P95), da área que desejo trabalhar, porém preciso fazer duas coisas:
1 - fazer a média somente dos casos acima do P95, ou seja, não consegui achar a lógica para que mostre os valores numéricos onde casos <- UV[UV>=resultado] e fazer a média somente destes casos acima do P95 para cada ponto de grade.
2- criar um arquivo .nc que contenha os valores de "casos".
Abaixo o link do meu script
https://www.dropbox.com/s/ph03e1b9l6ziylg/TESTE.R?dl=0
E aqui abaixo o link do dado .nc que estou utilizando https://www.dropbox.com/s/81k73nffpt1zwq7/310times_fillmissJAN.nc?dl=0
Desde já agradeço,
____________________________________________________________ ________________
MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET GRUPO DE ESTUDOS EM POLUIÇÃO ATMOSFÉRICA - GEPA UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 98112-5154 SKYPE: mateus.diasnunes ____________________________________________________________ ________________
Sent with Mailtrack <https://mailtrack.io/install?source=signature&lang=en&referral=nunes.mateusdias@gmail.com&idSignature=22>
_______________________________________________ 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.
participantes (2)
-
Cesar Rabak
-
Mateus Dias Nunes