# Testing if there is a difference in % proliferation between and 6h # and 10h BrdU labeling - Halisarca caerulea library(lmerTest) D = read.csv('prolif_6vs10.csv') D$hrs = as.factor(D$hrs) D$indiv = as.factor(D$indiv) M0 <-lmer( prolC ~ 1 + (1|indiv), data=D, REML=F ) Mh <-lmer( prolC ~ hrs + (1|indiv), data=D, REML=F ) summary(Mh) # The result from the lmer fit [ 'summary(Mh)' ] # shows that 'hours' is not a significant effect at all. # # Linear mixed model fit by maximum likelihood ['merModLmerTest'] # Formula: prolC ~ hrs + (1 | indiv) # Data: D # # AIC BIC logLik deviance # 137.4920 142.2043 -64.7460 129.4920 # # Random effects: # Groups Name Variance Std.Dev. # indiv (Intercept) 4.174 2.043 # Residual 9.810 3.132 # Number of obs: 24, groups: indiv, 8 # # Fixed effects: # Estimate Std. Error t value Pr(>|t|) # (Intercept) 17.6500 1.1139 15.846 2.52e-07 *** # hrs10 -0.7767 2.2278 -0.349 0.736 # --- # Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 # # Correlation of Fixed Effects: # (Intr) # hrs10 -0.500