Lucas ★ Brazil, 2014-02-19 15:27 (4082 d 02:20 ago) Posting: # 12447 Views: 27,042 |
|
Hello everybody. We've been asked to perform a qualitative analysis of the residuals of the ANOVA. Now, I don't know what that analysis is supposed to acomplish, so I don't know what is supposed to be done. Do you guys have any experience with that? Accept my apologies in advance if that is too much of beginner's question. Thank you. |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2014-02-19 16:22 (4082 d 01:25 ago) @ Lucas Posting: # 12450 Views: 25,285 |
|
Hi Lucas, ❝ We've been asked to perform a qualitative analysis of the residuals of the ANOVA. ANVISA? ❝ Now, I don't know what that analysis is supposed to acomplish, so I don't know what is supposed to be done. IMHO, assessing the studentized residuals makes sense. See also Canada’s Guidance, Section 2.3.4.
![]() — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Lucas ★ Brazil, 2014-02-19 16:58 (4082 d 00:49 ago) @ Helmut Posting: # 12453 Views: 25,019 |
|
Hi Helmut! Thanks for the fast response, as always. ❝ ANVISA? No, the sponsor wants us to perform that. The study would be submitted at ANVISA, but I think they want a study that can be (if they ever need) submitted at other agencys. I just love how they ask us that kind of analysis without even explaining their objective with that. So I thought that you guys had experience with that being required in other countries. The sponsor is from England, if that helps. ❝ IMHO, assessing the studentized residuals makes sense. See also Canada’s Guidance, Section 2.3.4. Doesn't seem to be a criteria for saying if the results are good or not, am I right? ❝ [...] and with a little bit of struggle – in Phoenix/WinNonlin. Can u tell me how it can be done at Phoenix? Thanks. |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2014-02-19 17:25 (4082 d 00:21 ago) @ Lucas Posting: # 12456 Views: 25,024 |
|
Hi Lucas, ❝ […] So I thought that you guys had experience with that being required in other countries. Limited experience. Once I had to deal with studies on a HVD (pre-scaling age) where the sponsor – based on previous experiences – stated in the protocols that based on assessing the residuals (boxplots and Shapiro-Wilk test) in case of non-normality the study will be evaluated by a nonparametric method. The Dutch agency didn’t like that. In a hearing in The Hague the told me that “We in The Netherlands don’t believe in nonparametric statistics”. ![]() ❝ Doesn't seem to be a criteria for saying if the results are good or not, am I right? In Canada you may exclude outliers if you have stated criteria in the protocol. ❝ ❝ [...] and with a little bit of struggle – in Phoenix/WinNonlin. ❝ Can u tell me how it can be done at Phoenix? I’m currently struggling with some PopPK/PD-models which send my CPUs’ load to 80%. I don’t want to start concurrently another PHX-session because it would exceed my RAM and Phoenix would start swapping memory to disk. If I don’t come up with something in the next two weeks, remind me. Some hints in my whitepaper (Section 3.6) about EMA’s ABEL. — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Lucas ★ Brazil, 2014-02-19 18:20 (4081 d 23:27 ago) @ Helmut Posting: # 12457 Views: 24,894 |
|
❝ I’m currently struggling with some PopPK/PD-models which send my CPUs’ load to 80%. I don’t want to start concurrently another PHX-session because it would exceed my RAM and Phoenix would start swapping memory to disk. If I don’t come up with something in the next two weeks, remind me. Some hints in my whitepaper (Section 3.6) about EMA’s ABEL. Thanks for the hints. If I can not perform that by reading the whitepaper I'll come back to you later. Best regards. |
Lucas ★ Brazil, 2014-02-20 20:07 (4080 d 21:40 ago) @ Helmut Posting: # 12464 Views: 24,723 |
|
Hi Helmut and ElMaestro! So, I got an answer from the sponsor about what would that qualitative analysis be. Both your's and ElMaestro's answers were right! This is how it would look like according to them: ![]() Now my problem is how I would do that in Phoenix. I haven't got the time to read your whitepaper but I will soon, I only want to show you guys how this analysis would look like. Maybe it helps someone of the forum in the future. Thank u. |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2014-02-20 21:11 (4080 d 20:36 ago) @ Lucas Posting: # 12465 Views: 25,802 |
|
Hi Lucas, ❝ This is how it would look like according to them: […] Interesting! I didn’t count, but ~200 subjects is not common in BE studies. It will not look that nice in the real world. The x-axis of the first plot is strange (negative means – I beg your pardon?). The distribution doesn’t look normal but skewed to the left* and the Q-Q plot looks extremely ugly. That’s actually a great example which EDA would suggest to use another transformation (e.g., ℯx or xn, where n>1). However, let’s see what you can get in PHX (I numbered the panels in your plot):
— Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2014-02-21 20:25 (4079 d 21:22 ago) @ Lucas Posting: # 12472 Views: 25,362 |
|
Hi Lucas, in the following a setup for PHX6.3+ (Clayton/Leslie’s data given in Chow/Liu’s table 6.9.1; validated against table 8.2.3); sequences coded as TR|RT. It seems complicated first, but works also for unbalanced sequences. This is not PHX’ standard model, but all effects fixed (like for EMA and in SAS Proc GLM). You have to code it only once and can use it as a template in other studies.
![]() ![]() ![]() ![]() ![]() In Phoenix Pre-Release V1.4 (Build 6.4.0.511, 2013-12-20) still no overlays in box plots and histograms. ![]() Eidt: I have to correct myself: In Phoenix (all versions) these are standardized (sometimes call internally studentized) residuals, not studentized (externally studentized) residuals. — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
BRN ☆ Turkey, 2014-03-04 16:15 (4069 d 01:32 ago) @ Helmut Posting: # 12554 Views: 24,489 |
|
Dear Helmut, I'm trying to run this process on WNL 6.3 but I'm getting confused in step 3 where I need to map sequence but there's no sequence since I exclude it in the previous step (2) because it's mentioned: ❝ Could you please kindly suggest? — Regards, BRN |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2014-03-04 17:52 (4068 d 23:55 ago) @ BRN Posting: # 12555 Views: 24,637 |
|
Hi BRN, ❝ […] I'm getting confused in step 3 where I need to map sequence but there's no sequence since I exclude it in the previous step (2) because it's mentioned: ❝ ❝ In Step 2 you have to keep Sequence mapped to None . I made a stupid error – sorry for the confusion caused.— Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Rosy ☆ Mexico, 2018-10-15 09:03 (2383 d 09:44 ago) @ Helmut Posting: # 19442 Views: 16,920 |
|
Dear Helmut: I do studentized residual test with R for the outliers detection, but now my regulation prohibid to use R. Can i do with phoenix 7 or another statistical software like ncss. than u Edit: Full quote removed. Please delete everything from the text of the original poster which is not necessary in understanding your answer; see also this post #5! [Helmut] |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2018-10-15 13:43 (2383 d 05:04 ago) @ Rosy Posting: # 19443 Views: 16,811 |
|
Hi Rosy, ❝ I do studentized residual test with R for the outliers detection, but now my regulation prohibid to use R. Bizarre. Consider presenting the documentand to ask the agency why it believes (!) that commercial software (SAS, Phoenix, …) is “better” than open source. Any software has to be validated* – which is easy for base ❝ Can i do with phoenix 7 or another statistical software like ncss. In Phoenix you get only standardized residuals (aka internally studentized residuals). For the setup see above. In SAS both internally and externally studentized residuals are possible. Haven’t used NCSS for 15+ years – no idea. * Examples where a commercial software screwed up:
— Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Rosy ☆ Mexico, 2018-10-18 02:38 (2380 d 16:08 ago) @ Helmut Posting: # 19462 Views: 16,602 |
|
Dear Helmut I know, what's wrong with my regulation, anyway... that's because I'm Mexican (if you remember Guillermo del Toro speach) LOL I will try to follow you directions on phoenix, but someone tell me that it's valid to make a linear regression between Ln_observed and predicted from "Residuals" output from Phoenix and obtain the Student residuals to find outliers, is that right? |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2018-10-18 19:14 (2379 d 23:33 ago) @ Rosy Posting: # 19470 Views: 16,528 |
|
Hi Rosy, ❝ I know, what's wrong with my regulation, anyway... that's because I'm Mexican (if you remember Guillermo del Toro speach) LOL ![]() ❝ I will try to follow you directions on phoenix, … Good luck. ❝ … but someone tell me that it's valid to make a linear regression between Ln_observed and predicted from "Residuals" output from Phoenix What? That’s not interesting because it is rather difficult to “see” anything useful.* You could plot the residuals (y) vs. Ln_predicted (x). Cave: In a 2×2 crossover for every subject you have two residuals with opposite signs. You need only one of them (commonly the first period is used). However, raw residuals are not very useful. Though you get an idea about trends (if they increase/decrease or show a strange shape the model is wrong) or homo-/heteroscedasticity (funnel shape) their magnitude is not telling. That’s why you should standardize/studentize them. ❝ … and obtain the Student residuals to find outliers, is that right? Yep. Though I would not call them outliers (regulators don’t like it) but “aberrant values”. ![]()
— Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Rosy ☆ Mexico, 2018-10-18 20:05 (2379 d 22:42 ago) @ Helmut Posting: # 19471 Views: 16,887 |
|
Dear Helmut, This video, i mean this became a meme for Mexicans, although everything it's black (wrong) you will see the light... Mexican Humor LOL ❝ Good luck. Thank u, but i have trouble, i finish with 4 residuals for each patients (from final variance parameters) i dont know how analyze, because in R appears a table with intra and inter subject residuals for each subject, if i wanna make the same table which values i have to pick up. ❝ ...That’s why you should standardize/studentize them. yep, i mean, i took "Ln_Observed" and "predicted" from Phoenix, then in ncss i run lineal regression and from that i do a regression diagnosis to obtain studentized residual, Do i have to do this with every sequence in every period (whole data) or just for the 1st period, of course the objetive it's to find aberrant values.... i know maybe this no make any sense but i have to comply my regulation thank u for your patience, i'm new in bioequivalence but i hope to learn so much Edit: youtube link is corrected. [Mittyri] |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2014-02-23 03:20 (4078 d 14:27 ago) @ Lucas Posting: # 12473 Views: 24,943 |
|
Hi Lucas, export the stuff to In the Phoenix-project navigate to Stud Res intra > Output Data > Result right-click > Export Filename: StudResIntra The default location is the desktop (which conveniently is also the default import site in R). Start R and paste this code into the console: res <- read.table("StudResIntra.txt", header=TRUE, sep="\t") ![]() Good luck! If you have a Connect license you even don’t have to leave PHX. — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Lucas ★ Brazil, 2014-02-25 20:28 (4075 d 21:18 ago) @ Helmut Posting: # 12497 Views: 24,868 |
|
Wow! Thanks for the help Helmut! I got tied up with other projects and only saw your responses now. Sorry about that. Thank god I got an e-mail from the sponsor a few hours ago saying that we no longer have to perform that analysis. ![]() But these hints will be very helpful for us. Thank you very much for your helpfulness! Best regards. |
ElMaestro ★★★ Denmark, 2014-02-19 16:32 (4082 d 01:15 ago) @ Lucas Posting: # 12451 Views: 25,079 |
|
Hi, it really depends on what your agency wants with that qualitative analysis of residuals. I would do a histogram of residuals, possibly a histogram of within-subject differences. Ideally that would amount to something that looks like a normal distribution, but realistically with a low number of subjects it may just look like the teeth of someone who pissed Mike Tyson seriously off. Looking for outliers (looking in this regards means looking, and not quantifying their potential presence) may be another aspect where this could come in handy. — Pass or fail! ElMaestro |
Helmut ★★★ ![]() ![]() Vienna, Austria, 2014-02-19 17:06 (4082 d 00:41 ago) @ ElMaestro Posting: # 12454 Views: 25,075 |
|
Hi ElMaestro, ❝ Ideally that would amount to something that looks like a normal distribution, but realistically with a low number of subjects it may just look like the teeth of someone who pissed Mike Tyson seriously off. With the data from above: ![]() 48 subjects. Agree with “tooth gaps” in small studies. Below ~30 you will see funny patterns. — Dif-tor heh smusma 🖖🏼 Довге життя Україна! ![]() Helmut Schütz ![]() The quality of responses received is directly proportional to the quality of the question asked. 🚮 Science Quotes |
Lucas ★ Brazil, 2014-02-19 17:13 (4082 d 00:34 ago) (edited on 2014-02-19 20:22) @ ElMaestro Posting: # 12455 Views: 25,116 |
|
Hello! ❝ [...] it may just look like the teeth of someone who pissed Mike Tyson seriously off. That is an extraordinary way to describe it! ![]() ❝ Looking for outliers (looking in this regards means looking, and not quantifying their potential presence) may be another aspect where this could come in handy. I don't know if I'm thinking too "statistically", but that wouldn't be very subjective? One can look and think that there are outliers and other can think that there aren't and both would not be wrong, since there is no criteria. Anyway, your's is a good way to attend to that requirement from the sponsor, but IMHO a qualitative analysis for the residuals would not add any value to the study, from a regulatory perspective. Am I wrong? |
ElMaestro ★★★ Denmark, 2014-02-19 21:40 (4081 d 20:07 ago) @ Lucas Posting: # 12458 Views: 24,921 |
|
Hi Lucas, ❝ I don't know if I'm thinking too "statistically", but that wouldn't be very subjective? One can look and think that there are outliers and other can think that there aren't and both would not be wrong, since there is no criteria. You're right, there is some degree of subjectivity in my proposal. I will argue that this is the case for many if not most qualitative assessments, like for instance a CoA for a injectable where you have a spec like "A colorless or slightly yellowish liquid". ❝ Anyway, your's is a good way to attend to that requirement from the sponsor, but IMHO a qualitative analysis for the residuals would not add any value to the study, from a regulatory perspective. Am I wrong? You're right, at least per principle. Regulatory decisions are as far as I know not based on residuals or their distributions in any territory. However, regulators may ask for anything they like ("Please submit all chromatograms for QCs and calibrators injected during validation, maintenance records from 2009-2011 for freezer ZX/Y353, and the phone number of the hot brunette from I.P. who attended the last pre-submission meeting" etc), and they may ultimately use that additional information against you as they please. No idea why a sponsor would ask you to do this, other than if some newly hired numbercruncher wants to impress her/his boss by asking something that sounds complicated. If Helmut's figure above applies to you data then just submit that, tell the sponsor you think this looks reasonably normal (as in Gaussian) without the presence of definitively aberrant values (you will deliberately not use the term outliers). — Pass or fail! ElMaestro |