LISTSERV at the University of Georgia
Menubar Imagemap
Home Browse Manage Request Manuals Register
Previous messageNext messagePrevious in topicNext in topicPrevious by same authorNext by same authorPrevious page (March 2009, week 1)Back to main SAS-L pageJoin or leave SAS-L (or change settings)ReplyPost a new messageSearchProportional fontNon-proportional font
Date:   Sun, 1 Mar 2009 16:03:00 -0800
Reply-To:   sej.in.md@GMAIL.COM
Sender:   "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:   sej.in.md@GMAIL.COM
Organization:   http://groups.google.com
Subject:   Re: PROC GLIMMIX: 2 questions: blocks, and covariance matrices
Comments:   To: sas-l@uga.edu
Content-Type:   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