Hi,
I have exposure data of PM2.5 on two different populations, lung function and demographic variables. I want to show the association between exposure, lung function and demographics. I ran 'xtgee FEV1pred_per pm25, family(gaussian) link(identity) corr(exchangeable) robust' but got an matsize too small. I increased the matsize but still did not get it to work. I am using Stata/MP 14.1 for mac and recently started working with Stata. Please kindly assist on how I can examine the association between these variables. Thank you in advance.
My code looks like this:
.
I have exposure data of PM2.5 on two different populations, lung function and demographic variables. I want to show the association between exposure, lung function and demographics. I ran 'xtgee FEV1pred_per pm25, family(gaussian) link(identity) corr(exchangeable) robust' but got an matsize too small. I increased the matsize but still did not get it to work. I am using Stata/MP 14.1 for mac and recently started working with Stata. Please kindly assist on how I can examine the association between these variables. Thank you in advance.
My code looks like this:
Code:
* Example generated by -dataex-. To install: ssc install dataex clear input long id float type byte age float FEFpred_per int(FEV1pred_per FEV1_FVCpred_per FVCpred_per PEFpred_per) byte(YearsDis MaritalStat Edu HHstove cookspace HHcooking) double datetime long n float(pm25 duration) 1 0 28 39 61 70 91 28 2 1 4 2 2 1 1937282760000 6 59.60243 2976 2 0 27 26 47 93 52 42 1 1 3 1 1 3 1938318720000 5 5 3124 3 0 24 32 58 92 65 39 1 2 2 1 1 1 1.9377129e+12 6 9.323024 3095 4 0 37 12 12 26 49 9 2 1 1 1 1 1 1.937472e+12 6 7.182249 2933 5 0 19 16 45 69 68 30 3 2 2 2 2 1 1938145440000 3 26.010786 2966 6 0 33 19 56 83 71 36 3 1 1 2 1 2 1938145440000 2 20.966175 2966 7 0 50 47 75 97 84 41 3 1 2 2 3 2 1938333360000 5 9.584315 2880 8 0 38 15 60 82 78 45 3 1 1 1 1 2 1937282760000 4 74.928185 2976 9 0 52 15 88 82 118 64 3 2 1 1 1 1 1937282760000 6 140.86449 2976 10 0 16 46 75 91 85 46 1 2 2 1 1 2 1937055420000 6 30.5924 3066 11 0 19 35 63 73 89 36 1 2 3 1 1 1 1937055420000 6 28.01949 3066 12 0 23 99 101 108 96 68 3 2 4 2 2 2 1937058480000 5 23.096144 3028 13 0 48 25 58 87 71 39 1 2 1 1 2 1 1937107560000 3 11.253942 2976 14 0 17 51 78 98 81 48 1 2 2 1 1 1 1937055420000 5 13.5291 3066 15 0 41 64 77 103 79 57 3 2 3 1 1 3 1937107560000 2 9.966941 2977 16 0 39 39 68 94 75 52 3 1 1 1 1 1 1937107560000 3 8.941095 2976 17 0 46 29 58 78 79 38 3 1 1 1 1 1 1937063640000 5 21.72209 2929 18 0 34 31 50 89 59 32 1 2 1 2 2 1 1.9377129e+12 5 5 3095 19 0 36 45 72 99 76 53 1 1 1 1 1 1 1938524580000 1 12.911551 2522 20 0 54 35 89 87 110 55 3 1 3 2 2 2 1938145440000 2 12.897917 2966 21 0 44 40 55 99 59 57 2 2 1 2 2 1 1937715360000 4 5 2967 22 0 40 20 57 72 83 38 3 1 1 1 1 1 1937715360000 4 6.77995 2967 23 0 19 71 95 102 95 69 1 2 1 1 1 2 1937715360000 4 11.296865 2967 24 0 17 27 61 79 80 40 3 2 3 1 1 2 1937904360000 2 75.117355 2767 25 0 26 38 66 73 93 32 1 1 1 1 1 1 1937908440000 6 20.14554 2777 26 0 41 20 39 46 88 17 3 1 1 1 1 2 1937904360000 1 290.9614 2845 27 0 23 21 51 81 65 37 2 2 3 1 1 2 1938333360000 5 12.80904 2880 28 0 27 55 83 97 89 63 2 1 1 2 2 1 1938318720000 6 10.407497 3124 29 0 26 13 33 57 59 18 2 1 2 1 1 1 1938318720000 4 6.911627 2917 30 0 34 18 50 79 66 36 3 1 2 1 1 1 1938319920000 1 54.62834 3117 31 0 22 45 71 88 84 39 3 2 3 1 1 2 1938318720000 6 22.844675 3137 32 0 38 29 68 83 86 41 1 1 1 1 1 1 1.9385049e+12 1 312.9952 4334 33 0 21 83 61 114 55 72 1 2 2 1 1 2 1929584880000 1 27.03701 3344 34 0 25 26 48 84 58 29 1 2 3 6 1 1 1.9297992e+12 6 29.351015 3247 35 0 30 39 50 112 46 32 3 1 1 1 1 2 1929799140000 1 26.44337 3458 36 0 35 61 78 98 84 54 1 2 3 2 1 3 1929584880000 3 5 3344 37 0 18 31 48 87 56 26 3 2 3 1 2 2 1930190040000 5 19.828266 3543 38 0 30 9 53 55 101 30 3 2 3 2 1 2 1933307760000 2 27.137356 3147 39 0 42 54 59 111 56 42 3 1 2 2 2 3 1.9314006e+12 2 5 3328 40 0 41 33 63 91 73 42 3 1 2 1 1 2 1930190040000 4 9.233819 3533 41 0 46 6 38 41 99 21 2 1 2 6 1 2 1.9314e+12 1 5 3339 42 0 18 42 64 92 71 46 3 2 2 1 1 1 1.9297992e+12 5 47.11475 3456 43 0 40 20 39 74 56 21 2 1 1 1 1 1 1.9314e+12 3 5 3360 44 0 29 8 36 41 92 22 2 1 3 1 3 3 1.9322733e+12 5 10.246486 3184 45 0 40 38 70 91 81 51 3 2 2 6 2 2 1930190820000 4 11.249393 3505 46 0 28 8 24 36 68 13 3 1 2 2 2 2 1931667180000 4 8.182312 3113 47 0 21 30 67 88 78 56 1 1 1 6 2 2 1931665380000 4 15.827188 3143 48 0 31 71 89 101 92 86 3 1 2 1 1 2 1.9314e+12 2 25.000776 3633 49 0 24 30 53 84 66 33 3 2 4 6 1 1 1931666340000 5 38.926 3141 50 0 27 21 38 89 45 30 2 1 3 2 1 2 1931402040000 6 12.05554 3598 51 0 28 34 62 88 74 46 3 2 3 6 1 1 1931402040000 6 13.79056 3328 52 0 31 5 21 36 62 14 3 2 3 1 1 2 1931400060000 6 11.596988 3338 53 0 33 40 87 88 101 81 2 1 3 1 3 3 1.9314006e+12 5 12.614867 3351 54 0 27 5 45 59 79 31 3 2 3 1 3 1 1.9314006e+12 4 5 3349 55 0 45 3 15 44 36 10 3 2 3 6 1 2 1931667180000 4 11.661978 3127 56 0 36 41 61 95 68 33 1 2 2 1 1 2 1931665380000 5 7.040653 3143 57 0 18 62 69 103 69 59 3 2 4 8 3 2 1.932618e+12 2 16.805027 3057 58 0 27 44 67 100 69 56 3 2 2 6 1 2 1932037020000 2 17.741343 3079 59 0 22 45 92 84 112 56 2 2 1 6 2 2 1932037020000 3 13.428303 3148 60 0 26 38 71 78 94 40 3 2 2 2 2 2 1931665380000 6 11.860586 857 61 0 34 63 75 103 76 57 2 1 1 6 1 2 1933304760000 3 28.274454 3163 62 0 28 54 63 108 61 40 2 1 3 6 2 1 1.932618e+12 3 8.281861 3057 63 0 32 43 81 92 92 54 3 1 4 1 1 2 1932042540000 1 5 3056 64 0 42 54 73 106 74 41 1 1 1 1 1 2 1932037020000 2 17.38264 3078 65 0 35 95 100 102 103 89 3 2 1 1 1 1 1.9322745e+12 3 48.21596 3178 66 0 30 76 76 102 77 62 3 2 3 1 1 1 1.9322733e+12 6 5.907212 3198 67 0 16 30 55 90 63 51 2 2 2 1 1 2 1.9322745e+12 6 64.7136 3164 68 0 27 37 70 92 79 47 1 2 1 6 1 1 1.932618e+12 4 63.99439 2992 69 0 37 21 52 84 66 46 2 1 2 6 2 2 1.932618e+12 3 8.0207615 2992 70 0 36 22 39 54 76 16 1 2 3 1 1 2 1933306260000 2 34.758698 3157 71 0 41 42 42 119 38 31 3 1 2 1 1 2 1933306980000 5 29.468744 3146 72 0 30 39 61 96 67 50 3 1 1 2 2 3 1934687820000 4 7.678277 3027 73 0 48 49 80 100 85 58 3 1 2 2 2 3 1933824840000 1 13.40274 3206 74 0 42 6 19 40 50 12 3 2 2 2 1 2 1.9338267e+12 1 5 3174 75 0 46 45 65 99 70 58 2 1 2 8 2 3 1.9338249e+12 4 38.80629 3585 76 0 23 59 84 98 89 61 1 2 3 8 1 2 1934687880000 3 7.925013 3072 77 0 31 22 65 78 88 39 3 2 2 2 2 2 1934425440000 5 5 3219 78 0 40 101 94 109 90 83 3 1 1 1 1 3 1934428020000 2 8.3684845 3144 79 0 33 25 66 86 80 60 1 1 1 1 1 3 1934425440000 3 5 3134 80 0 40 53 97 92 111 46 2 1 1 1 1 1 1934425440000 5 6.295738 3134 81 0 55 28 53 96 59 41 3 1 1 1 1 3 1934425380000 5 10.250114 3281 82 0 18 16 42 62 69 31 3 2 3 1 1 3 1934428020000 2 16.560717 3238 83 0 29 28 44 83 56 22 1 1 1 1 1 2 1934433540000 1 5 3052 84 0 35 44 44 117 39 36 1 1 1 1 1 3 1934433540000 2 12.050518 3052 85 0 26 18 37 74 53 20 1 1 1 6 1 3 1.9346889e+12 2 15.355927 3055 86 0 33 16 76 68 119 62 2 1 2 1 1 3 1.9346889e+12 1 30.472855 3009 87 0 33 35 56 96 61 59 1 1 4 2 1 2 1935730440000 1 15.674023 3165 88 0 28 37 64 91 73 35 2 2 1 1 1 2 1934694420000 4 17.564169 3179 89 0 19 43 69 88 81 42 3 2 2 6 2 1 1.9346889e+12 1 988.428 3055 90 0 25 45 68 94 75 52 3 2 3 6 1 2 1.935039e+12 4 16.089264 3071 91 0 31 52 81 96 88 69 3 2 3 2 1 3 1.935039e+12 3 13.08408 3127 92 0 44 19 42 78 57 19 3 2 2 2 3 2 1935736140000 2 5 3191 93 0 43 60 75 103 78 41 3 1 3 1 1 2 1935730440000 1 20.61944 3165 94 0 48 33 59 97 65 45 3 1 1 1 1 2 1.9357305e+12 5 24.95392 3164 95 0 52 21 42 91 49 35 3 1 2 2 1 2 1935239520000 6 5 3103 96 0 37 52 73 101 76 54 1 2 3 1 1 1 1.9357305e+12 6 10.727518 3164 97 0 26 47 56 105 55 43 2 2 3 1 1 2 1935736140000 6 7.124487 3191 98 1 28 34 68 84 84 37 2 1 2 2 2 1 1.9360746e+12 6 25.61064 2924 99 1 32 53 89 95 97 51 1 1 2 1 1 1 1.9360746e+12 3 26.497166 3591 100 1 26 61 86 97 92 60 3 2 2 2 2 1 1936243980000 4 25.342823 3134 end format %tcDD/NN/CCYY_HH:MM datetime label values id id label def id 1 "CK0001", modify label def id 2 "CK0002", modify label def id 3 "CK0003", modify label def id 4 "CK0004", modify label def id 5 "CK0005", modify label def id 6 "CK0006", modify label def id 7 "CK0007", modify label def id 8 "CK0008", modify label def id 9 "CK0009", modify label def id 10 "CK0010", modify label def id 11 "CK0011", modify label def id 12 "CK0012", modify label def id 13 "CK0013", modify label def id 14 "CK0014", modify label def id 15 "CK0015", modify label def id 16 "CK0016", modify label def id 17 "CK0017", modify label def id 18 "CK0018", modify label def id 19 "CK0019", modify label def id 20 "CK0020", modify label def id 21 "CK0021", modify label def id 22 "CK0022", modify label def id 23 "CK0023", modify label def id 24 "CK0024", modify label def id 25 "CK0025", modify label def id 26 "CK0026", modify label def id 27 "CK0027", modify label def id 28 "CK0028", modify label def id 29 "CK0029", modify label def id 30 "CK0030", modify label def id 31 "CK0031", modify label def id 32 "CK0032", modify label def id 33 "CT0001", modify label def id 34 "CT0002", modify label def id 35 "CT0003", modify label def id 36 "CT0004", modify label def id 37 "CT0005", modify label def id 38 "CT0006", modify label def id 39 "CT0007", modify label def id 40 "CT0008", modify label def id 41 "CT0009", modify label def id 42 "CT0010", modify label def id 43 "CT0011", modify label def id 44 "CT0012", modify label def id 45 "CT0013", modify label def id 46 "CT0014", modify label def id 47 "CT0015", modify label def id 48 "CT0016", modify label def id 49 "CT0017", modify label def id 50 "CT0018", modify label def id 51 "CT0019", modify label def id 52 "CT0020", modify label def id 53 "CT0021", modify label def id 54 "CT0022", modify label def id 55 "CT0023", modify label def id 56 "CT0024", modify label def id 57 "CT0025", modify label def id 58 "CT0026", modify label def id 59 "CT0027", modify label def id 60 "CT0028", modify label def id 61 "CT0029", modify label def id 62 "CT0030", modify label def id 63 "CT0031", modify label def id 64 "CT0032", modify label def id 65 "CT0033", modify label def id 66 "CT0034", modify label def id 67 "CT0035", modify label def id 68 "CT0036", modify label def id 69 "CT0037", modify label def id 70 "CT0038", modify label def id 71 "CT0039", modify label def id 72 "CT0040", modify label def id 73 "CT0041", modify label def id 74 "CT0042", modify label def id 75 "CT0043", modify label def id 76 "CT0044", modify label def id 77 "CT0045", modify label def id 78 "CT0046", modify label def id 79 "CT0047", modify label def id 80 "CT0048", modify label def id 81 "CT0049", modify label def id 82 "CT0050", modify label def id 83 "CT0051", modify label def id 84 "CT0052", modify label def id 85 "CT0053", modify label def id 86 "CT0054", modify label def id 87 "CT0055", modify label def id 88 "CT0056", modify label def id 89 "CT0057", modify label def id 90 "CT0058", modify label def id 91 "CT0059", modify label def id 92 "CT0060", modify label def id 93 "CT0061", modify label def id 94 "CT0062", modify label def id 95 "CT0063", modify label def id 96 "CT0064", modify label def id 97 "CT0065", modify label def id 98 "GWK0001", modify label def id 99 "GWK0002", modify label def id 100 "GWK0003", modify label values type type_label label def type_label 0 "Control", modify label def type_label 1 "Worker", modify label values YearsDis YearsDis label def YearsDis 1 "<5 years", modify label def YearsDis 2 "6-10 years", modify label def YearsDis 3 "Above 10 years", modify label values MaritalStat MaritalStat label def MaritalStat 1 "Married", modify label def MaritalStat 2 "Single", modify label values Edu Edu label def Edu 1 "None", modify label def Edu 2 "Primary School", modify label def Edu 3 "Middle/JHS", modify label def Edu 4 "Post JHS", modify label values HHstove HHstove label def HHstove 1 "Open mokyia (3-stone/mud/clay)", modify label def HHstove 2 "Closed mokyia (metal/clay)", modify label values cookspace cookspace label def cookspace 1 "Totally open outdoors", modify label def cookspace 2 "Partially enclosed", modify label def cookspace 3 "Fully enclosed", modify label values HHcooking HHcooking label def HHcooking 1 "1-hour", modify label def HHcooking 2 "2-hours", modify label def HHcooking 3 "3-hours", modify
