cajo debug

/Users/hqin/github/COVID19_transmission_MS

14:40  bug found. There are two dataframes, location and report, but with different rows. Correct the bug with using report consistently. 

bug found when cajo. Rt ~ g_workplace

Hardeman

Tennessee

47069

35.2077202

-88.991612

TRUE

TRUE

date::era5_t2m::era5_dewpoint::dailyCases::Rt::A_driving::g_retail::g_grocery::g_parks::g_transit::g_workplace::g_residential

2

There are only 2 variable, how did rank of 2 returned? 

Eigenvalues (lambda):

[1] 0.06519351 0.04638336 0.00000000

Values of teststatistic and critical values of test:

          test 10pct  5pct  1pct

r <= 1 | 21.51  7.52  9.24 12.97

r = 0  | 52.05 17.85 19.96 24.60

So, the return is correct numerically. 

report$order_cajo_gworkplace[i] = table( ratio_test_cval1Per > 1)['TRUE'] #return the number of TRUE, which 

                                                                                    # give 2 sometimes for 1 pair ca.jo test

HQ spent two hours on finding and fixing the cajo error. 

TODO: convert Rmd into R file.