16:00 lamba, 5% R code run. HQ could not find the output file.
16:25, 10% run on lambda.
This site is to serve as my note-book and to effectively communicate with my students and collaborators. Every now and then, a blog may be of interest to other researchers or teachers. Views in this blog are my own. All rights of research results and findings on this blog are reserved. See also http://youtube.com/c/hongqin @hongqin
I finished another cointegration analysis for 15 month year time window, 2020 May 1 - 2021 August 15. No rank of two is found for cajo Rt ~ T2m. So, either vaccination really alter the correlation, or second winter has to be included.
After filter adf p-value, NE cluster disappeared. Instead, Alaska and west coat are the new cluster.
tb <- tb %>% filter( adf_p_Rt>0.01 & adf_p_t2m >0.05)
(base) hqin@CS313BQin _run-20200501-20210215 % pwd
/Users/hqin/github/COVID19_transmission_MS/_run-20200501-20210215
(base) hqin@CS313BQin _run-20200501-20210215 % ll
total 1496
-rw-r--r-- 1 hqin staff 735K Mar 17 22:49 _autocajo_1-3151-2020-05-01-to-2021-02-15-cval=1percent-K=12trend-longrun.csv
-rw-r--r-- 1 hqin staff 11K Mar 17 22:34 auto_county_cajo_2020501-20210215-1percent-K12.R
Run can pause stochastically. Kill and restart R code three times, and the output is done.
pwd:
(base) hqin@CS313BQin _run-2020501-20200901 % pwd
/Users/hqin/github/COVID19_transmission_MS/_run-2020501-20200901
Rt ~ T2m is very strong!!!
In the summer, T2m are mostly non-stationary!!!
tryCatch fixed this running error.
i = 2677
Error in if (k < 0) stop("k negative") :
missing value where TRUE/FALSE needed
Calls: adf.test
In addition: There were 50 or more warnings (use warnings() to see the first 50)
Execution halted
get_adf_pvalue (tb_local$era5_t2m)
Loving, Texas, There are zero cases in the time window. So, adf.test is meaningless.
tried reinstall R and packages
https://stackoverflow.com/questions/49190251/caught-segfault-memory-not-mapped-error-in-r
Cater, TN
time window April 1 to September 1, 2020.
adf.test(Rt) --> p-value <=0.01. !!!!
This means Rt is stationary and leads to numerical errors in further analysis.
https://www.biorxiv.org/content/10.1101/2022.01.26.477967v1.abstract
3 forms of FadL channel in Vibrio
docking studies
energetic analysis for transportation
https://en.wikipedia.org/wiki/Molecular_dynamics
https://unsupervised.cs.princeton.edu/deeplearningtutorial.html