-
-
Notifications
You must be signed in to change notification settings - Fork 17
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
[Bug]: Imputation rule 50 % does not work #1208
Comments
The easiest solution and most flexible one would in my view would be:
|
@Melkiades I contacted gRED it looks like we may not have the same guidelines, I will update your asap. |
@Melkiades @ayogasekaram Please put this issue in the archive, it's not needed yet as this issue would tackle only a few studies with very specific data (1 %) and we may not we to be that precise for this summary table. We don't need to do the development of this one unless I come back later with more information that makes us think that we should do it. Thank you! |
What happened?
the function analyze_vars_in_cols.R would have to be refined.
It's quite a complex issue that we didn't think about when you develop the function I guess.
In this example, I was using:
imp_rule = "1/2",
avalcat_var = "AVALC",
na_str = "NE"
The thing is that the n NON BLQ >= 1/2 is not calculated correctly, in the example above, at 288H, there are 4 patients, only two have a pk sample, and 1 is BLQ, one is non BLQ.
The formula to report the statistics should be nNONBLQ >= 1/2, but taking in account all the patients of the treatment arm so: 1 / 4 >= 1/2, FALSE, so it the summary statistics should not appear for the timepoint 288H. (While now the calculation used is 1 / 2 >=1/2)
Another thing which would be more tricky to implement is:
In this one, it looks correct, but the thing is we actually don't have 3 patients in the treatment arm, but 6 patients! But 3 of them where considered as excluded (filter(ANL01FL == "Y")).
So the 1/2 rule should be applied on the n = 6, meaning at timepoint 4H, we have 3/6 nonBLQ values, fine, we have the statistics. But for the timepoint at 360H, we have actually 2/6 (<= 1/2) non BLQ so the summary statistics should not appear.
This issue comes from the fact the we didn't take in account when developing the functions that there could be excluded patients ANL01FL == "Y" (so this should be kept over the workflow to keep the 6 patients as a n() of reference for the calculation nNonBLQ / n >= 1/2). This is actually the first line of code I have in my script:
Would there be the possibility to fix that bugs? please
sessionInfo()
No response
Relevant log output
No response
Code of Conduct
Contribution Guidelines
Security Policy
The text was updated successfully, but these errors were encountered: