Update! [Software]
Hi Helmut and others,
Here is a strange one I just experience. 100% final plasma data from a 3-way (2 formulation vs reference) study. Both myself and the CRO used SAS Proc GLM to compute A vs C, B vs C (alpha=0.05 yes). Results? My lower 90% CI for formulation 1 was 80.24, CRO result showed 79.65 at the lower limit of the 90% CI. I checked my individual parameter values with theirs and they all matched! Okay their report only show 2 decimal places and I use 4 decimals. The values are in tenths (? I mean numbers like 20.50, etc, not into hundredths). The difference in arithmetic means are like .5 for some reason.
Geometric (least square) means are also off
Mines are 31.82 (T), 35.50 (R), Intra CV = 24.44
Theirs are 31.8 (T), 35.55 (R), Intrasub CV = 25.8
??? My suspicion is the number of decimal places they used. The plasma-time data is 3 decimal places. No time deviation (0-24 hrs sampling, all in house)
John
P.S. I tried running WinNonlin and the output confirmed my numbers as well.
❝ With the correct coding I never saw different results from SAS and PHX. Is this a hypothetical example or from the ‘real world’? On the other hand different results in higher-order Xovers are possible (again independent from software)...
Here is a strange one I just experience. 100% final plasma data from a 3-way (2 formulation vs reference) study. Both myself and the CRO used SAS Proc GLM to compute A vs C, B vs C (alpha=0.05 yes). Results? My lower 90% CI for formulation 1 was 80.24, CRO result showed 79.65 at the lower limit of the 90% CI. I checked my individual parameter values with theirs and they all matched! Okay their report only show 2 decimal places and I use 4 decimals. The values are in tenths (? I mean numbers like 20.50, etc, not into hundredths). The difference in arithmetic means are like .5 for some reason.
Geometric (least square) means are also off
Mines are 31.82 (T), 35.50 (R), Intra CV = 24.44
Theirs are 31.8 (T), 35.55 (R), Intrasub CV = 25.8
??? My suspicion is the number of decimal places they used. The plasma-time data is 3 decimal places. No time deviation (0-24 hrs sampling, all in house)
John
P.S. I tried running WinNonlin and the output confirmed my numbers as well.
Complete thread:
- SAS vs Winnonlin jag009 2013-07-24 17:30 [Software]
- SAS vs Winnonlin Helmut 2013-07-24 17:55
- SAS vs Winnonlin jag009 2013-07-24 22:56
- SAS vs Winnonlin Helmut 2013-07-24 23:53
- SAS vs Winnonlin jag009 2013-07-25 00:56
- EMA crippled approach d_labes 2013-07-25 09:13
- Hypothesis jag009 2013-07-25 15:38
- Two at a Time? Or All at Once? Helmut 2013-08-10 13:52
- Heteroscedasticity Helmut 2013-07-27 18:30
- Hypothesis jag009 2013-07-25 15:38
- EMA crippled approach d_labes 2013-07-25 09:13
- SAS vs Winnonlin jag009 2013-07-25 00:56
- SAS vs Winnonlin Helmut 2013-07-24 23:53
- Update!jag009 2013-07-27 05:59
- Rounding limbo? Helmut 2013-07-27 14:53
- Rounding limbo? jag009 2013-07-27 22:02
- Rounding limbo? Helmut 2013-07-27 14:53
- Winnonlin: exclude volunteers mittyri 2014-01-10 08:05
- Winnonlin: exclude incomplete data! Helmut 2014-01-10 13:41
- SAS vs Winnonlin jag009 2013-07-24 22:56
- SAS vs Winnonlin ElMaestro 2013-07-24 18:06
- SAS vs Winnonlin jag009 2013-07-24 22:59
- SAS vs Winnonlin Helmut 2013-07-24 17:55