Date: Sun, 1 Mar 2009 16:03:00 -0800 sej.in.md@GMAIL.COM "SAS(r) Discussion" sej.in.md@GMAIL.COM http://groups.google.com Re: PROC GLIMMIX: 2 questions: blocks, and covariance matrices To: sas-l@uga.edu text/plain; charset=ISO-8859-1

"Are you analysing the colonies as your outcome variable, or the termites in your colonies? If it is the termites, perhaps you need to specify colonies as an explanatory variable. Can you show us the structure of your file?"

Each colony is treated as an EU. The outcome variable is a count of this particular kind of termite in each colony. (I counted them each month, with the thought of doing a repeated-measures analysis, but decided just to look at the last month's count in the statistical analysis.) Below is the data entry step.

First, some explanation--this actually accounts for 3 different types of termites that appear within each colony: mrs, mn, fn. I count each of them through each of 6 months. I'm only going to analyze them at the 6 month mark--eg mrs6, mn6--and will analyze each type separately. The first column is the treatment group, which I turn into the dummy variable "trtmt." The "colony" column's numerical value is just the code name for each individual colony.

input group \$ block colony mrs1 mrs2 mrs3 mrs4 mrs5 mrs6 mn1 mn2 mn3 mn4 mn5 mn6 fn1 fn2 fn3 fn4 fn5 fn6; if group='c' then trtmt=0; else if group = 't' then trtmt=1;

cards;

c 1 804 0 0 0 0 0 0 0 0 1 1 1 1 0 0 0 0 0 0 c 1 806 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 c 1 811 0 1 1 1 1 1 0 0 1 1 1 1 0 0 0 0 0 0 c 1 814 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 c 1 833 0 1 1 1 1 1 0 1 1 0 0 0 0 0 0 0 0 0 c 1 838 0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 c 1 845 0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 c 1 902 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 c 2 823 0 0 0 1 0 1 0 0 0 0 1 1 0 0 0 0 0 0 c 2 848 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 c 2 850 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 c 2 871 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 c 2 927 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 c 2 929 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 c 2 930 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 t 1 805 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 t 1 808 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 4 5 5 t 1 817 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 t 1 853 0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 t 1 872 0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 t 1 895 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 t 1 903 0 0 0 0 0 0 0 0 1 1 1 1 0 0 0 0 0 0 t 1 938 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 t 2 827 0 0 0 0 0 0 0 0 1 1 1 1 0 0 2 2 2 2 t 2 857 0 0 0 0 0 0 0 0 0 0 2 2 0 0 0 0 0 0 t 2 876 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 t 2 905 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 t 2 908 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 t 2 913 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 t 2 928 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

Hope that helps. -SJ

Back to: Top of message | Previous page | Main SAS-L page