STATA slowness or problems in my data?

[Gianfranco Di Gennaro]

Dear STATA users, I am writing to ask for a general opinion on STATA.
I am running a logistic multilevel model with melogit.
This is a study with 900 patients from 23 hospitals each observed at different timepoints at which I assess the presence / absence of adverse drug events.
The model is also quite simple:

melogit events i.sex age time || center: || id:

Is it normal for this kind of model STATA takes 20 minutes to give me an output?

Or is it a problem with the structure of my data?
When I use only one of the two levels || id: OR || center: the output is immediate

Thank you very much.

Gianfranco

[Joseph Coveney]

With five observation intervals per patient and patients equally distributed among hospitals, I get about five seconds for the three-level hierarchical model and about one second when I move the top hierarchy over to the fixed effects equation. It's an ordinary Windows laptop, four cores, Stata MP 4. So, maybe something in your arrangement?

.ÿ
.ÿversionÿ17.0

.ÿ
.ÿclearÿ*

.ÿ
.ÿ//ÿseedem
.ÿsetÿseedÿ924341131

.ÿ
.ÿ//ÿHospitals
.ÿquietlyÿsetÿobsÿ23

.ÿgenerateÿbyteÿcidÿ=ÿ_n

.ÿgenerateÿdoubleÿcid_uÿ=ÿrnormal()

.ÿ
.ÿ//ÿPatients
.ÿquietlyÿexpandÿ`=ceil(900ÿ/ÿ_N)'

.ÿgenerateÿintÿpidÿ=ÿ_n

.ÿgenerateÿdoubleÿpid_uÿ=ÿrnormal()

.ÿgenerateÿintÿageÿ=ÿruniformint(1,ÿ99)

.ÿgenerateÿbyteÿsexÿ=ÿmod(_n,ÿ2)

.ÿ
.ÿ//ÿEvaluations
.ÿquietlyÿexpandÿ5

.ÿbysortÿpid:ÿgenerateÿbyteÿtimÿ=ÿ_n

.ÿ
.ÿ//ÿAdverseÿEvents
.ÿgenerateÿdoubleÿxbuÿ=ÿcid_uÿ+ÿpid_u

.ÿgenerateÿbyteÿoutÿ=ÿrbinomial(1,ÿinvlogit(xbu))

.ÿ
.ÿ*
.ÿ*ÿBeginÿhere
.ÿ*
.ÿtimerÿclear

.ÿtimerÿonÿ1

.ÿmelogitÿoutÿi.(sexÿtim)ÿc.ageÿ||ÿcid:ÿ||ÿpid:ÿ,ÿnogroupÿnolog

Mixed-effectsÿlogisticÿregressionÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿNumberÿofÿobsÿÿÿÿÿ=ÿÿÿÿÿÿ4,600

Integrationÿmethod:ÿmvaghermiteÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿIntegrationÿpts.ÿÿ=ÿÿÿÿÿÿÿÿÿÿ7

ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿWaldÿchi2(6)ÿÿÿÿÿÿ=ÿÿÿÿÿÿÿ8.03
Logÿlikelihoodÿ=ÿ-2765.9576ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿProbÿ>ÿchi2ÿÿÿÿÿÿÿ=ÿÿÿÿÿ0.2361
------------------------------------------------------------------------------
ÿÿÿÿÿÿÿÿÿoutÿ|ÿCoefficientÿÿStd.ÿerr.ÿÿÿÿÿÿzÿÿÿÿP>|z|ÿÿÿÿÿ[95%ÿconf.ÿinterval]
-------------+----------------------------------------------------------------
ÿÿÿÿÿÿÿ1.sexÿ|ÿÿÿ.0218186ÿÿÿ.0981237ÿÿÿÿÿ0.22ÿÿÿ0.824ÿÿÿÿ-.1705004ÿÿÿÿ.2141376
ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
ÿÿÿÿÿÿÿÿÿtimÿ|
ÿÿÿÿÿÿÿÿÿÿ2ÿÿ|ÿÿ-.2327173ÿÿÿ.1108187ÿÿÿÿ-2.10ÿÿÿ0.036ÿÿÿÿÿ-.449918ÿÿÿ-.0155167
ÿÿÿÿÿÿÿÿÿÿ3ÿÿ|ÿÿ-.1599272ÿÿÿ.1109845ÿÿÿÿ-1.44ÿÿÿ0.150ÿÿÿÿ-.3774528ÿÿÿÿ.0575984
ÿÿÿÿÿÿÿÿÿÿ4ÿÿ|ÿÿ-.0372126ÿÿÿ.1113779ÿÿÿÿ-0.33ÿÿÿ0.738ÿÿÿÿ-.2555093ÿÿÿÿ.1810842
ÿÿÿÿÿÿÿÿÿÿ5ÿÿ|ÿÿ-.2266723ÿÿÿ.1108305ÿÿÿÿ-2.05ÿÿÿ0.041ÿÿÿÿ-.4438962ÿÿÿ-.0094485
ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
ÿÿÿÿÿÿÿÿÿageÿ|ÿÿÿ.0012073ÿÿÿ.0017812ÿÿÿÿÿ0.68ÿÿÿ0.498ÿÿÿÿ-.0022838ÿÿÿÿ.0046984
ÿÿÿÿÿÿÿ_consÿ|ÿÿÿ.4653423ÿÿÿ.2416484ÿÿÿÿÿ1.93ÿÿÿ0.054ÿÿÿÿ-.0082799ÿÿÿÿ.9389645
-------------+----------------------------------------------------------------
cidÿÿÿÿÿÿÿÿÿÿ|
ÿÿÿvar(_cons)|ÿÿÿ.9161943ÿÿÿ.2903086ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.4923487ÿÿÿÿ1.704914
-------------+----------------------------------------------------------------
cid>pidÿÿÿÿÿÿ|
ÿÿÿvar(_cons)|ÿÿÿ1.028852ÿÿÿ.1283751ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.8056463ÿÿÿÿ1.313898
------------------------------------------------------------------------------
LRÿtestÿvs.ÿlogisticÿmodel:ÿchi2(2)ÿ=ÿ721.16ÿÿÿÿÿÿÿÿÿÿÿÿÿÿProbÿ>ÿchi2ÿ=ÿ0.0000

Note:ÿLRÿtestÿisÿconservativeÿandÿprovidedÿonlyÿforÿreference.

.ÿtimerÿoffÿ1

.ÿ
.ÿ//ÿ40ÿpidsÿ(200ÿobservations)ÿperÿcid
.ÿtimerÿonÿ2

.ÿmelogitÿoutÿi.(cidÿsexÿtim)ÿc.ageÿ||ÿpid:ÿ,ÿnolog

Mixed-effectsÿlogisticÿregressionÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿNumberÿofÿobsÿÿÿÿÿ=ÿÿÿÿÿÿ4,600
Groupÿvariable:ÿpidÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿNumberÿofÿgroupsÿÿ=ÿÿÿÿÿÿÿÿ920

ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿObsÿperÿgroup:
ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿminÿ=ÿÿÿÿÿÿÿÿÿÿ5
ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿavgÿ=ÿÿÿÿÿÿÿÿ5.0
ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿmaxÿ=ÿÿÿÿÿÿÿÿÿÿ5

Integrationÿmethod:ÿmvaghermiteÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿIntegrationÿpts.ÿÿ=ÿÿÿÿÿÿÿÿÿÿ7

ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿWaldÿchi2(28)ÿÿÿÿÿ=ÿÿÿÿÿ331.64
Logÿlikelihoodÿ=ÿ-2721.4333ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿProbÿ>ÿchi2ÿÿÿÿÿÿÿ=ÿÿÿÿÿ0.0000
------------------------------------------------------------------------------
ÿÿÿÿÿÿÿÿÿoutÿ|ÿCoefficientÿÿStd.ÿerr.ÿÿÿÿÿÿzÿÿÿÿP>|z|ÿÿÿÿÿ[95%ÿconf.ÿinterval]
-------------+----------------------------------------------------------------
ÿÿÿÿÿÿÿÿÿcidÿ|
ÿÿÿÿÿÿÿÿÿÿ2ÿÿ|ÿÿÿ1.537158ÿÿÿ.3220072ÿÿÿÿÿ4.77ÿÿÿ0.000ÿÿÿÿÿ.9060356ÿÿÿÿÿ2.16828
ÿÿÿÿÿÿÿÿÿÿ3ÿÿ|ÿÿ-.8018228ÿÿÿÿ.341756ÿÿÿÿ-2.35ÿÿÿ0.019ÿÿÿÿ-1.471652ÿÿÿ-.1319934
ÿÿÿÿÿÿÿÿÿÿ4ÿÿ|ÿÿÿ2.268876ÿÿÿ.3343387ÿÿÿÿÿ6.79ÿÿÿ0.000ÿÿÿÿÿ1.613584ÿÿÿÿ2.924167
ÿÿÿÿÿÿÿÿÿÿ5ÿÿ|ÿÿÿ1.026478ÿÿÿ.3187167ÿÿÿÿÿ3.22ÿÿÿ0.001ÿÿÿÿÿ.4018045ÿÿÿÿ1.651151
ÿÿÿÿÿÿÿÿÿÿ6ÿÿ|ÿÿÿ2.539473ÿÿÿ.3417587ÿÿÿÿÿ7.43ÿÿÿ0.000ÿÿÿÿÿ1.869638ÿÿÿÿ3.209307
ÿÿÿÿÿÿÿÿÿÿ7ÿÿ|ÿÿÿ1.620382ÿÿÿ.3230734ÿÿÿÿÿ5.02ÿÿÿ0.000ÿÿÿÿÿ.9871702ÿÿÿÿ2.253595
ÿÿÿÿÿÿÿÿÿÿ8ÿÿ|ÿÿÿ1.407583ÿÿÿ.3221164ÿÿÿÿÿ4.37ÿÿÿ0.000ÿÿÿÿÿ.7762466ÿÿÿÿ2.038919
ÿÿÿÿÿÿÿÿÿÿ9ÿÿ|ÿÿÿÿ2.61961ÿÿÿ.3445831ÿÿÿÿÿ7.60ÿÿÿ0.000ÿÿÿÿÿ1.944239ÿÿÿÿÿ3.29498
ÿÿÿÿÿÿÿÿÿ10ÿÿ|ÿÿÿ1.437929ÿÿÿ.3221672ÿÿÿÿÿ4.46ÿÿÿ0.000ÿÿÿÿÿ.8064933ÿÿÿÿ2.069365
ÿÿÿÿÿÿÿÿÿ11ÿÿ|ÿÿÿÿ2.20045ÿÿÿ.3325238ÿÿÿÿÿ6.62ÿÿÿ0.000ÿÿÿÿÿ1.548715ÿÿÿÿ2.852184
ÿÿÿÿÿÿÿÿÿ12ÿÿ|ÿÿÿÿ1.22788ÿÿÿ.3199982ÿÿÿÿÿ3.84ÿÿÿ0.000ÿÿÿÿÿ.6006951ÿÿÿÿ1.855065
ÿÿÿÿÿÿÿÿÿ13ÿÿ|ÿÿÿ1.267753ÿÿÿ.3205375ÿÿÿÿÿ3.96ÿÿÿ0.000ÿÿÿÿÿ.6395112ÿÿÿÿ1.895995
ÿÿÿÿÿÿÿÿÿ14ÿÿ|ÿÿÿ1.710861ÿÿÿ.3239335ÿÿÿÿÿ5.28ÿÿÿ0.000ÿÿÿÿÿ1.075963ÿÿÿÿ2.345759
ÿÿÿÿÿÿÿÿÿ15ÿÿ|ÿÿ-1.186634ÿÿÿ.3548775ÿÿÿÿ-3.34ÿÿÿ0.001ÿÿÿÿ-1.882181ÿÿÿ-.4910867
ÿÿÿÿÿÿÿÿÿ16ÿÿ|ÿÿÿ1.308686ÿÿÿ.3219386ÿÿÿÿÿ4.07ÿÿÿ0.000ÿÿÿÿÿ.6776978ÿÿÿÿ1.939674
ÿÿÿÿÿÿÿÿÿ17ÿÿ|ÿÿÿ1.891678ÿÿÿ.3264935ÿÿÿÿÿ5.79ÿÿÿ0.000ÿÿÿÿÿ1.251763ÿÿÿÿ2.531594
ÿÿÿÿÿÿÿÿÿ18ÿÿ|ÿÿÿ1.234973ÿÿÿÿ.319863ÿÿÿÿÿ3.86ÿÿÿ0.000ÿÿÿÿÿ.6080532ÿÿÿÿ1.861893
ÿÿÿÿÿÿÿÿÿ19ÿÿ|ÿÿÿ2.620398ÿÿÿÿ.344494ÿÿÿÿÿ7.61ÿÿÿ0.000ÿÿÿÿÿ1.945202ÿÿÿÿ3.295594
ÿÿÿÿÿÿÿÿÿ20ÿÿ|ÿÿÿ1.767052ÿÿÿ.3253705ÿÿÿÿÿ5.43ÿÿÿ0.000ÿÿÿÿÿ1.129338ÿÿÿÿ2.404767
ÿÿÿÿÿÿÿÿÿ21ÿÿ|ÿÿ-.1420833ÿÿÿ.3257514ÿÿÿÿ-0.44ÿÿÿ0.663ÿÿÿÿ-.7805443ÿÿÿÿ.4963776
ÿÿÿÿÿÿÿÿÿ22ÿÿ|ÿÿÿ1.371258ÿÿÿ.3221677ÿÿÿÿÿ4.26ÿÿÿ0.000ÿÿÿÿÿ.7398208ÿÿÿÿ2.002695
ÿÿÿÿÿÿÿÿÿ23ÿÿ|ÿÿÿ.9441795ÿÿÿ.3182537ÿÿÿÿÿ2.97ÿÿÿ0.003ÿÿÿÿÿ.3204137ÿÿÿÿ1.567945
ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
ÿÿÿÿÿÿÿ1.sexÿ|ÿÿÿ.0222675ÿÿÿ.0968704ÿÿÿÿÿ0.23ÿÿÿ0.818ÿÿÿÿ-.1675951ÿÿÿÿ.2121301
ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
ÿÿÿÿÿÿÿÿÿtimÿ|
ÿÿÿÿÿÿÿÿÿÿ2ÿÿ|ÿÿ-.2327679ÿÿÿÿ.110831ÿÿÿÿ-2.10ÿÿÿ0.036ÿÿÿÿ-.4499928ÿÿÿ-.0155431
ÿÿÿÿÿÿÿÿÿÿ3ÿÿ|ÿÿ-.1599793ÿÿÿ.1110027ÿÿÿÿ-1.44ÿÿÿ0.150ÿÿÿÿ-.3775407ÿÿÿÿÿ.057582
ÿÿÿÿÿÿÿÿÿÿ4ÿÿ|ÿÿ-.0372316ÿÿÿ.1114064ÿÿÿÿ-0.33ÿÿÿ0.738ÿÿÿÿ-.2555842ÿÿÿÿ.1811209
ÿÿÿÿÿÿÿÿÿÿ5ÿÿ|ÿÿ-.2267236ÿÿÿ.1108434ÿÿÿÿ-2.05ÿÿÿ0.041ÿÿÿÿ-.4439727ÿÿÿ-.0094746
ÿÿÿÿÿÿÿÿÿÿÿÿÿ|
ÿÿÿÿÿÿÿÿÿageÿ|ÿÿÿ.0012638ÿÿÿ.0017598ÿÿÿÿÿ0.72ÿÿÿ0.473ÿÿÿÿ-.0021853ÿÿÿÿ.0047128
ÿÿÿÿÿÿÿ_consÿ|ÿÿ-.8416388ÿÿÿ.2580759ÿÿÿÿ-3.26ÿÿÿ0.001ÿÿÿÿ-1.347458ÿÿÿ-.3358193
-------------+----------------------------------------------------------------
pidÿÿÿÿÿÿÿÿÿÿ|
ÿÿÿvar(_cons)|ÿÿÿ.9740969ÿÿÿ.1234807ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿ.7598015ÿÿÿÿ1.248832
------------------------------------------------------------------------------
LRÿtestÿvs.ÿlogisticÿmodel:ÿchibar2(01)ÿ=ÿ167.41ÿÿÿÿÿÿProbÿ>=ÿchibar2ÿ=ÿ0.0000

.ÿtimerÿoffÿ2

.ÿ
.ÿtimerÿlist
ÿÿÿ1:ÿÿÿÿÿÿ5.21ÿ/ÿÿÿÿÿÿÿÿ1ÿ=ÿÿÿÿÿÿÿ5.2060
ÿÿÿ2:ÿÿÿÿÿÿ0.98ÿ/ÿÿÿÿÿÿÿÿ1ÿ=ÿÿÿÿÿÿÿ0.9780

.ÿ
.ÿexit

endÿofÿdo-file


.