topic Re: Joint model in Statistical Procedures

Joint model

meleseenyew1919 — Mon, 17 May 2021 11:03:35 GMT

Hither! can I get the SAS code for The joint modeling of the "longitudinal ordinal and parametric proportional Hazards model" and I'm new to such code so please tell me with an example?

Re: Joint model

sbxkoenk — Mon, 17 May 2021 15:51:57 GMT

Hello @meleseenyew1919 ,

I think you should elaborate a bit more on the data you have at your disposal (or the data you will collect) and on the purpose of your research. Keep it concise but clear.

Also, why are you talking about joint modelling? Are there two targets (target variables) you want to model simultaneously?

I guess you will need PROC GLIMMIX or PROC PHREG (SAS/Stat) but to get a concrete code example, you will need to provide us with more detail. Otherwise the example may be completely wrong.

Cheers,

Koen

Re: Joint model

meleseenyew1919 — Tue, 18 May 2021 07:34:36 GMT

Hello! I am so sorry I'm too late! to clarify it I have two response variables one longitudinal ordinal and the other survival models. so when I run a separate analysis of each model using the SAS survival model fulfills the assumption of the parametric proportional model and the Weibull distribution is the best distribution and the Longitudinal model also as well. After this, I need a Joint model Analysis for these two so, please every one of you have a Joint SAS code help me.
you can send with out example also I can practice it!

Thank You!

Re: Joint model

meleseenyew1919 — Tue, 18 May 2021 08:02:52 GMT

Hello, sbxkoenk! sure your guess is correct I need PROC GLIMMIX (SAS) so please help me give me the code even without an example I will practice it!

If you have a question, please!
Thank you!

Re: Joint model

sbxkoenk — Tue, 18 May 2021 08:18:05 GMT

Hello @meleseenyew1919 ,

I will post an example in the afternoon (Brussels time). No time right now.

But I think (I am sure) you can find a lot on internet as well on joint modelling of longitudinal and survival time-to-event data with SAS PROC GLIMMIX.

If your data require 'special features' to be built in in the example code, let me know.

Koen

Re: Joint model

meleseenyew1919 — Tue, 18 May 2021 08:28:15 GMT

Sure I can get a lot of SAS cod but here is the thing the difficulty is that when the longitudinal model is categorical and ordinal to run as jointly!

Thanks for your good communication!

Re: Joint model

meleseenyew1919 — Tue, 18 May 2021 08:36:27 GMT

Simply the longitudinal submodel is the ordinal response!

Re: Joint model

meleseenyew1919 — Tue, 18 May 2021 17:12:28 GMT

Hello, sbxkoenk! I am waiting for you.....

Re: Joint model

sbxkoenk — Tue, 18 May 2021 17:21:51 GMT

Hello @meleseenyew1919 ,

I know. I should have said eve (Brussels time) instead of afternoon. Sorry for that. Have been busy until now. Will start looking for an appropriate example now. Check back in 30 minutes or so.

Koen

Re: Joint model

meleseenyew1919 — Tue, 18 May 2021 17:31:06 GMT

ok, thank you very much!

Re: Joint model

sbxkoenk — Tue, 18 May 2021 18:28:54 GMT

Hello @meleseenyew1919 ,

As you have probably figured out already you need the

dist=byobs(dist)

model option in the MODEL statement of PROC GLIMMIX, where dist is a variable in your data set (having two distinct values for your two targets).

But I don't find an appropriate distribution and link function for your ordinal outcome in PROC GLIMMIX. I guess that's exactly your problem.

I need to consult a colleague for this.

But are you sure you want to solve this via a 'real' JOINT MODEL?

Because, in principle, the joint modeling of longitudinal and time-to-event data could be considered a form of mixed outcome (outcomes of a different data type). And for mixed outcomes, joint models are only one of the three(?) possible approaches to achieve a 'combined' solution.

I'm sorry you have to wait till tomorrow.

You can also check:

Rizopoulos, D. (2012) Joint Models for Longitudinal and Time-to-Event Data. With Applications in R, Chapman & Hall/CRC,
Boca Raton, FL.

, but it will probably not contain SAS code 🤔😕.

Coming back on this tomorrow.

Good luck,

Koen

Re: Joint model

sbxkoenk — Tue, 18 May 2021 18:55:09 GMT

Hello @meleseenyew1919 ,

You marked my first "answer" as a solution. I think that was a bit too "flattered" as your actual question had not been answered at all. When community members see that a question has been solved, they don't look at it anymore. That way you have blocked possible good answers for yourself. But OK, we'll work it out. Unfortunately not today, but hopefully tomorrow.

Have a nice evening,

Koen

Re: Joint model

meleseenyew1919 — Wed, 19 May 2021 07:45:44 GMT

Hello @sbxkoenk,
I know that you are so busy! And as much as you can you try to help me thanks for that. sure I need it! I think it may be helpful but I can't get the cod there is one research on SAS:
Gilani N. Joint modeling of repeated ordinal measures and time to event data for CHD risk assessment. Biom Biostat Int J. 2019;8(6):204‒212.

DOI: 10.15406/bbij.2019.08.00290
If one can contact with authors or can search with it!

Thank you very much!

Re: Joint model

sbxkoenk — Wed, 19 May 2021 08:42:45 GMT

Hello @meleseenyew1919 ,

I will try to find some sample code today.

Is there a PROC GLIMMIX in the article:

Gilani N. Joint modeling of repeated ordinal measures and time to event data for CHD risk assessment. Biom Biostat Int J. 2019;8(6):204‒212.

? (I haven't opened up this article)

If there's some PROC GLIMMIX code in this article, does it meet your expectations or are there shortcomings? If the code there does not meet your expectations (case), explain why.

Thanks,

Koen

Re: Joint model

meleseenyew1919 — Wed, 19 May 2021 09:00:06 GMT

Hello Koen!
ok if possible in PROC GLIMMIX it fine but most of the
categorical longitudinal outcome uses PROC NLMIXED including this it has
a similar procedure to my analysis. fill free to use PROC GLIMMIX or PROC
NLMIXED I need the output only!

for the article alternative: try it to open
https://medcraveonline.com/BBIJ/BBIJ-08-00290.pdf

Thank you very much!

Re: Joint model

meleseenyew1919 — Thu, 20 May 2021 06:34:53 GMT

I can't open this, is there any attachment on View solution in original post?

Re: Joint model

meleseenyew1919 — Thu, 20 May 2021 07:15:55 GMT

Hello @sbxkoenk! Have you seen the article! And so far I can't get any solution.
I can't access this file/ your first solution:
https://communities.sas.com/t5/Statistical-Procedures/Joint-model/m-p/741898#M36084

Re: Joint model

sbxkoenk — Thu, 20 May 2021 09:16:49 GMT

Hello,

How do you mean: "You can't open / access it" ??

You are in it.

You are on the second page. Scroll to the bottom and click on "1" (page 1). There's a button bar showing a "1" and a "2" (page 1 and 2, respectively).

Click on 1 and you will see all the earlier posts.

Koen

Re: Joint model

sbxkoenk — Thu, 20 May 2021 09:27:19 GMT

Hello @meleseenyew1919 ,

I contacted a colleague meanwhile.

She says:

===========

GLIMMIX can estimate ordinal model using DIST=MULTINOMIAL LINK=CLOGIT.
However, GLIMMIX only supports G-side random effects for Multinomial models. In other words, you can not model R-side random effects (correlated errors for repeated measures) for Multinomial models in PROC GLIMMIX.

I do not have any examples of joint modeling of longitudinal and time-to-event (survival) data.

===========

[Later edit]

Also have a look at this entry:

Looking for SAS code for Joint Model between survival and longitudinal analysis model
https://communities.sas.com/t5/SAS-Health-and-Life-Sciences/Looking-for-SAS-code-for-Joint-Model-between-survival-and/m-p/263493#M2061

Koen

Re: Joint model

meleseenyew1919 — Thu, 20 May 2021 10:18:06 GMT

Hello Koen! sorry for that! on google, the question is put as a solved
question I know that is my mistake ( marked my first "answer" as a
solution, because this is my first question in SAS) so I would like to say
I can't rich the answer!

PROC NLMIXED is best for ordinal rather than PROC GLIMMIX as some research
states. so if you can give me the last chance with a great excuse I don't
get the code so far.

Thank you so much for your patience and all helps🙏!