Slightly off topic, but related :-) [Design Issues]

posted by PharmCat  – Russia, 2019-12-09 16:13 (1789 d 16:23 ago) – Posting: # 20966
Views: 8,466

❝ But would we really need a mixed model ??


Hello, ElMaestro!

Look at residuals:

resid(M1)
         1          2          3          4          5          6          7          8          9         10
-0.1540073  0.1540073 -0.4297153  0.4297153  0.1540073 -0.1540073 -1.1774560  1.1774560  1.6071713 -1.6071713
resid(M2)
            1             2             3             4             5             6             8             9            10
-1.540073e-01  1.540073e-01 -1.018443e+00  1.018443e+00  1.540073e-01 -1.540073e-01 -3.608225e-16  1.018443e+00 -1.018443e+00


This is how I understand:

For observation 8 we have -3.608225e-16, I think, that mean, that this subject affect on "intra-variation" and make it less, and you can see SE is smaller. "Intra-individual part" went to coefficient and began part of inter-individual. Brr.. I don't know how to say with сlever words.

In lm / glm we see each observation as at statistically independent. And with lm we make smart trick when "storing" inter-individual variation in model coefficients, and exclude it from model error. Missing observation violate this system.

In mixed model we have another situation - all observation of subject is one statistically independent observation - realization of multidimentional variable. In the case of the ML / REML estimation, such a significant transition of one component of the variation to another does not occur in the case when the data contains missing values. ML and REML variation estimates all biased, but less biased than lm with missing data.

Complete thread:

UA Flag
Activity
 Admin contact
23,279 posts in 4,889 threads, 1,657 registered users;
70 visitors (0 registered, 70 guests [including 10 identified bots]).
Forum time: 08:37 CET (Europe/Vienna)

Pharmacokinetics: one of the magic arts of divination
whereby needles are stuck into dummies in an attempt
to predict profits.    Stephen Senn

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5