[R-sig-ME] intercept and threshold values to calculate probability based on mixor package
K Imran M
drki@mu@@ @ending from gm@il@com
Tue Oct 30 07:49:57 CET 2018
Let me improve the question. On page 209, the authors wrote this for
Pr(Yij=1);
Pr(Yij=1) = 1/{1 + exp[7.31 + ....]} . Unfortunately. the values given for
the rest of the formula for Yij=2 and Yij=3.
My questions are:
1) Why v1 is 7.31 and not -7.31
2) What are the values for v2 and v3
Thanks David for your reply. Can you be more specific?
Best wishes,
Kamarul
On Tue, Oct 30, 2018 at 1:18 PM K Imran M <drki.musa using gmail.com> wrote:
> Hi everyone,
>
> I need help to identify the values for intercept and thresholds in the
> equations to calculate the probability (based on cumulative link models).
>
> On page 209 in Longitudinal Data Analysis book by Hedeker, the
> probabilities are written as:
>
> Pr(Yij=1) = 1/{1 + exp[-(v1 -(xijB+zijV))]}
> Pr(Yij=2) = 1/{1 + exp[-(v2 -(xijB+zijV))]} - Pr(Yij=1)
> Pr(Yij=3) = 1/{1 + exp[-(v3 -(xijB+zijV))]} - Pr(Yij=2)
> Pr(Yij=4) = 1 - 1/{1 + exp[-(v3 -(xijB+zijV))]}
>
> My question is what are the values for v1, v2 and v3?
>
> The R scripts below shows the codes and the results from running mixor.
>
> library(mixor)
> #> Loading required package: survival
> data("schizophrenia")
> mixor.mod <- mixor(imps79o ~ TxDrug + SqrtWeek + TxSWeek,
> data = schizophrenia, id = id, which.random.slope =
> 2, link = "logit")
> summary(mixor.mod)
>
> #>
> #> Estimate Std. Error z value P(>|z|)
> #> (Intercept) 7.318831 0.480778 15.2229 < 2.2e-16 ***
> #> SqrtWeek -0.882261 0.234568 -3.7612 0.0001691 ***
> #> TxDrug 0.057917 0.399102 0.1451 0.8846178
> #> TxSWeek -1.694861 0.268131 -6.3210 2.598e-10 ***
> #> (Intercept) (Intercept) 6.997646 1.369273 5.1105 3.213e-07 ***
> #> (Intercept) SqrtWeek -1.508514 0.536023 -2.8143 0.0048888 **
> #> SqrtWeek SqrtWeek 2.008916 0.453587 4.4290 9.469e-06 ***
> #> Threshold2 3.901260 0.213257 18.2937 < 2.2e-16 ***
> #> Threshold3 6.507172 0.289991 22.4392 < 2.2e-16 ***
> #> ---
> #> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
>
> Thank you very much.
>
> Best wishes,
>
> Kamarul Imran Musa
> MD PhD
> Universiti Sains Malaysia
> Malaysia
>
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list