Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Proc nlmixed type3

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 06-24-2019 04:32 AM
(920 views)

Hello,

I use PROC NLMIXED and I would like know if it's possible to have the p-value for fixe effect with the method of type3 in proc nlmixed ?

My code is :

proc nlmixed data=simul ;

parms sigsq0=0.2 sigsqe=0.3 beta0=1 beta1=0.15 beta2=0.55 beta3=0.04 beta4=-0.8 beta5=0.2;

bounds sigsq0 sigsqe >= 0;

pi=2*arsin(1);

mu=beta0+beta1*produit_bis+beta2*visite1+beta3*visite2+beta4*produit_bis*visite1+beta5*produit_bis*visite2 + a_i;

if Censure_aureus=0 then ll=(1/(sqrt(2*pi*sigsqe)))*exp(-(LB_qSTAanum_LOG_cop_cm-mu)**2/(2*sigsqe));

if Censure_aureus=1 then ll=probnorm((LB_qSTAanum_LOG_cop_cm-mu)/sqrt(sigsqe));

L=log(ll);

model LB_qSTAanum_LOG_cop_cm ~ general(L);

random a_i ~ normal(0,sigsq0) subject=id_bis;

by simul;

ods output ParameterEstimates=param_mle ;

run;

With this code, I obtain p-values for beta4 and beta5 which represent the product effect at visite1 and visite2 but, I would like obtain the global pvalue for the interaction. According you, it's possible with proc nlmixed ?

Thanks for your answer.

Clémence

5 REPLIES 5

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Ah, so you want a likelihood ratio test similar to the Type 3 Analysis described in the GENMOD documentation?

There is a SAS Usage note about comparing models by using a LR test.

SAS distributes a macro that compares two nested models by using the likelihood ratio test.

See the PROC NLMIXED example "Failure Time and Frailty Model" for an example.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Thank you. I'm going to read that.

Clémence

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Hello,

I have another question. Is it possible to add an option in proc nlmixed to adjust the alpha risk ?

Clémence

Build your skills. Make connections. Enjoy creative freedom. Maybe change the world. **Registration is now open through August 30th**. Visit the SAS Hackathon homepage.

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.