Good morning, I am trying to run a competing risks model, and Stata enters a loop with the following message: “Iteration 0: log pseudolikelihood = -94.22336 (not concave)”
Thanks!!
Code:
* Example generated by -dataex-. To install: ssc install dataex clear input str8 id byte sexo float(t_ev_combi ev_competi) "1619-290" 0 1.9329226 0 "1619-289" 1 2.0780287 0 "1619-310" 1 2.387406 0 "1619-311" 0 3.33744 0 "1619-312" 1 4.235455 0 "1321-54" 0 4.2956877 0 "1319-95" 0 4.87885 0 "1619-316" 1 5.234771 0 "1619-318" 0 6.028748 0 "1619-337" 1 6.036961 0 "1319-124" 0 6.110883 0 "1319-122" 0 7.597536 0 "1619-322" 0 8.197125 0 "1319-139" 1 9.28679 0 "1619-363" 1 9.371663 0 "1319-39" 0 9.568789 0 "1619-364" 1 9.837098 0 "1319-28" 1 10.01232 0 "1319-133" 1 10.80356 0 "1619-378" 1 11.143053 0 "1619-376" 1 11.288158 0 "1619-377" 1 11.35113 0 "1319-38" 0 11.36482 0 "1319-131" 1 11.920603 0 "1319-58" 0 12.303902 0 "1619-391" 0 12.36961 0 "1319-27" 0 13.483915 0 "1319-125" 1 13.50308 0 "1319-63" 0 14.439425 0 "1317-73" 0 15.334702 0 "1619-243" 1 15.340178 0 "1365-46" 1 15.50445 0 "1319-26" 0 15.570157 0 "1321-57" 0 15.709787 0 "1317-74" 0 15.854894 0 "1317-81" 1 15.86037 0 "1365-45" 0 15.95072 0 "1317-95" 0 16.095825 0 "1319-127" 1 16.109514 0 "1365-68" 0 16.172485 0 "1319-130" 0 16.49829 0 "1365-62" 0 16.53662 0 "1365-47" 1 16.629705 0 "1319-64" 1 16.673512 0 "1365-60" 1 16.758385 0 "1321-70" 1 16.837782 0 "1317-67" 1 17.100616 0 "1317-30" 0 17.779604 0 "1619-237" 1 17.979465 0 "1619-244" 0 18.168377 0 "1365-28" 1 18.236824 0 "1317-31" 0 18.35729 0 "1317-32" 0 18.710472 0 "1619-325" 0 18.932238 0 "1619-239" 0 18.970568 0 "1365-33" 0 19.071869 0 "1317-27" 0 19.318275 0 "1317-21" 0 19.40315 0 "1317-25" 0 19.44969 0 "1619-106" 1 19.4935 0 "1619-172" 1 19.65503 0 "1365-26" 1 19.83299 0 "1319-24" 0 20.06297 0 "1365-36" 0 20.093086 0 "1619-180" 0 20.238194 0 "1619-236" 1 20.528406 0 "1318-28" 1 20.55031 0 "1619-114" 1 20.87885 0 "1321-36" 0 21.065023 0 "1619-238" 0 21.49213 0 "1619-175" 0 21.549623 0 "1645-55" 0 21.61807 0 "1645-42" 1 21.98494 0 "1318-22" 1 22.058865 0 "1365-40" 1 22.069815 0 "1619-105" 1 22.146475 0 "1319-36" 1 22.16564 0 "1619-166" 0 22.349075 0 "1619-185" 0 22.450377 0 "1319-66" 0 22.71595 0 "1365-38" 0 23.255304 0 "1645-38" 1 23.27173 0 "1619-159" 1 23.28268 0 "1619-110" 0 23.54004 0 "1645-51" 1 23.627653 0 "1365-32" 0 23.737165 0 "1619-146" 1 23.84668 0 "1619-128" 1 23.917864 0 "1319-135" 0 23.93703 0 "1619-158" 0 24.323067 0 "1619-153" 1 24.372347 0 "1645-34" 1 24.807667 0 "1317-80" 1 24.974674 0 "1619-147" 0 25.13073 0 "1317-50" 0 25.141684 0 "1619-171" 0 25.330595 0 "1619-176" 0 25.33607 0 "1619-266" 0 25.462013 0 "1619-294" 1 26.05339 0 "1321-39" 0 26.209446 0 "1317-92" 1 26.54346 0 "1645-37" 0 26.67214 0 "1619-269" 1 26.759754 0 "1619-252" 1 26.85558 0 "1619-162" 0 26.910336 0 "1619-186" 0 27.5154 0 "1317-53" 1 27.701574 0 "1319-60" 0 27.917864 0 "1317-85" 0 28.17522 0 "1645-44" 0 28.399727 0 "1619-303" 0 28.58864 0 "1317-77" 0 28.71184 0 "1319-23" 0 28.71732 0 "1619-291" 0 29.42916 0 "1645-41" 0 29.749487 0 "1619-178" 0 29.92197 0 "1619-299" 0 29.954826 0 "1317-58" 0 30.10541 0 "1321-41" 1 30.89117 0 "1619-111" 1 31.12115 0 "1319-128" 1 31.18412 0 "1619-112" 0 31.28816 0 "1317-24" 1 31.43874 0 "1619-136" 0 31.52635 0 "1619-154" 1 31.676933 0 "1619-168" 0 32.060234 0 "1619-148" 1 32.109512 0 "1317-69" 1 32.29569 0 "1318-16" 1 32.457222 0 "1619-384" 0 32.687202 0 "1319-74" 0 32.821354 0 "1645-35" 0 32.960987 0 "1318-21" 0 33.52772 0 "1365-39" 1 33.535934 0 "1317-89" 1 34.12731 0 "1619-282" 1 34.351814 0 "1619-109" 0 34.707733 0 "1319-129" 1 34.836414 0 "1619-271" 0 34.91855 0 "4667-5" 0 35.268993 0 "1619-188" 1 35.44148 0 "1317-83" 0 35.775497 0 "1319-144" 0 36.01369 0 "1319-123" 1 36.095825 0 "1619-115" 1 36.566734 0 "1619-119" 1 37.1718 0 "1365-77" 1 37.18275 0 "1619-163" 1 37.546886 0 "1319-30" 0 37.59343 0 "1619-249" 0 37.73032 0 "1319-69" 1 37.77139 0 "1317-64" 1 38.07803 0 "1317-84" 1 38.28337 0 "1318-26" 1 38.30801 0 "1619-313" 1 38.41205 0 "1619-354" 1 38.56536 0 "1619-274" 1 38.70226 0 "1619-256" 1 38.86653 0 "1619-352" 0 38.88296 0 "1619-160" 1 39.06639 0 "1365-49" 1 39.12663 0 "1319-138" 1 40.13689 0 "1319-136" 1 40.38877 0 "1619-253" 0 40.47365 0 "1619-360" 0 40.58316 0 "1319-73" 0 40.80219 0 "1619-117" 1 40.81041 0 "1619-118" 0 41.31691 0 "1365-48" 0 41.32786 0 "1619-326" 0 41.61533 0 "1619-116" 0 41.88912 0 "1317-34" 1 42.06434 0 "1365-59" 1 42.25873 0 "1321-56" 1 42.37098 0 "1619-277" 1 42.37098 0 "1619-143" 1 42.69678 0 "1619-273" 0 43.18412 0 "1619-240" 1 43.23888 0 "1619-284" 0 43.67693 0 "1619-113" 0 43.70431 0 "1619-169" 1 43.76181 0 "1319-119" 1 43.84121 0 "1319-101" 1 44.25188 0 "1619-247" 0 44.75291 0 "1619-122" 1 44.96646 0 "1318-25" 1 45.1143 0 "1619-355" 0 45.24572 0 "1319-141" 1 45.82615 0 "1619-123" 1 45.90554 0 "1365-31" 1 45.96578 0 "1619-179" 1 46.10267 0 "1619-353" 0 46.15195 0 "1619-174" 1 46.16564 0 "1321-48" 1 46.2423 0 "1619-357" 1 46.69678 0 "1619-107" 1 46.95688 0 "1317-88" 1 47.5154 0 "1365-37" 0 47.55373 0 "1321-45" 0 47.58659 0 "2993-3" 0 47.61944 0 "1619-138" 0 47.66325 0 "1619-165" 0 47.71526 0 "1319-72" 0 47.87954 0 "1321-59" 0 48.26831 0 "1619-141" 1 48.27652 0 "1619-374" 0 48.36961 0 "1321-55" 1 48.38056 0 "1619-124" 1 48.71458 0 "1619-125" 1 48.75291 0 "1365-75" 1 48.84873 0 "1619-369" 0 49.0267 0 "1619-151" 1 49.27858 0 "1318-10" 0 49.53046 0 "1619-182" 1 49.9603 0 "1319-115" 1 50.14374 0 "1619-301" 1 50.24778 0 "1319-126" 0 50.25599 0 "1317-22" 1 50.51609 0 "1319-140" 1 50.75428 0 "1317-17" 1 50.94593 0 "1317-54" 1 50.97878 0 "1619-127" 1 51.08282 0 "1619-126" 1 51.35387 0 "1619-358" 1 51.57837 0 "1619-373" 1 51.75359 0 "1619-120" 0 51.88227 0 "1619-261" 1 51.8987 0 "1619-156" 1 52.03012 0 "1619-133" 1 52.89528 0 "1619-144" 1 53.10336 0 "1317-75" 1 53.32786 0 "1321-43" 1 53.3744 0 "1619-300" 1 53.39631 0 "1619-272" 1 53.60986 0 "1365-53" 1 53.83162 0 "1619-258" 0 54.01232 0 "1318-24" 0 54.7269 0 "1619-229" 0 55.63586 0 "1319-132" 0 55.81109 0 "1619-267" 1 55.96167 0 "1619-121" 1 56.63244 0 "1619-250" 1 56.68446 0 "1645-49" 1 56.77208 0 "1619-292" 1 57.30869 0 "1619-152" 1 57.80151 0 "1619-350" 0 58.5681 0 "1365-72" 1 58.80082 0 "1619-142" 1 59.436 0 "1619-328" 1 59.83847 0 "1317-78" 1 59.8768 0 "1317-23" 0 60.20808 0 "1321-51" 1 60.40794 0 "1365-44" 0 60.55305 0 "2993-6" 1 60.80219 0 "1321-60" 1 60.84326 0 "1619-264" 1 61.51951 0 "1619-263" 1 62.10815 0 "1619-259" 1 62.56536 0 "1319-25" 0 63.70979 0 "1619-145" 1 63.93977 0 "1619-286" 0 64.10951 0 "1619-155" 1 64.402466 0 "1319-57" 1 64.94182 0 "1619-161" 0 66.13005 0 "1619-351" 0 66.231346 0 "1619-298" 1 66.97878 0 "1619-260" 1 67.10746 0 "1619-324" 1 67.52909 0 "1619-365" 1 68.67625 0 "1619-356" 0 68.69815 0 "1321-47" 1 71.7974 0 "1619-359" 1 71.99452 0 "1619-132" 1 73.48392 0 "1319-145" 1 73.801506 0 "1317-82" 0 75.46064 0 "1619-139" 0 81.72211 0 "1317-14" 0 12.5 3 "1318-15" 0 12.8 2 "1318-23" 0 14.4 2 "1317-57" 0 15.3 2 "1317-26" 0 17.4 2 "1619-283" 0 17.5 2 "1619-287" 0 17.7 2 "1318-14" 0 18.5 2 "1619-257" 0 18.338125 0 "1317-66" 0 19.8 3 "1619-85" 0 20 3 "1619-98" 1 20.4 2 "1321-28" 0 20.6 2 "1619-346" 1 21 2 "1321-50" 0 21.7 3 "1619-305" 0 21.8 2 "1619-225" 0 22.1 2 "1619-278" 0 22.1 2 "1317-72" 1 22.4 3 "1321-44" 1 22.7 2 "1619-309" 1 22.8 2 "1319-53" 0 22.8 3 "1675-27" 0 23 2 "1317-71" 0 23 3 "1619-293" 0 23 2 "1317-15" 0 23.1 2 "1365-29" 0 23.6 2 "1619-297" 0 24.8 2 "1321-30" 0 25.2 3 "1317-33" 1 25.7 2 "1619-280" 0 25.4538 0 "1319-65" 0 25.9 4 "1317-65" 0 26 3 "1619-327" 0 26.5 2 "1319-37" 0 26.5 2 "2994-2" 0 26.7 3 "1321-38" 0 26.9 2 "1319-35" 0 26.896646 0 "1619-75" 0 27 3 "1619-88" 1 27.4 2 "1321-37" 0 27.4 3 "1619-349" 0 27.9 1 "1317-87" 1 28 2 "1619-285" 0 28.3 2 "1619-223" 1 28.6 2 "1321-33" 0 28.8 3 "1619-219" 0 29.4 2 "1317-18" 1 29.4 2 "1319-137" 1 29.7 3 "1321-46" 0 29.8 2 "1317-19" 1 29.9 4 "1619-95" 0 30.5 3 "1619-103" 1 30.7 4 "2993-4" 0 30.9 2 "1319-103" 0 30.9 2 "1319-102" 0 31.2 4 "1619-296" 1 31.2 2 "1619-380" 1 31.3 2 "1619-315" 0 31.4 2 "1319-67" 1 31.7 4 "1317-59" 1 31.7 2 "1619-368" 1 31.7 2 "1317-70" 1 31.8 2 "1619-92" 0 31.8 3 "1619-295" 1 32 2 "1319-70" 1 32.2 2 "2994-1" 0 32.2 3 "1317-51" 0 32.5 2 "1619-233" 0 33 2 "1317-49" 1 33.2 2 "1321-49" 0 33.2 4 "1619-276" 0 33.3 2 "1619-254" 0 33.3 3 "1317-37" 0 33.7 2 "1619-340" 0 33.7 2 "1619-255" 1 33.7 2 "1619-383" 1 29.8 2 "1619-212" 1 33.8 2 "1619-246" 1 33.8 2 "1619-279" 0 34.2 2 "1675-22" 0 34.2 2 "1318-11" 0 34.4 2 "1365-25" 0 34.5 2 "1319-34" 1 34.7 2 "1619-99" 1 35 3 "1619-268" 0 35.1 2 "1619-77" 1 35.1 3 "1321-53" 0 35.2 2 "1619-366" 1 34.4 2 "1619-304" 0 35.7 2 "1619-220" 0 35.8 2 "1619-81" 0 35.8 3 "1319-59" 0 35.9 2 "1675-28" 0 35.9 2 "1319-29" 0 36.18891 0 "1675-24" 1 36.2 3 "1619-275" 1 36.309376 0 "1619-336" 1 36.4 2 "1317-13" 0 36.6 3 "1619-314" 0 36.7 2 "1317-35" 0 36.9 2 "1645-58" 1 36.462696 0 "1619-262" 0 37.4 2 "1619-224" 0 37.5 2 "1675-29" 1 37.9 2 "1619-101" 1 37.9 2 "1619-215" 1 38.5 4 "1619-94" 0 38.7 3 "1318-13" 0 38.7 2 "1645-62" 0 39.2 2 "1675-26" 1 39.2 3 "1365-64" 0 39.6 2 "1321-35" 0 40 2 "1619-78" 1 40 3 "1619-370" 0 36.5 2 "1619-248" 0 40.5 3 "1645-43" 0 40.8 2 "1619-102" 0 40.9 4 "1619-80" 1 41.1 4 "1317-28" 0 41.1 2 "1619-76" 0 41.1 3 "1319-68" 0 41.26489 0 "1319-54" 0 41.6 3 "1619-375" 1 41.7 2 "1675-25" 0 42.4 3 "1317-16" 1 42.4 3 "1619-339" 1 42.7 2 "1317-56" 1 43.02533 0 "1619-281" 0 43.3 2 "1675-31" 0 43.6 4 "1619-83" 0 43.6 3 "1319-62" 1 43.7 2 "1619-235" 0 43.8 2 "1619-231" 0 44.2 2 "1619-214" 1 44.5 4 "1317-52" 1 44.7 2 "1619-79" 1 44.9 4 "1619-302" 1 45.2 2 "1619-211" 0 45.2 2 "1619-332" 0 45.4 1 "1619-89" 1 45.9 4 "1645-36" 0 46 4 "1619-93" 1 46 3 "1645-39" 0 46.3 2 "1619-86" 0 46.4 2 "1619-137" 0 46.4 3 "1645-56" 0 45.4 2 "1619-251" 0 47.3 2 "1645-61" 1 47.4 3 "1318-17" 1 47.7 2 "1619-288" 0 48.5 2 "1619-91" 1 48.7 4 "1619-306" 0 49.3 2 "1619-245" 1 49.3 2 "1319-56" 1 49.3 2 "1321-27" 1 49.4 3 "1619-270" 0 49.4 2 "1319-55" 0 49.4 2 "1619-329" 1 49.9 2 "1619-104" 0 50 2 "1645-40" 0 50.1 3 "1365-27" 0 50.4 2 "1321-52" 0 50.4 2 "1619-379" 1 42.6 2 "1619-371" 0 48.7 2 "2993-7" 1 51.7 2 "1619-227" 1 52.1 2 "1645-1" 0 52.7 2 "1619-164" 1 53 3 "1619-386" 1 43.3 2 "1317-60" 1 53.2 3 "1619-74" 0 53.4 3 "1619-367" 0 44.1 2 "1619-232" 0 54.9 2 "1319-22" 0 55.1 2 "1619-226" 0 55.2 2 "1645-45" 0 55.9 2 "1619-96" 0 56 3 "1619-97" 1 56 3 "1619-230" 0 56.4 2 "1619-335" 1 56.58316 0 "1619-82" 0 56.7 3 "1321-40" 0 57.5 2 "1619-84" 1 59.1 4 "1619-381" 0 48.6 2 "1619-216" 1 60.7 4 "1619-90" 0 60.9 2 "1619-342" 1 61.2 3 "1645-63" 1 38.8 2 "1619-385" 0 62.1 2 "1321-24" 0 62.2 2 "1619-228" 1 62.4 3 "1619-222" 1 62.4 2 "1619-343" 1 64.1 2 "1321-34" 0 64.5 2 "1321-58" 1 -20.25462 0 "1619-221" 0 69.1 2 "1619-265" 1 69.3 3 "2993-8" 1 69.6 2 "1619-331" 1 69.9 2 "1619-338" 0 71.37577 0 "1619-345" 0 72.6 1 "1318-18" 1 74.6 2 "1619-234" 0 75 1 "1619-344" 1 78 1 "1619-348" 0 83.3 1 "1619-333" 1 85.1937 0 "1619-382" 1 39.2 2 "1619-347" 0 47.9 2 "1619-87" 1 64 1 end label values sexo sexo_ label def sexo_ 0 "Hombre", modify label def sexo_ 1 "Mujer", modify label values ev_competi d_evento label def d_evento 0 "Alive", modify label def d_evento 1 "Death", modify label def d_evento 2 "Elective Ao surgery", modify label def d_evento 3 "AAS type A", modify label def d_evento 4 "AAS type B", modify
stset t_ev_combi, failure(ev_competi == 4)
stcrreg sexo, compete(ev_competi== 1 2 3)