[R-sig-ME] c++ exception (unknown reason) when using an offset of the slope with glmer
xavier.paoletti at curie.fr
xavier.paoletti at curie.fr
Tue Jan 27 14:24:01 CET 2015
Thanks for your attention.
Here there are. 180 rows, 4 columns + obs number.
Obs patid cycle dose DLT
1 1 1 4.1 0
2 2 1 4.8 0
3 3 1 5.3 0
4 4 1 5.7 0
5 5 1 6.0 1
6 6 1 5.7 1
7 7 1 5.7 0
8 8 1 5.3 0
9 9 1 5.3 0
10 10 1 5.7 0
11 11 1 5.3 0
12 12 1 5.3 0
13 13 1 5.3 0
14 14 1 5.3 0
15 15 1 5.7 0
16 16 1 5.3 0
17 17 1 5.3 0
18 18 1 5.3 0
19 19 1 5.3 0
20 20 1 5.7 0
21 21 1 5.7 1
22 22 1 5.7 1
23 23 1 5.7 0
24 24 1 5.3 0
25 25 1 5.3 0
26 26 1 5.3 0
27 27 1 5.3 1
28 28 1 5.7 0
29 29 1 5.7 1
30 30 1 5.7 0
31 1 2 4.1 0
32 2 2 4.8 0
33 3 2 5.3 0
34 4 2 5.7 0
35 5 2 6.0 0
36 6 2 5.7 0
37 7 2 5.7 1
38 8 2 5.3 0
39 9 2 5.3 0
40 10 2 5.7 0
41 11 2 5.3 0
42 12 2 5.3 0
43 13 2 5.3 0
44 14 2 5.3 0
45 15 2 5.7 0
46 16 2 5.3 0
47 17 2 5.3 0
48 18 2 5.3 0
49 19 2 5.3 0
50 20 2 5.7 0
51 21 2 5.7 0
52 22 2 5.7 0
53 23 2 5.7 1
54 24 2 5.3 0
55 25 2 5.3 0
56 26 2 5.3 0
57 27 2 5.3 0
58 28 2 5.7 0
59 29 2 5.7 0
60 30 2 5.7 1
61 1 3 4.1 0
62 2 3 4.8 0
63 3 3 5.3 0
64 4 3 5.7 0
65 5 3 6.0 1
66 6 3 5.7 1
67 7 3 5.7 1
68 8 3 5.3 0
69 9 3 5.3 0
70 10 3 5.7 0
71 11 3 5.3 0
72 12 3 5.3 0
73 13 3 5.3 0
74 14 3 5.3 0
75 15 3 5.7 0
76 16 3 5.3 0
77 17 3 5.3 0
78 18 3 5.3 0
79 19 3 5.3 1
80 20 3 5.7 0
81 21 3 5.7 0
82 22 3 5.7 0
83 23 3 5.7 1
84 24 3 5.3 0
85 25 3 5.3 0
86 26 3 5.3 1
87 27 3 5.3 0
88 28 3 5.7 0
89 29 3 5.7 0
90 30 3 5.7 1
91 1 4 4.1 0
92 2 4 4.8 0
93 3 4 5.3 0
94 4 4 5.7 0
95 5 4 6.0 0
96 6 4 5.7 0
97 7 4 5.7 0
98 8 4 5.3 0
99 9 4 5.3 0
100 10 4 5.7 0
101 11 4 5.3 1
102 12 4 5.3 1
103 13 4 5.3 0
104 14 4 5.3 0
105 15 4 5.7 1
106 16 4 5.3 0
107 17 4 5.3 0
108 18 4 5.3 0
109 19 4 5.3 0
110 20 4 5.7 0
111 21 4 5.7 0
112 22 4 5.7 0
113 23 4 5.7 0
114 24 4 5.3 0
115 25 4 5.3 0
116 26 4 5.3 0
117 27 4 5.3 0
118 28 4 5.7 0
119 29 4 5.7 0
120 30 4 5.7 0
121 1 5 4.1 0
122 2 5 4.8 0
123 3 5 5.3 0
124 4 5 5.7 0
125 5 5 6.0 0
126 6 5 5.7 1
127 7 5 5.7 0
128 8 5 5.3 1
129 9 5 5.3 0
130 10 5 5.7 0
131 11 5 5.3 0
132 12 5 5.3 0
133 13 5 5.3 1
134 14 5 5.3 0
135 15 5 5.7 0
136 16 5 5.3 0
137 17 5 5.3 0
138 18 5 5.3 0
139 19 5 5.3 0
140 20 5 5.7 0
141 21 5 5.7 1
142 22 5 5.7 0
143 23 5 5.7 0
144 24 5 5.3 0
145 25 5 5.3 0
146 26 5 5.3 0
147 27 5 5.3 0
148 28 5 5.7 0
149 29 5 5.7 0
150 30 5 5.7 0
151 1 6 4.1 0
152 2 6 4.8 0
153 3 6 5.3 0
154 4 6 5.7 1
155 5 6 6.0 1
156 6 6 5.7 0
157 7 6 5.7 0
158 8 6 5.3 0
159 9 6 5.3 0
160 10 6 5.7 1
161 11 6 5.3 0
162 12 6 5.3 0
163 13 6 5.3 0
164 14 6 5.3 0
165 15 6 5.7 0
166 16 6 5.3 0
167 17 6 5.3 0
168 18 6 5.3 0
169 19 6 5.3 0
170 20 6 5.7 1
171 21 6 5.7 1
172 22 6 5.7 0
173 23 6 5.7 0
174 24 6 5.3 0
175 25 6 5.3 0
176 26 6 5.3 1
177 27 6 5.3 0
178 28 6 5.7 0
179 29 6 5.7 1
180 30 6 5.7 0
Ben Bolker <bbolker at gmail.com>
27/01/2015 14:14
A
<xavier.paoletti at curie.fr>
cc
"r-sig-mixed-models at r-project.org" <r-sig-mixed-models at r-project.org>
Objet
Re: [R-sig-ME] c++ exception (unknown reason) when using an offset of the
slope with glmer
Your posted data set got removed by the mailing list machinery. Can
you post it somewhere publicly accessible?
On Tue, Jan 27, 2015 at 3:45 AM, <xavier.paoletti at curie.fr> wrote:
>
> Dear all,
>
> I am new on this forum and I hope my request follows the right format.
>
> I use R 3.1.2 and lme1.7 on a mac OS X (snow leopard) or Windows OS.
>
> I try to fit a longitudinal logistic mixed effect model on dose-time
> response data where the response is measured several times at the same
> dose. The probability of response increases with the dose.
> There is a set of discrete doses (let's say 6) but most of the data are
> measured at 1 or 2 doses.
>
> I use a very simple logstic model with a random intercept and the dose
> effect.
> In the simplest case, there is not time effect.
> Furthermore, I would like to set the slope of the dose to some value
using
> offset slope=1.5.
>
> The command line,
> glmer(DLTb ~ offset(slope*dose) + (1 | patid),
> family=binomial,data=dataAllCRM,nAGQ=10)
> gives the following error:
> Under Mac OSX: (function (fr, X, reTrms, family, nAGQ = 1L, verbose =
0L,
> control = glmerControl(), :
> c++ exception (unknown reason)
>
> Under Windows: Error: (maxstephalfit) PIRLS step-halvings failed to
reduce
> deviance in pwrssUpdate
>
> Whatever the value of the offset, I get the same error.
> If I remove the offset or if I remove the variable, some estimates are
> obtained.
>
> Please find attached an example of dataset as an illustration; I get the
> same error for all tested datasets.
>
>
> Thank you very much for your help.
>
> Best regards,
>
> Xavier
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
-------------- section suivante --------------
Un texte encapsulé et encodé dans un jeu de caractères inconnu a été nettoyé...
Nom : dataEx_Glmer.txt
URL : <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20150127/3652cc70/attachment.txt>
More information about the R-sig-mixed-models
mailing list