Unequal variances on NLMIXED

[marcioadg]

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

[SteveDenham]

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

[marcioadg]

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

[SteveDenham]

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

[lvm]

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