intercept and threshold values to calculate probability based on mixor package
Hi Kamarul, I'm not familiar with this specific package so I don't know how to interpret this output. I can see the standard line for "(Intercept)", but the there is another line for "(Intercept)(Intercept)" which means nothing to me. It makes no sense to me to add/subtract intercepts and threshold values because intercept is a value on the y-axis and the thresholds are values on the x-axis. I'm putting the list back in the cc to potentially reach others who could help. All the best, Roi
On Thu, 1 Nov 2018 at 03:29, K Imran M <drki.musa at gmail.com> wrote:
Dear Roi, Really appreciate your response. Perhaps the first v1 equals 0 minus the intercept (7.31). That gives rise to v1= -7.31. How about v2 and v3. Can I say v2 = 3.90 and v3 = 6.51? Best wishes, Kamarul On Tue, Oct 30, 2018 at 5:40 PM roee maor <roeemaor at gmail.com> wrote:
Hi Kamarul,
I'm not familiar with the mixor package, but looking at your model's
output
it seems that the threshold values you're looking for are Threshold2 and
Threshold3.
The first and last thresholds in any threshold model are minus infinity
and
infinity, respectively. In your current model, that leaves 3 thresholds to
be recovered. To simplify calculations, one of these is assumed to be 0,
which leaves only two thresholds to estimate from the data. These two are
the last two parameters listed in your model's output.
However, I'm not familiar with package mixor, so I can't say which one of
the three thresholds is the one set to 0. In MCMCglmm it is the second
threshold, but I don't think that has to be so by necessity. Finding out
which of the thresholds is 0 may (or may not) clarify why you get a
positive value where you expect a negative one.
Best,
Roi Maor
[[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
Roi Maor PhD candidate School of Zoology, Tel Aviv University Centre for Biodiversity and Environment Research, UCL [[alternative HTML version deleted]]