Monday, March 21, 2022

K18 run

 16:00 lamba, 5% R code run. HQ could not find the output file. 

16:25, 10% run on lambda. 







K16 conintegration

It seems K=16 days is the optimal lag between Rt and T2m. 










 




K14 cointegration


auto_county_cajo_20200501-20210215-01percent-K14.R 













Sunday, March 20, 2022

15 month time window, no Rt-T2m rank of 2

 


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. 




Nat. Biotechnol. | 使用深度学习来注释蛋白质宇宙

 

https://mp.weixin.qq.com/s/iJcS3AP2VXeJSAp_7ftNyA


Friday, March 18, 2022

unit root test

 https://en.wikipedia.org/wiki/Dickey%E2%80%93Fuller_test






cajo Rt - T2m pattern

 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)



Thursday, March 17, 2022

cajo 2020 May 1 to 2021 Feb 15

 

(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. 




2020May1-Sep1 covid19 peak analysis

 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!!!


















adf.test error

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. 

re-install R and Rstudio to remove segfaulty memory not mapped error.

tried reinstall R and packages

https://stackoverflow.com/questions/49190251/caught-segfault-memory-not-mapped-error-in-r


cajo, Carter, TN,

 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. 



Turgeson defense

 

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