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

posted by yjlee168 Homepage – Kaohsiung, Taiwan, 2010-05-05 22:13 (5879 d 06:25 ago) – Posting: # 5302
Views: 30,561

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.6:- 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,653 posts in 4,991 threads, 1,570 registered users;
122 visitors (0 registered, 122 guests [including 22 identified bots]).
Forum time: 04:38 CEST (Europe/Vienna)

In theory, there is no difference between theory and practice.
But, in practice, there is.    Jan L.A. van de Snepscheut

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