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
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
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
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
It is not clear what you are trying to do. I suspect you have two groups and want a different variance for each group?
Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.
Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.
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.