################################################################################################################# #Model to explicitly test for differences in dominance between West and bonsall creek measures (Ucrit in this example) ################################################################################################################# ucrit<-read.csv(file.choose(),header=TRUE,stringsAsFactors = FALSE, na.strings="") #Read in Data set ucrit$crosstype <-as.factor(ucrit$crosstype)#set types of data ucrit$family <-as.factor(ucrit$family) ucrit$indv<-as.factor(ucrit$indv) ucrit$additive<-as.numeric(ucrit$additive) ucrit$dom<-as.numeric(ucrit$dom) ucrit$location<-as.factor(ucrit$location) library(nlme) #Load the nlme package- should treat data as nested z<-lme(ucrit~location+additive, random=~1|family, data=ucrit) summary(z) # parameter estimates and overall model fit y<-lme(ucrit~location+additive+dom, random=~1|family, data=ucrit) summary(y) # parameter estimates and overall model fit w<-lme(ucrit~location+additive+dom+location*dom, random=~1|family, data=ucrit) summary(w) t<-lme(ucrit~location+additive+dom+location*dom+additive*location, random=~1|family, data=ucrit) summary(t)