[R-sig-ME] intercept and threshold values to calculate probability based on mixor package

[K Imran M]

Dear Roi,

Really appreciate it.

Regards,
Kamarul


On Thu, Nov 1, 2018 at 6:06 PM roee maor <roeemaor using gmail.com> wrote:

> 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 using 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 using 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 using 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]]