yjlee168
★★★
avatar
Homepage
Kaohsiung, Taiwan,
2008-05-08 07:57
(5803 d 11:28 ago)

Posting: # 1831
Views: 13,049
 

 stab for R [🇷 for BE/BA]

Dear All,

   Again, we like to announce and share another tool for the stability data analysis using R.
It is called stab for R.
We follow the ICH guideline 'Q1E Evaluation for Stability Data' to design this tool. This guideline describes when extrapolation should be considered as proposing a retest period for a drug substance or a shelf-life of a drug product that extends beyond the period covered by available data from the stability study under the long-term storage condition.

Thanks.

All the best,
--- Yung-jin Lee & Hsin-ya Lee (KMU, Taiwan)
Helmut
★★★
avatar
Homepage
Vienna, Austria,
2008-05-08 15:36
(5803 d 03:49 ago)

@ yjlee168
Posting: # 1836
Views: 11,601
 

 stab for R

Dear Yung-jin Lee & Hsin-ya,

again a fine piece of work!

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
kegoosse
☆    

Belgium,
2017-08-08 17:01
(2424 d 02:24 ago)

@ yjlee168
Posting: # 17683
Views: 9,160
 

 stab for R

Indeed a nice tool.

I do have 1 question. My csv file is deemed as an incorrect dataset.
It contains 3 columns with batch, time and assay information (including column titles!)
This, however, doesn't work. How should the csv input file look like?

(PS I already tried to work without column titles. Alas.)
yjlee168
★★★
avatar
Homepage
Kaohsiung, Taiwan,
2017-08-08 21:19
(2423 d 22:05 ago)

@ kegoosse
Posting: # 17684
Views: 9,154
 

 A single-batch drug stability study?

Dear kegoosse,

It happens only if your dataset was obtained from a single-batch drug stability study (or you selected batch type as single-batch from GUI). If this is the case, there is no need to include "batch" in your dataset. Just delete the "batch" column from your dataset and it should work. There only requires two columns of "time" and "assay" for the dataset of a single-batch drug stability. Sorry about this error. I should put more clear information about this in the window of error pop-up dialogue for next release. Finally, the column title is required. Thanks for bringing this error up here.

[edited] if there is still an error, please try to use period(.) as the decimal point instead of comma (,). See this thread for detailed information.

❝ ...

❝ (PS I already tried to work without column titles. Alas.)


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

2017-08-16 03:18
(2416 d 16:06 ago)

@ yjlee168
Posting: # 17688
Views: 8,984
 

 stab for R

Dear Yung-jin Lee & Hsin-ya Lee,
Chapeau for the stab package in its new GUI look, a real addition to your group of apps through R.
We all are waiting your update for bear2.8.3 to include the missing pieces of the puzzle (2x3x3 designs, HVD) and (NTI Drugs).
Good luck in all the coming.
Greetings.

Mahmoud Teaima, PhD.
yjlee168
★★★
avatar
Homepage
Kaohsiung, Taiwan,
2017-08-16 21:11
(2415 d 22:13 ago)

(edited by yjlee168 on 2017-08-16 21:54)
@ mahmoud-teaima
Posting: # 17696
Views: 9,016
 

 bear is still going...

Dear Dr. Teaima,

I am working on bear right now with partial replicate BE. So I have already started studying the package replicateBE (v1.0.3 and still not available from CRAN yet) developed by Helmut Schütz, Michael Tomashevskiy, and Detlew Labes. Regarding HVD and NTI, I still don't any plan/clue yet. Any suggestion/comment is welcomed. Once I have finished it, I will first announce the update here. So please stay tuned and thanks for reminding me.

❝ We all are waiting your update for bear2.8.3 to include the missing pieces of the puzzle (2x3x3 designs, HVD) and (NTI Drugs).


All the best,
-- Yung-jin Lee
bear v2.9.1:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here
Helmut
★★★
avatar
Homepage
Vienna, Austria,
2017-08-16 21:30
(2415 d 21:55 ago)

@ yjlee168
Posting: # 17697
Views: 9,070
 

 package replicateBE

Hi Yung-jin,

❝ I have already started studying the package replicateBE v1.0.3 …


Please clone or download the current version (1.05.9000) from GitHub.

❝ … and still not available from CRAN yet


We haven’t decided yet whether we will submit it to CRAN. ;-)

❝ Regarding HVD and NTI, I still don't any plan/clue yet.


I reckon a complete implementation of the FDA’s methods is not possible in R.

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
ElMaestro
★★★

Denmark,
2017-08-16 22:52
(2415 d 20:32 ago)

@ Helmut
Posting: # 17698
Views: 8,953
 

 package replicateBE

Hi Hötzi,

❝ I reckon a complete implementation of the FDA’s methods is not possible in R.


Satterthwaite's denominator df? :-)

Pass or fail!
ElMaestro
mittyri
★★  

Russia,
2017-08-16 23:40
(2415 d 19:44 ago)

@ ElMaestro
Posting: # 17699
Views: 9,010
 

 FDA RSABE in R is still a Riemann hypothesis

Hi ElMaestro,

❝ Satterthwaite's denominator df? :-)

Satterthwaite could be solved using lmerTest package.

AFAIK most difficult thing is
REPEATED/GRP=TRT SUB=SUBJ;

I do not give up (yet) with FDA code

PS: will send you some mail

Kind regards,
Mittyri
yjlee168
★★★
avatar
Homepage
Kaohsiung, Taiwan,
2017-08-16 23:53
(2415 d 19:31 ago)

@ Helmut
Posting: # 17700
Views: 8,959
 

 package replicateBE

Dear Helmut,

Thank you so much for your quick response. Now I can see more detailed about the package - replicateBE from GIT. That really helps. I was considering about the output format generated by replicateBE. Will you try to generate SAS-like (or similar at least) outputs in the future (no hurry at all)? I knew that you have done cross-validation with some other software, including SAS. So it seems possible. Let me know if I am wrong.

❝ ...

❝ I reckon a complete implementation of the FDA’s methods is not possible in R.


OK. got it.

All the best,
-- Yung-jin Lee
bear v2.9.1:- created by Hsin-ya Lee & Yung-jin Lee
Kaohsiung, Taiwan https://www.pkpd168.com/bear
Download link (updated) -> here
Helmut
★★★
avatar
Homepage
Vienna, Austria,
2017-08-17 01:49
(2415 d 17:35 ago)

@ yjlee168
Posting: # 17701
Views: 9,003
 

 package replicateBE

Hi Yung-jin,

❝ I was considering about the output format generated by replicateBE. Will you try to generate SAS-like (or similar at least) outputs in the future (no hurry at all)?


The current output is a hobby-project. However, I added goodies inspectors love and commercial software don’t provide out of the box: Which data was used, who performed the analysis, when, on which system, which OS, which version (of the OS, R, package/functions used).

❝ I knew that you have done cross-validation with some other software, including SAS. So it seems possible.


The main purpose was cross-validation. For that we need only the CVwR and the 100(1–2α) CI. What I use in the function calls is foo(..., details=TRUE). I included more human-readable stuff because beta-testers asked for it. I don’t want to send the model tables to the output. However, if one wants them: foo(..., verbose=TRUE).
Couldn’t resist looking at the Type I Error. Try

method.A(path.out = "~/", ola = TRUE, adj = TRUE, data = rds03)

and have a look at the end of the results-file in your home-folder.

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
UA Flag
Activity
 Admin contact
22,957 posts in 4,819 threads, 1,638 registered users;
76 visitors (0 registered, 76 guests [including 7 identified bots]).
Forum time: 18:25 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