Yes, but.... [🇷 for BE/BA]

posted by ElMaestro  – Denmark, 2024-02-23 18:41 (441 d 01:58 ago) – Posting: # 23878
Views: 3,051

Hi Mittyri,

❝ I have nothing to add to this: that's a price for identifiability. The only solution I see besides you already gave is to turn off the contrasts as shown there:

f <- sample(gl(3, 4, labels = letters[1:3]))

❝ # [1] c a a b b a c b c b a c

❝ #Levels: a b c


❝ g <- sample(gl(3, 4, labels = LETTERS[1:3]))

❝ # [1] A B A B C B C A C C A B

❝ #Levels: A B C

❝ X0 <- model.matrix(~ f + g, contrasts.arg = list(

❝                    f = contr.treatment(n = 3, contrasts = FALSE),
❝                    g = contr.treatment(n = 3, contrasts = FALSE)))

:-D:-D
Nice, but are you aware that this model matrix is also not full rank, but rank-deficient? A cardinal question was how to get the (or a) full rank matrix.
The rank of your X0 is 5 but there are 7 columns. Seems you went to the opposite extreme of full rank, i.e. one column for every level of factors, plus intercept. Deficiency here = 2 redundant columns.
Anyways, I guess I will just live with it and do my stuff in a dumb but necessary fashion. Keep me posted if you find a nifty way of generating a full rank matrix model, please.:-)

Pass or fail!
ElMaestro

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,670 registered users;
117 visitors (0 registered, 117 guests [including 8 identified bots]).
Forum time: 21:39 CEST (Europe/Vienna)

Freedom is always and exclusively
freedom for the one
who thinks differently.    Rosa Luxemburg

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