jabriges
☆    

US,
2020-06-24 21:46
(1372 d 19:16 ago)

Posting: # 21577
Views: 2,861
 

 random effects in BEAR (lme_lm.mod.r) [🇷 for BE/BA]

I have a question about getting subject level random effect. The call of lme() in lme_lm.mod.r uses

random=~drug - 1|subj

I'm not seeing how this gives you a random effect for each subject. I understand the -1 removes the intercept, but doesn't the inclusion of drug there make the random effect differ by formulation, and hence period (for the same subject)?
Any insight appreciated.

-Jay
ElMaestro
★★★

Denmark,
2020-06-25 11:40
(1372 d 05:22 ago)

@ jabriges
Posting: # 21578
Views: 2,252
 

 random effects in BEAR (lme_lm.mod.r)

Hi all,

random=~drug - 1|subj


I do not offer an answer to the question at all, only I wish to say that I have no ideas, generally, how syntactically R works for lme and lmer when specifying random effects.
All the stuff that involves "~", "|" for random effects is entirely obscure to me, and I tried to read the various online resources for it, and help files for packages. I even ordered the book by Pinheiro & Bates to learn from it, but that book assumes that the reader knows everything there is to about mixed models and R syntax before you open it on page 1. It is also not clear to me how covariance matrix structures are actually set using the R functions.

So, I think the question you asked jabriges is also highly relevant for me in a very general sense. I hope someone will chime in an educate me about the syntax.

Pass or fail!
ElMaestro
jabriges
☆    

US,
2020-06-25 20:12
(1371 d 20:50 ago)

@ ElMaestro
Posting: # 21579
Views: 2,231
 

 R formulas

Hi, Maestro.

❝ So, I think the question you asked jabriges is also highly relevant for me in a very general sense. I hope someone will chime in an educate me about the syntax.


There are places that you can read about R's formula syntax in general, like here.

My question is specifically about the formula for random (see my OP) which seems to make the random coefficient a slope on the drug term, rather than a random intercept term (grouped by subject, of course), which was my understanding of the intended model.

-Jay
UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,638 registered users;
69 visitors (0 registered, 69 guests [including 9 identified bots]).
Forum time: 16:02 CET (Europe/Vienna)

Nothing shows a lack of mathematical education more
than an overly precise calculation.    Carl Friedrich Gauß

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