Model is not identifiable...
summary from lme4 is returned
some computational error has occurred in lmerTest


Model is not identifiable...
Error in base::chol2inv(x, ...) : 'a' must be a numeric matrix

> Rstats <- read.csv("~/Dropbox/R stats.csv")
> Rstats
    Treatment  Sig Rep  Values
1         luc  low   1 3081182
2         luc  low   1 3067997
3         luc  low   2 4517539
4         luc  low   2 4862481
5         luc  low   3 4836874
6         luc  low   3 4687831
7         utr  low   1 2110208
8         utr  low   1 2147653
9         utr  low   2 2692765
10        utr  low   2 3222273
11        utr  low   3 1599556
12        utr  low   3 1601848
13        luc high   1 2046070
14        luc high   1 2046958
15        luc high   2 3053723
16        luc high   2 3054667
17        luc high   3 1969402
18        luc high   3 1950767
19        utr high   1 1178373
20        utr high   1 1100908
21        utr high   2 1297684
22        utr high   2 1379638
23        utr high   3 1136914
24        utr high   3 1131103
25         m1  low   1 2278695
26         m1  low   1 2406247
27         m1  low   1 2649573
28         m1  low   2 2578858
29         m1  low   2 2568366
30         m1  low   3 4955344
31         m1  low   3 4986120
32         m1 high   1 1280173
33         m1 high   1 1571484
34         m1 high   1 1822083
35         m1 high   2 1628416
36         m1 high   2 1573733
37         m1 high   3 1919091
38         m1 high   3 1724770
39         m2  low   1 1330070
40         m2  low   1 1174064
41         m2  low   1 1314733
42         m2  low   2 2369844
43         m2  low   2 2348595
44         m2  low   3 4477456
45         m2  low   3 4386348
46         m2 high   1  851196
47         m2 high   1  789824
48         m2 high   1  797485
49         m2 high   2  872469
50         m2 high   2  885174
51         m2 high   3 2137142
52         m2 high   3 2026491
53         m3  low   1 3633856
54         m3  low   1 4008909
55         m3  low   2 4066817
56         m3  low   2 4626628
57         m3  low   3 4063864
58         m3  low   3 4200093
59         m3  low   3 4131145
60         m3 high   1 1392872
61         m3 high   1 1775025
62         m3 high   2 2634198
63         m3 high   2 2646698
64         m3 high   3 2548206
65         m3 high   3 2666324
66         m4  low   1 3756921
67         m4  low   1 4285239
68         m4  low   1 3802821
69         m4  low   2 4297518
70         m4  low   2 4685541
71         m4  low   3 5510624
72         m4  low   3 5755805
73         m4 high   1 1202740
74         m4 high   1 1519717
75         m4 high   1 1069819
76         m4 high   2 2674748
77         m4 high   2 2658240
78         m4 high   3 3029938
79         m4 high   3 2893930
80         m5  low   1 1667185
81         m5  low   1 1794860
82         m5  low   2 2282347
83         m5  low   2 2481697
84         m5  low   3 4319992
85         m5  low   3 4314514
86         m5 high   1 1168027
87         m5 high   1 1104565
88         m5 high   2 1121619
89         m5 high   2 1721151
90         m5 high   3 1300998
91         m5 high   3 1320560
92         m6  low   1 3635387
93         m6  low   1 4013003
94         m6  low   2 5218195
95         m6  low   2 5271826
96         m6  low   3 4282878
97         m6  low   3 4381196
98         m6 high   1 2079993
99         m6 high   1 2264893
100        m6 high   1 1831285
101        m6 high   2 2989859
102        m6 high   2 3264232
103        m6 high   3 1959999
104        m6 high   3 1976519

> Anova1=lmer(Values~Treatment*Sig+(1|Rep), data=Rstats)
> summary(Anova1)
Model is not identifiable...
summary from lme4 is returned
some computational error has occurred in lmerTest
Linear mixed model fit by REML ['lmerMod']
Formula: Values ~ Treatment * Sig + (1 | Rep)
   Data: Rstats

REML criterion at convergence: 2630.4

Scaled residuals:
    Min      1Q  Median      3Q     Max
-1.7070 -0.7135 -0.0029  0.6699  2.4632

Random effects:
 Groups   Name        Variance  Std.Dev.
 Rep      (Intercept) 2.729e+11 522409
 Residual             3.707e+11 608835
Number of obs: 104, groups:  Rep, 3

Fixed effects:
                    Estimate Std. Error t value
(Intercept)          2353598     390833   6.022
Treatmentm1          -625863     338923  -1.847
Treatmentm2         -1077287     338923  -3.179
Treatmentm3           -76377     351511  -0.217
Treatmentm4          -121665     338923  -0.359
Treatmentm5         -1064111     351511  -3.027
Treatmentm6            66570     338923   0.196
Treatmentutr        -1149494     351511  -3.270
Siglow               1822053     351511   5.183
Treatmentm1:Siglow   -264417     479029  -0.552
Treatmentm2:Siglow   -530434     479029  -1.107
Treatmentm3:Siglow    -53345     488302  -0.109
Treatmentm4:Siglow    612995     479029   1.280
Treatmentm5:Siglow   -301440     497112  -0.606
Treatmentm6:Siglow    224860     488292   0.461
Treatmentutr:Siglow  -797106     497112  -1.603

Correlation matrix not shown by default, as p = 16 > 12.
Use print(x, correlation=TRUE)  or
     vcov(x)     if you need it

> lsmeans(Anova1, test.effs="Treatment")
Model is not identifiable...
Error in base::chol2inv(x, ...) : 'a' must be a numeric matrix