1

我在理解 glmmTMB 中权重的行为时遇到了一些麻烦,因为它导致与 lmer 无法比较的结果。

这个例子可以帮助理解问题:

m9a <- lmer(wi_gamma ~ level + (1 | TAG_IDc), weights=sub$pesi, data=sub)  
m9b <- lmer(wi_gamma ~ level + (1 | TAG_IDc), weights=sub$pesi/max(sub$pesi), data=sub)

这两个模型给出了相同的结果,因此我认为 lmer 将权重设为 1。

惠特TMB:

m9c <- glmmTMB(wi_gamma ~ level + (1 | TAG_IDc), weights=sub$pesi, data=sub)

输出类似于 m9a,但我怀疑这可能是偶然的,因为我们尝试了其他模型,但并非总是如此

m9d <- glmmTMB(wi_gamma ~ level + (1 | TAG_IDc), weights=sub$pesi/max(sub$pesi), data=sub)

这与 m9c 和 m9b 完全不同。我们无法弄清楚为什么...

我附上一个数据集以防有人可以提供帮助。

   wi_gamma level           TAG_IDc
1   1.0697674     3       i_16334_ALT
2   1.3043478     3       i_16336_ALT
3   1.4754098     3       i_16344_ALT
4   1.2592593     3       i_16350_ALT
5   1.0344828     3       i_16351_ALT
6   0.6422018     3       i_16355_ALT
7   0.9583333     3       i_16356_ALT
8   1.3129771     3       i_16363_ALT
9   0.4203152     3       i_16373_ALT
10  1.0714286     3      i_16438A_ALT
11  1.0576923     3       i_16490_ALT
12  0.9236234     3       i_16509_ALT
13  1.0210210     3       i_16523_ALT
14  0.9677419     3       i_16527_ALT
15  0.9181637     3       i_16564_ALT
16  1.4261745     3       i_16566_ALT
17  1.0242588     3       i_16579_ALT
18  1.3157895     3       i_16588_ALT
19  1.2761506     3       i_16680_ALT
20  1.5660685     3       i_16684_ALT
21  0.8371385     3       i_16685_ALT
22  3.3333333     3       i_17004_ALT
23  0.5693950     3 i_16318_CAPODARSO
24  0.5573770     3 i_16333_CAPODARSO
25  0.3292181     3 i_16343_CAPODARSO
26  0.1683502     3 i_16353_CAPODARSO
27  0.7380074     3 i_16357_CAPODARSO
28  0.6269592     3 i_16358_CAPODARSO
29  0.3125000     3 i_16364_CAPODARSO
30  0.3278689     3 i_16379_CAPODARSO
31  0.4575163     3      i_pic12_gela
32  0.6463878     3      i_pic13_gela
33  0.3048780     3      i_pic15_gela
34  0.2222222     3      i_pic16_gela
35  0.2811245     3      i_pic18_gela
36  0.3083700     3      i_pic19_gela
37  1.3274336     3       i_16191_GRA
38  1.0884354     3       i_16328_GRA
39  0.9302326     3       i_16340_GRA
40  1.1111111     3       i_16360_GRA
41  1.6853933     3       i_16361_GRA
42  1.3227513     3       i_16362_GRA
43  2.2660099     3       i_16365_GRA
44  0.6621005     3       i_16371_GRA
45  1.0185185     3       i_16375_GRA
46  1.5909091     3       i_16377_GRA
47  0.8737864     3       i_16438_GRA
48  2.3090278     3       i_16479_GRA
49  0.6951872     3       i_16520_GRA
50  1.1654135     3       i_16550_GRA
51  0.6250000     3       i_16558_GRA
52  1.0619469     3      i_16673A_GRA
53  2.0068027     3       i_17001_GRA
54  1.3461538     3       i_40401_MAT
55  1.1764706     3       i_40432_MAT
56  0.0010000     3       i_40441_MAT
57  1.8181818     3       i_40442_MAT
58  0.5494505     3       i_40449_MAT
59  0.4651163     3       i_40450_MAT
60  1.3636364     3       i_40468_MAT
61  3.0136986     3       i_40472_MAT
62  1.8918919     3       i_40475_MAT
63  0.2127660     3       i_40508_MAT
64  2.0000000     3       i_40514_MAT
65  2.7058824     3 i_41380_TARQUINIA
66  1.9780220     3 i_41397_TARQUINIA
67  0.9876543     3 i_41399_TARQUINIA
68  1.2941176     3 i_41429_TARQUINIA
69  1.2745098     3 i_41436_TARQUINIA
70  0.4444444     3 i_41469_TARQUINIA
71  2.9026193     2       i_16191_GRA
72  0.8282671     2 i_16318_CAPODARSO
73  1.4331228     2       i_16328_GRA
74  0.8658397     2 i_16333_CAPODARSO
75  1.6631572     2       i_16334_ALT
76  2.0049406     2       i_16336_ALT
77  0.7523895     2       i_16340_GRA
78  0.4074540     2 i_16343_CAPODARSO
79  1.0833450     2       i_16344_ALT
80  1.8238924     2       i_16350_ALT
81  1.7721126     2       i_16351_ALT
82  0.2646174     2 i_16353_CAPODARSO
83  0.6909334     2       i_16355_ALT
84  1.5339254     2       i_16356_ALT
85  1.1015470     2 i_16357_CAPODARSO
86  1.0525894     2 i_16358_CAPODARSO
87  1.9734806     2       i_16360_GRA
88  0.9523918     2       i_16361_GRA
89  1.4169294     2       i_16362_GRA
90  1.6450794     2       i_16363_ALT
91  0.1894661     2 i_16364_CAPODARSO
92  5.1850062     2       i_16365_GRA
93  1.1319997     2       i_16371_GRA
94  0.6967395     2       i_16373_ALT
95  1.1435281     2       i_16375_GRA
96  0.7470534     2       i_16377_GRA
97  0.4281720     2 i_16379_CAPODARSO
98  1.7035234     2       i_16438_GRA
99  1.7459570     2      i_16438A_ALT
100 5.5133775     2       i_16479_GRA
101 1.5966948     2       i_16490_ALT
102 1.4920488     2       i_16509_ALT
103 0.7903890     2       i_16520_GRA
104 0.7804563     2       i_16523_ALT
105 1.3696917     2       i_16527_ALT
106 3.0779570     2       i_16550_GRA
107 0.3343953     2       i_16558_GRA
108 1.3276080     2       i_16564_ALT
109 2.0664791     2       i_16566_ALT
110 1.5499096     2       i_16579_ALT
111 2.0303196     2       i_16588_ALT
112 1.7798461     2      i_16673A_GRA
113 1.3379980     2       i_16680_ALT
114 2.6767394     2       i_16684_ALT
115 1.1750567     2       i_16685_ALT
116 3.5022469     2       i_17001_GRA
117 2.3362075     2       i_17004_ALT
118 1.4023931     2       i_40401_MAT
119 1.1088690     2       i_40432_MAT
120 0.0010000     2       i_40441_MAT
121 2.6989452     2       i_40442_MAT
122 0.9312767     2       i_40449_MAT
123 0.3017808     2       i_40450_MAT
124 1.5893789     2       i_40468_MAT
125 2.9138612     2       i_40472_MAT
126 3.8217888     2       i_40475_MAT
127 0.1440525     2       i_40508_MAT
128 5.3833800     2       i_40514_MAT
129 3.7664419     2 i_41380_TARQUINIA
130 3.0025754     2 i_41397_TARQUINIA
131 2.6042539     2 i_41399_TARQUINIA
132 2.2046273     2 i_41429_TARQUINIA
133 1.8770580     2 i_41436_TARQUINIA
134 0.7418886     2 i_41469_TARQUINIA
135 0.6431548     2      i_pic12_gela
136 0.9083064     2      i_pic13_gela
137 0.6377614     2      i_pic15_gela
138 0.0704814     2      i_pic16_gela
139 1.5867438     2      i_pic18_gela
140 1.8477214     2      i_pic19_gela
4

0 回答 0