Bioequivalence and Bioavailability Forum

keep it simple! [🇷 for BE/BA]

posted by ElMaestro  – Denmark, 2015-04-20 17:39 (3658 d 16:18 ago) – Posting: # 14716
Views: 27,816

Hi Helmut,

❝ Statler & Waldorf-day today?

Every day is a Statler & Waldorf day Not good for my blood pressure. Today's major trouble is bacterial DNA, though. Why do people publish sensational papers about having sequenced a new bacterial species when they are not publishing the genome?

❝ cnames <- c("subj", "drug", "seq", "prd", "Cmax", MeanR),

❝ (...) "\nT:", sprintf("%.6f", MeanT), "\n")

Ohdearohgoodnessmelordhavemercy. Yuck.

❝ I’m not promoting SAS.

Don't worry. I know that for a long time already.

There is still something that really bothers me and that is the log Likelihood difference of R versus WNL. With REML I think optimisation switches back and forth be in terms of the estimates. You fiddle with the covariance matrix and read the logLikelihood. Then you fiddle with the fixed effects and read the logLikelihood. Then you fiddle with the covariance matrix, and so forth. Repeat until some convergence criterion is met. ML would only be plain and simple covariance matrix fiddling, I think (?). Thus I have a vague and completely unsubstantiated feeling this could possibly be fixed by proper inputs and that you'd then see agreement about the likelihoods and probably the treatment effects too. It is a guess, and undertunately I do not know my way around with R or mixed models to an extent where I can do it.

Thread locked

Complete thread:

• Reproducing Bear results in Phoenix mittyri 2015-04-14 08:35 [🇷 for BE/BA] 
  • Bug in bear? Helmut 2015-04-14 14:07
    • Bug in bear? mittyri 2015-04-14 20:43
      • Bug in bear? Helmut 2015-04-15 02:17
    • use the same dataset? yjlee168 2015-04-20 10:13
  • Bugs in bear with replicated demo data set? yjlee168 2015-04-15 09:38
  • Bugs fixed in bear? yjlee168 2015-04-17 01:21
    • Bugs fixed in bear? Helmut 2015-04-17 14:25
      • Bugs fixed in bear? yjlee168 2015-04-17 18:50
        • Bugs fixed in bear? ElMaestro 2015-04-17 23:14
          • lme() in bear yjlee168 2015-04-18 11:12
            • lme() in bear ElMaestro 2015-04-18 11:28
              • lme() in bear yjlee168 2015-04-18 12:02
                • lme() in bear ElMaestro 2015-04-18 13:42
                  • Mean means Helmut 2015-04-18 14:08
                    • Mean means ElMaestro 2015-04-18 16:01
                      • Mean means Helmut 2015-04-19 00:59
                        • Mean means ElMaestro 2015-04-19 09:14
                          • LL and AIC Helmut 2015-04-19 11:08
                            • Confused ElMaestro 2015-04-19 11:42
                        • lsmeans for mixed model in R yjlee168 2015-04-19 23:36
                          • lsmeans() & lme() Helmut 2015-04-20 01:28
                            • lsmeans() & lme() yjlee168 2015-04-20 10:23
                            • lsmeans() & lme() ElMaestro 2015-04-20 10:35
                              • keep it simple! Helmut 2015-04-20 14:34
                                • keep it simple!ElMaestro 2015-04-20 15:39
                                  • ML vs. REML Helmut 2015-04-20 16:35
            • lme() with NA in R yjlee168 2015-04-19 23:44
        • Dataset Helmut 2015-04-18 13:30
          • Dataset and lme() in bear yjlee168 2015-04-18 23:14