首先要去SAS的官网下载GLMMIX的安装程序,要不你9.13里面没有,9.2以上才有这个过程,所以要独立安装,地址如下,
http://www.sas.com/apps/demosdownloads/setupcat.jsp;jsessionid=81A2F5677BD9B44787C52F1D46C77B79.tomcat1?cat=SAS%2FSTAT+Software

你的自变量的类型你说的不是很明确,不知道那些是分类的那些是连续的,是否有等级的,我按我的理解给你写一小段,然后告诉你都是什么东西,根据你自己的变量类型在自己写好了,

proc glimmix data=work  method=RMPL ;    /*work是你的数据文件名,method是你用的计算方法*/
class fid  precipitation;         /*你的家庭id或家庭的分类的只是变量,做随机效应,如果你还

                                    有其他的分类变量都写在后面我不知道你这个变量是什么类

                                    型的我假设它是分类的了precipitation,

                                     所以也放这个语句后面,如果不是就自行删除*/
model Time= temperature daylength cloudcover，precipitation/solution  chisq dist=logit ;
                                   /*这个是固定效应的语句,所有固定效应都放这里,你上面说其他变

                                    量都是混合效应,是不对的fixed,effects是固定效应的意思不是混

                                     合效应,混合效应是说含有随机效应和固定效应两种情况下称为混

                                    合效应,dist是定义广义线性混合模型的link function,如果你的

                                     反应变量是连续的则定义为normal,如果是分类的用logit,我推断

                                      你是分类的所以就用了logit*/
random intercept /subject=fid  type=vc ;     

                                  /*定义随机效应fid有关,type是定义协方差矩阵的类型,

                                       有很多种,默认为VC,你可根据你的情况去自己看下用什么协方

                                       差阵,我曾经做双生子研究估计加性遗传效应的时候用过

                                        teop(1)的矩阵,你的数据具体是什么样我不清楚不好帮

                                          你做选择*/
freq count;                        /*这个是要强调你的repeated,每个家庭的重复数count,

                                       这个过程需要你在你的数据集work中先做好*/
run;

就这样了,不知道对你有没什么帮助~~~黑体是语句,/**/是说明,这个不知道怎么空格显示的更清楚

论坛上有SAS9.13的glimmix试用程序 自己去找找