[R] help with binary outcome

markus eichel Wed, 04 Feb 2009 04:36:13 -0800

conducted a training with 2 training groups and one control-group (TG: Training 
Groups 1-3 in my data). After that training I checked every week within 4 weeks 
(Zeit= 1-4 in my data), if the participants used the training skills (FLEXPV). 
unfortunately I didnt use any real scale for this item, so it´s only binary 
coded 0/1 0=not put in practice 1=put in practice.


I´m interessted if there are any variances in the slopes of the persons within 
the training groups. 

First I tried do do it with the nlme-package function:

>model.1=gls(FLEXPV~ZEIT,data=data,na.action=na.omit)

>model.2=lme(FLEXPV~ZEIT,random=~1|ID,data=data,na.action=na.omit)
>model.3=lme(FLEXPV~ZEIT,random=~1+ZEIT|ID,data=data,na.action=na.omit)
 anova(model.2, model.3)

Then I was wondering if this is the right function for my binary datas. 
Therefore I tried the following thing with the MASS package

>mod.2<-glmmPQL(FLEXPV~ZEIT+TG,random=~1|ID, 
>family="binomial",data=data,na.action=na.omit)
 
>mod.3<-glmmPQL(FLEXPV~ZEIT*TG,random=~1|ID, 
>family="binomial",data=data,na.action=na.omit)

But I couldn´t compare mod.2 vs. mod.3 with the anova. And how can I build a 
level-1 model ??

I do attache the text belwo, so that you can get an impression of my data.

 

I really appriciate your help and thank you in advance for helping me.

Regards Markus Eichel

 


> data
 ID TG X ZEIT X.1 FLEXPV
1 1 1 NA 0 NA 0
2 1 1 NA 1 NA 1
3 1 1 NA 2 NA 1
4 1 1 NA 3 NA 1
5 2 1 NA 0 NA 0
6 2 1 NA 1 NA 0
7 2 1 NA 2 NA 0
8 2 1 NA 3 NA 0
9 3 1 NA 0 NA 0
10 3 1 NA 1 NA 1
11 3 1 NA 2 NA 1
12 3 1 NA 3 NA 1
13 4 1 NA 0 NA 0
14 4 1 NA 1 NA 0
15 4 1 NA 2 NA 0
16 4 1 NA 3 NA 0
17 5 1 NA 0 NA 0
18 5 1 NA 1 NA 1
19 5 1 NA 2 NA 0
20 5 1 NA 3 NA 0
21 6 1 NA 0 NA 0
22 6 1 NA 1 NA 0
23 6 1 NA 2 NA 0
24 6 1 NA 3 NA 0
25 7 1 NA 0 NA 0
26 7 1 NA 1 NA 0
27 7 1 NA 2 NA 0
28 7 1 NA 3 NA 0
29 8 1 NA 0 NA 0
30 8 1 NA 1 NA 0
31 8 1 NA 2 NA 1
32 8 1 NA 3 NA 0
33 9 1 NA 0 NA 0
34 9 1 NA 1 NA 1
35 9 1 NA 2 NA 1
36 9 1 NA 3 NA 0
37 10 1 NA 0 NA 0
38 10 1 NA 1 NA 0
39 10 1 NA 2 NA 1
40 10 1 NA 3 NA 1
41 11 1 NA 0 NA 0
42 11 1 NA 1 NA 1
43 11 1 NA 2 NA 1
44 11 1 NA 3 NA 1
45 12 1 NA 0 NA 1
46 12 1 NA 1 NA 0
47 12 1 NA 2 NA 0
48 12 1 NA 3 NA 0
49 13 1 NA 0 NA 0
50 13 1 NA 1 NA 1
51 13 1 NA 2 NA 1
52 13 1 NA 3 NA 1
53 14 2 NA 0 NA 0
54 14 2 NA 1 NA 0
55 14 2 NA 2 NA 0
56 14 2 NA 3 NA 0
57 15 2 NA 0 NA 0
58 15 2 NA 1 NA 0
59 15 2 NA 2 NA 0
60 15 2 NA 3 NA 0
61 16 2 NA 0 NA 0
62 16 2 NA 1 NA 1
63 16 2 NA 2 NA 0
64 16 2 NA 3 NA 0
65 17 2 NA 0 NA 0
66 17 2 NA 1 NA 0
67 17 2 NA 2 NA 0
68 17 2 NA 3 NA 0
69 18 2 NA 0 NA 0
70 18 2 NA 1 NA 1
71 18 2 NA 2 NA 0
72 18 2 NA 3 NA 0
73 19 2 NA 0 NA 1
74 19 2 NA 1 NA 1
75 19 2 NA 2 NA 0
76 19 2 NA 3 NA 1
77 20 2 NA 0 NA 1
78 20 2 NA 1 NA 0
79 20 2 NA 2 NA 1
80 20 2 NA 3 NA 0
81 21 2 NA 0 NA 0
82 21 2 NA 1 NA 0
83 21 2 NA 2 NA 0
84 21 2 NA 3 NA 0
85 22 2 NA 0 NA 0
86 22 2 NA 1 NA 1
87 22 2 NA 2 NA 0
88 22 2 NA 3 NA 0
89 23 2 NA 0 NA 0
90 23 2 NA 1 NA 0
91 23 2 NA 2 NA 0
92 23 2 NA 3 NA 0
93 24 2 NA 0 NA 1
94 24 2 NA 1 NA 1
95 24 2 NA 2 NA 1
96 24 2 NA 3 NA 1
97 25 2 NA 0 NA 0
98 25 2 NA 1 NA 0
99 25 2 NA 2 NA 0
100 25 2 NA 3 NA 1
101 26 2 NA 0 NA 0
102 26 2 NA 1 NA 0
103 26 2 NA 2 NA 0
104 26 2 NA 3 NA 0
105 27 3 NA 0 NA 0
106 27 3 NA 1 NA 0
107 27 3 NA 2 NA 0
108 27 3 NA 3 NA 0
109 28 3 NA 0 NA 1
110 28 3 NA 1 NA 0
111 28 3 NA 2 NA 0
112 28 3 NA 3 NA 0
113 29 3 NA 0 NA 1
114 29 3 NA 1 NA 1
115 29 3 NA 2 NA 1
116 29 3 NA 3 NA 1
117 30 3 NA 0 NA 1
118 30 3 NA 1 NA 1
119 30 3 NA 2 NA 1
120 30 3 NA 3 NA 1
121 31 3 NA 0 NA 1
122 31 3 NA 1 NA 1
123 31 3 NA 2 NA 0
124 31 3 NA 3 NA 0
125 32 3 NA 0 NA 1
126 32 3 NA 1 NA 0
127 32 3 NA 2 NA 1
128 32 3 NA 3 NA 1
129 33 3 NA 0 NA 0
130 33 3 NA 1 NA 1
131 33 3 NA 2 NA 1
132 33 3 NA 3 NA 1
133 34 3 NA 0 NA 0
134 34 3 NA 1 NA 1
135 34 3 NA 2 NA 1
136 34 3 NA 3 NA 1
137 35 3 NA 0 NA 1
138 35 3 NA 1 NA 0
139 35 3 NA 2 NA 0
140 35 3 NA 3 NA 1
141 36 3 NA 0 NA 1
142 36 3 NA 1 NA 1
143 36 3 NA 2 NA 0
144 36 3 NA 3 NA 0
145 37 3 NA 0 NA 1
146 37 3 NA 1 NA 1
147 37 3 NA 2 NA 0
148 37 3 NA 3 NA 1
149 38 3 NA 0 NA 1
150 38 3 NA 1 NA 0
151 38 3 NA 2 NA 0
152 38 3 NA 3 NA 0


        
Deutschlands größte Online-Videothek schenkt Ihnen 12.000 Videos!* 
*http://entertainment.web.de/de/entertainment/maxdome/index.html* 
[http://entertainment.web.de/de/entertainment/maxdome/index.html]

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] help with binary outcome

Reply via email to