log to 3 decimals guess work [🇷 for BE/BA]

posted by yjlee168 Homepage – Kaohsiung, Taiwan, 2010-05-05 22:13 (5479 d 04:54 ago) – Posting: # 5302
Views: 24,958

Dear D, Labes,

Sorry to confuse you again.:-D Actually it is not rounding, but directly transforming in bear.

❝ I guess that you compared the logs calculated to full precision

TotalData$lnCmax <- log(TotalData$Cmax)

❝ lm(lnCmax ~ seq + subj:seq + prd + drug, data=TotalData)

❝ versus my suggestion

lm(log(Cmax) ~ seq + subj:seq + prd + drug, data=TotalData).


Yes, this is exactly the way that we do in bear with the code as follow:
TotalData<-data.frame(subj = as.factor(Total$subj), drug = as.factor(Total$drug), Cmax = Total$Cmax, AUC0t = Total$AUC0t, AUC0INF = Total$AUC0INF,lnCmax = log(Total$Cmax),lnAUC0t = log(Total$AUC0t),lnAUC0INF = log(Total$AUC0INF))

❝ To get Helmut's concern change

TotalData$lnCmax <- round(log(TotalData$Cmax),3)

❝ or even more drastically round to 2 decimals.


No! In bear, only "T/R ratio (%)" and "Power (%)" in "Sample size estimation" use round() function right now. Later (v2.5.0), in Welch t test, we will also use round() to present the point estimate and 90%CIs.

All the best,
-- Yung-jin Lee
bear v2.9.2:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,672 registered users;
78 visitors (0 registered, 78 guests [including 4 identified bots]).
Forum time: 03:08 CEST (Europe/Vienna)

It’s always fun to have your models validated,
but is way more fun to have them trashed.
Finding out you are completely wrong
is a great part of science.    G. Randall Gladstone

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