Load the simulated nutrition dataset from this package
Source:R/data_io.R, R/datasets.R
load_simulated_nutrition_data.RdThe full simulated nutrition dataset used in the tutorial is shipped as an RDS
under inst/extdata/ (to avoid loading it automatically).
Loads the file inst/extdata/simulated_nutrition_data.rds.
Examples
nut <- load_simulated_nutrition_data()
str(nut)
#> List of 14
#> $ Y_mat : num [1:15021, 1:42] 2.652 -0.472 -1.917 3.265 6.645 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : chr [1:15021] "1" "2" "3" "4" ...
#> .. ..$ : NULL
#> $ Y_list :List of 12
#> ..$ : num [1:1362, 1:42] 2.652 -0.472 -1.917 3.265 6.645 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:1362] "1" "2" "3" "4" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:217, 1:42] -1.759 0.613 -6.417 3.336 -1.319 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:217] "1363" "1364" "1365" "1366" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:417, 1:42] 2.671 0.288 3.039 1.862 7.672 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:417] "1580" "1581" "1582" "1583" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:1012, 1:42] 0.517 8.507 5.274 0.841 5.233 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:1012] "1997" "1998" "1999" "2000" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:2241, 1:42] 0.577 0.624 -0.172 6.048 0.719 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:2241] "3009" "3010" "3011" "3012" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:205, 1:42] 5.4242 2.3453 0.0947 -0.7641 1.0229 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:205] "5250" "5251" "5252" "5253" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:2403, 1:42] 2.026 2.434 -0.268 1.051 -1.021 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:2403] "5455" "5456" "5457" "5458" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:3775, 1:42] 2.25 3.59 3.53 4.06 2.64 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:3775] "7858" "7859" "7860" "7861" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:1790, 1:42] 5.95 5.97 7.04 5.54 4.89 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:1790] "11633" "11634" "11635" "11636" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:761, 1:42] 3.129 6.97 0.472 0.462 4.993 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:761] "13423" "13424" "13425" "13426" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:373, 1:42] 2.8 4.66 3.03 7.31 9.65 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:373] "14184" "14185" "14186" "14187" ...
#> .. .. ..$ : NULL
#> ..$ : num [1:465, 1:42] 4.246 3.899 6.168 1.545 0.268 ...
#> .. ..- attr(*, "dimnames")=List of 2
#> .. .. ..$ : chr [1:465] "14557" "14558" "14559" "14560" ...
#> .. .. ..$ : NULL
#> $ N_s : num [1:12] 1362 217 417 1012 2241 ...
#> $ A : num [1:15021, 1:12] 1 1 1 1 1 1 1 1 1 1 ...
#> $ X : num [1:15021, 1:32] 1 0 1 0 0 0 1 1 1 0 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : chr [1:15021] "1" "2" "3" "4" ...
#> .. ..$ : chr [1:32] "Age_Group18-44" "Age_Group45-55" "Age_Group55-74" "SexFemale" ...
#> $ Beta : num [1:42, 1:32] 0.997 0.714 0.208 0.23 0.933 ...
#> $ Phi : num [1:42, 1:4] -0.94 0.692 0 -0.862 0 ...
#> $ SigmaPhi : num [1:42, 1:42] 2.1699 -1.0863 0 1.2546 0.0267 ...
#> $ LambdaList :List of 12
#> ..$ : num [1:42, 1:2] 0.96 0 -0.997 0 0 ...
#> ..$ : num [1:42, 1:3] 0.96 0 -0.997 0 0 ...
#> ..$ : num [1:42, 1:3] -0.972 -0.738 0 0.653 0 ...
#> ..$ : num [1:42, 1:3] -0.972 -0.738 0 0.653 0 ...
#> ..$ : num [1:42, 1:2] -0.864 0.825 0 0 0 ...
#> ..$ : num [1:42, 1:3] 0.96 0 -0.997 0 0 ...
#> ..$ : num [1:42, 1:3] -0.664 -0.836 0 0 0.739 ...
#> ..$ : num [1:42, 1:2] -0.911 0 -0.857 0 -0.997 ...
#> ..$ : num [1:42, 1:2] 0.608 0.72 -0.781 0 0 ...
#> ..$ : num [1:42, 1:3] -0.664 -0.836 0 0 0.739 ...
#> ..$ : num [1:42, 1:3] -0.664 -0.836 0 0 0.739 ...
#> ..$ : num [1:42, 1:3] -0.864 0.825 0 0 0 ...
#> $ SigmaLambdaList:List of 12
#> ..$ : num [1:42, 1:42] 1.327 0 -0.957 0 0.487 ...
#> ..$ : num [1:42, 1:42] 2.8304 0.0665 -0.9575 0.1537 0 ...
#> ..$ : num [1:42, 1:42] 1.385 1.272 0 -0.635 -0.491 ...
#> ..$ : num [1:42, 1:42] 1.69125 0.00459 0 -0.63507 0 ...
#> ..$ : num [1:42, 1:42] 0.747 -0.713 0 0 0 ...
#> ..$ : num [1:42, 1:42] 2.432 0 -0.849 -0.525 1.715 ...
#> ..$ : num [1:42, 1:42] 1.27 0.555 0.781 0 0.417 ...
#> ..$ : num [1:42, 1:42] 0.829 0 0.781 0 0.908 ...
#> ..$ : num [1:42, 1:42] 0.369 0.437 -0.474 0 0 ...
#> ..$ : num [1:42, 1:42] 0.81 0.992 -0.474 0 -0.491 ...
#> ..$ : num [1:42, 1:42] 2.226 -0.205 0 -0.588 0.384 ...
#> ..$ : num [1:42, 1:42] 1.537 -1.473 0 -0.588 0 ...
#> $ SigmaMarginal :List of 12
#> ..$ : num [1:42, 1:42] 3.927 -1.086 -0.957 1.255 0.514 ...
#> ..$ : num [1:42, 1:42] 5.4487 -1.0198 -0.9575 1.4083 0.0267 ...
#> ..$ : num [1:42, 1:42] 4.52 0.186 0 0.619 -0.464 ...
#> ..$ : num [1:42, 1:42] 4.0678 -1.0817 0 0.6195 0.0267 ...
#> ..$ : num [1:42, 1:42] 3.082 -1.7991 0 1.2546 0.0267 ...
#> ..$ : num [1:42, 1:42] 4.937 -1.086 -0.849 0.73 1.741 ...
#> ..$ : num [1:42, 1:42] 4.413 -0.532 0.781 1.255 0.444 ...
#> ..$ : num [1:42, 1:42] 3.122 -1.086 0.781 1.255 0.934 ...
#> ..$ : num [1:42, 1:42] 3.1697 -0.6489 -0.4744 1.2546 0.0267 ...
#> ..$ : num [1:42, 1:42] 3.5691 -0.0943 -0.4744 1.2546 -0.464 ...
#> ..$ : num [1:42, 1:42] 4.709 -1.292 0 0.666 0.411 ...
#> ..$ : num [1:42, 1:42] 4.2798 -2.5591 0 0.6664 0.0267 ...
#> $ Psi_list :List of 12
#> ..$ : num [1:42, 1:42] 0.431 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.448 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.965 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.207 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.166 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.335 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.973 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.123 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.631 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.589 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.313 0 0 0 0 ...
#> ..$ : num [1:42, 1:42] 0.573 0 0 0 0 ...
#> $ Psi_mat : num [1:42, 1:12] 0.431 0.684 0.639 0.587 0.132 ...
#> $ big_T : num [1:12, 1:23] 1 1 1 1 1 1 1 1 1 1 ...