DATA JessicaPCR; options formdlim = '*';
INPUT ID DAY TRT BLOCK AMPK IGF MHCI MHCIIA MHCIIX RAGA AKT EIF;
CARDS;
219	0	0	2	0.201	0.015	0.416	0.310	0.190	0.846	0.543	0.897
223	0	0	1	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000
287	0	0	1	0.543	1.097	0.784	0.338	0.561	0.756	0.816	0.634
296	0	0	2	0.455	0.920	0.494	0.213	0.246	1.253	0.887	0.765
654	0	0	2	0.657	1.209	.	0.002	0.185	0.930	0.890	0.969
668	0	0	1	0.845	0.712	0.898	0.497	0.584	.	0.907	1.081
681	0	0	2	0.243	0.397	0.613	0.433	0.332	1.167	0.859	0.751
692	0	0	1	0.446	0.600	0.634	0.371	0.401	0.975	0.678	0.959
723	0	0	1	0.980	1.498	0.773	0.503	0.282	1.267	.	.
920	0	0	2	1.002	1.778	1.663	0.631	0.444	1.792	0.787	0.731
256	0	8	1	0.908	0.892	0.947	0.312	0.645	1.290	0.810	0.756
653	0	8	1	0.634	0.784	1.088	0.587	0.372	1.268	0.693	0.916
664	0	8	2	0.417	0.739	0.598	0.250	0.313	1.102	0.535	0.682
674	0	8	1	0.554	0.987	1.212	0.919	0.280	1.458	0.660	0.981
709	0	8	1	2.139	5.094	0.885	0.025	0.146	.	0.852	0.251
710	0	8	2	1.083	1.154	1.207	0.785	1.007	1.050	0.496	0.861
715	0	8	2	0.697	0.820	0.704	0.317	0.652	1.137	0.678	0.731
745	0	8	1	0.477	0.790	1.207	0.686	0.424	0.691	0.578	0.757
924	0	8	2	0.612	0.685	1.264	0.461	0.540	1.318	0.624	.
934	0	8	2	0.464	1.135	0.923	0.431	0.482	0.925	0.438	0.718
219	14	0	2	0.564	0.523	1.077	0.666	0.265	1.230	0.713	0.611
223	14	0	1	0.851	0.881	0.745	0.461	0.688	0.901	0.973	0.343
287	14	0	1	0.603	1.821	1.099	0.597	0.571	0.974	0.428	0.466
296	14	0	2	0.437	0.317	0.461	0.296	0.287	1.060	0.973	0.678
654	14	0	2	0.314	1.594	1.031	0.675	0.418	0.585	0.481	0.786
668	14	0	1	0.640	0.641	1.139	0.576	0.442	0.897	.	.
681	14	0	2	0.876	0.437	0.600	0.342	0.362	1.384	1.115	0.827
692	14	0	1	0.367	0.501	0.559	0.323	0.392	1.219	0.624	0.747
723	14	0	1	0.918	0.996	0.981	0.452	0.374	1.733	0.721	0.317
920	14	0	2	0.663	0.729	0.904	0.552	0.393	1.338	0.821	0.516
256	14	8	1	1.105	0.783	0.942	0.252	0.558	1.038	0.791	0.444
653	14	8	1	0.555	0.381	1.304	0.645	0.247	0.724	0.347	0.376
664	14	8	2	0.360	0.825	0.968	0.357	0.321	0.766	0.403	.
674	14	8	1	0.622	0.422	1.031	0.849	0.273	0.844	0.557	0.502
709	14	8	1	0.635	1.213	0.827	0.557	0.627	1.151	0.639	0.883
710	14	8	2	0.554	1.476	1.018	0.602	0.689	1.754	0.506	1.040
715	14	8	2	0.594	0.072	0.896	0.449	0.368	0.960	.	.
745	14	8	1	0.896	1.424	1.407	0.827	0.687	1.223	1.020	0.549
924	14	8	2	0.842	0.783	1.039	0.402	0.551	1.563	.	1.637
934	14	8	2	0.530	0.772	0.492	0.244	0.334	1.104	0.617	0.892
219	28	0	2	0.837	0.949	0.489	0.647	0.238	1.247	0.808	0.844
223	28	0	1	1.128	0.874	0.696	0.562	0.535	1.344	0.692	0.831
287	28	0	1	0.763	2.910	0.381	0.197	0.264	1.240	0.503	0.382
296	28	0	2	0.200	1.969	0.443	0.341	0.257	1.146	0.013	0.023
654	28	0	2	0.398	0.722	0.642	0.257	0.371	0.625	0.729	0.598
668	28	0	1	0.584	1.403	0.876	0.324	0.431	1.137	0.720	0.649
681	28	0	2	0.235	0.455	0.319	0.162	0.106	0.908	0.529	0.561
692	28	0	1	.	.	.	.	.	.	0.528	0.676
723	28	0	1	0.815	0.548	0.975	0.873	0.219	1.373	0.451	1.335
920	28	0	2	0.410	0.092	0.956	0.555	0.339	0.960	0.252	0.418
256	28	8	1	.	.	.	.	.	.	0.709	0.192
653	28	8	1	0.735	1.398	0.990	0.456	0.239	1.042	0.772	0.965
664	28	8	2	0.306	1.586	0.710	0.206	0.264	0.655	0.675	1.447
674	28	8	1	0.617	1.122	0.936	0.921	0.522	1.104	0.993	0.894
709	28	8	1	0.556	1.094	0.678	0.315	0.250	1.122	0.455	0.531
710	28	8	2	0.416	1.313	0.696	0.296	0.430	0.886	0.183	0.802
715	28	8	2	0.050	0.580	0.644	0.786	0.266	0.207	1.093	0.969
745	28	8	1	0.460	0.996	0.831	0.400	0.264	0.874	0.485	0.515
924	28	8	2	.	.	1.175	1.243	0.316	3.198	0.014	.
934	28	8	2	0.465	1.032	0.595	0.245	0.381	0.814	1.017	.
219	42	0	2	2.890	2.807	1.647	2.235	0.380	1.425	0.937	0.377
223	42	0	1	1.825	9.135	1.948	1.468	1.018	1.639	1.621	0.620
287	42	0	1	1.988	8.114	2.402	1.059	1.080	1.606	0.914	0.278
296	42	0	2	1.862	5.359	3.145	1.851	1.195	2.027	1.361	0.662
654	42	0	2	1.545	1.721	1.695	1.780	0.290	1.199	0.905	0.662
668	42	0	1	2.529	9.389	3.571	2.240	0.690	1.994	1.210	0.394
681	42	0	2	2.827	7.804	3.683	2.480	1.366	2.369	.	0.809
692	42	0	1	1.938	3.998	4.075	1.805	1.568	1.843	1.061	0.590
723	42	0	1	1.842	3.616	2.118	2.266	0.451	1.848	0.763	0.338
920	42	0	2	2.765	7.593	.	3.460	1.591	3.305	1.183	0.647
256	42	8	1	2.932	.	3.331	1.110	1.062	2.360	1.038	0.454
653	42	8	1	2.131	3.518	1.521	0.877	0.359	1.332	.	.
664	42	8	2	2.000	10.341	3.064	1.202	1.031	2.486	0.803	0.610
674	42	8	1	2.245	.	3.435	.	0.526	2.283	1.575	0.923
709	42	8	1	1.069	10.767	.	2.207	0.725	1.106	1.021	0.546
710	42	8	2	2.368	7.210	5.365	2.946	.	2.517	.	0.824
715	42	8	2	1.996	4.831	4.770	1.789	1.679	1.777	0.822	0.447
745	42	8	1	2.341	4.224	3.387	3.063	1.488	2.236	1.438	0.455
924	42	8	2	2.076	6.155	5.813	1.455	1.535	2.751	0.954	0.863
934	42	8	2	0.466	7.911	0.884	0.235	0.005	1.088	.	.
219	56	0	2	1.983	6.742	4.048	3.116	0.427	2.701	1.258	0.960
223	56	0	1	2.840	8.211	4.184	.	.	3.247	1.492	0.738
287	56	0	1	1.915	5.872	3.018	1.590	1.245	2.693	0.880	0.364
296	56	0	2	1.914	5.125	4.039	1.164	1.035	2.352	0.999	1.444
654	56	0	2	2.859	6.062	3.410	2.565	1.753	2.894	1.269	0.426
668	56	0	1	2.315	3.428	2.739	1.993	0.906	2.476	.	.
681	56	0	2	1.504	4.409	2.506	2.050	0.851	3.118	0.591	0.540
692	56	0	1	1.784	10.243	1.396	1.076	0.609	3.441	.	.
723	56	0	1	2.751	4.972	3.638	2.383	1.186	3.784	1.005	0.384
920	56	0	2	.	.	.	.	.	.	1.115	0.598
934	56	8	2	1.900	4.658	2.231	0.953	0.766	2.950	0.587	0.510
256	56	8	1	2.315	3.816	3.167	0.884	1.067	2.060	0.898	0.410
653	56	8	1	2.084	5.695	4.437	3.609	0.741	0.826	0.851	0.471
664	56	8	2	1.356	2.970	5.130	1.782	0.788	1.603	0.696	0.405
674	56	8	1	1.259	2.308	2.318	1.954	0.561	1.434	0.884	0.636
709	56	8	1	1.656	7.138	1.875	1.547	0.787	2.099	0.959	0.645
710	56	8	2	1.194	2.724	2.813	1.496	0.977	2.334	0.001	.
715	56	8	2	1.373	4.080	2.411	2.056	0.937	.	0.739	0.828
745	56	8	1	2.614	4.007	3.309	2.298	1.202	3.032	1.204	0.318
924	56	8	2	2.358	5.706	3.569	1.737	1.537	3.740	0.918	0.382
;
dm "output; clear; log; clear";
Proc glimmix;
Class ID TRT BLOCK DAY;
Model AMPK = TRT|DAY/ddfm = KR;
random intercept / subject=BLOCK;
random intercept / subject=ID*TRT;
random _residual_/subject = ID*TRT type=arh(1);
LSMEANS TRT|DAY/ pdiff lines;
nloptions tech=nrridg;
Run;
Proc glimmix;
Class ID TRT BLOCK DAY;
Model IGF = TRT|DAY/ddfm = KR;
random intercept / subject=BLOCK;
random intercept / subject=ID*TRT;
random _residual_/subject = ID*TRT type=csh;
LSMEANS TRT|DAY/ pdiff lines;
nloptions tech=nrridg;
Run;
Proc glimmix;
Class ID TRT BLOCK DAY;
Model MHCI = TRT|DAY/ddfm = KR;
random intercept / subject=BLOCK;
random intercept / subject=ID*TRT;
random _residual_/subject = ID*TRT type=csh;
LSMEANS TRT|DAY/ pdiff lines;
nloptions tech=nrridg;
Run;
Proc glimmix;
Class ID TRT BLOCK DAY;
Model MHCIIA = TRT|DAY/ddfm = KR;
random intercept / subject=BLOCK;
random intercept / subject=ID*TRT;
random _residual_/subject = ID*TRT type=csh;
LSMEANS TRT|DAY/ pdiff lines;
nloptions tech=nrridg;
Run;
Proc glimmix;
Class ID TRT BLOCK DAY;
Model MHCIIX = TRT|DAY/ddfm = KR;
random intercept / subject=BLOCK;
random intercept / subject=ID*TRT;
random _residual_/subject = ID*TRT type=ante(1);
LSMEANS TRT|DAY/ pdiff lines;
nloptions tech=nrridg;
Run;
Proc glimmix;
Class ID TRT BLOCK DAY;
Model RAGA = TRT|DAY/ddfm = KR;
random intercept / subject=BLOCK;
random intercept / subject=ID*TRT;
random _residual_/subject = ID*TRT type=csh;
LSMEANS TRT|DAY/ pdiff lines;
nloptions tech=nrridg;
Run;
Proc glimmix;
Class ID TRT BLOCK DAY;
Model AKT = TRT|DAY/ddfm = KR;
random intercept / subject=BLOCK;
random intercept / subject=ID*TRT;
random _residual_/subject = ID*TRT type=csh;
LSMEANS TRT|DAY/ pdiff lines;
nloptions tech=nrridg;
Run;
Proc glimmix;
Class ID TRT BLOCK DAY;
Model EIF = TRT|DAY/ddfm = KR;
random intercept / subject=BLOCK;
random intercept / subject=ID*TRT;
random _residual_/subject = ID*TRT type=ante(1);
LSMEANS TRT|DAY/ pdiff lines;
nloptions tech=nrridg;
Run;