Hi Simone,
Thanks again for your reply.
simone77 wrote:There are a lot of things that can go wrong. It's hard to help you without knowing more about your study. Perhaps it would help if you could copy your GEPAT here, how are the data formatted? Headed format? inp? What is the range of values of your individual covariate? are they numerical values? what are the GEMACO sentences you are using for the other parameters?
In this study, I am trying to estimate some vital rates and pair-bond disruption rates. In the transition part, Step 1 models survival, Step 2 models widowhood, Step 3 models breeding probability and Step 4 models divorce (marked in red color with questions).
In Step 4, I am trying to test whether the divorce rate is associated with individual personality. So I use personality score as individual covariates.
All other steps work well since this is the best model after a step down model selection. The problem happens when I tried to test the individual covariates in step 4 of transition.
I pasted my GEPAT and GEMACO code here and mark the session with questions in red color.
-------- INIT -------------
1
= STEP 1 =
1
to(1,2,3,4:7)
0
1 7
p p p p p p *
--------- TRANSITION --------
4
= STEP 1 =
1
from(1 4, 2 3 5, 6 7).t(1:9,10:19,20:29,30:39,40:45)
0
8 8
s - - - - - - *
- s - - - - - *
- - s - - - - *
- - - s - - - *
- - - - s - - *
- - - - - s - *
- - - - - - s *
- - - - - - - *
= STEP 2 =
1
from(1,2,3, 4 5).t(1:9,10:19,20:29,30:39,40:45)
0
8 8
* - - - - - w -
- * - - - - w -
- - * - - - w -
- - - * - - w -
- - - - * - w -
- - - - - * - -
- - - - - - * -
- - - - - - - *
= STEP 3 =
1
from.t(1:19,20_45)
0
8 12
b - - - - * - - - - - -
- b - - - - * - - - - -
- - b - - - * - - - - -
- - - b - * - - - - - -
- - - - b - * - - - - -
- - - - - - - b - * - -
- - - - - - - - b - * -
- - - - - - - - - - - *
= STEP 4 =
1
from.t(1:45).g(1)+from.t(1:45).g(2).[i+xind]0
12 8
* d - - - - - -
* d - - - - - -
* d - - - - - -
* d - - - - - -
* d - - - - - -
- - - * - d - -
- - - - * d - -
- * - - - - - -
- - * - - - - -
- - - - - * - -
- - - - - - * -
- - - - - - - *
---------- CAPTURE -----------
1
= STEP 1 =
1
firste.[from(4 5).to(9)+to(1 12 13, 2 3 7,4,5,6,10,11)+from(6,7).to(9)+from(1,2,3).to(8)]+nexte.[from(4 5,6,7).to(9)+from(1,2 3).to(7)+to(2,3,4,6,10,11,12,13)+from(1,2).to(5)+from(1,2,3).to(8)]+nexte.[from(4 5,7).to(9)+from(1,2 3).to(7)+to(2,3,4,6,10,11,12,13)+from(1,2).to(5)+from(1,2,3).to(8)].t(2_46)*x(1)+nexte.[from(4 5,6,7).to(9)+from(1,2 3).to(7)+to(2,3,4,6,10,11,12,13)+from(1,2).to(5)+from(1,2,3).to(8)].t(2_46)*x(2)
0
8 13
* p - - p - p p - - - - -
* - p - p p p p - - - - -
* - - p - - p p - - - - -
* - - - - - - - p - - - -
* - - - - - - - p - - - -
* - - - - - - - p p - p -
* - - - - - - - p - p - p
* - - - - - - - - - - - -
simone77 wrote:It's not clear to me if you have empty spaces in your dataset for a group (which would create problems).
In general, it is recommended that the values of the individual covariates be standardized (0 mean, unit variance).
I am using the head format for data input and I pasted part of my input data here.
- Code: Select all
H: H: H: H: H: H: S: $COV:per COV:personality
0 1 4 6 3 0 1 2 6.648117
0 1 0 1 0 0 1 2 4.684725
2 0 0 1 0 0 1 1 0
1 1 0 1 0 1 1 2 4.519807
1 0 1 0 1 0 1 2 4.519807
0 1 0 1 0 1 1 1 0
$COV:per is used to define groups. And for individuals in group 1, we don't have personality scores, so I put 0 in COV:personality. I tried to standerdize my personality score, and it did not help. Would be appreciate if you can help me with that.
Looking forward to the reply. Many thanks in advance