Miettinen-Nurminen (MN) method has increasingly been requested by regulatory agencies, rather than the traditional Wald method that is based on the asymptotic normal distribution. This is particularly relevant for non-inferiority trials where it's appropriate for the variance to be estimated under the null hypothesis. Calculate MN Test Statistics within PROC FREQ.
Learning Bisection Method
When I learned the Miettinen and Nurminen Test algorithm from the rate-compare
article (Unstratified and Stratified Miettinen and Nurminen Test), I found that its CI is given by the roots of an equation. In order to reproduce its algorithm, I'd like to learn more about the Bisection method first.
Shift Value for MCMC imputation in Tipping Point Analysis
This is a continuation of the previous article, Tipping Point Analysis in Multiple Imputation Using SAS. In the last post, we talked about the tipping point analysis in monotone imputation, but how to implement TPA in MCMC imputation since the MNAR
statement can only support the shift adjustment in monotone and FCS.
Tipping Point Analysis in Multiple Imputation using SAS
The tipping point analysis has been a useful sensitivity analysis for multiple imputation to assess the robustness of the deviations from the MCAR or MAR assumptions. It aims to find out how severe departures from MAR will overturn the conclusions from primary analysis. If the departures are considered unlikely, this can give strong evidence supporting the treatment effect found in the primary analysis under the MAR assumptions.
Simply Understanding Log Rank Test
The logrank test is the most commonly used statistical test in clinical trials to compare the survival distributions in different treatment groups. We usually just use the logrank results to test whether there is a difference between two survival curves. But what does this difference mean?
Quality control of SDTM using the sdtmchecks package
如果想对SDTM有个快速全面的check,Pinnacle 21毫无疑问是首选,其能帮助我们对SDTM进行CDISC标准验证,发现data quality issue以及mapping不合理的之处。
Estimated LS-means from Multiple Imputation of mice using emmeans package
Symbols count in article: 4.7k Reading time ≈ 4 mins.
Continue with the question in the previous article (Multiple Imputaton - Linear Regression in R), where we just discussed how to compute the pooled coefficients of ANCOVA using mice
package but left out the Ls-means and hypothesis test. Luckly I find out that emmeans
package have wrapped this process inside so we can use it to obtain the pooled Ls-means estimation and p-value straightforward wihout pool
function of mice
.
Multiple Imputaton - Linear Regression in R
Symbols count in article: 5.6k Reading time ≈ 5 mins.
We have discussed the multiple imputation in non-monotone pattern of missingness in the article of Understanding Multiple Imputation in SAS, and sort out how to implement it in SAS. While here, I would like to learn how to use linear regression in multiple imputation to deal with monotone pattern data in R.
使用rtables生成Time-To-Event汇总表
在临床试验中,通常使用SAS来完成统计分析和生成图表,但我们不应该只局限于一种编程方法,而且这个所用的编程语言SAS并不是开源的。毫无疑问SAS能完成的事情,R和Python同样能做;但有些R和Python能做的,SAS却很难完成,我想这就是开源和不开源的区别。
Common Survival analysis of Oncology trials in R
Symbols count in article: 8.7k Reading time ≈ 8 mins.
Time-to-event endpoints are widely used in oncology trials, such as OS and PFS. And survival analysis is a common method for estimating time-to-event endpoints. In this blog, I’d like to make a note of how to summarize the essential results for survival analysis in oncology trials in R and also compare them with SAS.
Evaluation of Best Overall Response per RECIST in R
The Best Overall Response (BOR) is a very common evaluation of efficacy in oncology trials. Usually, it is defined as the best response among all time-point responses from the treatment start until the first disease progression, in the order of CR, PR, SD, PD, and NE per RECIST 1.1. For non-randomized trials, BOR is not only the best among all responses but also requires confirmation for CR and PR to ensure the result is not a measurement error. More details can be found in the RECIST 1.1 document, which I will not expand on here.
Response Rate and Odd Ratio in R and SAS
Symbols count in article: 7k Reading time ≈ 6 mins.
As we know, the objective response rate (ORR) is used as a key endpoint to demonstrate the efficacy of a treatment in oncology and is also valuable for clinical decision making in phase I-II trials, especially in single-arm trials.
Hypothesis testing of MMRM
Contrasts and Hypothesis Tests of emmeans
In the article Definition of least-squares means (LS means), we have known how to compute the LS mean step by step and how to implement it in the emmeans
package that will calculate the estimated mean value for different factor variables and assume the mean value for continuous variables.
Hexo迁移 - 更换ECS服务器
趁着最近阿里云双11的优惠活动,我计划更换下博客所在的ECS服务器(其实为了响应消费降级~),咨询了下售前和售后,最终顺利完成迁移,记录一下迁移过程以备后续所需。
Understanding Mixed Model Repeated Measures (MMRM) in SAS and R
Symbols count in article: 7k Reading time ≈ 6 mins.
Mixed models for repeated measures (MMRM) is widely used for analyzing longitdinal continuous outcomes in randomized clinical trials. Repeated measures refer to multiple measures taken from the same experimental unit, such as a couple of tests over time on the same subject. And the advantage of this model is that it can avoid model misspcification and provide unbiased estimation for data that is missing completely at random (MCAR) or missing at random (MAR).
mcradds R Package
I'm tickled pink to announce the release of mcradds
(version 1.0.1) helps with designing, analyzing and visualization in In Vitro Diagnostic trials.
Releasing R Package to CRAN
Recently, I've been developing my R package - mcradds, which will be my first package released to CRAN. To be honest, finishing coding is just the first step for R package development, whereas I feel like the submission to CRAN is the most challenging for me. This blog is to keep track of something I came across during the submission process to help giving me a reminder when I would develop other packages in next steps. If you are a beginner like me, this blog will be beneficial to you as well.
Convert Plots to Editable Format in R
推荐一个R包(officer
)可以用于生成editable图片在PPT中。这里的editable是指图片中每个元素包括散点、X/Y轴、标签都能修改,常用于图片的再修饰
Multiple Imputation in Non-inferiority and Superiority Trials
In the previous article (Understanding Multiple Imputation in SAS), we talked about how to implement multiple imputation in the SAS procedure to compare the difference between the treatment and placebo groups. Let's look at how to do it in non-inferiority and superiority trials, which differ from common use.