cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
The Discovery Summit 2025 Call for Content is open! Submit an abstract today to present at our premier analytics conference.
Choose Language Hide Translation Bar
View Original Published Thread

How to calculate REML Estimates of Varaince Componenets based on USP <1033> Biological Assay Validation?

Luk_s
Level I

Hello, 

I have encountered a problem with calculating the results from REML analysis based on the data and results presented in USP <1033>. I have trouble with computing the data in Table. 4 and getting the same results contained in Table 9. Is there a person who can handle this data and explain how to proceed in making the REML analysis using this example?

Luk_s_0-1729071408615.png

Thanks!

 

1 ACCEPTED SOLUTION

Accepted Solutions
MRB3855
Super User


Re: How to calculate REML Estimates of Varaince Componenets based on USP <1033> Biological Assay Validation?

Hi @Luk_s , and welcome. I'm able to reproduce tables 5 and 6. But, not 9.  It appears the results in table 9 are an example using different data. That can be seen by the difference in Var(error). In table 5 it is 0.000766. In table 9 it is 0.0022. They should be the same if both using the same data. So, table 9 is not based on the data shown in table 4. It is an example from something else, as implied by the statement " Variance component estimates obtained from such an analysis are presented in Table 9" (italics mine); notice it does not refer back to the data in table 4. And, as they correctly point under table 5, " Analyst and media lot have not been included because of the small number of levels (2 levels) for each factor." , so it would be ill-advised to carry out the type of analysis shown in table 9.

 

FWIW, and as ill-advised as this is, I get the following (Table 4, Level 0.5):

MRB3855_1-1729080729836.png

 

 

 

 

View solution in original post

3 REPLIES 3
Luk_s
Level I


Re: How to calculate REML Estimates of Varaince Componenets based on USP <1033> Biological Assay Validation?

For more information, I provide the sources: 1033 BIOLOGICAL ASSAY VALIDATION 

MRB3855
Super User


Re: How to calculate REML Estimates of Varaince Componenets based on USP <1033> Biological Assay Validation?

Hi @Luk_s , and welcome. I'm able to reproduce tables 5 and 6. But, not 9.  It appears the results in table 9 are an example using different data. That can be seen by the difference in Var(error). In table 5 it is 0.000766. In table 9 it is 0.0022. They should be the same if both using the same data. So, table 9 is not based on the data shown in table 4. It is an example from something else, as implied by the statement " Variance component estimates obtained from such an analysis are presented in Table 9" (italics mine); notice it does not refer back to the data in table 4. And, as they correctly point under table 5, " Analyst and media lot have not been included because of the small number of levels (2 levels) for each factor." , so it would be ill-advised to carry out the type of analysis shown in table 9.

 

FWIW, and as ill-advised as this is, I get the following (Table 4, Level 0.5):

MRB3855_1-1729080729836.png

 

 

 

 

Luk_s
Level I


Re: How to calculate REML Estimates of Varaince Componenets based on USP <1033> Biological Assay Validation?

Hi @MRB3855 thank you kindly for your attention to the problem which I have presented. I was convinced that I was doing something wrong in developing the data for the model. I also did a check and we come out with the same results, which I'm really pleased about. I have dealt with the previous data table ones as well, however, I was surprised by the results obtained in this example. 

Luk_s_0-1729082568200.png

Thanks!