CI 95% Graph mixed model prediction

[Augusto Sao]

I can't put IC95 on this graph. Can you help me, please? Thank you very much

twoway connected ava temps if grupos_ava_new==1, c(L) ms(i) ytitle(" ") color("242 183 153") ///
|| connected ava temps if grupos_ava_new==2, c(L) ms(i) ytitle(" ") color("242 183 153") ///
||connected ava temps if grupos_ava_new==3, c(L) ms(i) ytitle(" ") color("242 183 153") ///
|| line predinter temps if grupos_ava_new==1 , ytitle(" ") color(red) ///
|| line predinter temps if grupos_ava_new==2 , ytitle(" ") color(blue) ///
|| line predinter temps if grupos_ava_new==3 , ytitle(" ") color(green) ///
||, title (Evolución de la predicción con interacción ) xsc(range(0,100)) xlabel(0(10)100)



----------------------- copy starting from the next line -----------------------

* Example generated by -dataex-. To install: ssc install dataex
clear
input float id byte time float(grupos_ava_new temps predinter)
  6 1 2         0 .8922561
  6 2 2  44.89256 .7242961
  7 1 2         0 .8922561
  7 2 2  49.95041 .7053728
 11 1 2         0 .8922561
 11 2 2 15.438017 .8344966
 13 1 2         0 .8922561
 13 2 2 15.768595 .8332598
 14 1 2         0 .8922561
 14 2 2 33.421486 .7672137
 17 1 2         0 .8922561
 17 2 2 17.818182 .8255915
 18 1 2         0 .8922561
 18 2 2  24.46281 .8007315
 19 1 2         0 .8922561
 19 2 2  75.96694 .6080351
 20 1 2         0 .8922561
 20 2 2  77.95042 .6006142
 21 1 2         0 .8922561
 21 2 2 17.057852 .8284362
 23 1 2         0 .8922561
 23 2 2 108.26447 .4871979
 25 1 2         0 .8922561
 25 2 2 18.049587 .8247257
 29 1 3         0 .7664383
 29 2 3 13.785124 .7415267
 30 1 1         0 .7521105
 30 2 1  49.32232 .6182449
 31 1 1         0 .7521105
 31 2 1  64.79339 .5762549
 32 1 3         0 .7664383
 32 2 3 24.264463 .7225891
 33 1 2         0 .8922561
 33 2 2 29.652893 .7813134
 36 1 1         0 .7521105
 36 2 1 32.363636 .6642724
 37 1 2         0 .8922561
 37 2 2  48.33058 .7114332
 39 1 2         0 .8922561
 39 2 2 17.785124 .8257152
 40 1 3         0 .7664383
 40 2 3  50.64463 .6749166
 41 1 2         0 .8922561
 41 2 2 15.438017 .8344966
 42 1 2         0 .8922561
 42 2 2  25.68595 .7961553
 44 1 2         0 .8922561
 44 2 2 35.471073 .7595454
 46 1 1         0 .7521105
 46 2 1  37.38843 .6506346
 48 1 2         0 .8922561
 48 2 2 31.768595 .7733978
 49 1 2         0 .8922561
 49 2 2  46.34711 .7188541
 51 1 2         0 .8922561
 51 2 2 28.099174 .7871265
 52 1 2         0 .8922561
 52 2 2  37.05785 .7536087
 54 1 1         0 .7521105
 54 2 1 13.024794 .7167599
 55 1 2         0 .8922561
 55 2 2 12.264462 .8463701
 57 1 2         0 .8922561
 57 2 2  15.83471 .8330125
 60 1 2         0 .8922561
 60 2 2 32.859505 .7693163
 62 1 2         0 .8922561
 62 2 2 30.710743 .7773556
 63 1 2         0 .8922561
 63 2 2  56.76033 .6798943
 64 1 1         0 .7521105
 64 2 1 14.247934 .7134402
 65 1 1         0 .7521105
 65 2 1 32.694214 .6633752
 71 1 1         0 .7521105
 71 2 1  41.05785 .6406755
 72 1 1         0 .7521105
 72 2 1  49.19008 .6186038
 74 1 1         0 .7521105
 74 2 1  24.03306 .6868824
 76 1 3         0 .7664383
 76 2 3 19.438017 .7313112
 82 1 1         0 .7521105
 82 2 1 32.727272 .6632855
 86 1 2         0 .8922561
 86 2 2  24.33058 .8012262
 90 1 3         0 .7664383
 90 2 3  27.50413 .7167346
 92 1 3         0 .7664383
 92 2 3  44.06612 .6868049
 95 1 1         0 .7521105
 95 2 1 15.735538 .7094027
 96 1 1         0 .7521105
 96 2 1  14.31405 .7132608
 97 1 2         0 .8922561
 97 2 2 37.619835 .7515061
 99 1 3         0 .7664383
 99 2 3  85.12397 .6126078
100 1 1         0 .7521105
100 2 1  31.23967  .667323
end

[Andrew Musau]

Your data example lacks the variable "ava". You should present a dataex of the data just prior to graphing (in case you change the data in the intermediate), including the missing variable.

[Augusto Sao]

Hi Andrés, sorry. I am sending you the data and my script. Thank you very much

* Example generated by -dataex-. To install: ssc install dataex
clear
input float id byte time double ava float grupos_ava_new byte(grupos_ava_new_1 grupos_ava_new_2 grupos_ava_new_3) float temps
  6 1   .608291504740883 2 0 1 0         0
  6 2   .608291504740883 2 0 1 0  44.89256
  7 1  .9640962872164949 2 0 1 0         0
  7 2  .8847874817675939 2 0 1 0  49.95041
 11 1  .9347099274725276 2 0 1 0         0
 11 2         .593760888 2 0 1 0 15.438017
 13 1  .9379896114285715 2 0 1 0         0
 13 2          .69115024 2 0 1 0 15.768595
 14 1  .9823652287640451 2 0 1 0         0
 14 2   .733529154625789 2 0 1 0 33.421486
 17 1  .8497082362352942 2 0 1 0         0
 17 2        .8263172358 2 0 1 0 17.818182
 18 1  .7759268039084843 2 0 1 0         0
 18 2  .5890944297076025 2 0 1 0  24.46281
 19 1   .803663069767442 2 0 1 0         0
 19 2 .36863203689578716 2 0 1 0  75.96694
 20 1 1.0083939493670886 2 0 1 0         0
 20 2  .7027977145892352 2 0 1 0  77.95042
 21 1  .7441464057672504 2 0 1 0         0
 21 2  .6256894527196653 2 0 1 0 17.057852
 23 1  .8975977142857144 2 0 1 0         0
 23 2  .5031995456396396 2 0 1 0 108.26447
 25 1  .9626872314634146 2 0 1 0         0
 25 2  .7758475603200002 2 0 1 0 18.049587
 29 1  .7475865419532325 3 0 0 1         0
 29 2  .5236507660364842 3 0 0 1 13.785124
 30 1  .9729205994005995 1 1 0 0         0
 30 2  .5971402792262406 1 1 0 0  49.32232
 31 1  .9252739575671141 1 1 0 0         0
 31 2  .6850608105726873 1 1 0 0  64.79339
 32 1  .8740445673667714 3 0 0 1         0
 32 2  .6314419713957935 3 0 0 1 24.264463
 33 1  .7048584781549817 2 0 1 0         0
 33 2  .8147887196776795 2 0 1 0 29.652893
 36 1  .6816420061480076 1 1 0 0         0
 36 2  .6215593810671257 1 1 0 0 32.363636
 37 1   .856163253426043 2 0 1 0         0
 37 2  .7775259856027555 2 0 1 0  48.33058
 39 1  .9167635017793595 2 0 1 0         0
 39 2  .6920691700182815 2 0 1 0 17.785124
 40 1  .8302465357941178 3 0 0 1         0
 40 2   .881187195327103 3 0 0 1  50.64463
 41 1  .9072617574730355 2 0 1 0         0
 41 2  .8097242534899328 2 0 1 0 15.438017
 42 1  .7587562767978291 2 0 1 0         0
 42 2  .6069040326562501 2 0 1 0  25.68595
 44 1  .9590712951894272 2 0 1 0         0
 44 2   .815150804096802 2 0 1 0 35.471073
 46 1  .9706954922155688 1 1 0 0         0
 46 2  .5976263179500431 1 1 0 0  37.38843
 48 1  .9785942292193311 2 0 1 0         0
 48 2  .9552726454072792 2 0 1 0 31.768595
 49 1  .8288768073811931 2 0 1 0         0
 49 2  .7522327016241299 2 0 1 0  46.34711
 51 1  .9817475000000001 2 0 1 0         0
 51 2  .5924144914285715 2 0 1 0 28.099174
 52 1  .8550156787200001 2 0 1 0         0
 52 2  .8357190753982301 2 0 1 0  37.05785
 54 1  .8695999716279069 1 1 0 0         0
 54 2  .8555849704579026 1 1 0 0 13.024794
 55 1  .8435253753089534 2 0 1 0         0
 55 2  .6400372687627909 2 0 1 0 12.264462
 57 1  .9009804519402986 2 0 1 0         0
 57 2  .7987385815348101 2 0 1 0  15.83471
 60 1  .8865410795072788 2 0 1 0         0
 60 2   .908294147224457 2 0 1 0 32.859505
 62 1  .8791286255144032 2 0 1 0         0
 62 2  .8755889521487604 2 0 1 0 30.710743
 63 1  .9559690993865031 2 0 1 0         0
 63 2  .3805422892125318 2 0 1 0  56.76033
 64 1                  1 1 1 0 0         0
 64 2  .5264874868965518 1 1 0 0 14.247934
 65 1  .8803710615831135 1 1 0 0         0
 65 2  .7663112426338895 1 1 0 0 32.694214
 71 1  .8669959869026551 1 1 0 0         0
 71 2  .4760065129972376 1 1 0 0  41.05785
 72 1  .9914749064220184 1 1 0 0         0
 72 2  .6418337156979406 1 1 0 0  49.19008
 74 1  .4897313529090909 1 1 0 0         0
 74 2  .4435697328440368 1 1 0 0  24.03306
 76 1 .49576903890410967 3 0 0 1         0
 76 2 .47123880000000007 3 0 0 1 19.438017
 82 1  .6413148669047619 1 1 0 0         0
 82 2  .4440519461538462 1 1 0 0 32.727272
 86 1  .6576163928089889 2 0 1 0         0
 86 2  .3398297961089495 2 0 1 0  24.33058
 90 1  .6693000648717949 3 0 0 1         0
 90 2  .5785765266666667 3 0 0 1  27.50413
 92 1  .7093125687499999 3 0 0 1         0
 92 2  .6908008044491526 3 0 0 1  44.06612
 95 1  .8168139200000001 1 1 0 0         0
 95 2   .478718780952381 1 1 0 0 15.735538
 96 1  .7595625394736844 1 1 0 0         0
 96 2  .7324576022773974 1 1 0 0  14.31405
 97 1  .8943993123037977 2 0 1 0         0
 97 2 .45224769942428034 2 0 1 0 37.619835
 99 1  .5131483520165461 3 0 0 1         0
 99 2 .22917484002668448 3 0 0 1  85.12397
100 1 .21114843750000004 1 1 0 0         0
100 2 .34576603376097237 1 1 0 0  31.23967
end

------------------ copy up to and including the previous line ------------------


gen interAVA1=temps*grupos_ava_new_1
gen interAVA2=temps*grupos_ava_new_2
gen interAVA3=temps*grupos_ava_new_3


xi:xtmixed ava temps grupos_ava_new_2 grupos_ava_new_3 interAVA2 interAVA3 || id:temps, var mle





twoway connected ava temps if grupos_ava_new==1, c(L) ms(i) ytitle(" ") color(red) ///
|| connected ava temps if grupos_ava_new==2, c(L) ms(i) ytitle(" ") color(blue) ///
||connected ava temps if grupos_ava_new==3, c(L) ms(i) ytitle(" ") color(green) ///
|| , title (Evolución Individual )


xi: xtmixed ava temps grupos_ava_new_2 grupos_ava_new_3 || id: ,mle iterate(40)


estimates store mod1


xi: xtmixed ava temps grupos_ava_new_2 grupos_ava_new_3 || id:temps ,mle iterate(40)


estimates store mod2


lrtest mod1 mod2


gen interAVA1=temps*grupos_ava_new_1
gen interAVA2=temps*grupos_ava_new_2
gen interAVA3=temps*grupos_ava_new_3


xi:xtmixed ava temps grupos_ava_new_2 grupos_ava_new_3 interAVA2 interAVA3 || id:temps, var mle


estimates store mod3


lrtest mod2 mod3


xi:xtmixed ava temps grupos_ava_new_2 grupos_ava_new_3 interAVA2 interAVA3 || id:temps, var mle



predict predinter , xb

sort id temps



twoway connected ava temps if grupos_ava_new==1, c(L) ms(i) ytitle(" ") color("242 183 153") ///
|| connected ava temps if grupos_ava_new==2, c(L) ms(i) ytitle(" ") color("242 183 153") ///
||connected ava temps if grupos_ava_new==3, c(L) ms(i) ytitle(" ") color("242 183 153") ///
|| line predinter temps if grupos_ava_new==1 , ytitle(" ") color(red) ///
|| line predinter temps if grupos_ava_new==2 , ytitle(" ") color(blue) ///
|| line predinter temps if grupos_ava_new==3 , ytitle(" ") color(green) ///
||, title (Evolución de la predicción con interacción ) xsc(range(0,100)) xlabel(0(10)100)

[Andrew Musau]

I do not see the relationship between your estimated models and the graphs. You can create twoway linear prediction plots with 95% CIs as below:

* Example generated by -dataex-. To install: ssc install dataex
clear
input float id byte time double ava float grupos_ava_new byte(grupos_ava_new_1 grupos_ava_new_2 grupos_ava_new_3) float temps
  6 1   .608291504740883 2 0 1 0         0
  6 2   .608291504740883 2 0 1 0  44.89256
  7 1  .9640962872164949 2 0 1 0         0
  7 2  .8847874817675939 2 0 1 0  49.95041
 11 1  .9347099274725276 2 0 1 0         0
 11 2         .593760888 2 0 1 0 15.438017
 13 1  .9379896114285715 2 0 1 0         0
 13 2          .69115024 2 0 1 0 15.768595
 14 1  .9823652287640451 2 0 1 0         0
 14 2   .733529154625789 2 0 1 0 33.421486
 17 1  .8497082362352942 2 0 1 0         0
 17 2        .8263172358 2 0 1 0 17.818182
 18 1  .7759268039084843 2 0 1 0         0
 18 2  .5890944297076025 2 0 1 0  24.46281
 19 1   .803663069767442 2 0 1 0         0
 19 2 .36863203689578716 2 0 1 0  75.96694
 20 1 1.0083939493670886 2 0 1 0         0
 20 2  .7027977145892352 2 0 1 0  77.95042
 21 1  .7441464057672504 2 0 1 0         0
 21 2  .6256894527196653 2 0 1 0 17.057852
 23 1  .8975977142857144 2 0 1 0         0
 23 2  .5031995456396396 2 0 1 0 108.26447
 25 1  .9626872314634146 2 0 1 0         0
 25 2  .7758475603200002 2 0 1 0 18.049587
 29 1  .7475865419532325 3 0 0 1         0
 29 2  .5236507660364842 3 0 0 1 13.785124
 30 1  .9729205994005995 1 1 0 0         0
 30 2  .5971402792262406 1 1 0 0  49.32232
 31 1  .9252739575671141 1 1 0 0         0
 31 2  .6850608105726873 1 1 0 0  64.79339
 32 1  .8740445673667714 3 0 0 1         0
 32 2  .6314419713957935 3 0 0 1 24.264463
 33 1  .7048584781549817 2 0 1 0         0
 33 2  .8147887196776795 2 0 1 0 29.652893
 36 1  .6816420061480076 1 1 0 0         0
 36 2  .6215593810671257 1 1 0 0 32.363636
 37 1   .856163253426043 2 0 1 0         0
 37 2  .7775259856027555 2 0 1 0  48.33058
 39 1  .9167635017793595 2 0 1 0         0
 39 2  .6920691700182815 2 0 1 0 17.785124
 40 1  .8302465357941178 3 0 0 1         0
 40 2   .881187195327103 3 0 0 1  50.64463
 41 1  .9072617574730355 2 0 1 0         0
 41 2  .8097242534899328 2 0 1 0 15.438017
 42 1  .7587562767978291 2 0 1 0         0
 42 2  .6069040326562501 2 0 1 0  25.68595
 44 1  .9590712951894272 2 0 1 0         0
 44 2   .815150804096802 2 0 1 0 35.471073
 46 1  .9706954922155688 1 1 0 0         0
 46 2  .5976263179500431 1 1 0 0  37.38843
 48 1  .9785942292193311 2 0 1 0         0
 48 2  .9552726454072792 2 0 1 0 31.768595
 49 1  .8288768073811931 2 0 1 0         0
 49 2  .7522327016241299 2 0 1 0  46.34711
 51 1  .9817475000000001 2 0 1 0         0
 51 2  .5924144914285715 2 0 1 0 28.099174
 52 1  .8550156787200001 2 0 1 0         0
 52 2  .8357190753982301 2 0 1 0  37.05785
 54 1  .8695999716279069 1 1 0 0         0
 54 2  .8555849704579026 1 1 0 0 13.024794
 55 1  .8435253753089534 2 0 1 0         0
 55 2  .6400372687627909 2 0 1 0 12.264462
 57 1  .9009804519402986 2 0 1 0         0
 57 2  .7987385815348101 2 0 1 0  15.83471
 60 1  .8865410795072788 2 0 1 0         0
 60 2   .908294147224457 2 0 1 0 32.859505
 62 1  .8791286255144032 2 0 1 0         0
 62 2  .8755889521487604 2 0 1 0 30.710743
 63 1  .9559690993865031 2 0 1 0         0
 63 2  .3805422892125318 2 0 1 0  56.76033
 64 1                  1 1 1 0 0         0
 64 2  .5264874868965518 1 1 0 0 14.247934
 65 1  .8803710615831135 1 1 0 0         0
 65 2  .7663112426338895 1 1 0 0 32.694214
 71 1  .8669959869026551 1 1 0 0         0
 71 2  .4760065129972376 1 1 0 0  41.05785
 72 1  .9914749064220184 1 1 0 0         0
 72 2  .6418337156979406 1 1 0 0  49.19008
 74 1  .4897313529090909 1 1 0 0         0
 74 2  .4435697328440368 1 1 0 0  24.03306
 76 1 .49576903890410967 3 0 0 1         0
 76 2 .47123880000000007 3 0 0 1 19.438017
 82 1  .6413148669047619 1 1 0 0         0
 82 2  .4440519461538462 1 1 0 0 32.727272
 86 1  .6576163928089889 2 0 1 0         0
 86 2  .3398297961089495 2 0 1 0  24.33058
 90 1  .6693000648717949 3 0 0 1         0
 90 2  .5785765266666667 3 0 0 1  27.50413
 92 1  .7093125687499999 3 0 0 1         0
 92 2  .6908008044491526 3 0 0 1  44.06612
 95 1  .8168139200000001 1 1 0 0         0
 95 2   .478718780952381 1 1 0 0 15.735538
 96 1  .7595625394736844 1 1 0 0         0
 96 2  .7324576022773974 1 1 0 0  14.31405
 97 1  .8943993123037977 2 0 1 0         0
 97 2 .45224769942428034 2 0 1 0 37.619835
 99 1  .5131483520165461 3 0 0 1         0
 99 2 .22917484002668448 3 0 0 1  85.12397
100 1 .21114843750000004 1 1 0 0         0
100 2 .34576603376097237 1 1 0 0  31.23967
end

set scheme s1color
twoway lfitci ava temps if grupos_ava_new==1, clcolor(red) clwidth(medthick)  fcolor(red%30) || ///
lfitci ava temps if grupos_ava_new==2, clcolor(blue) clwidth(medthick) fcolor(blue%30)|| ///
lfitci ava temps if grupos_ava_new==3,  clcolor(green) clwidth(medthick) fcolor(green%30) ||, ///
title (Evolución Individual) legend(order(1 2 "Predicted (grupos 1)" 3 4 "Predicted (grupos 2)" 5 6 "Predicted  (grupos 3)"))

[Augusto Sao]

Thank you very much, Andrés, for your help; it has benefited me. Thanks again.