[R-sig-ME] Analyzing similarity scores between subjects

Han Zhang h@nzh @ending from umich@edu
Wed Aug 8 19:22:59 CEST 2018

Hi all,

I have a modeling problem involving similarity scores between subjects.
During 4 time points in my experiment, I sampled eye movements of my
subjects. At each time point, subjects had either one of two different
states, Y or N. I have no control of the state, it is purely observational.
My data produces 4 similarity matrices - for each sampling, every subject
was compared to every other subject on some similarity measure of eye
movements (self-comparisons excluded). Each matrix contains three types of
comparison: N-N, N-Y, and Y-Y. My hypothesis is that the eye movements of
those in state N were more similar to each other, compared to N-Y, or Y-Y.
So N-N > N-Y or Y-Y.

I came up with a model like this:

lmer(dist ~ type + (1|sub_i) + (1|sub_i:type) + (1|segment) +
(1|segment:type) + (1|sub_i: segment) + (1|sub_i: segment:type), data,

where dist is the similarity score, type is a 3-level factor (n-n, n-y,
y-y), sub_i is subject ID, segment is sample ID. I was
trying to build a model with a "maximal" random structure.

Have I correctly specified my model? I have two concerns:
(1) because any given data point in the matrix belongs to two subjects, i
and j, should I include random effects for both subject i and subject j?

(2) Becuase each matrix is symmetrical, I am duplicating my data in the
above model. Should I use only the unique pairwise comparisons and do
something like this:

lmer(dist ~ type + (1|segment) + (1|segment:type), half_data, REML=F)


Han Zhang
Graduate Student
Combined Program in Education and Psychology
University of Michigan, Ann Arbor
Email: hanzh using umich.edu
Phone: 1-734-680-6031

	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list