Bioequivalence and Bioavailability Forum • A small point for the code

A small point for the code [GxP / QC / QA]

posted by sameep – India, 2022-06-01 13:17 (36 d 05:02 ago) – Posting: # 23037
Views: 392

Dear Helmut,

» for (i in 3:n)

Works for all crossover datasets, whereas

» for (i in subj:n) {

Works only if the sequence for the first two subjects is the same. If the first two subject's sequences are different (i.e. TR and RT), then it gives the below error, as calculations are not feasible due to zero degree of freedom.

Error in power.TOST(CV = res$CV[i - 1], theta0 = res$PE[i - 1], n = length(unique(tmp$Subject))) : 

  n too small. Degrees of freedom <1!

In addition: Warning messages:

1: In qt(a, object$df.residual) : NaNs produced

2: In anova.lm(m) :

  ANOVA F-tests on an essentially perfect fit are unreliable

» I was only talking about importing the data from xls(x).

This is what I meant.

» I’m still working on the script (in the meantime a monster with 525 lines of code), allowing not only 2×2×2 but also higher-order crossovers. Not fail-safe.

Thanks in advance, looking forward to the goodie.

Regards,
Sameep.

Complete thread:

Another two... Ohlbe 2021-09-16 23:41 [GxP / QC / QA]
- Another two... ElMaestro 2021-09-17 08:42
  - Dates Ohlbe 2021-09-17 10:15
  - Some more info Ohlbe 2021-09-17 20:28
    - EU Article 31 referral started Ohlbe 2022-01-24 17:47
- Cheating Helmut 2021-09-17 16:50
  - Increased variability Ohlbe 2021-09-17 17:56
    - Increased variability ElMaestro 2021-09-17 18:51
      - Visualization Helmut 2021-09-18 15:26
- Another two... jag009 2021-09-19 09:18
- Another two... ElMaestro 2022-05-24 12:00
  - Blind monitors or greedy sponsors? Helmut 2022-05-24 12:51
    - Blind monitors or greedy sponsors? Ohlbe 2022-05-24 14:21
      - Blind monitors or greedy sponsors? Helmut 2022-05-24 16:35
    - Blind monitors or greedy sponsors? ElMaestro 2022-05-25 08:11
      - (Cumulative) T/R-ratio vs. time Helmut 2022-05-25 09:19
        
        Sponsors and CRO selection Ohlbe 2022-05-25 10:53
        
        I still think that T/R-ratios are useful Helmut 2022-05-25 12:04
        
        T/R-ratios are useful Ohlbe 2022-05-25 14:40
        
        T/R-ratios are useful Helmut 2022-05-25 15:35
        Bust the Buster Helmut 2022-05-26 15:32
        
        complicate the assessor's life mittyri 2022-05-26 18:28
        
        complicate the assessor's life PharmCat 2022-05-27 09:57
        complicate the assessor's life Helmut 2022-05-27 10:11
        
        Thanks for Busting the Buster sameep 2022-05-31 08:34
        
        Thanks for Busting the Buster Helmut 2022-05-31 15:22
        
        A small point for the codesameep 2022-06-01 13:17
        
        A small point for the code Helmut 2022-06-01 14:11