Programming the statistical procedures from SAS

Unequal variances on NLMIXED

Reply
Contributor
Posts: 45

Unequal variances on NLMIXED

This is the part of the code I used for equal variances and work fine.

model y ~ normal(L + U*(z1) + lvar, s2e);

random lvar ~ normal(0,blockvar) subject=block out=blups;

For unequal variances, I tried the one below and it outputted errors in the log.

model y ~ normal(L + U*(z1) + lvar, s2a s2b);

random lvar ~ normal(0,blockvar) subject=block out=blups;

Could you please help me identify what I am missing,

Thank you,

Marcio

Respected Advisor
Posts: 2,655

Re: Unequal variances on NLMIXED

The normal argument to the MODEL statement takes only a mean and variance.  I think you will have to code up the log likelihood function with the two variances, and use the general argument.

Steve Denham

Contributor
Posts: 45

Re: Unequal variances on NLMIXED

Thank you Steve,

But I am not sure how to do that. I looked into some papers (http://www.education.umd.edu/EDMS/fac/Harring/MIsc/WUSS-Paper.pdf) and websites (SAS/STAT(R) 9.22 User's Guide) but I am not being able to do it.

Would you recommend any resource on how to do that?

Thank you,

Marcio

Respected Advisor
Posts: 2,655

Re: Unequal variances on NLMIXED

I have never done it, so I can't offer much beyond searching this site and the SAS-L archives.  I'm sorry.

Steve Denham

Valued Guide
Valued Guide
Posts: 684

Re: Unequal variances on NLMIXED

It is not clear what you are trying to do. I suspect you have two groups and want a different variance for each group?

Ask a Question
Discussion stats
  • 4 replies
  • 281 views
  • 6 likes
  • 3 in conversation