Monday, October 5, 2020

how to cite ERS5-land

 


How to cite ERA5-Land

(1) Please acknowledge the use of ERA5-Land as stated in the Copernicus C3S/CAMS License agreement:

  • "5.1.2 Where the Licensee communicates or distributes Copernicus Products to the public, the Licensee shall inform the recipients of the source by using the following or any similar notice: 'Generated using Copernicus Climate Change Service Information [Year]'.

  • 5.1.3 Where the Licensee makes or contributes to a publication or distribution containing adapted or modified Copernicus Products, the Licensee shall provide the following or any similar notice: 'Contains modified Copernicus Climate Change Service Information [Year]';

Any such publication or distribution covered by clauses 5.1.1 and 5.1.2 shall state that neither the European Commission nor ECMWF is responsible for any use that may be made of the Copernicus Information or Data it contains."

(2) cite the ERA5-Land dataset (as part of the bibliography) as follows:

Muñoz Sabater, J., (2019): ERA5-Land hourly data from 1981 to present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). (<date of access>), 10.24381/cds.e2161bac

Muñoz Sabater, J., (2019): ERA5-Land monthly averaged data from 1981 to present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). (<date of access>), 10.24381/cds.68d2bb30


Thursday, October 1, 2020

cross corelation

https://youtu.be/6ao9-39zw40

absolute humidity

Bukhari, 2020, Enviromental Resaerch and Public Health, MDPI

cross correlation between two time series

correlation is a linear measure of similarity between two signals. Cross-correlation is somewhat a generalization of the correlation measure as it takes into account the lag of one signal relative to the other. If lag == 0, then correlation = cross-correlation. Cross-correlation is particularly important to assess the causal relationship between two signals in time. If you suspect that there is a non-linear relationship between the two signals, then you should consider measures such as mutual information and partial mutual information, which are the information-theoretic equivalent of correlation and cross-correlation. 


You can refer to the comment above for information on how to compute correlation and cross-correlation. Check this paper instead (and references therein) for details on mutual information and partial mutual information http://www.sciencedirect.com/science/article/pii/S0022169400003462 https://youtu.be/6ao9-39zw40 https://youtu.be/L6YJqhbsuFY  

 ccf in R Ziwei Ma:
I checked the ccf function which calculates the correlation for x_{t+h} and y_t, so in our case, the ccf( dailyCases (x), dew_points (y) ) report they are positively relative, and the peak happen at arround 7 which says dailyCases lags dew_points or dew_points predict dailyCases in 7 days. 

The following website have more details there. https://online.stat.psu.edu/stat510/lesson/8/8.2

Wednesday, September 30, 2020

Genome-wide analysis of SARS-CoV-2 virus strains circulating worldwide implicates heterogeneity

https://www.nature.com/articles/s41598-020-70812-6#Sec2 Article Open Access Published: 19 August 2020 Genome-wide analysis of SARS-CoV-2 virus strains circulating worldwide implicates heterogeneity
Oct 8, Hamilton county to allow big party in restaurant? https://newschannel9.com/news/local/gov-lee-lifts-restrictions-bradley-co-restaurant-calls-executive-order-a-god-send Marl 11, 2020 WHO pandemic declaired.

ARHGAP11B

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

computational genomics with R

http://compgenomr.github.io/book/

worldmet

Functions to import data from more than 30,000 surface meteorological sites around the world managed by the National Oceanic and Atmospheric Administration (NOAA) Integrated Surface Database (ISD, see ). https://cran.r-project.org/web/packages/worldmet/index.html https://github.com/davidcarslaw/worldmet Hong note: it seems worldmet only provide pre-defined sites. In Bukhari et al, 2020, effects of weather on coronavirus pandmic, USA weather were provided at the state level. For many other countries, weather were estiamted from the busiest airports.

Tuesday, September 29, 2020

EpiEstim cran pacakge

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

EpiNow2 package

https://github.com/epiforecasts/EpiNow2 https://epiforecasts.io/EpiNow2/reference/estimate_infections.html#examples https://github.com/epiforecasts/covid-rt-estimates


These worked
reporting_delay <- bootstrapped_dist_fit(rlnorm(100, log(4), 1), max_value = 30)
generation_time <- get_generation_time(disease = "SARS-CoV-2", source = "ganyani")
incubation_period <- get_incubation_period(disease = "SARS-CoV-2", source = "lauer")

estimates <- epinow(reported_cases = local_cases, samples=100, output='samples',
                    generation_time = generation_time,
                    delays = list(incubation_period, reporting_delay))

localRtTb = estimates$estimates$summarised[ estimates$estimates$summarised$variable=='R' , ]



It seems the default setting is recommended. 
# Note: all examples below have been tuned to reduce the runtimes of examples # these settings are not suggesed for real world use. # run model with default setting


It seems that by default, Rt varies. 
# run model with stationary Rt assumption (likely to provide biased real-time estimates) stat <- estimate_infections(reported_cases, generation_time = generation_time, delays = list(incubation_period, reporting_delay), stan_args = list(warmup = 200, cores = ifelse(interactive(), 4, 1), control = list(adapt_delta = 0.9)), stationary = TRUE)

Rt D3 package

https://epiforecasts.io/RtD3/

Monday, September 28, 2020

possible criteria of a head.

Ability and experience evaluating faculty and staff work performance in an equitable and clearly documented way Ability to communicate with faculty in an open and transparent way Ability to foster an environment in which everyone is comfortable contributing, and in which all voices are considered Ability to unify and support a large, complex, multidisciplinary department that maintains and promotes its disciplinary identities, strengths, and resources Ability to resolve conflict to the advantage of the department as a whole Accessibility and approachability to faculty, staff, and students Commitment and ability to attract and support a diverse and inclusive faculty, staff, and student body Effective advocacy of department needs and promotion of department strengths to upper administration Experience in an administrative role that includes program, budget, and personnel management Motivation and dedication to develop and maintain positive relationships with donors, alumni, and the professional community Support of faculty excellence in teaching across all modalities Support of faculty excellence in research including dissemination and funding Understanding that the department head position is complex, time consuming, and requires total commitment and dedication to the role Other (Please describe any other qualification(s), ability(s), and/or skill(s) that you think should be prioritized in our Department Chair search. Please only describe priorities here that would be among your top five priorities. In other words, priorities that you describe here plus those that you choose from the priorities listed above should amount to no more than five total.)

Friday, September 25, 2020

down load ERA5 land hourly data

 

ERA5 land hours data was last updated on July 12. 

I selected all 2020, all day, and all hours, it is 45.8 GB, which will take 12 hours to download at home.


I ssh into my simcenter workstation, 

wget url, 

which runs with an estimated time of a little over 2 hours.