code: def test_rank_test_output_against_R_no_censorship():
      """
      > time <- c(10,20,30,10,20,50)
      > status <- c(1,1,1,1,1,1)
      > treatment <- c(1,1,1,0,0,0)
      > survdiff(Surv(time, status) ~ treatment)
      """
      result = stats.multivariate_logrank_test(
          [10, 20, 30, 10, 20, 50], [1, 1, 1, 0, 0, 0]
      )
      r_p_value = 0.614107
      r_stat = 0.254237
      assert abs(result.p_value - r_p_value) < 10e-6
      assert abs(result.test_statistic - r_stat) < 10e-6
  
imports: import numpy
  import pandas
  import numpy.testing
  import pytest
  from pandas.testing import assert_frame_equal
  from lifelines import statistics
  from lifelines import CoxPHFitter, KaplanMeierFitter, WeibullFitter
  from lifelines.exceptions import StatisticalWarning
  from lifelines.datasets import load_waltons, load_g3, load_lymphoma, load_dd, load_regression_dataset, load_leukemia
label_desc: Use image node_id 7 for calling a global function with example usage: test_rank_test_output_against_R_no_censorship() without return types