[R-sig-ME] intercept and threshold values to calculate probability based on mixor package
roee maor
roeem@or @ending from gm@il@com
Thu Nov 1 11:07:08 CET 2018
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]]
More information about the R-sig-mixed-models
mailing list